本人第三个点的 python 程序算 $s_4$ 如下：

n = 1000000000000000
mod = 2 ** 64
s4 = n * (n + 1) * (2 * n + 1) * (3 * (n**2) + 3 * n - 1) / 30
print("%d"%(s4 % mod))
# 0.0

有什么错误的地方吗？谢谢指出