479. Roots on the Rise
https://projecteuler.net/problem=479
令 a_k, b_k, c_k 为方程式 1/x = (k/x)^2 (k+x^2) - kx 的三个解 (含实数及复数)。
例如当 k = 5 时，{a_k, b_k, c_k} 的近似值为 {5.727244, -0.363622+2.057397i,
-0.363622-2.057397i}。
令 S(n) = Σ(a_k+b_k)^p (b_k+c_k)^p (c_k+a_k)^p,
其中 p, k 为整数, 满足 1≦p,k≦n。
有趣的是，S(n) 永远是整数，如 S(4) = 51160。
求 S(10^6) mod 1 000 000 007。

看完题目感想: 这什么高中程度的题目...

看起来不难，但是不是要跑很久啊...

这真的好简单啊...早点看到说不定可以抢进前几名XDi大,10^6已经是PE近来少见的小了,跑很久可能是方法不够精简喔

基本上这题只要小改真的可以出成给高中生的考题例如固定某个 p 变动 k 或固定某个 k 变动 p