Skip to Content
Concrete Math

Concrete Math

2018-08-10legacy assets may be incomplete

Concrete Math

Divisibility

因子求和恒等式

m/nam=km>0am[n=mk]\sum_{m/n}a_m = \sum_k\sum_{m>0}a_m[n=mk]

Prime

算数基本定理

正整数的数系表示<n2,n3,n5,...><n_2,n_3,n_5,...>

Prime Examples

素数

6k+1/6k-1
素数平方 8n+1

Euclid number

en=e1e2...en1+1=en12en1+1e_n=e_1*e_2*...*e_{n-1}+1=e_{n-1}^2-e_{n-1}+1

Mersenne number

2p12^p-1

素数分布

Pnnln(n)P_n\sim nln(n)
P1P2...Pn+2P_1P_2...P_n + 2P1P2...Pn+Pn+11P_1P_2...P_n+P_{n+1}-1 间都是合数 (可生成任意多合数)

Factorial factor

Stirling formula

n!2πn(ne)nn!\sim \sqrt{2\pi n}(\frac{n}{e})^n

相对误差约112n\frac{1}{12n}

整除阶乘的最大幂

ϵp(n!)=k>0[npk]<np1\epsilon_p(n!)=\sum_{k>0}[\frac{n}{p^k}] < \frac{n}{p-1} ϵ2(n!)=nv2(n)\epsilon_2(n!)=n-v_2(n)

Relative Primality

记号

gcd(a,b)=1=defabgcd(a,b)=1 \overset{def}{=} a\perp b
ab,acabca\perp b, a\perp c\rightarrow a\perp bc

Stern-Brocot树

mediant of m/n and m'/n' is m+mn+n\frac{m+m'}{n+n'}
从(0/1,1/0)构造可得此树;从(0/1,1/0,0/-1,-1/0,0/1)可得Stern-Brocot环
性质:任何阶段两个相邻分数有mnnm=1m'n-n'm=1

表示有理数的数系

LR字符串->有理数

M(S)=(n,n;m,m)M(S) = (n,n';m,m') M(I)=(1,0;0,1)M(I) = (1,0;0,1) L=(1,1;0,1)R=(1,0;1,1)L = (1,1;0,1) R = (1,0;1,1) Lk=(1,k;0,1)Rk=(1,0;k,1)L^k = (1,k;0,1) R^k = (1,0;k,1)

有理数->LR字符串

二分搜索

Farey serires

FnF_n是0与1间分母不超过N的所有最简分数的集合

The Congruence relation

规定

amod0=a;amod1=[a]a \bmod 0 = a; a \bmod 1 = [a]

平方剩余解的个数

x2=1(modp)kx^2 = 1 \pmod p^k p2x=1\1p\not=2,x=1\backslash-1 p=2,k>=3x=1\1\2k11\2k+1+1p=2,k>=3,x=1\backslash-1\backslash2^{k-1}-1\backslash2^{k+1}+1 若m有r个不同素因子,个数为2r+[8m]+[4m][2m]2^{r+[8|m]+[4|m]-[2|m]}

Wilson 定理

(n-1)!=-1(mod n) 等价于 n>1,n是素数

欧拉函数与莫比乌斯函数

性质

ϕ(x)=x(11p)\phi(x) = x*\prod(1-\frac{1}{p}) dmμ(d)=[m=1]\sum_{d|m}\mu(d)=[m=1]

反演定理

g=fI<>f=μgg = f*I <-> f = \mu * g g(x)=d>=1f(x/d)<>d>=1μ(d)g(x/d)g(x) = \sum_{d>=1}f(x/d) <-> \sum_{d>=1}\mu(d)g(x/d)