π \pi π有一种简单的积分求法:
π 4 = a r c t a n ( 1 ) = ∫ 0 1 1 1 + x 2 d x = ∫ 0 1 1 − x 2 + x 4 − . . . d x = 1 − 1 3 + 1 5 − 1 7 + . . . \frac{\pi}{4}=arctan(1)=\int_0^1\frac{1}{1+x^2}dx =\int_0^1 1-x^2+x^4-...dx=1-\frac{1}{3}+\frac{1}{5}-\frac{1}{7}+... 4π=arctan(1)=∫011+x21dx=∫011−x2+x4−...dx=1−31+51−71+...
现在我们有另一个问题:求满足 x 2 + y 2 = r x^2+y^2=r x2+y2=r的 ( x , y ) (x,y) (x,y)个数,即半径为 r \sqrt r r的圆上格点数,
我们要借助这个问题来推导 π \pi π
任意一个数可以分解为若干个素数,且分解方案唯一(除了乘 − 1 -1 −1),
相似地,素数又可以分解为高斯素数,且分解方案唯一(除了乘 − 1 , i , − i -1,i,-i −1,i,−i)
2 2 2比较特殊,等于 ( 1 + i ) ( 1 − i ) (1+i)(1-i) (1+i)(1−i)
对于其他的素数,只有 4 k + 1 4k+1 4k+1都能被分解, 4 k + 3 4k+3 4k+3都不能( x 2 x^2%4 x2只有 0 , 1 0,1 0,1两种取值)
考虑 5 5 5这个素数,它等于 ( 2 + i ) ( 2 − i ) (2+i)(2-i) (2+i)(2−i),那么 25 = ( 2 + i ) ( 2 − i ) ( 2 + i ) ( 2 − i ) 25=(2+i)(2-i)(2+i)(2-i) 25=(2+i)(2−i)(2+i)(2−i)
我们要把它分成两组,共有三种分法:
( 2 + i ) ( 2 − i ) , ( 2 + i ) ( 2 − i ) (2+i)(2-i),(2+i)(2-i) (2+i)(2−i),(2+i)(2−i);
( 2 + i ) ( 2 + i ) , ( 2 − i ) ( 2 − i ) (2+i)(2+i),(2-i)(2-i) (2+i)(2+i),(2−i)(2−i);
( 2 − i ) ( 2 − i ) , ( 2 + i ) ( 2 + i ) (2-i)(2-i),(2+i)(2+i) (2−i)(2−i),(2+i)(2+i);
这其实就是 25 25 25中 5 5 5的幂次 + 1 +1 +1
对于那些不能分解的素数,如果是奇数个,那分配到两边后,两边的积一定不是共轭复数
本来如果是分解方法的话这是算的,但我们的问题是圆上格点数,那就不能算
现在考虑 2 2 2,如果我们把 ( 1 + i ) (1+i) (1+i)分在左边, ( 1 − i ) (1-i) (1−i)分在右边,
交换后左边相当于除了 i i i,右边乘了 − i -i −i,可以省略
假如我们有一个数 450 450 450,它等于 2 ∗ 3 2 ∗ 5 2 2*3^2*5^2 2∗32∗52
2 2 2对答案不造成贡献, 3 3 3也不对答案造成贡献, 5 5 5对答案造成的贡献是 2 + 1 = 3 2+1=3 2+1=3
最后不要忘了乘以 4 4 4
现在我们定义一个函数 χ ( n ) \chi(n) χ(n),
χ ( n ) = { 1 , n = 4 k + 1 − 1 , n = 4 k + 3 0 , n = 2 k \chi(n)= \begin{cases} 1,&n=4k+1\\ -1,&n=4k+3\\ 0,&n=2k \end{cases} χ(n)=⎩⎪⎨⎪⎧1,−1,0,n=4k+1n=4k+3n=2k
而且, χ ( n ) \chi(n) χ(n)是完全积性函数,这样,我们就可以把上面的东西用 χ ( n ) \chi(n) χ(n)来表示:
( χ ( 2 0 ) + χ ( 2 1 ) ) ( χ ( 3 0 ) + χ ( 3 1 ) + χ ( 3 2 ) ) ( χ ( 5 0 ) + χ ( 5 1 ) + χ ( 5 2 ) ) (\chi(2^0)+\chi(2^1))(\chi(3^0)+\chi(3^1)+\chi(3^2))(\chi(5^0)+\chi(5^1)+\chi(5^2)) (χ(20)+χ(21))(χ(30)+χ(31)+χ(32))(χ(50)+χ(51)+χ(52))
由于除了 χ ( 2 0 ) \chi(2^0) χ(20)外, χ ( 2 k ) = 0 \chi(2^k)=0 χ(2k)=0,所以第一部分 = 1 =1 =1
由于 ( 4 k + 3 ) 2 p (4k+3)^{2p}%4=1 (4k+3)2p, ( 4 k + 3 ) 2 p + 1 (4k+3)^{2p+1}%4=3 (4k+3)2p+1,所以第二部分= 1 − 1 + 1 − 1 + . . . 1-1+1-1+... 1−1+1−1+...,满足偶数为 1 1 1,奇数为 − 1 -1 −1
由于 ( 4 k + 1 ) p (4k+1)^p%4=1 (4k+1)p,所以第三部分= 1 + 1 + 1 + . . . = 1+1+1+...= 1+1+1+...=次数 + 1 +1 +1
由此,我们发现:半径为 n n n的圆上格点数 = 4 ∑ d ∣ n χ ( d ) =4\sum_{d|n}\chi(d) =4∑d∣nχ(d),即每个 χ ( d ) \chi(d) χ(d)都会对其倍数造成贡献
那么, lim r → + ∞ = 4 ∑ r 2 i χ ( i ) \lim_{r\to+\infty}=4\sum\frac{r^2}{i}\chi(i) limr→+∞=4∑ir2χ(i),即 lim r → + ∞ π r 2 = 4 r 2 ∑ χ ( i ) i \lim_{r\to+\infty}\pi r^2=4r^2\sum\frac{\chi(i)}{i} limr→+∞πr2=4r2∑iχ(i)
π 4 = 1 − 1 3 + 1 5 − 1 7 + . . . \frac{\pi}{4}=1-\frac{1}{3}+\frac{1}{5}-\frac{1}{7}+... 4π=1−31+51−71+...