이항 계수는 수학의 많은 영역, 특히 조합론에서 발생합니다.기호 k {는 보통 "nchoose k"로 읽힙니다.이는 고정된n 요소 집합에서 k 요소의 (순서 없는) 서브셋을 선택하는 이(k 있기 때문입니다.를 , {1, 2, 3, {{,3 {1, {에서 2개의 요소를 선택하는 방법은 (입니다
이항 계수는 임의의 복소수 z와 정수 k ≤ 0에 대해 (k) { {로일반화할 수 있으며, 이들의 성질은 이 보다 일반적인 형태로 계속 유지됩니다.
안드레아스 폰 에팅스하우젠은 1826년에 [1]표기법(을 도입했다. 비록 그 숫자들이 수세기 전에 알려져 있었지만 말이다(파스칼의 삼각형 참조).이항계수에 대한 가장 초기의 상세한 논의는 고대 산스크리트어 문헌인 핑갈라의 찬다샤스트라에 대한 할라유다의 10세기 주석에서 볼 수 있다.이항 계수에 대한 두 번째 초기 설명은 Al-Karaji에 의해 제공됩니다.약 1150년, 인도의 수학자 바스카라차랴는 그의 책 Lālavatī에서 이항 계수를 설명했습니다.[2]
대체 표기법으로는 C(n,k),Ck, Ck, Ckn, Cnk, C, C, C가n,k있습니다.C는 모두 조합 또는 선택지를 나타냅니다.대부분의 계산기는 C 표기법을 한 줄 표시로 나타낼 수 있기 때문에 C 표기법의 변형을 사용합니다.이 형식에서 이항 계수는 P(n,k)로 표기된 n의 k-변환과 쉽게 비교된다.
정의와 해석
k
n
0
1
2
3
4
⋯
0
1
0
0
0
0
⋯
1
1
1
0
0
0
⋯
2
1
2
1
0
0
⋯
3
1
3
3
1
0
⋯
4
1
4
6
4
1
⋯
⋮
⋮
⋮
⋮
⋮
⋮
⋱
처음 몇 개의 이항 계수 파스칼의 왼쪽 삼각형에서
자연수(0 포함) n및 k의 경우, 이항 계수)(\ {는 (1+ X)n의 확장 시 단항k X의 계수로 정의할 수 있다.이항식에서도 동일한 계수가 발생합니다(kn n인 경우).
(표준)
(가환환의 모든 요소 x, y에 대해 유효)를 지정하면, 「변환 계수」라고 하는 이름이 됩니다.
이 숫자의 또 다른 발생은 조합론에서 k개의 객체가 n개의 객체 중에서 선택될 수 있는 방법(또는 n개의 요소 집합의 k개의 조합)을 제공하는 것이다.이 숫자를 쓴 첫번째 정의의 동등한 상대로, 독립적으로 어떤 공식을 아래의 계산을 수행해야 만약 각 전력의 n요소 의 한(1+X)n 일시적으로 나는(1n까지 달리기) 다음 k지수의 각 부분 집합의 확장 기여 Xk게 해 주는 지수, m이상인 계수를 'X이라고 볼 수 있o결과의 nomial은 그러한 서브셋의 수가 됩니다.이는 특히 ( k { {이가) 자연수 n과 k에 대한 자연수임을 나타냅니다.이항계수(이항계수식에 의해 답이 주어지는 문제)에는 다른 많은 조합적 해석이 있다. 예를 들어, 합계가 k인 n비트(자릿수 0 또는 1)로 구성된 단어의 수는 ( kk에 의해주어지며, 반면 wr에 대한 방법은 다음과 같다. k + 2 + + _ { } + _ {2 + \ +_ { } (여기서i 모든 a는 음이 아닌 정수)는 ( - n -\ { + { )로쉽게 알 수 있습니다.
공식은 집합 {1, 2, 3, ...,n}을 고려하고 (a) 모든 그룹에서 특정 집합 요소(예를 들어 "i")를 포함하는 k-요소 그룹을 개별적으로 계수하는 것으로부터 뒤따른다 (i"는 모든 그룹에서 한 자리를 채우도록 이미 선택되었기 때문에, 우리는 "i"를 포함하지 않는 나머지n - 1에서 k-group만을 선택하면 된다).는 n개의 요소의 가능한 모든 k 조합을 열거한 것입니다.또한 (1 + X)(n−11 + X)의 X에 대한k 기여도를 추적합니다.(1 + X)n에는 X 또는−1 X가 0이므로n+1 k > n또는k < 0일 때 ( k (\}=으로 정의를 확장할 수 있습니다.이 재귀적 공식은 파스칼의 삼각형을 만들 수 있으며, 이 삼각형은 0, 즉 사소한 계수가 있을 수 있는 공백으로 둘러싸여 있습니다.
곱셈 공식
개별 이항 계수를 계산하는 보다 효율적인 방법은 다음 공식에 의해 제공됩니다.
여기서 첫 번째 k {\ n {k의 분자는하강 요인 검정력으로 표현된다.이 공식은 이항 계수의 조합적 해석을 위해 가장 이해하기 쉽습니다.분자는 n개의 개체 집합에서 k개의 개별 개체 순서를 선택하는 방법을 제공합니다.분모는 순서가 무시될 때 동일한 k-결합을 정의하는 구별되는 시퀀스의 수를 카운트합니다.
k와 n- k에 대한 이항계수의 대칭성으로 인해 위의 제품의 상한을 k와 n - k 중 작은 값으로 설정하여 계산을 최적화할 수 있다.
요인 공식
마지막으로, 계산상 부적절하지만, 흔히 증명과 도출에 사용되는 콤팩트 형식이 있어 익숙한 요인 함수를 반복적으로 사용한다.
여기서 n!은 n의 계수를 나타냅니다.이 공식은 분자와 분모에 (n- k!)를 곱한 위의 곱셈 공식에서 따랐습니다. 그 결과 분자와 분모에 공통되는 많은 인자를 포함합니다.공통 요인이 먼저 취소되지 않는 한(특히 요인 값이 매우 빠르게 증가하므로) 명시적 계산(k가 작고 n이 큰 경우)에는 덜 실용적입니다.공식은 곱셈 공식에서 덜 명확한 대칭을 나타낸다(정의에서이다).
곱셈 공식은 n을 임의의 수α(음수, 실수, 복소수)로 대체하거나 모든 양의 정수가 반전되는 모든 교환환의 요소로 대체함으로써 이항 계수의 정의를 확장할[3] 수 있도록 한다.
이 정의에서는 (변수 중 하나를 1로 설정한 상태에서) 이항식을 일반화하면 ( k) \이항 계수를 호출할 수 있습니다.
(2)
이 공식은 X < 1의 모든 복소수α와 X에 유효하다.그것은 또한 X에서 형식멱급수의 동일성으로 해석될 수 있으며, 여기서 그것은 실제로 상수 계수가 1인 임의의 멱급수의 정의로 작용할 수 있다; 요점은 이 정의로 모든 멱급수가 지수에 대해 기대하는 것을 유지한다는 것이다, 특히
α가 음이 아닌 정수n이면 k > n인 항은 모두 0이고 무한 급수는 유한합이 되어 이항식을 회복한다.그러나 음의 정수와 유리수를 포함한 α의 다른 값들의 경우, 급수는 정말로 무한하다.
파스칼 삼각형
파스칼 삼각형의 1000번째 행은 수직으로 배열되며 계수의 십진수 자릿수를 회색 스케일로 표현하여 오른쪽 정렬됩니다.이미지의 왼쪽 경계는 대략 이항 계수의 로그 그래프에 해당하며 로그 오목형 시퀀스를 형성함을 나타냅니다.
행 번호 n에는 k = 0, …,n에 대한 k { {k})가 포함됩니다.가장 바깥쪽 위치에 1s를 먼저 배치한 다음 각 내부 위치를 바로 위의 두 숫자의 합으로 채우는 방식으로 구성됩니다.이 방법을 사용하면 분수나 곱셈 없이 이항 계수를 빠르게 계산할 수 있습니다.예를 들어, 삼각형의 5번째 행을 보면, 그것을 빠르게 읽을 수 있다.
조합과 통계
이항 계수는 특정 빈번한 계수 문제에 대해 준비된 공식을 제공하기 때문에 조합론에서 중요합니다.
n개의 요소 집합에서 k개의 요소를 선택하는 방법은( {가지가 있습니다.「조합」을 참조해 주세요.
반복이 허용되는 경우 n개의 요소 집합에서 k개의 요소를 선택하는 방법은( + - k) { {n+ { ) 있습니다.멀티셋을 참조해 주세요.
각 다항식k k})은 정수 값입니다. 모든 정수 t{ t에 정수 값이 있습니다. (이를 증명하는 한 가지 방법은 파스칼의 항등식을 사용하여 k에 유도하는 것입니다.)따라서 이항 계수 다항식의 모든 정수 선형 조합도 정수 값입니다.반대로 (4)는 모든 정수값 다항식이 이러한 이항 계수 다항식의 정수 선형 조합임을 보여준다.보다 일반적으로 특성 0 필드 K의 하위 링 R에 대해, K[t]의 다항식은 이항 계수 다항식의 R-선형 조합인 경우에만 모든 정수의 값을 취한다.
예
정수값 다항식 3t(3t + 1) / 2는 다음과 같이 다시 쓸 수 있습니다.
이항 계수와 관련된 동일성
요인 공식을 사용하면 인근의 이항 계수를 쉽게 연결할 수 있습니다.예를 들어 k가 양의 정수이고 n이 임의의 경우
(5)
그리고 조금만 더 하면
우리는 또한
또한 다음 사항이 유용할 수 있습니다.
상수 n에 대해 다음과 같은 반복이 있습니다.
요약하자면, 우리는
이항 계수의 합
공식
( ∗ ( )
파스칼 삼각형의 n번째 행에 있는 원소는 항상 2의 n제곱을 더한다고 합니다.이것은 x = 1과 y = 1을 설정하여 이항 정리(정리)에서 구한다.이 공식은 또한 자연 조합 해석을 가지고 있다: 왼쪽은 k = 0, 1, …, n 크기의 {1, …, n}의 서브셋의 수를 합하여 총 서브셋 수를 나타낸다.(즉, 왼쪽은 {1, …, n}의 전원 세트를 카운트합니다).단, 이들 서브셋은 각 요소 1, …, n을 순차적으로 선택하거나 제외하여 생성할 수도 있습니다.n개의 독립된 바이너리 선택지(비트 문자열)는 총 의\ 2개의를 허용합니다.왼쪽과 오른쪽은 동일한 부분 집합 집합을 세는 두 가지 방법이기 때문에 같습니다.
(1 + x)(m1 + n−mx) = (1 + nx)의 팽창에서 x({ x의 계수를 식 (2)를 이용하여 조사하면 알 수 있다.m = 1이면 방정식(7)은 방정식(3)으로 줄어든다.특수한 경우 n = 2m,k = m에서 (1)을 사용하면 팽창(7)은 (오른쪽 파스칼 삼각형에서 볼 수 있음)이 된다.
파스칼의 삼각형, 0행에서 7행까지. m = 3에 대한방정식 8은 3행과 6행에 2 + 2 + + 으로 나타나 {\ 1+ + }
여기서 m과 d는 복소수입니다.그 후 (10)을xx) : + x에즉시 하고Qx) {displaystyle PP(m)} {displaystyleQ(x)} {displaystyle Q(의 차수가 n 이하임을 관찰한다.
k - k j 1 ( + ) ( x- ) { { {{k - \ _ { j=}^{ \ } { \ 1 { \ + } { }{1 {은(는) k÷ 2에 대해 수렴됩니다.이 공식은 독일 탱크 문제의 분석에 사용됩니다. 값은 k- j M ( + k ) ( -) - ( + k - ) ( { - 1 )( { k} \ { j= 0} { \ {} { \ + x + x x + x} { k } { k } { k } { k } { k} { k } ) k 。에 유도하여 증명하는frac k-1{Mk-1}}}.
조합 증명과 동일성
이항 계수와 관련된 많은 동일성은 조합 평균에 의해 입증될 수 있다.예를 들어 음이 아닌 n{의 경우 ID는 다음과 같습니다.
(q = 1)은 다음과 같이 이중 계수 증명을 제공할 수 있다.왼쪽은 최소 q개의 요소가 있는 [n] = {1, 2, …, n}의 하위 집합을 선택하고 선택한 요소 중에서 q개의 요소를 표시하는 방법을 카운트합니다.오른쪽도 마찬가지입니다.마킹할 q 요소 세트를 선택하는 방법은 (q) { {n} [n]의 나머지 요소 중 어느 것이 서브셋에 속하는지 선택하는 방법은 - { 2개입니다.
파스칼의 정체성 속에서
양쪽 모두 [n]의 k-소수 부분 집합의 수를 계산합니다. 오른쪽에 있는 두 개의 항은 요소를 포함하는 항과 포함하지 않는 항으로 그룹화합니다.
빈 정사각형이 스타일 씩 일렬로 되어 있고 그 중 n개를 표시(선택)한다고 가정합니다.( ( ( 방법이 있습니다한편, 나머지 n개의 정사각형에서 첫 번째 n과 -k(\ 정사각형 중에서 k개의 정사각형을 선택하여 n개의 정사각형을 선택할 수 있습니다.0 ~ n의 임의의 k가 동작합니다.이것으로 알 수 있다.
(1)을 신청하면 결과가 나옵니다.
F(0) = F(1) = 1이 되도록 지수화된 피보나치 숫자의 순서를 F(i)로 나타내면, 동일성
에는 다음과 같은 조합 [10]증거가 있습니다.유도를 통해 F(n)가 n × 1 정사각형 조각이 2 × 1및 1 × 1 타일로 덮일 수 있는 방법을 계산한다는 것을 알 수 있다.한편, 이러한 타일이 2 × 1 타일 중 정확히 k개를 사용하는 경우, 1 × 1 타일중 n - 2k를 사용하므로n - k 타일 합계가 사용됩니다.이러한 타일을 주문하는 방법은( ){ {가지가 있습니다.따라서 k의 가능한 모든 값에 대해 이 계수를 합하면 동일성이 부여됩니다.
1852년 Kummer는 m과 n이 음이 아닌 정수이고 p가 소수인 경우 p 나누기 m의 최대 제곱(\은 p가 된다는c 것을 증명했습니다.여기서 c는 m과 n을 기저 p에 더했을 때의 반송파 수입니다.마찬가지로 ( k의 소수 p in ({의 지수는 k/p의j소수 부분이 n/p의j소수 부분보다 크도록 음이 아닌 정수 j의 수와 같다.이를 통해 (k) { {은(는) n/gcd(n,k)로 나눌 수 있음을추론할 수 있습니다.따라서 p는 모든r 양의 정수 r 및 s에 대해 (s {{}}{를나눕니다.단, p의 고승수에는 해당되지 않습니다.예를 들어 는되지 않습니다tbinom {9
David Singmaster(1974)의 다소 놀라운 결과는 모든 정수가 거의 모든 이항 계수를 나눈다는 것입니다.보다 정확하게는 정수 d를 고정하고 f(N)는 n < N으로 이항계수k의 수(\displaystyle\tbinomn}{k})를 나타내며 d는 n< k) \ {k
n < N인 이항계수k) {\ {k의 수는 N(N + 1) / 2이므로 d로 나눌 수 있는 이항계수의 밀도가 1이 된다는 것을 의미한다.
이항 계수에는 연속 정수의 최소 공배수와 관련된 약수 특성이 있습니다.예를 [12]들어 다음과 같습니다.
( +)({ {\ {+1, +}) n {을를)분할합니다.
+) { {+ } { k는 lcm ( , n + , , + ) k ) ,( )、 ...,k )\ { lcm } ( n , n , n + \ )의 배수입니다
또 다른 사실:모든 중간 이항 계수가 다음과 같은 경우에만 정수 n ≤ 2가 소수이다.
는 n으로 나눌 수 있습니다.
증명: p가 소수일 경우 p는 나눗셈
(k ) p ( - ) ( - + 1)( k- 1 ) {\ ( \ { \ { ={ p \ cdot ( p - ) \ ( - k + 1) } { \cdots}
){ {은 자연수이고 p는 분자를 나누지만 분모는 나누지 않기 때문입니다.n이 합성인 경우 p를 n의 최소 소수 인자로 하고 k = n/p로 한다.다음으로 0 < p < n 및
그렇지 않으면 분자 k(n- 1)(n- 2)θ(n- p + 1)는 n = k×p로 나누어져야 하며, 이는 (n- 1)θ(n- p + 1)가 p로 나누어질 때만 해당될 수 있다.그러나 n은 p로 나눌 수 있으므로 p는 n - 1, n - 2, …, n -p + 1을 나누지 않습니다. p는 소수이기 때문에 p는 (n - 1)(n - 2)δ(n- p+ 1)를 나누지 않으므로 분자는 n으로 나눌 수 없습니다.
p ( _를 양의 정수 n의 p-adic 가치로 하고, n을 나누는소수 p의 최대 거듭제곱의 지수라고 하자.그럼.
{ \ { x \} denotes、 x { x의 부분을 나타냅니다.
한계 및 점근 공식
k)의 다음 경계는 1µkµn이 되도록 n과 의 모든 값에 대해
첫 번째 불평등은 ...라는 사실에서 비롯된다.
그리고 이 제품의은 두 번째 을 나타내려면 유사한 인수를 사용할 수 있다최종 엄밀한 부등식은 e > / e > 에 해당하며, 이는 RHS가 지수 e j 0 k /j! { e {의 항이므로
여기서 연결 계수는 다항 계수입니다.라벨 부착 조합 객체의 관점에서 연결 계수는 첫 번째 k개의 라벨이 식별되거나 무게 m + n - k의 새로운 라벨 부착 조합 객체를 얻기 위해 함께 접착된 체중 m과 n의 라벨 부착 조합 객체 쌍에 m +n - k 라벨을 할당하는 방법의 수를 나타낸다(즉, sepa).라벨을 접착된 부분, 첫 번째 개체의 접착되지 않은 부분 및 두 번째 개체의 접착되지 않은 부분에 적용할 세 부분으로 등급을 매깁니다.)이와 관련하여 이항 계수는 지수 생성 급수에 대한 계수이며 하강 인자는 일반 생성 급수에 대한 계수입니다.
이항 계수는 지정된 집합에서 지정된 크기의 부분 집합을 카운트합니다.관련된 조합 문제는 소정의 크기의 멀티셋을 소정의 집합에서 추출한 요소로 카운트하는 것이다.즉, 같은 요소를 반복적으로 선택할 수 있는 소정의 집합에서 특정 수의 요소를 선택하는 방법의 수를 카운트하는 것이다.결과값을 [17]멀티세트계수라고 부릅니다.n개의 요소 집합에서 k개의 항목을 "멀티초즈"(즉, 치환으로 선택)하는 방법은(k)\로 표시됩니다
이 기사에서n의 주요 표현과의 모호성 및 혼동을 피하기 위해 f = n= r+ (k- 1) 및 r= f- (k- 1)로 합니다.
다중 집합 계수는 규칙에 따라 이항 계수로 표시할 수 있습니다.
One possible alternative characterization of this identity is as follows: We may define the falling factorial as
and the corresponding rising factorial as
so, for example,
Then the binomial coefficients may be written as
while the corresponding multiset coefficient is defined by replacing the falling with the rising factorial:
음의 정수에 대한 일반화 n
음수 및 분수 n에 대해 확장된 이항 계수C(n,k)이며, 단순 이항으로 표시됩니다.파스칼의 삼각형이 회전하고 다른 항이 부정되는 것을 관찰할 수 있다.케이스 = -1은 Grandi의 시리즈를 나타냅니다.
임의의 n에 대해서
특히, 음의 정수 n에서 평가된 이항 계수는 부호화된 다중 집합 계수에 의해 주어진다.특수한 n - { n=-의 경우, 이 값은(-(- ) (-) .{\(-로 감소합니다
예를 들어, n = -4이고 k = 7이면 r = 4이고 f = 10입니다.
2개의 실제 또는 복잡한 가치 있는 인수
이항 계수는 감마 함수또는 베타 함수를 사용하여 두 개의 실제 또는 복잡한 값 인수로 일반화된다.
이 정의는 \ \ 에서 다음 추가 속성을 상속합니다.
게다가.
그 결과 나온 함수는 거의 연구되지 않았으며, 최초의 그래프는 (Fowler 1996)에 나타나 있다. 많은 이항식별이 실패한다 ( m ) ( - ) ( m )= ( - ( - n- m ) -n -m)\ - } { - { - so\ ( )이 동작은 매우 복잡하고 다양한 8진수(즉, x축과 y축과 y)에 대해)에서 현저하게 다릅니다. 음의 x에 대한 동작은 음의 정수 값에 특이점을 가지며 양의 영역과 음의 영역의 체커보드를 가집니다.
in the octant it is a smoothly interpolated form of the usual binomial, with a ridge ("Pascal's ridge").
in the octant and in the quadrant the function is close to zero.
in the quadrant the function is alternatingly very large positive and negative on the parallelograms with vertices
in the octant the behavior is again alternatingly very large positive and negative, but on a square grid.
in the octant it is close to zero, except for near the singularities.
The definition of the binomial coefficient can be generalized to infinite cardinals by defining:
where A is some set with cardinality. One can show that the generalized binomial coefficient is well-defined, in the sense that no matter what set we choose to represent the cardinal number , will remain the same. For finite cardinals, this definition coincides with the standard definition of the binomial coefficient.
Assuming the Axiom of Choice, one can show that for any infinite cardinal .
are very slow and are useless for calculating factorials of very high numbers (in languages such as C or Java they suffer from overflow errors because of this reason). A direct implementation of the multiplicative formula works well:
defbinomial_coefficient(n:int,k:int)->int:ifk<0ork>n:return0ifk==0ork==n:return1k=min(k,n-k)# Take advantage of symmetryc=1foriinrange(k):c=c*(n-i)/(i+1)returnc
(In Python, range(k) produces a list from 0 to k−1.)
Pascal's rule provides a recursive definition which can also be implemented in Python, although it is less efficient:
defbinomial_coefficient(n:int,k:int)->int:ifk<0ork>n:return0ifk>n-k:# Take advantage of symmetryk=n-kifk==0orn<=1:return1returnbinomial_coefficient(n-1,k)+binomial_coefficient(n-1,k-1)
The example mentioned above can be also written in functional style. The following Scheme example uses the recursive definition
Rational arithmetic can be easily avoided using integer division
The following implementation uses all these ideas
(define (binomialnk);; Helper function to compute C(n,k) via forward recursion(define (binomial-iternkiprev)(if (>= ik)prev(binomial-iternk(+ i1)(/ (* (- ni)prev)(+ i1)))));; Use symmetry property C(n,k)=C(n, n-k)(if (< k(- nk))(binomial-iternk01)(binomial-itern(- nk)01)))
When computing in a language with fixed-length integers, the multiplication by may overflow even when the result would fit. The overflow can be avoided by dividing first and fixing the result using the remainder:
Implementation in the C language:
#include<limits.h>unsignedlongbinomial(unsignedlongn,unsignedlongk){unsignedlongc=1,i;if(k>n-k)// take advantage of symmetryk=n-k;for(i=1;i<=k;i++,n--){if(c/i>=ULONG_MAX/n)// return 0 on potential overflowreturn0;c=c/i*n+c%i*n/i;// split c * n / i into (c / i * i + c % i) * n / i}returnc;}
Another way to compute the binomial coefficient when using large numbers is to recognize that
where denotes the natural logarithm of the gamma function at . It is a special function that is easily computed and is standard in some programming languages such as using log_gamma in Maxima, LogGamma in Mathematica, gammaln in MATLAB and Python's SciPy module, lngamma in PARI/GP or lgamma in C, R,[19] and Julia. Roundoff error may cause the returned value to not be an integer.
^See (Graham, Knuth & Patashnik 1994), which also defines for . Alternative generalizations, such as to two real or complex valued arguments using the Gamma function assign nonzero values to for , but this causes most binomial coefficient identities to fail, and thus is not widely used by the majority of definitions. One such choice of nonzero values leads to the aesthetically pleasing "Pascal windmill" in Hilton, Holton and Pedersen, Mathematical reflections: in a room with many mirrors, Springer, 1997, but causes even Pascal's identity to fail (at the origin).
^Boardman, Michael (2004), "The Egg-Drop Numbers", Mathematics Magazine, 77 (5): 368–372, doi:10.2307/3219201, JSTOR3219201, MR1573776, it is well known that there is no closed form (that is, direct formula) for the partial sum of binomial coefficients.
Grinshpan, A. Z. (2010), "Weighted inequalities and negative binomials", Advances in Applied Mathematics, 45 (4): 564–606, doi:10.1016/j.aam.2010.04.004
Knuth, Donald E. (1997). The Art of Computer Programming, Volume 1: Fundamental Algorithms (Third ed.). Addison-Wesley. pp. 52–74. ISBN0-201-89683-4.
Singmaster, David (1974). "Notes on binomial coefficients. III. Any integer divides almost all binomial coefficients". Journal of the London Mathematical Society. 8 (3): 555–560. doi:10.1112/jlms/s2-8.3.555.