主题
Search

幂和


通常考虑两种幂和。第一种是一组 p 次幂,即一组 n 个变量 x_k 的和,

 S_p(x_1,...,x_n)=sum_(k=1)^nx_k^p,
(1)

第二种是特殊情况 x_k=k,即,

 S_p(n)=sum_(k=1)^nk^p.
(2)

一般幂和在统计学中很常见。例如,k 统计量最常根据幂和定义。幂和通过 牛顿-吉拉德公式对称多项式 相关。

x 的 k 次幂乘以 k 的和可以通过解析方式给出:

 sum_(k=0)^nkx^k=(x-(n+1)x^(n+1)+nx^(n+2))/((x-1)^2).
(3)

其他解析和包括:

(sum_(k=0)^(infty)x^k)^p=(1-x)^(-p)
(4)
=1/((p-1)!)sum_(n=0)^(infty)((n+p-1)!)/(n!)x^n
(5)
=1/((p-1)!)sum_(n=0)^(infty)(n+1)_(p-1)x^n
(6)

对于 |x|<1,其中 (n)_p 是一个 波赫哈默尔符号。有限版本具有简洁的闭合形式

 (sum_(k=0)^nx^k)^p=1/((p-1)!)sum_(k=0)^(np)((n-|n-k|+p-1)!)/((n-|n-k|)!)x^k
(7)

对于 p=1 和 2。另一个和由下式给出:

 (sum_(n=0)^inftya_nx^n)^2=sum_(n=0)^inftya_n^2x^(2n)+2sum_(n=1; i+j=n; i<j)^inftya_ia_jx^n.
(8)

整数 和的解析解是

S_p(n)=sum_(k=1)^(n)k^p
(9)
=zeta(-p)-zeta(-p,1+n)
(10)
=H_n^((-p)),
(11)

其中 zeta(z)黎曼 zeta 函数zeta(z;a)赫尔维茨 zeta 函数,而 H_n^((k)) 是广义 调和数。对于 p正整数 的特殊情况,福尔哈伯公式 明确给出了

 S_p(n)=1/(p+1)sum_(k=1)^(p+1)(-1)^(delta_(kp))(p+1; k)B_(p+1-k)n^k,
(12)

其中 delta_(kp)克罗内克 delta(n; k)二项式系数,而 B_k伯努利数。同样正确的是,正如伯努利 (Boyer 1943) 所述,这种展开式中各项的 系数 之和为 1。

伯努利使用了 图形数三角形 的性质:

 sum_(i=0)^na_(ij)=((n+1)a_(nj))/(j+1),
(13)

以及他通过归纳推导出的 a_(nj) 的形式,以计算高达 n=10 的和(Boyer 1968,第 85 页)。对于 p in Z>0,和由下式给出:

 sum_(k=1)^nk^p=((B+n+1)^([p+1])-B^([p+1]))/(p+1),
(14)

其中 符号 B^([k]) 表示所讨论的量被提升到适当的 k,并且所有 形式为 B^m 的项都用相应的伯努利数 B_m 替换。用 和显式表示为:

sum_(k=1)^(n)k^p=n^p+sum_(k=0)^(p)(B_kp!)/(k!(p-k+1)!)n^(p-k+1)
(15)
=sum_(k=1)^(p+1)((-1)^(p-k+1)B_(p-k+1)p!)/(k!(p-k+1)!)n^k
(16)
=sum_(k=1)^(p+1)b_(pk)n^k,
(17)

其中

 b_(pk)=((-1)^(p-k+1)B_(p-k+1)p!)/(k!(p-k+1)!).
(18)

同样正确的是,各项 b_(pk)系数 之和为 1,

 sum_(k=1)^(p+1)b_(pk)=1,
(19)

伯努利在没有证明的情况下陈述了这一点。

S_p(n) 的 双级数 解由下式给出:

 S_p(n)=sum_(i=1)^psum_(j=0)^(i-1)(-1)^j(i-j)^p(n+p-i+1; n-i)(p+1; j).
(20)

计算 p=1, ..., 10 的和得到

sum_(k=1)^(n)k=1/2(n^2+n)
(21)
sum_(k=1)^(n)k^2=1/6(2n^3+3n^2+n)
(22)
sum_(k=1)^(n)k^3=1/4(n^4+2n^3+n^2)
(23)
sum_(k=1)^(n)k^4=1/(30)(6n^5+15n^4+10n^3-n)
(24)
sum_(k=1)^(n)k^5=1/(12)(2n^6+6n^5+5n^4-n^2)
(25)
sum_(k=1)^(n)k^6=1/(42)(6n^7+21n^6+21n^5-7n^3+n)
(26)
sum_(k=1)^(n)k^7=1/(24)(3n^8+12n^7+14n^6-7n^4+2n^2)
(27)
sum_(k=1)^(n)k^8=1/(90)(10n^9+45n^8+60n^7-42n^5+20n^3-3n)
(28)
sum_(k=1)^(n)k^9=1/(20)(2n^(10)+10n^9+15n^8-14n^6+10n^4-3n^2)
(29)
sum_(k=1)^(n)k^(10)=1/(66)(6n^(11)+33n^(10)+55n^9-66n^7+66n^5-33n^3+5n)
(30)

或因式分解形式为,

