Projecteuler (288) An enormous factorial - 拼圖

Projecteuler (288) An enormous factorial

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

對於任意質數p, 函數N(p,q)被定義成 N(p,q)=Σ T_n*p^n
n=0 to q
T_n的產生是藉由下列的亂數產生器

S_0=290797
S_n+1=(S_n)^2 mod 50515093
T_n=S_n mod p

令Nfac(p,q)為N(p,q)的階乘 (譯注:意即 (N(p,q))!)

令NF(p,q)為Nfac(p,q)中含有之因數p的個數
(譯注: 即(N(p,q))!因數分解後，p的指數即為Nfac(p,q))

現已知 NF(3,10000) mod 3^20 = 624955285

求NF(61,10^7) mod 61^10 = ?

--