在平面直角坐标系中,画一个边长为R的正方形:
用 表示该正方形内格点(lattice)的数量,再用 表示该正方形内素格点(横纵坐标为互素的整数)的数量,求当 时素格点密度 的极限。下面我们将求解这个有趣的谜题。
很明显,这是一个数论问题,而且还是一个数格点问题。根据图像的性质,易得 ,而 的求法就比较特殊了。现在我们画一条对角线:
可以发现只要能求出橙色线下方素格点的数量,就能得出 的表达式。
由图可知,若蓝线的横坐标为n,则蓝线上素格点的数量就是比n小且与n互素数的数量,即 。因此 ,所以原问题就变成了求解如下极限:
对于分子,我们可以利用欧拉函数的狄利克雷卷积性质 ,得到:
于是代入回原来的表达式,得:
事实上,我们可以利用Dirichlet级数的乘法来转换问题。根据莫比乌斯反演,可知:
所以 。综上所述当正方形的边长不断增大时,素格点的密度会逐渐向 靠拢。
又一个意想不到的圆周率
估算一个上界。思路是每一轮都寻求一条最短线段,将当前包含天使的多边形,按面积等分成两个新的子多边形。再假设天使的运气足够好,每次都瞬移到等分效率较低的子多边形。
直观看出,取平行于正三角形一条边的线段来等分其面积,等分效率最高。令此线段长度 ,三角形边长 ,则:
这样,初始正三角形被分成一个新的小正三角形和一个等腰梯形,易见等腰梯形的等分效率远高于新的小正三角形,于是根据假设,天使将瞬移到新的小正三角形当中。如此循环,至于无穷,天使将被锁定在初始正三角形的一个顶点。计算魔鬼走过的耗时路程:
记魔鬼速度 ,则捉住天使的时间:
这个题目如此离散,不借助于数值离散优化不易得到全局最优解,建议大家来改进这个上界吧。
按照 @yyx 说的圆弧线等分正三角形以及后续的扇形,上界可以改进为: