Project Euler 291（4）

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

C++でMiller-Rabin法を実装してみた。除算はdoubleで概算して微調整という手抜き。
これで15分。最初の方法では1分だが、ミラー・ラビン法のコードがあれば、ほとんどコードを書くことなく答えを出すことができる。