Квант
Основные статьи

Арифметика биномиальных коэффициентов

Фукс Д. Б., Фукс М. Б.

17

АРИФМЕТИКА
БИНОМИАЛЬНЫХ
КОЭФФИЦИЕНТОВ

Д. Б. Фукс, М. Б. Фукс

(1+x)2=1+2x+x2(1+x)^2=1+2x+x^2
(1+x)3=1+3x+3x2+x3(1+x)^3=1+3x+3x^2+x^3

Эти формулы знает любой школьник. Числа 1, 2, 1; 1, 3, 3, 1, а также числа, получающие аналогичным образом при возведении в степени 4, 5, …, называются биомиальные коэффициентами. В этой статье расссказывается о различных свойствах биомиальных коэффициентов. В первом параграфе излагается «общая теория»: большинство из доказываемых здесь теорем входило равше в школьную программу. Во втором параграфе мы показем, как можно совсем легко найти остаток от деления любого из биомиальных коэффициентов на любое простое число. Некоторые утиверждения этого параграфа сформулированы в виде гипотез. Можно быть, читатели «Кванта» попробуют их доказать.

18

§ 1. Определение и простейшие свойства биномiales коэффициентов

Если двучлен (бином) 1+x1+x возвеси в некоторую натуральную степень nn, то, очевидно, получится многочлен nn-й степени (т. е. наибольшая степень, в котором xx войдет в этот многочлен, будет равна nn). Например:

(1+x)0=1,(1+x)^0 = 1,
(1+x)1=1+x,(1+x)^1 = 1+x,
(1+x)2=1+2x+x2,(1+x)^2 = 1+2x+x^2,
(1+x)3=1+3x+3x2+x3,(1+x)^3 = 1+3x+3x^2+x^3,
(1+x)4=1+4x+6x2+4x3+x4,(1+x)^4 = 1+4x+6x^2+4x^3+x^4,
(1+x)5=1+5x+10x2+10x3+(1+x)^5 = 1+5x+10x^2+10x^3+
+5x4+x5,+5x^4+x^5,

Коэффициенты этих многочленов называются биномиальным. Для них есть специальные обозначения: коэффициент при xmx^m в (1+x)n(1+x)^n обозначает через CnmC_n^m. Например, C21=2C_2^1=2, C42=6C_4^2=6, C53=10C_5^3=10. В старых учебниках по алгебре число CnmC_n^m называется «числом сочетаний из nn по mm». На то есть свои причины, но мы не

Отсюда легко получить, что и

(a+x)n=Cn0an+Cn1an−1x…(a+x)^n = C_n^0a^n + C_n^1a^{n-1}x\ldots
…+Cnnxn\ldots + C_n^nx^n

(достаточно применить формулу (1) к вычислению (1+x)an(1+x)a^n и умножить обе части полученного равенства на ana^n). Последняя формула называется биномом Ньютона. Отсюда и называние коэффициентов — биномиальные.

Ясно, что CnmC_n^m — целые неотрицательные числа и что Cnm=0C_n^m=0 при m>nm>n (многочлен (1+x)n(1+x)^n имеет степень nn, и xmx^m с m>nm>n в него не входит). Легко убедиться, далее, в том, что Cn0=Cn1=1C_n^0=C_n^1=1. Другие биномиальные коэффициенты CnmC_n^m, где 0<m<n0<m<n, можно находить, возвода 1+x1+x в различных степени. Их значения для n⩽10n\leqslant10 приведены в таблице.

Мы видим, что биномиальные коэффициенты — довольно-tаки быстро растущее числа. Наблюдательный читатель заметит в расположении этих чисел ряд закономерностей. Как правило, такие закономерности легко выводятся из определения биномиальных коэффициентов. Мы доказем здесь только главные из них. Начнем с главного из главных — с так называемого тожества Паскала.

Будем выделять в них: для нас это не важно. Итак,

(1+x)n=Cn0+Cn1x+Cn2x2+…(1+x)^n = C_n^0+C_n^1x+C_n^2x^2+\ldots
…+Cnnxn.\ldots +C_n^nx^n.

(1)

Теорема 1 (тоождество Паскала).

Cnm=Cn−1m+Cn−1m−1C_n^m = C_{n-1}^m+C_{n-1}^{m-1}

(2)

для любых натуральных n,mn, m.

19

