首先有结论:半径Rsqrt{R}R的圆经过的整点数是4∑d∣Rχ(d)4sum_{d|R}chi(d)4∑d∣Rχ(d),其中χ(d)=1 (d​​mod  4=1), −1 (d​​mod  4=3), 0 (d​​mod  2=0).chi(d)=1~(d!!mod 4=1),~-1~(d !!mod 4 =3),~0~(d !!mod 2=0).χ(d)=1 (dmod4=1), −1 (dmod4=3), 0 (dmod2=0).
所以答案是
4∑i=1Ri∑d∣iχ(d)4sum_{i=1}^{R}isum_{d|i}chi(d)4i=1∑Rid∣i∑χ(d)
=4∑d=1Rd×χ(d)∑i=1⌊nd⌋i=4sum_{d=1}^{R}d imes chi(d)sum_{i=1}^{lfloor{frac{n}{d}} floor}i=4d=1∑Rd×χ(d)i=1∑⌊dn⌋i
然后就可以O(R)O(sqrt{R})O(R)计算了。