OFFSET
2,13
COMMENTS
Repeated prime factors are counted.
LINKS
Seppo Mustonen, Table of n, a(n) for n = 2..250
S. Mustonen, On prime factors of numbers m^n+-1
Seppo Mustonen, On prime factors of numbers m^n+-1 [Local copy]
EXAMPLE
For n=14, 2^n+1=16385=5*29*113 has two prime factors of form, namely 29=2n+1, 113=8n+1. Thus a(14)=2.
MATHEMATICA
m = 2; n = 2; nmax = 250;
While[n <= nmax, {l = FactorInteger[m^n + 1]; s = 0;
For[i = 1, i <= Length[l],
i++, {p = l[[i, 1]];
If[IntegerQ[(p - 1)/n] == True, s = s + l[[i, 2]]]; }];
a[n] = s; } n++; ];
Table[a[n], {n, 2, nmax}]
Table[{p, e}=Transpose[FactorInteger[2^n+1]]; Sum[If[Mod[p[[i]], n] == 1, e[[i]], 0], {i, Length[p]}], {n, 2, 50}]
CROSSREFS
KEYWORD
nonn
AUTHOR
Seppo Mustonen, Nov 24 2010
STATUS
approved