As hinted by
Abigail-II, exponentiating mod a number is better accomplished by keeping the numbers reduced at every step. The straightforward algorithm would be
# modexp ( base, exp, n ) = base**exp mod n
sub modexp {
my ($base, $exp, $n) = @_;
die "Negative exponent" if $exp<0;
my $res = 1;
while ($exp>0) {
$res = ( $res * $base ) % $n;
$exp--;
}
return $res;
}
This has the disadvantage of being of linear complexity in the exponent. An algorithm linear in log($exp) is as follows:
sub binmodexp {
my ($base, $exp, $n) = @_;
die "Negative exponent" if $exp<0;
my $res = 1;
my $mul = $base % $n;
while ($exp) {
if ($exp & 1) {
$res = ( $res * $mul ) % $n;
}
$exp>>=1;
$mul = ($mul*$mul) % $n;
}
return $res;
}
It is based on the simple (but powerful) observation that
if you have to take, say, the 37-th power of a number a, as 37 = 2^5 + 2^2 + 1, it is enough to multiply the factors
a,
a^(2^2) = (a^2) ^ 2
and
a^(2^5) = ((((a^2) ^ 2) ^ 2) ^ 2) ^ 2
Best regards
Antonio Bellezza
The stupider the astronaut, the easier it is to win the trip to Vega - A. Tucket
Posts are HTML formatted. Put <p> </p> tags around your paragraphs. Put <code> </code> tags around your code and data!
Titles consisting of a single word are discouraged, and in most cases are disallowed outright.
Read Where should I post X? if you're not absolutely sure you're posting in the right place.
Please read these before you post! —
Posts may use any of the Perl Monks Approved HTML tags:
- a, abbr, b, big, blockquote, br, caption, center, col, colgroup, dd, del, details, div, dl, dt, em, font, h1, h2, h3, h4, h5, h6, hr, i, ins, li, ol, p, pre, readmore, small, span, spoiler, strike, strong, sub, summary, sup, table, tbody, td, tfoot, th, thead, tr, tt, u, ul, wbr
You may need to use entities for some characters, as follows. (Exception: Within code tags, you can put the characters literally.)
| |
For: |
|
Use: |
| & | | & |
| < | | < |
| > | | > |
| [ | | [ |
| ] | | ] |
Link using PerlMonks shortcuts! What shortcuts can I use for linking?
See Writeup Formatting Tips and other pages linked from there for more info.