跳转至

数论入门

因数

整除

bb 能整除 aa,则记为 aba\mid b,如 2122\mid 12. 若 bb 不能整除 aa,则记为 aba\nmid b,如 5125\nmid 12.

aba\nmid b,则 b÷ab\div a 存在余数 rr0<r<a0<r<a,记 r=a mod br=a\ \mathrm{mod}\ b. 例如,3 mod 2=13\ \mathrm{mod}\ 2=1.

整除具有以下性质:

  1. aba\mid baca\mid c,则 x,y\forall x,y,有 axb+yca\mid xb+yc.
  2. aba\mid bbcb\mid c,则 aca\mid c.
  3. aba\mid bbab\mid a,则 a=±ba=\pm b.
  4. m0m\neq 0,则 aba\mid b,当且仅当 mambma\mid mb.

最大公因数与最小公倍数

a,ba,b 是两个不为 00 的整数,能使 dad\mid adbd\mid b 成立的最大整数 dd,称为 a,ba,b 的最大公因数,记作 gcd(a,b)\gcd(a,b).

a,ba,b 是两个不为 00 的整数,能使 ada\mid dbdb\mid d 成立的最小整数 dd,称为 a,ba,b 的最小公倍数,记作 lcm(a,b)\mathrm{lcm}(a,b).

  • a,bN, gcd(a,b)lcm(a,b)=ab\forall a,b\in\N,\ \gcd(a,b)\cdot\mathrm{lcm}(a,b)=ab

证明:设 d=gcd(a,b),a0=ad,b0=bd\displaystyle d=\gcd(a,b),a_0=\frac{a}{d},b_0=\frac{b}{d}. 根据最大公因数的定义,有 gcd(a0,b0)=1\gcd(a_0,b_0)=1         \\\ \ \ \ \ \ \ \ \ \text{} 再根据最小公倍数的定义,有 lcm(a,b)=a0b0\mathrm{lcm}(a,b)=a_0\cdot b_0.          \\\ \ \ \ \ \ \ \ \ \text{} 于是 lcm(a,b)=lcm(a0d,b0d)=lcm(a0,b0)d=a0b0d=abd\displaystyle\mathrm{lcm}(a,b)=\mathrm{lcm}(a_0\cdot d,b_0\cdot d)=\mathrm{lcm}(a_0,b_0)\cdot d=a_0\cdot b_0\cdot d=\frac{a\cdot b}{d},原命题得证。

  • 九章算术 更相减损术: gcd(a,b)=gcd(b,ab)=gcd(a,ab),gcd(2a,2b)=2gcd(a,b)\gcd(a,b)=\gcd(b,a-b)=\gcd(a,a-b),\gcd(2a,2b)=2\gcd(a,b).

证明:da,db    d(ab)d\mid a,d\mid b\implies d\mid(a-b).

  • 欧几里得算法:a,bN,b0,gcd(a,b)=gcd(b,a mod b)\forall a,b\in\N,b\neq 0,\gcd(a,b)=\gcd(b,a\ \mathrm{mod}\ b).

证明:分类讨论。

  1. a<ba<b,则有 gcd(b,a mod b)=gcd(b,a)=gcd(a,b)\gcd(b,a\ \mathrm{mod}\ b)=\gcd(b,a)=\gcd(a,b).
  2. aba\geq b,不妨设 a=qb+r (0r<b)a=q\cdot b+r\ (0\leq r<b)。显然 r=a mod br=a\ \mathrm{mod}\ b. 对于 a,ba,b 的任意公因数 dd\\ 因为 da,dqbd\mid a,d\mid q\cdot b,故 d(aqb)d\mid (a-q\cdot b),即 drd\mid r. 因此 dd 也是 b,rb,r 的公因数,反之亦成立。\\a,ba,b 的公因数集合与 b,a mod bb,a\ \mathrm{mod}\ b 的公因数集合相同。于是它们的最大公因数也相等。

裴蜀定理

a,bZ,ab0a,b\in\Z,ab\neq 0,则 x,yZ\exist x,y\in\Z 使 ax+by=gcd(a,b)=a(x+bu)+b(yau)ax+by=\gcd(a,b)=a(x+bu)+b(y-au)

由上易推出两个整数 a,ba,b 互素的充分必要条件是 x,yZ\exist x,y\in\Z 使得 ax+by=1ax+by=1

比如要证明 21n+421n+414n+314n+3 互素,知道 3(14n+3)2(21n+4)=13(14n+3)-2(21n+4)=1 即可。

例 1:证明 n!+1n!+1(n+1)!+1(n+1)!+1 互素。

(n+1)(n!+1)[(n+1)!+1]=n(n+1)(n!+1)-[(n+1)!+1]=n,于是 [d=gcd(n!+1,(n+1)!+1)]n[d=\gcd(n!+1,(n+1)!+1)]|n,又因为 dn!d|n!,结合 d(n!+1)d|(n!+1) 得到 d=1d=1,原命题得证。


例 2:记费马数 Fk=22k+1,k0F_k=2^{2^k}+1,k\geq 0,证明若 mnm\neq ngcd(Fm,Fn)=1\gcd(F_m,F_n)=1.

m>nm>n,只要利用 Fn(Fm2)F_n|(F_m-2) 证明 xZ\exist x\in\Z 使得 Fm+xFn=2F_m+xF_n=2【基本技巧例 2】

