Project Euler 361(2)

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

あれからいろいろ考えて、なんとかAnの計算量をO(log n)にしたかったのだが、結局断念した。
それでもだいぶ速くなった。

1018まで 0.06s
1050まで 0.65s
10100まで 0.41s
10200まで3s
10400まで23s