Polinomios para sumas de potencias de progresiones aritméticas

Son polinomios en función de una variable que, cuando la variable coincide con el número de sumandos, calculan la suma de potencias con bases en progresión aritmética y exponente igual al número anterior a su grado. El problema es por tanto encontrar polinomios tales que:

$S_n^p(h,d)=\sum _{k=0}^{n-1}(h+kd{)}^p=h^p+(h+d{)}^p+\cdots +(h+(n-1)d{)}^p,$

con variable $n$ y parámetros $p,h,d$ del plinomio $S_n^p(h,d)$; $n$ y $p$ enteros no negativos, $h$ primer término de una progresión aritmética y $d\ne 0$ diferencia de la misma progresión, siendo $h$ y $d$ cualquier número real o complejo

$S_n^p(1,1)=\sum _{k=0}^{n-1}(1+k{)}^p=\sum _{k=1}^n{k}^p=1^p+2^p+\dots +n^p$ son los polinomios identificados por la fórmula de Faulhaber presentada póstumamente por Jacob Bernoulli en 1713;

$S_n^p(0,1)=\sum _{k=0}^{n-1}{k}^p=0^p+1^p+\dots +(n-1{)}^p$ son los polinomios diferenciándose de los anteriores sólo en el signo de un monomio de grado p;

$S_n^p(1,2)=\sum _{k=0}^{n-1}(1+2k{)}^p=1^p+3^p+\dots +(2n-1{)}^p$ son los polinomios por sumas de potencias de números impares sucesivos.

Método con matrices

Para cualquier $m$ entero positivo, el caso general se resuelve mediante la siguiente fórmula:

$\overrightarrow{S_n}(h,d)=T(h,d)A^{-1}\overrightarrow{N_n},$ donde

$[\overrightarrow{S_n}(h,d){]}_r=S_n^{r-1}(h,d),[\overrightarrow{N_n}{]}_r=n^r,$

$[T(h,d){]}_{r,c}=\{\begin{array}{cc}0,& \text{si }c>r,\\ (\genfrac{}{}{0ex}{}{r-1}{c-1})h^{r-c}{d}^{c-1}& \text{si }c\le r.\end{array},[A{]}_{r,c}=\{\begin{array}{cc}0,& \text{si }c>r,\\ (\genfrac{}{}{0ex}{}{r}{c-1}),& \text{si }c\le r,\end{array},$

con $1\le r\le m$ y $1\le c\le m$ donde $r$ (fila), $c$ (columna) y $m$ (orden de la matriz) son enteros.

Ejemplo

La fórmula en el caso particular  $m=5(p=0,1,...,m-1)$ se convierte en :  

$\left(\begin{array}{c}{S}_n^0(h,d)\\ S_n^1(h,d)\\ S_n^2(h,d)\\ S_n^3(h,d)\\ S_n^4(h,d)\end{array}\right)=\left(\begin{array}{ccccc}1& 0& 0& 0& 0\\ h& d& 0& 0& 0\\ h^2& 2hd& d^2& 0& 0\\ h^3& 3h^{2}d& 3h{d}^2& d^3& 0\\ h^4& 4h^{3}d& 6h^2{d}^2& 4h{d}^3& d^4\end{array}\right){\left(\begin{array}{ccccc}1& 0& 0& 0& 0\\ 1& 2& 0& 0& 0\\ 1& 3& 3& 0& 0\\ 1& 4& 6& 4& 0\\ 1& 5& 10& 10& 5\end{array}\right)}^{-1}\left(\begin{array}{c}n\\ n^2\\ n^3\\ n^4\\ n^5\end{array}\right)$

Y en el caso especial $m=5,h=1,d=2$, calcula la suma de los $n$ primeros números impares consecutivos

Calculando la matriz T(h,d), cuyos elementos siguen el teorema del binomio con los valores asignados, es decir, T(1,2), y hallando la matriz inversa de la matriz triangular inferior A obtenida a partir del triángulo de Pascal privado del último elemento de cada fila (matriz formada a partir del números de Bernoulli, mostrada en rojo), tenemos :

$T(1,2)=\left(\begin{array}{ccccc}1& 0& 0& 0& 0\\ 1& 2& 0& 0& 0\\ 1& 4& 4& 0& 0\\ 1& 6& 12& 8& 0\\ 1& 8& 24& 32& 16\end{array}\right),A^{-1}=\left(\begin{array}{ccccc}1& 0& 0& 0& 0\\ -\frac{1}{2}& \frac{1}{2}& 0& 0& 0\\ \frac{1}{6}& -\frac{1}{2}& \frac{1}{3}& 0& 0\\ 0& \frac{1}{4}& -\frac{1}{2}& \frac{1}{4}& 0\\ -\frac{1}{30}& 0& \frac{1}{3}& -\frac{1}{2}& \frac{1}{5}\end{array}\right)$

multiplicando las filas por las columnas de las dos matrices se obtiene

$\left(\begin{array}{c}{S}_n^0(1,2)\\ S_n^1(1,2)\\ S_n^2(1,2)\\ S_n^3(1,2)\\ S_n^4(1,2)\end{array}\right)=\left(\begin{array}{ccccc}1& 0& 0& 0& 0\\ 0& 1& 0& 0& 0\\ -\frac{1}{3}& 0& \frac{4}{3}& 0& 0\\ 0& -1& 0& 2& 0\\ \frac{7}{15}& 0& -\frac{8}{3}& 0& \frac{16}{5}\end{array}\right)\left(\begin{array}{c}n\\ n^2\\ n^3\\ n^4\\ n^5\end{array}\right)=\left(\begin{array}{c}n\\ n^2\\ -\frac{1}{3}n+\frac{4}{3}{n}^3\\ -n^2+2n^4\\ \frac{7}{15}n-\frac{8}{3}{n}^3+\frac{16}{5}{n}^5\end{array}\right).$

y por tanto:

$S_n^0(1,2)=n,S_n^1(1,2)=n^2,S_n^2(1,2)=-\frac{1}{3}n+\frac{4}{3}{n}^3,S_n^3(1,2)=-n^2+2n^4,S_n^4(1,2)=\frac{7}{15}n-\frac{8}{3}{n}^3+\frac{16}{5}{n}^5.$. Por último, si interesan las sumas de los tres primeros sumandos

$S_3^0(1,2)=3,S_3^1(1,2)=9,S_3^2(1,2)=35,S_3^3(1,2)=153,S_3^4(1,2)=707=1^4+3^4+5^4.$

Método con polinomios de Bernoulli

La siguiente fórmula resuelve el problema de forma implícita utilizando polinomios de Bernoulli: $S_n^p(h,d)=\frac{d^p}{p+1}(B_{p+1}(n+\frac{h}{d})-B_{p+1}(\frac{h}{d}))$

En particular:

$S_n^p(1,1)=\frac{B_{p+1}(n+1)-B_{p+1}(1)}{p+1}$

$S_n^p(0,1)=\frac{B_{p+1}(n)-B_{p+1}(0)}{p+1}$

$S_n^p(1,2)=2^p\frac{B_{p+1}(n+\frac{1}{2})-B_{p+1}(\frac{1}{2})}{p+1}$