[d=gcd(Fm,Fn)]2[d=\gcd(F_m,F_n)]|2. FnF_n 显然为奇数故 d=1d=1. 因此费马数两两互素。


例 3:设 a>1,m,n>0a>1,m,n>0,证明 gcd(am1,an1)=agcd(m,n)1\boxed{\gcd(a^m-1,a^n-1)=a^{\gcd(m,n)}-1}

D=gcd(am1,an1)D=\gcd(a^m-1,a^n-1),只要证明 [agcd(m,n)1]D[a^{\gcd(m,n)}-1]|DD[agcd(m,n)1]D|[a^{\gcd(m,n)}-1] 即可。

显然 [agcd(m,n)1](am1)[a^{\gcd(m,n)}-1]|(a^m-1)[agcd(m,n)1](an1)[a^{\gcd(m,n)}-1]|(a^n-1),可以证明 [agcd(m,n)1]D[a^{\gcd(m,n)}-1]|D.

d=gcd(m,n)d=\gcd(m,n),可选 u,v>0u,v>0,使 munv=dmu-nv=d. 因此 D(am1),D(amu1)D|(a^m-1),D|(a^{mu}-1),扩展开来 D(amuanv)=anv(ad1)D|(a^{mu}-a^{nv})=a^{nv}(a^d-1),故 D(ad1)D|(a^d-1). 综合上述原命题得证。


例 4:设 m,n>0,mn(m2+n2)m,n>0,mn|(m^2+n^2),证明 m=nm=n.

gcd(m,n)=d,m=m1d,n=n1d,gcd(m1,n1)=1\gcd(m,n)=d,m=m_1d,n=n_1d,\gcd(m_1,n_1)=1①. 已知化为 m1n1(m12+n12)m_1n_1|(m_1^2+n_1^2),从而 m1n12,n1m12m_1|n_1^2,n_1|m_1^2. 结合①可得 m1=n1=1,m=nm_1=n_1=1,m=n.


例 5:设 kk 为正奇数,证明 i=1ni\displaystyle\sum_{i=1}^ni 整除 i=1nik\displaystyle\sum_{i=1}^ni^k.

等价于证明 n2i=1nikn|2\displaystyle\sum_{i=1}^ni^k(n+1)2i=1nik(n+1)|2\displaystyle\sum_{i=1}^ni^k,显然

2i=1nik=[1k+(n1)k]+[2k+(n2)k]++[(n1)k+1k]+2nk=[1k+nk]+[2k+(n1)k]++[nk+1k]\begin{aligned}2\sum_{i=1}^ni^k&=[1^k+(n-1)^k]+[2^k+(n-2)^k]+\dots+[(n-1)^k+1^k]+2n^k\\&=[1^k+n^k]+[2^k+(n-1)^k]+\dots+[n^k+1^k]\end{aligned}

nnn+1n+1 的倍数。

由上可知,为了证明 bab|a,只需将 bb 分解成若干个两两互素的整数 b1,b2,,bnb_1,b_2,\dots,b_n 之积,证明 biab_i|a 即可。

质数

  • 定义:一个正整数无法被除了 11 和它自身之外的任何自然数整除,则成为质数,否则为合数。
  • 数量:对于一个足够大的整数 NN,不超过 NN 的质数有 π(x)NlnN\displaystyle\pi(x)\approx\frac{N}{\ln N} 个,所以第 nn 个质数约为 nlnnn\ln n.
  • 判定:试除法,若 NN 为合数,则存在一个能整除 NN 的数 TT2TN2 \leq T \leq \sqrt{N}.

  • 算数基本定理:N=p1c1p2c2pncnN = p_1^{c_1}p_2^{c_2}\dots p_n^{c_n}NN 的正约数个数 =i=1n(ci+1)=\red{\boxed{\displaystyle\prod_{i=1}^{n}(c_i+1)}}NN 的所有正约数的和 =i=1n(j=0ci(pi)j)=i=1npici+11pi1=\displaystyle\prod_{i=1}^{n}(\displaystyle\sum_{j=0}^{c_i}(p_i)^j)=\red{\boxed{\prod_{i=1}^n\frac{p_i^{c_i+1}-1}{p_i-1}}}.

一个正整数 NN 的约数个数上界为 2N2\sqrt{N}11 ~ NN 每个数的约数个数的总和大约为 NlogNN\log N.

  • bacb|ac,若 gcd(b,c)=1\gcd(b,c)=1,则 bab|a.
  • a,bNa,b\in\N^*a,ba,b 之积是一个整数的 k(k2)k(k\geq 2) 次幂,若 gcd(a,b)=1\gcd(a,b)=1,则 a,ba,b 都是整数的 kk 次幂。一般地,设正整数 a,b,,ca,b,\dots,c 之积是一个整数的 kk 次幂,且 a,b,,ca,b,\dots,c 两两互质,则 a,b,,ca,b,\dots,c 都是整数的 kk 次幂。

  • 对任意正整数 mm 及质数 pp,记 pcmp^c\Vert m 表示 pcmp^c|mpc+1mp^{c+1}\nmid m,即 pcp^cmm 的标准分解中出现的 pp 的幂。

  • n>1,pn>1,p 为质数,pcpn!p^{c_p}\Vert n!,则 cp=l=1+nplc_p=\displaystyle\sum_{l=1}^{+\infty}\lfloor\frac{n}{p^l}\rfloor

例 1:证明无穷数列 10001,100010001,10001,100010001,\dots 中无质数。

