Project Euler 428

http://projecteuler.net/index.php?section=problems&id=428

25着。
とにかく数式変形が大変で。Mathematica的なものがあれば楽だったんだろうけど。昔自作したんだけど、途中で挫折したんだっけ?
まだPyPyで40分以上かかっているのでなんとかしたい。帰りの電車の中で全然違う方法を考えていたけど、なんとかなるような気も。