17
算术
二项式
ODDS
D.B.福克斯,M.B.福克斯
(1+x)2=1+2x+x2
(1+x)3=1+3x+3x2+x3
任何小学生都知道这些公式。数字 1、2、1; 1, 3, 3, 1,以及以类似方式获得的数字的 4, 5, … 次方,称为生物项系数。本文讨论生物系数的各种属性。第一段阐述了“一般理论”:这里证明的大部分定理也包含在学校课程中。在第二段中,我们将展示如何在将任何生物系数除以任何素数时轻松找到余数。本段中的一些陈述是作为假设提出的。也许读者«Kvant»会尝试证明它们。
18
§ 1. 二项式系数的定义和最简单的性质
如果二项式(二项式)1+x用某种自然幂n进行加权,那么显然,结果是一个n次数的多项式(即x出现在该多项式中的最高次数将等于n)。例如:
(1+x)0=1,
(1+x)1=1+x,
(1+x)2=1+2x+x2,
(1+x)3=1+3x+3x2+x3,
(1+x)4=1+4x+6x2+4x3+x4,
(1+x)5=1+5x+10x2+10x3+
+5x4+x5,
这些多项式的系数称为二项式。它们有特殊的符号:(1+x)n中xm的系数用Cnm表示。例如,C21=2、C42=6、C53=10。在旧的代数教科书中,数字Cnm被称为“n到m的组合数”。这是有原因的,但我们没有
从这里很容易得到
(a+x)n=Cn0an+Cn1an−1x…
…+Cnnxn
(将公式(1)应用于(1+x)an的计算并将所得等式两边乘以an就足够了)。最后一个公式称为牛顿二项式。因此,系数被称为“二项式”。
显然,Cnm是非负整数,Cnm=0对应m>n(多项式(1+x)n的次数为n,xm和m>n不包含在其中)。很容易进一步验证Cn0=Cn1=1。其他二项式系数 Cnm(其中 0<m<n)可以通过 1+x 的各种幂得到。表中给出了 n⩽10 的值。
我们看到二项式系数增长得相当快。细心的读者会注意到这些数字的排列有多种模式。通常,这种模式很容易从二项式系数的定义中导出。我们在这里仅证明主要的。让我们从最重要的开始——所谓的帕斯卡恒等式。
我们将在其中强调:对我们来说这并不重要。所以,
(1+x)n=Cn0+Cn1x+Cn2x2+…
…+Cnnxn.
(1)
定理 1(帕斯卡恒等式)。
Cnm=Cn−1m+Cn−1m−1
(2)
对于任何自然n,m。
19
证明。根据定义,Cnm是多项式(1+x)n中xm的系数。为了找到这个系数,您需要乘以等于 1+x 的 n 多项式。将它们相乘n−1,我们得到多项式1+Cn−11x+Cn−12x2+⋯+xn−1。执行最后一个乘法,我们得到
(1+x)(1+Cn−11x+Cn−12x2+⋯+xn−1)=1+Cn−11x+Cn−12x2+⋯+xn−1+x+Cn−11x2+Cn−12x3+⋯+xn=1+(Cn−11+1)x+(Cn−12+Cn−11)x2+⋯+xn.
在结果表达式中,xm的系数等于Cn−1m+Cn−1m−1。因此,
Cnm=Cn−1m+Cn−1m−1.
帕斯卡恒等式可以方便地用于计算生物系数。例如,如果我们想在表中再添加一行(第十一行),那么我们只需将上一行(第十行)中相邻的数字成对添加即可:
C110=1C111=C101+C100=11C112=C102+C101=55C113=C103+C102=165C114=C104+C103=330C115=C105+C104=462
这说明将生物项系数写成三角表的形式是很方便的(见
3*
20
因此,公式(3)对于n=1成立,并且从其对于n−1的有效性可知它对于n也成立。这证明了对于任意n的公式(3)。
我们建议初次熟悉公式(3)的读者从中推导出已知的Cn0=Cnn=1和Cnm=0对于m>n的等式。
公式(3)很有趣,因为它右边的分数等于一个整数,即它的分母中的所有数字都会与分子中的数字相消。
将进一步使用以下定理。
定理3. 如果数字n和m互质(即数字n和m的最大公约数为1),则Cnm能被n整除。
证明:
Cnm=1⋅2…mn(n−1)…(n−m+1)=mn⋅1⋅2…(m−1)(n−1)…(n−m+1)=mn⋅Cn−1m
因此,
mCnm=nCn−1m
即数字mCnm可以被n整除。但由于 m 和 n 互质,即 m 不能被 n 的任何素因数整除,因此 Cnm 可以被 n 整除。
例如,C94=126可以被9整除,C103=120可以被10整除。
让我们指出生物系数的更多属性:
- Cnm=Cn−1m
- Cnm+Cn+1m+…+Cm+km=Cm+k+1m
- Cn0+Cn1+…+Cnn=2n
- Cn0+Cn2+Cn4+…=Cn1+Cn3+Cn5+…
我们将这些性质的证明留给读者。
§ 2. 生物系数除以素数的余数
在本段(以及下一段)中,我们经常需要使用“数字 a 和 b 除以 p 时余数相等”的形式。通常,这种形式被缩写符号“a≡b(modp)”代替。换句话说,公式“a≡b (mod p)”意味着a−b可以被p整除。例如,4≡1 (mod 3)、999999≡222222 (mod 7) 等。公式“a≡b (mod p)”有时读作如下:“数字a 与数字b 模p 相当”(但是,我们也不会使用表达式)。
以下是符号 ≡ 的两个特性:
- 如果 a≡b (mod p) 和 k 是整数,则 ka≡kb (mod p)。
事实上,如果a−b能被p整除,那么ka−kb=k(a−b)也能被p整除。
- 如果 a≡b (mod p) 和 b≡c (mod p),则a≡c (mod p)。
事实上,如果a−b能被p整除并且b−c能被p整除,那么a−c=(a−b)+(b−c)也能被p整除。
让我们回想一下,任何自然数a都可以除以自然数p“有余数”,即a可以唯一地写成a=bp+c的形式,其中b,c是整数,使得0≤c<p。
本节的主要目标是证明以下陈述。
定理4.设p为素数,m和n为自然数。此外,令k和l为数字m和n除以p的部分,s和t为余数(即m=kp+s、n=lp+t,其中k、 l、s、t 为整数,0≤s<p、0≤t<p)。然后
Cn∗=Cl∗⋅Cs∗(modp).
正如我们将在下面看到的,该定理使我们能够找到生物系数除以素数的余数,几乎无需执行任何计算。
23
这样C11933≡0 (mod 3),即C11933可以被3整除而没有余数。
请注意,如果我们将定理 4 应用于 Cnm,其中 n≥m,那么通过写出 m=kp+s,n=lp+t,我们当然会得到 l≥k;然而,无法预测s、t 哪个数字更大。如果有s>t,那么根据我们的定理Cnm≡ClkCst=0 (mod p),即Cnm能被p整除。正如我们所看到的,为了确定数字Cnm除以p的余数,通常需要多次使用定理4;每当类似上述的现象发生时,每当这种情况发生时,这意味着类似的数字Cnm可以被p整除。同样,我们得到C11933可以被3整除。
我们看到,n越大,Cnm越有可能被p整除。很容易证明以下更精确的说法:Cnm和0≤n≤pr、0≤m≤n的数字总数还有2pr(pr+1),其中不能被p整除的正好是2rpr(p+1)r(这里p很简单,r很自然;证明仅基于定理3,我们留下给读者)。需要强调的是,对于大r,数字2rpr(p+1)r比数字2pr(pr+1)小许多倍。
例如,在数字 Cnm 和 0≤n≤35、0≤m≤n 中,大约有 26.2% 不能被 3 整除;如果0≤n≤310,则约为3.6%,如果0≤n≤315,则约为0.45%。
最后,我们来说说定理 4 的视觉解释的权重,它可以通过考虑所谓的“根据 p 模型的帕斯卡三角形”来获得。这是从帕斯卡三角形中得到的三角表的名称,如果其中的每个数字被除以 p 时的余数替换。我们不证明关于这个三角形的任何定理,但我们邀请读者看两张描绘帕斯卡三角形模 2 和 3 的图片(封面上和第 22 页)。想想这些三角形在 图 上看不到的部分是如何排列的。尝试将定理 4 表述为关于帕斯卡三角形模 p 结构的定理。
§ 3. 关于生物系数除以素数幂的余数的一些知识
我们不以任何一般方式考虑将生物系数除以合数的余数问题。 (不过,读者可以自己思考一下;例如,数字C11933除以4的余数是多少:一或trem?)我们将仅限于谈论一个我们尚未完全解释的惊人现象。
让我们从一些计算开始。使用生物系数公式(参见§1,公式(2)),我们得到
C21=2,
C42=6,
C84=70,
C168=12870,
C3216=601080390.
(当然,每个人都知道 1、2、4、8、16、32…是 2 的连续幂。)结果数字本身不应该以任何特殊方式使用;然而,它们的连续差异揭示了惊人的特性。让我们看看它们:
6−2=4=22,
70−6=64=26,
12870−70=12800=29⋅25,
24
running_head:
folio:
601 080 390−12 870=601 067 520=212⋅146 745.
我们看到这些差异被分成了二的高次方——大到这种现象不会造成伤害。有可能证明一个至少部分解释它的定理。
定理 5. 如果 n>1,则数
αn=C2n+12n−C2n2n−1
能被22n+2整除。
扫地。
1°。假设 n>1 显着,因为 α1 等于 4 并且不能被 22n+2=24=16 整除。
2°。人们可能认为n>1的αn甚至可以被23n整除:对于n=2,3,4也是如此。但我们不知道如何证明这一点。
证明。让我们从一般观察开始,如果数字 r 是奇数,则 C2nr 可以被 2n 整除。事实上,由于 r 是奇数,并且 2n 除了 2 之外没有素因数,因此数字 r 和 2n 是互质的,我们需要的陈述来自定理 3。
让我们把理论
P(x)=(1+x)2n+1−(1−x2)2n.
多项式P(x)包含x2n,系数为C2n+12n−C2n2n−1=αn(这里我们使用n>1这一事实:当(1−x2)=1+(−x2)出现2n次方时,我们得到系数C2n2n−1不是x2n,而是(-x^2)^{2n-1} = (-1)^{2n-1}x^2^n,而数字(−1)2n−1为1 n>1 且在 n=1 处等于 −1。
另一边,
P(x)=(1+x)2n+1−(1+x)2n=(1+x)2n[(1+x)2n−(1−x)2n].
同时
(1+x)2n−(1−x)2n=(1+x)2n−(1+(−x))2n=1+C2n1x+C2n2x2+C2n3x3+⋯+x2n−1−C2n1(−x)−C2n2(−x)2−C2n3(−x)3−⋯−(−x)2n=
(因为 (−x)k 等于 xk 对于甚至 k 且相等 −xk 对于奇数 k)
=2(C2n1x+C2n3x3+C2n5x5+⋯+C2n2n−1x2n−1).
我们强调,所得多项式 x 仅以奇数次幂出现。
我们想知道 x2n 以什么系数进入 P(x),即进入乘积
(1+x)2n[(1+x)2n−(1−x)2n]=2(1+C2n1x+C2n2x2+C2n3x3+⋯+x2n)×(C2n1x+C2n3x3+C2n5x5+⋯+C2n2n−1x2n−1).
显然,第一个因子乘以x2n−1,第二个因子乘以x,第二个因子乘以x2n−3,第二个因子乘以x3,等等即可得到x2n。因此,x2n的P(x)中的系数等于αn,也等于
2(C2n1C2n2n−1+C2n3C2n2n−3+⋯+C2n2n−1C2n1).
众所周知,每个数字C2n1,C2n3,…,C2n2n−1都可以被2n整除。因此,括号中的每一项都除以2n⋅2n=22n;此外,两者位于整个表达式之前,
25
此外,括号中的每个术语出现两次。因此,根据需要,αn 可以被 22n+2整除。
这样,数αn能被2的高次方整除的奇怪现象就在一定程度上得到了解释。但如果你用三、五、七代替二,就不会观察到类似的情况。真的吗,
C93−C31=84−3=81=34;
C279−C93=4686825−84=4686741=37⋅2143;
C8127−C279=2306279447501851002720−4686825=2306279447501846315895=310⋅39057044954221855;
C255−C51=43130−5=43125=55⋅69;
C497−C71=85900584−7=85900577=75⋅5111.
简而言之,对于素数 p 来说,这个数字似乎是正确的
Cpn+1pn−Cpnpn−1
可以被数字 p 的高次幂整除。但我们不知道如何证明这一点。
顺便说一句,如果 p 不是素数,则不会观察到任何此类情况。例如,
C164−C41=1820−4=1816
甚至不能被42整除;
C366−C61=1947792−6=1947786
甚至不能被62整除。
我们相信«Kvant»的读者之一能够理解生物系数算术这个复杂的问题。