Re: Small prime divisors of very large numbers
- --- In firstname.lastname@example.org, "David Broadhurst" <d.broadhurst@...> wrote:
>Feature number one, where some primes stay forever as divisors, is a natural consequence. Specifically, if a factor appears twice in a row at two successive powerings, it will remain from then on. Why some prime divisors appear only once while others flit in and out and others never make an appearance, those are details that could be worked out I'm sure.
> --- In email@example.com,
> "Mark Underwood" <mark.underwood@> wrote:
> > Here's a hint: (n-1) is a factor of (n^n-n).
> The stakes have risen, Mark.
> You have to explain 3 features of the data that
> Richard Heylen has exposed:
> 1) some primes stay as divisors at successive higher powerings;
> 2) some disappear for ever after a single division;
> 3) some pop in, pop out, then pob back, then stay put.
> To how many of those features of the data does your
> gnomic "hint" apply? If not to all 3, then why not?
> Thanks, in any case, for your interest
- Dear all,
I come back to this topic, looking for a working "pow_mod".
but 2^2^2 % 5 = 16 % 5 = 1.
(and idem for 2^2^2^2 - of course.).
David gave some other code, in
Below is a Pari-GP procedure "pmod(a,m)" to compute
a^(a^(a^ ... ^(a[k-1]^a[k]) ... ) modulo m
where the modulus "m" need not be prime.
Although this may work for prime m
(at least, pmod([2,2,2],5) = Mod(1,5) as should),
= Mod(1, 4)
which is most certainly wrong.
So, to put it short, has anyone a working pmod() in his "library" ?
Thanks in advance!
--- In firstname.lastname@example.org, "David Broadhurst"
I consider that his code should not be trusted.