13/02/11 07:24:09.83
>>157 Io言語
powerMod := method(a, b, m,
r := 1
for(i, 0, b log2 floor,
if(b at(i) == 1, r = r * a % m)
a = a * a % m
)
r
)
Io> powerMod(2,57885161,10**10)-1
==> 1724285951
13/02/11 07:24:09.83
>>157 Io言語
powerMod := method(a, b, m,
r := 1
for(i, 0, b log2 floor,
if(b at(i) == 1, r = r * a % m)
a = a * a % m
)
r
)
Io> powerMod(2,57885161,10**10)-1
==> 1724285951