sum_(k=1)^(n)k=1/2n(n+1)
(31)
sum_(k=1)^(n)k^2=1/6n(n+1)(2n+1)
(32)
sum_(k=1)^(n)k^3=1/4n^2(n+1)^2
(33)
sum_(k=1)^(n)k^4=1/(30)n(n+1)(2n+1)(3n^2+3n-1)
(34)
sum_(k=1)^(n)k^5=1/(12)n^2(n+1)^2(2n^2+2n-1)
(35)
sum_(k=1)^(n)k^6=1/(42)n(n+1)(2n+1)(3n^4+6n^3-3n+1)
(36)
sum_(k=1)^(n)k^7=1/(24)n^2(n+1)^2(3n^4+6n^3-n^2-4n+2)
(37)
sum_(k=1)^(n)k^8=1/(90)n(n+1)(2n+1)(5n^6+15n^5+5n^4-15n^3-n^2+9n-3)
(38)
sum_(k=1)^(n)k^9=1/(20)n^2(n+1)^2(n^2+n-1)(2n^4+4n^3-n^2-3n+3)
(39)
sum_(k=1)^(n)k^(10)=1/(66)n(n+1)(2n+1)(n^2+n-1)(3n^6+9n^5+2n^4-11n^3+3n^2+10n-5)
(40)

(OEIS A064538A079618)。

Sumk

特殊情况 S_1(n)=n(n+1)/2 的一个简单的图形证明也可以通过构建一系列盒子堆栈来给出,每个盒子宽度为 1 单位,高度为 k 单位,其中 k=1、2、...、n。现在在顶部添加一个旋转副本,如上图所示。请注意,所得图形的 宽度n高度n+1,因此 面积n(n+1)。所需的和是它的一半,因此和中盒子的 面积n(n+1)/2。由于盒子是单位宽度,因此这也是和的值。

S_1(n)=n(n+1)/2 也可以使用第一个 欧拉-麦克劳林积分公式 计算

 sum_(k=1)^nf(k)=int_1^nf(x)dx+1/2f(1)+1/2f(n)+1/(2!)B_2[f^'(n)-f^'(1)]+...
(41)

其中 f(k)=k。那么

sum_(k=1)^(n)k=int_1^nxdx+1/2·1+1/2·n+1/6(1-1)+...
(42)
=1/2(n^2-1)-1/2+h+1/2n
(43)
=1/2n(n+1).
(44)
CubeSquare

令人惊讶的恒等式

S_3(n)=sum_(k=1)^(n)k^3
(45)
=(sum_(k=1)^(n)k)^2,
(46)

被称为 尼科马科斯定理,也可以用图形方式说明(Wells 1991,第 198-199 页)。

Schultz (1980) 表明,可以通过写出以下公式找到和 S_p(n)

 S_p(n)=c_(p+1)n^(p+1)+...+c_1n
(47)

并求解 p+1 个方程组:

 sum_(i=j+1)^(p+1)(-1)^(i-j+1)(i; j)c_i=delta_(j,p)
(48)

对于 j=0, 1, ..., p (Guo 和 Qi 1999) 获得,其中 delta_(j,p)克罗内克 delta。例如,对于 p=2 要解的三个方程是

0=c_1-c_2+c_3
(49)
0=2c_2-3c_3
(50)
1=3c_3,
(51)

得到 c_1=1/6c_2=1/2c_3=1/3,或

 S_2(n)=1/6n+1/2n^2+1/3n^3,
(52)

正如预期的那样。

S_i(n) 通过以下公式与 二项式定理 相关:

 (1+n)^(k+1)=1+sum_(i=0)^k(k+1; i)S_i(n)
(53)

(Guo 和 Qi 1999)。


另请参阅

丢番图方程福尔哈伯公式多重次数方程尼科马科斯定理级数

使用 Wolfram|Alpha 探索

参考文献

Boyer, C. B. "Pascal's Formula for the Sums of Powers of the Integers." Scripta Math. 9, 237-244, 1943.Boyer, C. B. A History of Mathematics. New York: Wiley, 1968.Brualdi, R. A. Introductory Combinatorics, 3rd ed. New York: Elsevier, p. 119, 1997.Cao, J.-T. "A Method of Summing Series and Some Corollaries" [Chinese]. Math. Pract. Th. 20, 77-84, 1990.Conway, J. H. and Guy, R. K. The Book of Numbers. New York: Springer-Verlag, p. 106, 1996.Guo, S.-L. and Qi, F. "Recursion Formulae for sum_(m=1)^(n)m^k." Z. Anal. Anwendungen 18, 1123-1130, 1999.Havil, J. Gamma: Exploring Euler's Constant. Princeton, NJ: Princeton University Press, pp. 81-82, 2003.Schultz, H. J. "The Sums of the kth Powers of the First n Integers." Amer. Math. Monthly 87, 478-481, 1980.Sloane, N. J. A. Sequences A064538 and A079618 in "The On-Line Encyclopedia of Integer Sequences."Struik, D. A Source Book in Mathematics, 1200-1800. Cambridge, MA: Harvard University Press, 1969.Wells, D. The Penguin Dictionary of Curious and Interesting Geometry. London: Penguin, pp. 198-199, 1991.Yang, B.-C. "Formulae Related to Bernoulli Number and for Sums of the Same Power of Natural Numbers" [Chinese]. Math. Pract. Th. 24, 52-56 and 74, 1994.Zhang, N.-Y. "Euler's Number and Some Sums Related to Zeta Function" [Chinese]. Math. Pract. Th. 20, 62-70, 1990.

在 Wolfram|Alpha 中引用

幂和

请这样引用

Weisstein, Eric W. "幂和。" 来自 MathWorld——Wolfram 网络资源。 https://mathworld.net.cn/PowerSum.html

主题分类