an=1+104++104(n1)=104n11041\displaystyle a_n=1+10^4+\dots+10^{4(n-1)}=\frac{10^{4n}-1}{10^4-1},接下来分 nn 为奇偶讨论即可。

a2k=108k11041=108k1108110811041a_{2k}=\frac{10^{8k}-1}{10^4-1}=\frac{10^{8k}-1}{10^8-1}\cdot\frac{10^8-1}{10^4-1}
a2k+1=104(2k+1)11041=102(2k+1)11021102(2k+1)+1102+1a_{2k+1}=\frac{10^{4(2k+1)}-1}{10^4-1}=\frac{10^{2(2k+1)}-1}{10^2-1}\cdot\frac{10^{2(2k+1)}+1}{10^2+1}

例 2:证明 nN,n>1,n4+4n\forall n\in\N^*,n>1,n^4+4^n 不是质数。

nn 为偶数显然。nn 为奇数时,设 n=2k+1n=2k+1

n4+4n=n4+442k=n4+4(2k)4=n4+4n2(2k)2+4(2k)44n2(2k)2=(n2+222k)2(2n2k)2=(n2+2k+1n+22k+1)(n22k+1n+22k+1)\begin{aligned}n^4+4^n&=n^4+4\cdot 4^{2k}=n^4+4\cdot (2k)^4=n^4+4n^2\cdot(2^k)^2+4\cdot(2k)^4-4n^2\cdot(2^k)^2\\&=(n^2+2\cdot 2^{2k})^2-(2n\cdot 2^k)^2=(n^2+2^{k+1}n+2^{2k+1})(n^2-2^{k+1}n+2^{2k+1})\end{aligned}

即把 4n4^n 看成 4y44y^4,有 x4+4y4=(x2+2y2+2xy)(x2+2y22xy)\red{\boxed{x^4+4y^4=(x^2+2y^2+2xy)(x^2+2y^2-2xy)}}

练 1:设 a,b,c,dN,ab=cda,b,c,d\in\N^*,ab=cd,证明 a+b+c+da+b+c+d 不是质数。

例 3:证明:若 a,bZ,2a2+a=3b2+ba,b\in\Z,2a^2+a=3b^2+b,则 aba-b2a+2b+12a+2b+1 都为完全平方数。

已知化为 (ab)(2a+2b+1)=b2(a-b)(2a+2b+1)=b^2,设 d=gcd(ab,2a+2b+1)d=\gcd(a-b,2a+2b+1),若 d>1d>1,则 dd 有质因子 pppb2    pbp|b^2\implies p|bp(ab)    pa    p1p|(a-b)\implies p|a\implies p|1,这不可能,因此 d=1d=1. 因为 b2b^2 是完全平方数,由此可得 ab,2a+2b+1|a-b|,|2a+2b+1| 也是完全平方。

只要证 aba-b2a+2b+12a+2b+1 不能同时 <0<0 即可。设 ab<0a-b<0,则 ba=r2b-a=r^2. 显然 rb,rar|b,r|a,令 b=b1r,a=a1rb=b_1r,a=a_1r,代入题目 a12+6a1r+3r2+1=0a_1^2+6a_1r+3r^2+1=0,得 a1=3r±6r21a_1=-3r\pm\sqrt{6r^2-1},显然 6r216r^2-1 应为完全平方数,而 6r212 (mod 3)6r^2-1\equiv 2\ (\mathrm{mod}\ 3),矛盾,原命题得证。

不定方程

例 1:证明两个连续正整数之积不能是完全平方或完全立方。

反证法,即设 x(x+1)=y2x(x+1)=y^2 有解,则 (2x+1)2=4y2+1(2x+1)^2=4y^2+1,分解为

