Método general para la suma indefinida $\def\nn{\mathbb{N}}$ $\def\zz{\mathbb{Z}}$
Definir el operador de diferencia hacia delante:
$D = ( \text{function $ f $ on $ \zz $} \mapsto ( \text{int $ n $} \mapsto f(n+1) - f(n) ) )$
Es decir, para cualquier función $f$ en $\zz$ y $n \in \zz$ , $D(f)(n) = f(n+1) - f(n)$ .
Si pensamos en las funciones como secuencias (infinitas en ambas direcciones), tomar la diferencia hacia delante significa sustituir cada término por el valor del término siguiente menos él mismo. Por ejemplo, lo siguiente muestra lo que sucede cuando se toma repetidamente la diferencia hacia adelante de la secuencia de cubos:
...,-27,-8,-1, 0, 1, 8,27,...
..., 19, 7, 1, 1, 7,19,37,...
...,-12,-6, 0, 6,12,18,24,...
..., 6, 6, 6, 6, 6, 6, 6,...
..., 0, 0, 0, 0, 0, 0, 0,...
..., 0, 0, 0, 0, 0, 0, 0,...
Entonces tenemos:
$D\left( \text{int $ n $} \mapsto \binom{n}{k+1} \right) = \left( \text{int $ n $} \mapsto \binom{n}{k} \right)$ para cualquier $k \in \zz$ .
Esto es de esperar porque se deduce directamente del triángulo de Pascal, especialmente si definimos $\binom{n}{k}$ utilizando el triángulo.
Esto significa que si tenemos cualquier función $f$ en $\zz$ tal que $f(n) = \sum_{k=0}^\infty a_k \binom{n}{k}$ para cualquier $n \in \zz$ , entonces obtenemos:
$D(f)(n) = \sum_{k=0}^\infty a_{k+1} \binom{n}{k}$ para cualquier $n \in \zz$ .
Desde una perspectiva de alto nivel, ésta es la versión discreta de la serie de Taylor, y de hecho para tal función vemos fácilmente que $f(n) = \sum_{k=0}^\infty D^k(f)(0) \binom{n}{k}$ para cualquier $n \in \zz$ porque $\binom{0}{0} = 1$ mientras que $\binom{0}{k} = 0$ para cualquier $k \in \nn^+$ .
Esto funciona para cualquier función polinómica $f$ en $\zz$ ya que $D^k(f)$ es la función cero una vez $k$ es mayor que el grado de $f$ por lo que podemos utilizarlo para encontrar inmediatamente la serie de $(\text{int n} \mapsto n^3)$ y luego simplemente tomar la antidiferencia desplazando los coeficientes de la serie hacia el otro lado. La constante indeterminada que aparece desaparecerá una vez que realicemos una suma definida como si quisiéramos la suma de la primera $m$ cubos.
Suma de $p$ poderes
Por ejemplo, si queremos $\sum_{k=1}^{n-1} k^2$ primero hallamos las diferencias hacia adelante iteradas de la secuencia de cuadrados $( n^2 )_{n \in \zz}$ :
...,0,1,4,9,...
...,1,3,5,...
...,2,2,...
...,0,...
Así que inmediatamente obtenemos $n^2 = 0 \binom{n}{0} + 1 \binom{n}{1} + 2 \binom{n}{2}$ y por lo tanto $\sum_{k=0}^{n-1} = 0 \binom{n}{1} + 1 \binom{n}{2} + 2 \binom{n}{3}$ . Es trivial simplificarlo a $\sum_{k=0}^{n-1} = \frac{1}{6}(n-1)n(2n-1)$ .
Eficacia de cálculo
Esto es mucho más eficiente que el método dado por ncmathsadist porque la serie que utiliza coeficientes binomiales es fácil de manipular y de calcular. Para la suma de $p$ -poderes que sólo necesitamos $O(p^2)$ operaciones aritméticas para hallar las diferencias hacia adelante y luego $O(p^2)$ más para simplificar la forma en serie a una forma polinómica estándar. En cambio, el otro método requiere $O(p^3)$ operaciones aritméticas.
Suma indefinida de no polinomios
Además, para una amplia clase de funciones no polinómicas, todavía podemos calcular la suma indefinida sin utilizar la serie, utilizando el análogo discreto de la integración por partes, aquí llamado suma por partes.
Para derivarlo, basta con comprobar que $D(f \times g)(n) = f(n+1) g(n+1) - f(n) g(n) = f(n+1) D(g)(n) - D(f)(n) g(n)$ y así obtenemos la regla del producto:
$D(f \times g) = R(f) \times D(g) + D(f) \times g$
donde $R$ es el operador de desplazamiento a la derecha definido como
$R = ( \text{function $ f $ on $ \zz $} \mapsto ( \text{int $ n $} \mapsto f(n+1) ) )$
Es decir, para cualquier función $f$ en $\zz$ y $n \in Z$ , $R(f)(n) = f(n+1)$ .
Por comodidad definimos también el operador de suma:
$S = ( \text{function $ f $ on $ \zz $} \mapsto ( \text{int $ n $} \mapsto \sum_{k=0}^{n-1} f(k) ) )$
Entonces tenemos la importante propiedad de que $DS(f) = f$ para cualquier función $f$ en $\zz$ análogo al teorema fundamental del cálculo.
Ahora, sustituyendo $f$ con $S(f)$ en la regla del producto y sumando a ambos lados obtenemos la suma por partes:
$S( f \times g ) = S(f) \times g - S( R(S(f)) \times D(g) ) + c$ para alguna función constante $c$ en $\zz$ .
Ejemplo de suma indefinida
Así podemos calcular fácilmente cosas como $\sum_{k=1}^n k^3 3^k$ aplicándolo tres veces, reduciendo cada vez el grado de la parte polinómica. Hay otras formas de conseguirlo utilizando la diferenciación, pero este método es puramente discreto.
5 votos
Consideremos el examen del Teorema de Faulhaber: theoremoftheday.org/Binomial/Faulhaber/TotDFaulhaber.pdf
2 votos
Creo que es una mala idea ver la fórmula de Faulhaber por sí misma... no mejorarás ningún conocimiento o comprensión, sólo verás una receta. El punto clave es entender los polinomios de Bernoulli que son la solución para la ecuación $\Delta f_n (k)=nk^{n-1}$ donde $f_n(k)$ son polinomios de grado n (los polinomios de Bernoulli) y después ver la relación con esta pregunta. Puedes probar las primeras páginas de este texto por ejemplo.
0 votos
Quizá también habría que preguntarse si se quiere (en Cálculo 2). Para calcular una suma de Riemann, sólo necesitamos conocer el coeficiente del primer término de la fórmula (ya que los demás términos se reducen a cero a medida que $n$ se hace arbitrariamente grande). Por ejemplo $$\sum^{n}_{i=1}i^2 = \frac{1}{3}n^3+an^2+bn+c$$ es suficiente para calcular el área bajo una parábola.