質問します〜

Rabin法を用いた1億以下の素数の高速判定アルゴリズムが欲しいです
C言語でお願いします