(2x+1+2y)(2x+12y)=1    {2x+1+2y=12x+12y=1    x=y=0(2x+1+2y)(2x+1-2y)=1\implies\begin{cases}2x+1+2y=1\\2x+1-2y=1\end{cases}\implies x=y=0

显然不可能。类似的,对于完全立方,由于 xxx+1x+1 互质,可得它们都是立方数。设 x=u3,x+1=v3,y=uv,v3u3=1=(vu)(v2+uv+u2)x=u^3,x+1=v^3,y=uv,v^3-u^3=1=(v-u)(v^2+uv+u^2),显然不可能。

类似的,可证明连续两个正整数之积不能是整数的 k(k2)k(k\geq 2) 次幂。

练 1:设 kN,k2k\in\N^*,k\geq 2,证明:连续三个正整数之积不能是整数的 kk 次幂。

提示:x(x21)=yk    akbk=(ab)k,1=a2kbk=(a2)kbkx(x^2-1)=y^k\implies a^kb^k=(ab)^k,1=a^{2k}-b^k=(a^2)^k-b^k,导出矛盾。


例 2:证明方程 y2+y=x+x2+x3y^2+y=x+x^2+x^3 没有 x0x\neq 0 的整数解。

转化为 (yx)(y+x+1)=x3(y-x)(y+x+1)=x^3,可以证明 gcd(yx,y+x+1)=1\gcd(y-x,y+x+1)=1,设 yx=a3,y+x+1=b3,x=aby-x=a^3,y+x+1=b^3,x=ab,可得 b3a3=(ba)(b2+ab+a2)=2ab+1b^3-a^3=(b-a)(b^2+ab+a^2)=2ab+1③,证明该方程无解即可。ab>0ab>0 时,ba1b-a\geq 1,③的左边 3ab>\geq 3ab> 右边;ab<0ab<0 时,③的左边的绝对值 2(a2+b2ab)>2ab>\geq 2(a^2+b^2-|ab|)>2|ab|> ③的右边的绝对值。因此原命题成立。

练 2:求方程 (x2y2)2=1+16y(x^2-y^2)^2=1+16y 的所有整数解。

提示:y0,xy\geq 0,|x|\geqy+1\leq y+1,均可得出 (x2y2)2(2y1)2(x^2-y^2)^2\geq(2y-1)^2,得 (2y1)21+16y(2y-1)^2\leq 1+16y,答案 (x,y)=(±1,0),(±4,3)(±4,5)(x,y)=(\pm 1,0),(\pm 4,3)(\pm4,5).


例 3:设 x,y,zN,2xx=yy+zzx,y,z\in\N^*,2x^x=y^y+z^z,证明 x=y=zx=y=z.

(x+1)x+1>xx+1+(x+1)xx>2xx    y,zx(x+1)^{x+1}>x^{x+1}+(x+1)x^x>2x^x\implies y,z\leq x

反之,设 yx+1y\geq x+1,则 yy+zz>(x+1)x+1>2xxy^y+z^z>(x+1)^{x+1}>2x^x,矛盾。而 yy+zzxx+xx=2xxy^y+z^z\leq x^x+x^x=2x^x,所以 x=y=zx=y=z.

欧拉函数

11 ~ NN 中与 NN 互质的数的个数被称为欧拉函数,记为 φ(N)\varphi(N). 用数学符号表示即为 φ(N)=i=1N[gcd(i,N)=1]\varphi(N)=\displaystyle\sum_{i=1}^{N}[\gcd(i,N)=1].

N=p1c1p2c2pncnN = p_1^{c_1}p_2^{c_2}\dots p_n^{c_n},则 φ(N)=N×p11p1×p21p2××pn1pn=N质数pN(11p)\displaystyle\varphi(N)=N\times\frac{p_1-1}{p_1}\times\frac{p_2-1}{p_2}\times\dots\times\frac{p_n-1}{p_n}=N\cdot\displaystyle\prod_{质数p|N}\left(1-\frac{1}{p}\right)

证明:设 p,qp,qNN 的不同质因子,11 ~ NNpp 的倍数有 Np\displaystyle \frac{N}{p} 个,qq 的倍数有 Nq\displaystyle \frac{N}{q} 个。若把 Np+Nq\displaystyle \frac{N}{p}+\frac{N}{q} 个数去掉,则 Npq\displaystyle\frac{N}{pq} 被计算了 22 次(容斥原理)。因此, 11 ~ NN 中不与 NN 含有相同质因子的 ppqq 数量为 NNpNq+Npq=N(11p)(11q)\displaystyle N-\frac{N}{p}-\frac{N}{q}+\frac{N}{pq}=N\left(1-\frac{1}{p}\right)\left(1-\frac{1}{q}\right),对 NN 的全部质因子继续容斥即可得到公式。

nn112233445566778899101011111212131314141515
φ(n)\varphi(n)1111222244226644664410104412126688

欧拉函数的性质:

  1. n>1\forall n>111 ~ nn 中与 nn 互质的数的和为 n×φ(n)2\displaystyle\frac{n\times\varphi(n)}{2}.
  2. a,ba,b 互质,则 φ(ab)=φ(a)φ(b)\varphi(ab)=\varphi(a)\varphi(b).
  3. pp 为质数,φ(p)=p1\varphi(p)=p-1.
  4. pp 为质数,kN,φ(pk)=pk1×(p1)=pkpk1k\in\N^*,\varphi(p^k)=p^{k-1}\times(p-1)=p^k-p^{k-1}.
  5. pp 为质数,若 pnp\mid np2np^2\mid n,则 φ(n)=φ(np)×p\displaystyle\varphi(n)=\varphi\left(\frac{n}{p}\right)\times p.
  6. pp 为质数,若 pnp\mid np2np^2\nmid n,则 φ(n)=φ(np)×(p1)\displaystyle\varphi(n)=\varphi\left(\frac{n}{p}\right)\times (p-1).
  7. nn 为奇数,则 φ(2n)=φ(n)\varphi(2n)=\varphi(n).
  8. n3, nN, 2φ(n)\forall n\geq 3,\ n\in\N^*,\ 2\mid\varphi(n).
  9. 欧拉反演:dnφ(d)=n\displaystyle\sum_{d\mid n}\varphi(d)=n.

证明:

  1. 因为 gcd(n,x)=gcd(n,nx)\gcd(n,x)=\gcd(n,n-x),所以与 nn 不互质的数 x,nxx,n-x 成对出现,平均值为 n2\displaystyle\frac{n}{2}. \\ 因此与 nn 互质的数的平均值也是 n2\displaystyle\frac{n}{2},进而得到性质 1。
  2. 根据欧拉函数的计算式可直接获得性质 2。
  3. 根据欧拉函数的定义可直接获得性质 3。
  4. 11 ~ pkp^{k} 中的所有数,除了 pk1p^{k-1}pp 的倍数外都与 pkp^k 互素。
  5. pnp\mid np2np^2\mid n,则 n,np\displaystyle n,\frac{n}{p} 包含相同的质因子,只是 pp 的指数不同。\\ 按照欧拉函数的计算公式,φ(n)φ(np)=nnp=p\displaystyle \frac{\varphi(n)}{\varphi(\frac{n}{p})}=\frac{n}{\frac{n}{p}}=p,得到性质 5。
  6. pnp\mid np2np^2\mid n,则 n,np\displaystyle n,\frac{n}{p} 互质。因为 φ(n)=φ(np)φ(p)\displaystyle\varphi(n)=\varphi\left(\frac{n}{p}\right)\varphi(p),而 φ(p)=p1\varphi(p)=p-1,得到性质 6。
  7. m<2n,mN\forall m<2n,m\in\N^*,若 2m2\mid m,则 gcd(m,2n)1\gcd(m,2n)\neq 1,也就是对欧拉函数的值没有贡献;\\2m2\nmid m,则 gcd(m,2n)=1\gcd(m,2n)=1,欧拉函数的值也就与 2n2n 中的偶质因子无关。
  8. n3n\geq 3 时,与 nn 互质的数不可能为 n2    x<n,gcd(x,n)=gcd(nx,n)\displaystyle\frac{n}{2}\implies\forall x<n,\gcd(x,n)=\gcd(n-x,n),也就是存在一一对应关系。
  9. f(n)=dnφ(d)f(n)=\displaystyle\sum_{d\mid n}\varphi(d),利用 φ\varphi 是积性函数,得到:\\n,mn,m 互质,则 f(nm)=dnmφ(d)=dnφ(d)dmφ(d)=f(n)f(m)f(nm)=\displaystyle\sum_{d\mid nm}\varphi(d)=\displaystyle\sum_{d\mid n}\varphi(d)\cdot\displaystyle\sum_{d\mid m}\varphi(d)=f(n)f(m),即 f(n)f(n) 是积性函数。\\ 对于单个质因子有:f(pm)=dpmφ(d)=i=0mφ(pi)=φ(1)+φ(p)+φ(p2)+φ(p3)++φ(pm)=1+(p1)+(p1)p+(p1)p2++(p1)pm1=1+(p1)+(p2p)+(p3p2)++(pmpm1)=pm\begin{aligned}f(p^m)&=\displaystyle\sum_{d\mid p^m}\varphi(d)=\displaystyle\sum_{i=0}^{m}\varphi(p^i)=\varphi(1)+\varphi(p)+\varphi(p^2)+\varphi(p^3)+\dots+\varphi(p^m)\\&= 1+(p-1)+(p-1)p+(p-1)p^2+\dots+(p-1)p^{m-1}\\&=1+(p-1)+(p^2-p)+(p^3-p^2)+\dots+(p^m-p^{m-1})=p^m\end{aligned} \\ 所以 f(n)=i=1mf(pici)=i=1mpici=nf(n)=\displaystyle\prod_{i=1}^{m}f(p_i^{c_i})=\displaystyle\prod_{i=1}^{m}p_i^{c_i}=n

积性函数与完全积性函数

若函数 f(x)f(x) 满足 f(1)=1f(1)=1x,yN,gcd(x,y)=1\forall x,y\in\N^{*},\gcd(x,y)=1 都有 f(xy)=f(x)f(y)f(xy)=f(x)f(y),则 f(x)f(x) 是积性函数。

若函数 f(x)f(x) 满足 f(1)=1f(1)=1x,yN\forall x,y\in\N^{*} 都有 f(xy)=f(x)f(y)f(xy)=f(x)f(y),则 f(x)f(x) 是完全积性函数。

性质:

  1. f(x)f(x)g(x)g(x) 均为积性函数,则以下函数也为积性函数: $h(x)=f(xp)h(x)=fp(x)h(x)=f(x)g(x)h(x)=dxf(d)g(xd)\begin{aligned}h(x)&=f(x^p) \\ h(x)&=f^p(x) \\ h(x)&=f(x)g(x) \\ h(x)&=\displaystyle\sum_{d\mid x}f(d)g\left(\frac{x}{d}\right)\end{aligned}$
  2. x=picix=\displaystyle\prod p_i^{c_i}

f(x)f(x) 为积性函数,则 f(x)=f(pici)f(x)=\displaystyle\prod f(p_i^{c_i})

f(x)f(x) 为完全积性函数,则 f(x)=fci(pi)f(x)=\displaystyle\prod f^{c_i}(p_i)

例子:

积性函数:φ(n),σk(n)=dndk,σ0(n)\varphi(n),\sigma _k(n)=\displaystyle\sum_{d\mid n}d^k,\sigma _0(n) 通常简记为 d(n)d(n)τ(n)\tau(n)σ1(n)\sigma _1(n) 通常简记为 σ(n)\sigma(n)

完全积性函数:ε(n)=[n=1],idk(n)=nk,id1(n)\varepsilon(n)=[n=1],\text{id}_k(n)=n^k,\text{id}_1(n) 通常简记为 id(n),f(n)=1\text{id}(n),f(n)=1.

同余

  • 平方数 mod 40\mathrm{mod}\ 4\equiv 01,mod 801,\mathrm{mod}\ 8\equiv0114,mod 304,\mathrm{mod}\ 3\equiv 01,mod 501,\mathrm{mod}\ 5\equiv0±1\pm 1.
  • 立方数 mod 90\mathrm{mod}\ 9\equiv 0±1\pm 1.
  • 四次幂 mod 160\mathrm{mod}\ 16\equiv 011.

例如可证明 a,b,c,dN,a4b+da4c+da,b,c,d\in\N^*,a^{4b+d}-a^{4c+d} 能被 240=24×3×5240=2^4\times 3\times 5 整除。

例 1:设 a,b,cZ,a+b+c=0,d=a1999+b1999+c1999a,b,c\in\Z,a+b+c=0,d=a^{1999}+b^{1999}+c^{1999},证明 dd 可被 66 整除。

u1999u(mod 2)    da+b+c0(mod 2)    2du^{1999}\equiv u(\mathrm{mod}\ 2)\implies d\equiv a+b+c\equiv 0(\mathrm{mod}\ 2)\implies2|d
u3u(mod 3)    u1999uu1998uu666uu222u75u25u9u(mod 3)u^3\equiv u(\mathrm{mod}\ 3)\implies u^{1999}\equiv u\cdot u^{1998}\equiv u\cdot u^{666}\equiv u\cdot u^{222}\equiv u^{75}\equiv u^{25}\equiv u^9\equiv u(\mathrm{mod}\ 3)

注意 u3u(mod 3)u^3\equiv u(\mathrm{mod}\ 3) 是费马小定理的特殊情形。

练 1:设 x,y,zZ,(xy)(yz)(zx)=x+y+zx,y,z\in\Z,(x-y)(y-z)(z-x)=x+y+z,证明 27(x+y+z)27|(x+y+z).

只要证明 x,y,zx,y,z 两两模 33 同余即可。

练 2:证明 11,111,1111,11,111,1111,\dots 不是完全平方数。只要知道 11 mod 4011\ \mathrm{mod}\ 4\neq 0 即可。


例 2:数列 {xn}\set{x_n}1,3,5,11,1,3,5,11,\dots,满足 xn+1=xn+2xn1x_{n+1}=x_n+2x_{n-1}.

数列 {yn}\set{y_n}7,17,55,161,7,17,55,161,\dots,满足 yn+1=2yn+3yn1y_{n+1}=2y_n+3y_{n-1}. 证明:两个数列无相同项。

考虑以 88 为模,因为 x23,x35x_2\equiv 3,x_3\equiv 5. 若 xn13,xn5x_{n-1}\equiv 3,x_n\equiv 5,则 xn+13,xn+25x_{n+1}\equiv 3,x_{n+2}\equiv 5.

显然 {xn}\set{x_n} 是模 88 后的周期数列,同理 {yn}\set{y_n} 也是( 7,1,7,1,7,1,7,1,\dots ),于是命题得证。

费马小定理与欧拉定理

若整数 aa 和整数 bb 除以正整数 mm 的余数相等,则称 a,ba,bmm 同余,记为 ab (mod m)a\equiv b\ (\mathrm{mod}\ m).

并且注意到 k mod i=kkii\displaystyle k\ \mathrm{mod}\ i=k-\left\lfloor\frac{k}{i}\right\rfloor\cdot i,且同余满足同加性、同乘性、同幂性,但不满足同除性。

对于 a[0,m1]\forall a \in [0, m - 1],集合 {a+km (kZ)}\set{a+km\ (k\in\Z)} 的所有数模 mm 同余,余数都是 aa. 该集合称为一个模 mm 的同余类,简记为 a\overline{a}.

mm 的同余类一共有 mm 个,分别为 0,1,2,,m1\overline{0},\overline{1},\overline{2},\dots,\overline{m-1}. 它们构成 mm 的完全剩余系。

11 ~ mm 中与 mm 互质的数代表的同余类共有 φ(m)\varphi(m) 个,它们构成 mm 的简化剩余系。\\ 例如,模 1010 的简化剩余系为 {1,3,7,9}\set{\overline{1},\overline{3},\overline{7},\overline{9}}

若从某个非空数集中任选两个元素(同一元素可重复选出),选出的这两个元素通过某种(或几种)运算后的得数仍是该数集中的元素,那么,就说该集合对于这种(或几种)运算是封闭的。\\ 例如若一个集合中的元素,如果能够做到做加法运算的结果还在这个集合中,就说这个集合对加法运算封闭。\\ 例如 N\N 对加法、乘法运算是封闭的;Z\Z 对加、减、乘法运算是封闭的;Q,C\mathbb{Q}, \mathbb{C} 对四则运算是封闭的。

简化剩余系关于模 mm 乘法封闭。这是因为若 a,b (1a,bm)a,b\ (1\leq a,b\leq m)mm 互质,则 abab 也与 mm 互质。\\ 由余数的定义得 ab mod mab\ \mathrm{mod}\ m 也与 mm 互质,即 ab mod mab\ \mathrm{mod}\ m 也属于 mm 的简化剩余系。

  • 费马小定理:若 pp 是质数,则对于任意整数 aa,有 apa (mod p)a^p\equiv a\ (\mathrm{mod}\ p).
  • 欧拉定理:若正整数 a,na,n 互质,则 aφ(n)1 (mod n)a^{\varphi(n)}\equiv 1\ (\mathrm{mod}\ n).

证明:设 nn 的简化剩余系为 {a1,a2,,aφ(n)}\set{\overline{a_1},\overline{a_2},\dots,\overline{a_{\varphi(n)}}}. ai,aj\forall a_i,a_j,若 aaiaaj (mod n)a\cdot a_i\equiv a\cdot a_j\ (\mathrm{mod}\ n),则 a(aiaj)0.a\cdot(a_i-a_j)\equiv 0. 因为 a,na,n 互质,所以 aiaja_i\equiv a_j. 故当 aiaja_i\neq a_j 时,aai,aajaa_i,aa_j 也代表不同的同余类。

又因为简化剩余系关于模 mm 乘法封闭,故 aa1\overline{aa_1} 也在简化剩余系集合中。因此,集合 {a1,a2,,aφ(n)}\set{\overline{a_1},\overline{a_2},\dots,\overline{a_{\varphi(n)}}} 与集合 {aa1,aa2,,aaφ(n)}\set{\overline{aa_1},\overline{aa_2},\dots,\overline{aa_{\varphi(n)}}} 都能表示 nn 的简化剩余系。综上所述:

aφ(n)a1a2aφ(n)(aa1)(aa2)(aaφ(n))a1a2aφ(n) (mod n)a^{\varphi(n)}a_1a_2\dots a_{\varphi(n)}\equiv(aa_1)(aa_2)\dots(aa_{\varphi(n)})\equiv a_1a_2\dots a_{\varphi(n)}\ (\mathrm{mod}\ n)

因此 aφ(n)1 (mod n)a^{\varphi(n)}\equiv 1\ (\mathrm{mod}\ n)。当 pp 为质数时,满足 φ(p)=p1\varphi(p)=p-1,两边同乘 aa 即可得到费马小定理。

另外,当 aapp 的倍数,费马小定理显然成立。


  • 欧拉定理的推论:若正整数 a,na,n 互质,则对于任意正整数 bb,有 abab mod φ(n) (mod n)a^b\equiv a^{b\ \mathrm{mod}\ \varphi(n)}\ (\mathrm{mod}\ n)

证明:设 b=qφ(n)+rb=q\cdot \varphi(n)+r,其中 0r<φ(n)0\leq r<\varphi(n),即 r=b mod φ(n)r=b\ \mathrm{mod}\ \varphi(n). 于是有:

abaqφ(n)+r(aφ(n))qar1qararab mod φ(n) (mod n)a^b\equiv a^{q\cdot\varphi(n)+r}\equiv (a^{\varphi(n)})^q\cdot a^r\equiv 1^q\cdot a^r\equiv a^r\equiv a^{b\ \mathrm{mod}\ \varphi(n)}\ (\mathrm{mod}\ n)

证毕。面对 a+b,ab,aba+b,a-b,a\cdot b 这样的算式,可以在计算前先把 a,ba,bpp 取模。面对 aba^b 这样的乘方算式,可以先把底数对 pp 取模、指数对 φ(p)\varphi(p) 取模,再计算乘方。

ab(a mod p)b mod φ(p) (mod p)a^b\equiv (a\ \mathrm{mod}\ p)^{b\ \mathrm{mod}\ \varphi(p)}\ (\mathrm{mod}\ p)


  • 扩展欧拉定理

$ab{ab mod φ(m)if gcd(a,m)=1abif gcd(a,m)1,b<φ(m)ab mod φ(m)+φ(m)if gcd(a,m)1,bφ(m) (mod m)a^b\equiv\begin{cases}a^{b\ \mathrm{mod}\ \varphi(m)}&\text{if }\gcd(a,m)=1\\ a^b&\text{if }\gcd(a,m)\neq 1,b<\varphi(m)\\ a^{b\ \mathrm{mod}\ \varphi(m)+\varphi(m)}&\text{if }\gcd(a,m)\neq 1,b\geq\varphi(m)\end{cases}\ (\mathrm{mod}\ m)$

证明十分显然,略。


  • 若正整数 a,na,n 互质,则满足 ax1 (mod n)a^x\equiv 1\ (\mathrm{mod}\ n) 的最小正整数 x0x_0φ(n)\varphi(n) 的约数。

反证法,假设满足 ax1 (mod n)a^x\equiv 1\ (\mathrm{mod}\ n) 的最小正整数 x0x_0 不能整除 φ(n)\varphi(n).

φ(n)=qx0+r (0<r<x0)\varphi(n)=qx_0+r\ (0<r<x_0),因为 ax01 (mod n)a^{x_0}\equiv 1\ (\mathrm{mod}\ n),所以 aqx01 (mod n)a^{qx_0}\equiv 1\ (\mathrm{mod}\ n). 根据欧拉定理 aφ(n)1 (mod n)a^{\varphi(n)}\equiv 1\ (\mathrm{mod}\ n),所以 ar1 (mod n)a^r\equiv 1\ (\mathrm{mod}\ n). 这与 x0x_0 最小矛盾。故假设不成立,原命题成立。

中国剩余定理

m1,m2,,mnm_1,m_2,\dots,m_n 是两两互质的整数,m=i=1nmi,Mi=mmi,tim=\displaystyle\prod_{i=1}^n m_i,M_i=\frac{m}{m_i},t_i 是线性同余方程 Miti1 (mod mi)M_it_i\equiv 1\ (\mathrm{mod}\ m_i) 的一个解。对于任意的 nn 个整数,方程组

{xa1 (mod m1)xa2 (mod m2)            xan (mod mn)\begin{cases}x\equiv a_1\ (\mathrm{mod}\ m_1) \\ x\equiv a_2\ (\mathrm{mod}\ m_2)\\ \ \ \ \ \ \ \ \ \ \ \ \ \vdots\\ x\equiv a_n\ (\mathrm{mod}\ m_n)\end{cases}

有整数解,解为 x=i=1naiMitix=\displaystyle\sum_{i=1}^{n}a_iM_it_i,通解可以表示为 x+km(kZ)x+km(k\in\Z),最小非负整数解为 x mod mx\ \mathrm{mod}\ m.

证明:因为 Mi=mmi\displaystyle M_i=\frac{m}{m_i} 是除了 mim_i 之外所有模数的倍数,所以 ki,aiMiti0 (mod mi)\forall k\neq i,a_iM_it_i\equiv 0\ (\mathrm{mod}\ m_i)          \\ \ \ \ \ \ \ \ \ \ \ \text{} 所以代入 x=i=1naiMitix=\displaystyle\sum_{i=1}^{n}a_iM_it_i,原方程组成立。

威尔逊定理

  • pp 为素数 (p1)!1 (mod p)\xLeftrightarrow{}(p-1)!\equiv -1\ (\mathrm{mod}\ p)

证明:

  1. 充分性

对于 pp 不是素数的情况,有 {p=1(11)!0 (mod 1)p=4(41)!2 (mod 4)p>4分类讨论得出 (p1)!0 (mod p)\begin{cases} p=1 & (1-1)!\equiv 0\ (\mathrm{mod}\ 1) \\ p=4 & (4-1)!\equiv 2\ (\mathrm{mod}\ 4) \\ p>4 & 分类讨论得出\ (p-1)!\equiv 0\ (\mathrm{mod}\ p) \end{cases}

(a) 当 pp 为完全平方数。

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
则 $p=k^2$ 且 $k>2$,用相减法比较 $2k$ 与 $p$ 的大小得 $2k<p$,于是有

$$\begin{aligned}(p-1)!&=1\times 2\times\dots\times k\times\dots\times 2k\times\dots\times (p-1)\\ &=2nk^2\\&=2np\end{aligned}$$

所以 $(p-1)!\equiv 0\ (\mathrm{mod}\ p)$ 且 $p$ 为完全平方数。

(b) 当 $p$ 不为完全平方数。

则 $p$ 可以表示为两个不相等的数 $a$ 和 $b$ 的乘积,设 $a<b$,则有 $p=ab,1<a<b<p$

$$\begin{aligned}(p-1)!&=1\times 2\times\dots\times a\times\dots\times b\times\dots\times (p-1)\\&=a\times b\times n\\ &=np\end{aligned}$$

所以 $(p-1)!\equiv 0\ (\mathrm{mod}\ p)$ 且 $p$ 不为完全平方数。
  1. 必要性

pp 为素数时,考虑二次剩余式 x21 (mod p)x^2\equiv 1\ (\mathrm{mod}\ p),化简得 (x1)(x+1)0 (mod p)(x-1)(x+1)\equiv 0\ (\mathrm{mod}\ p)

于是 x1 (mod p)x\equiv 1\ (\mathrm{mod}\ p)xp1 (mod p)x\equiv p-1\ (\mathrm{mod}\ p)。现在先抛开 11p1p-1 不管,

a[2,p2]\forall a\in [2,p-2],必然存在一个和它不相等的逆元 a1[2,p2]a^{-1}\in[2,p-2],满足 aa11 (mod p)aa^{-1}\equiv 1\ (\mathrm{mod}\ p)

所以必然有 p32\displaystyle\frac{p-3}{2} 对数相乘的乘积为 11,即 (p2)!1 (mod p)(p-2)!\equiv 1\ (\mathrm{mod}\ p)

等式两边同时乘 p1p-1 就得到威尔逊定理。

例题:nNn\in\N^*n106n\leq 10^6,求 SnS_n.

Sn=k=1n(3k+6)!+13k+7(3k+6)!3k+7S_n=\displaystyle\sum_{k=1}^{n}\left\lfloor\frac{(3k+6)!+1}{3k+7}-\left\lfloor\frac{(3k+6)!}{3k+7}\right\rfloor\right\rfloor

d=3k+7d=3k+7,原式化简为

Sn=k=1n(d1)!+1d(d1)!dS_n=\displaystyle\sum_{k=1}^{n}\left\lfloor\frac{(d-1)!+1}{d}-\left\lfloor\frac{(d-1)!}{d}\right\rfloor\right\rfloor

由威尔逊定理,当 dd 为素数时,(d1)!+1d\displaystyle\frac{(d-1)!+1}{d} 必然是整数,而 (d1)!d\displaystyle\left\lfloor\frac{(d-1)!}{d}\right\rfloor 必然比 (d1)!+1d\displaystyle\frac{(d-1)!+1}{d}11,于是有:

{p 为素数(d1)!+1d(d1)!d=1p 为合数(d1)!+1d(d1)!d=0\begin{cases}p\ 为素数 & \displaystyle\left\lfloor\frac{(d-1)!+1}{d}-\left\lfloor\frac{(d-1)!}{d}\right\rfloor\right\rfloor=1 \\ p\ 为合数 & \displaystyle\left\lfloor\frac{(d-1)!+1}{d}-\left\lfloor\frac{(d-1)!}{d}\right\rfloor\right\rfloor=0\end{cases}

所以只需统计 [1,3×106+7][1,3\times 10^6+7] 中的素数即可得出答案。

拉格朗日定理

pp 是质数,则对于模 pp 意义下的 nn 次整系数多项式 f(x)=anxn+an1xn1++a0(pan)f(x)=a_nx^n+a_{n-1}x^{n-1}+\dots+a_0(p\nmid a_n) 同余方程 f(x)0(mod p)f(x)\equiv 0(\mathrm{mod}\ p) 至多有 nn 个不同的解。