Доказательство. По определению, CnmC_n^m есть коэффициент при xmx^m в многочлене (1+x)n(1+x)^n. Для того чтобы найти этот коэффициент, нужно перемножить nn многочленов, равных 1+x1+x. Перемножив n−1n-1 из них, получим многочлен 1+Cn−11x+Cn−12x2+⋯+xn−11+C_{n-1}^1x + C_{n-1}^2x^2 + \dots + x^{n-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.(1+x)(1+C_{n-1}^1x + C_{n-1}^2x^2 + \dots + x^{n-1}) = 1 + C_{n-1}^1x + C_{n-1}^2x^2 + \dots + x^{n-1} + x + C_{n-1}^1x^2 + C_{n-1}^2x^3 + \dots + x^n = 1 + (C_{n-1}^1 + 1)x + (C_{n-1}^2 + C_{n-1}^1)x^2 + \dots + x^n.

В полученном выражении коэффициент при xmx^m равен Cn−1m+Cn−1m−1C_{n-1}^m + C_{n-1}^{m-1}. Таким образом,

Cnm=Cn−1m+Cn−1m−1.C_n^m = C_{n-1}^m + C_{n-1}^{m-1}.

Тождество Паскала удобно использовать для вычисления биомиальные коэффициентов. Например, если мы захотим дополить нашу таблицу еще одной строккой (одиннадцатой), то нам достаточно сложить попарно числа, стоящие рядом друг с другом в предыдущей (десятой) строке:

C110=1C111=C101+C100=11C112=C102+C101=55C113=C103+C102=165C114=C104+C103=330C115=C105+C104=462C_{11}^0 = 1 \quad C_{11}^1 = C_{10}^1 + C_{10}^0 = 11 \quad C_{11}^2 = C_{10}^2 + C_{10}^1 = 55 \quad C_{11}^3 = C_{10}^3 + C_{10}^2 = 165 \quad C_{11}^4 = C_{10}^4 + C_{10}^3 = 330 \quad C_{11}^5 = C_{10}^5 + C_{10}^4 = 462

Отсюда видно, что биомиальные коэффициенты удобно записать в виде треугольной таблицы (см. застав-

3*

20

Таким образом, формула (3) имеет место при n=1n=1, и из ее справедливости для n−1n-1 вытекает, что она верна для nn. Этим формула (3) доказана при любом nn.

Мы рекомендуем читателю, впервые познакомившись с формулой (3), вывести из нее уже известные нам равенства Cn0=Cnn=1C_n^0 = C_n^n = 1 и Cnm=0C_n^m = 0 при m>nm>n.

Формула (3) интересна уже тем, что дробь, стоящая в ее правой части, равняется целому числу, т. е. все числа, стоящие в ее знаменателе, сократятся с числами, стоящими в числитель.

Следующая теорема будет использована в дальнейшем.

Теорема 3. Если числа nn и mm взаимно просты (т. е. наибольший общий делитель чисел nn и mm равен 1), то CnmC_n^m делится на nn.

Доказательство:

Cnm=n(n−1)…(n−m+1)1⋅2…m=nm⋅(n−1)…(n−m+1)1⋅2…(m−1)=nm⋅Cn−1mC_n^m = \frac{n(n-1)\ldots(n-m+1)}{1\cdot2\ldots m} = \frac{n}{m} \cdot \frac{(n-1)\ldots(n-m+1)}{1\cdot2\ldots(m-1)} = \frac{n}{m} \cdot C_{n-1}^m

Таким образом,

mCnm=nCn−1mmC_n^m = nC_{n-1}^m

т. е. число mCnmmC_n^m делится на nn. Но так как mm и nn взаимно просты, т. е. mm не делится ни на один простой делитель числа nn, то CnmC_n^m делится на nn.

Например, C94=126C_9^4 = 126 делится на 9, C103=120C_{10}^3 = 120 делится на 10.

Укажем еще несколько свойств биомиальных коэффициентов:

  1. Cnm=Cn−1mC_n^m = C_{n-1}^m
  2. Cnm+Cn+1m+…+Cm+km=Cm+k+1mC_n^m + C_{n+1}^m + \ldots + C_{m+k}^m = C_{m+k+1}^m
  3. Cn0+Cn1+…+Cnn=2nC_n^0 + C_n^1 + \ldots + C_n^n = 2^n
  4. Cn0+Cn2+Cn4+…=Cn1+Cn3+Cn5+…C_n^0 + C_n^2 + C_n^4 + \ldots = C_n^1 + C_n^3 + C_n^5 + \ldots

Доказательства этих свойств мы оставляем читателю.

§ 2. Остатки от деления биомиальных коэффициентов на простые числа

В этом (да и в следующем) параграфе нам придется часто употреблять форму «числа aa и bb имеют равные остатки от деления на pp». Общий эту форму заменяет сокращенной записью «a≡ba \equiv b (mod pp)». Иначе говоря, формула «a≡ba \equiv b (mod pp)» означает, что a−ba-b делится на pp. Например, 4≡14 \equiv 1 (mod 3), 999999≡222222999999 \equiv 222222 (mod 7) и т. д. Формула «a≡ba \equiv b (mod pp)» читается иногда так: «число aa сравним с числом bb по модulo pp» (впрочем, мы также выражения употребляются не будем).

Вот два очень свойства символа ≡:

  1. Если a≡ba \equiv b (mod pp) и kk — целое число, то ka≡kbka \equiv kb (mod pp).

В самом деле, если a−ba-b делится на pp, то и ka−kb=k(a−b)ka-kb = k(a-b) делится на pp.

  1. Если a≡ba \equiv b (mod pp) и b≡cb \equiv c (mod pp), то a≡ca \equiv c (mod pp).

В самом деле, если a−ba-b делится на pp и b−cb-c делится на pp, то и a−c=(a−b)+(b−c)a-c = (a-b)+(b-c) делится на pp.

Напомним, что всякое натуральное число aa можно поделить на натуральное число pp «с остатком», т. е. aa можно единственным образом записать в виде a=bp+ca=bp+c, где bb, cc — такие целые числа, что 0≤c<p0 \leq c < p.

Главной целью этого параграфа является доказательство следующего утверждения.

Теорема 4. Пусть pp — простое число, mm и nn — натуральные числа. Пусть, далее, kk и ll — частиные от деления чисел mm и nn на pp, а ss и tt — остатки (т. е. m=kp+sm=kp+s, n=lp+tn=lp+t, где kk, ll, ss, tt — целые числа и 0≤s<p0 \leq s < p, 0≤t<p0 \leq t < p). Тогда

Cn∗=Cl∗⋅Cs∗(modp).C_n^* = C_l^* \cdot C_s^* \pmod p.

Как мы увидим ниже, эта теорема позволяет находить остатки от деления биомиальных коэффициентов на простые числа, почти не производа вычислений.

23

таким образом, что C11933≡0C_{119}^{33} \equiv 0 (mod 3), т. е. C11933C_{119}^{33} делится без остатка на 3.

Отметим, что если мы применяем нашу теорему 4 к CnmC_n^m, где n≥mn \geq m, то, записьвая m=kp+sm = kp + s, n=lp+tn = lp + t, мы, конечно, будем иметь l≥kl \geq k; однако невозможно предвидеть, какое из чисел ss, tt оказывается большим. Если будет s>ts > t, то согласно нашей теореме Cnm≡ClkCst=0C_n^m \equiv C_l^k C_s^t = 0 (mod pp), т. е. CnmC_n^m делится на pp. Как мы видели, для определения остатка от деления числа CnmC_n^m на pp приходит, как правило, использовать теорему 4 несколько раз; и каждый раз может возникнуть явление, подобное описанному выше, причем, когда бы это ни случилось, это означает, что и схо-дное число CnmC_n^m делится на pp. Подобным образом мы получили, что C11933C_{119}^{33} делится на 3.

Мы видим, что чем больше nn, тем более верятна делимость числа CnmC_n^m на pp. Легко доказать следующее, более точное утверждение: всего чисел CnmC_n^m с 0≤n≤pr0 \leq n \leq p^r, 0≤m≤n0 \leq m \leq n имеется pr(pr+1)2\frac{p^r (p^r + 1)}{2}, из них не делится на ppровно pr(p+1)r2r\frac{p^r (p + 1)^r}{2^r} (здесь pp просто, rr натуральное; доказательство опирается только на теорему 3, мы оставляем его читателю). Подчерком, что при больших rr числа pr(p+1)r2r\frac{p^r (p + 1)^r}{2^r} во много раз меньше числа pr(pr+1)2\frac{p^r (p^r + 1)}{2}.

Так, например, из чисел CnmC_n^m с 0≤n≤350 \leq n \leq 3^5, 0≤m≤n0 \leq m \leq n не делятся на 3 примерно 26,2%; если 0≤n≤3100 \leq n \leq 3^{10}, то примерно 3,6%, а если 0≤n≤3150 \leq n \leq 3^{15}, то примерно 0,45%.

В заключение скажем несколько слов о весыта наглядной интерпретации теоремы 4, которое можно получить, рассматривая так называемый «треугольник Паскала по модели pp». Так называется треугольная таблица, которая получится из треугольника Паскала, если в нем заменить каждое число его остатком от деления на pp. Мы не доказываем про этот треугольник никаких теорем, но предлагаем читателю посмотреть две картины (на обложке и на стр. 22), на которых изображены треугольники Паскала по модулям 2 и 3. Подумайте, как устроены эти треугольники в той части, которая не видна на рисунке. Постраайтесь формулировать теорему 4 таким образом, чтобы она стала теоремой об устройстве треугольника Паскала по модуляю pp.

§ 3. Немного об остатках от деления биомиальных коэффициентов на степени простых чисел

Мы не рассматриваем в сколькоибудь общей постановке вопрос об остатках от деления биомиальных коэффициентов на составные числа. (Впрочем, читатель может подумать об этом сам; чему, например, равен остаток от деления числа C11933C_{119}^{33} на 4: единице или trem?) Мы ограничимся тем, что расскажем об одном удивительном, до конца не обьясенном нами явлении.

Начнем с некоторых вычислений. Пользовясь формулой биомиальных коэффициентов (см. § 1, формула (2)), получаем

C21=2,C_2^1 = 2,
C42=6,C_4^2 = 6,
C84=70,C_8^4 = 70,
C168=12870,C_{16}^8 = 12870,
C3216=601080390.C_{32}^16 = 601080390.

(Все знают, конечно, что 1, 2, 4, 8, 16, 32, … — это последовательные степени двойки.) Сами полученные числа исчем особенным не применять; их последовательные разности, однако, обнаруживают порAZительные свойства. Посмотрим на них:

6−2=4=22,6 - 2 = 4 = 2^2,
70−6=64=26,70 - 6 = 64 = 2^6,
12870−70=12800=29⋅25,12870 - 70 = 12800 = 2^9 \cdot 25,