OFFSET
1,3
LINKS
Robert Israel, Table of n, a(n) for n = 1..10000
MAPLE
N:= 100: # for a(1) .. a(N)
V:= Vector(N):
p:= 2:
do
p:= nextprime(p);
if p > 2*N-1 then break fi;
J:= [seq((p^i+1)/2, i = 1 .. ilog[p](2*N-1))];
V[J]:= V[J] +~ 1;
od:
ListTools:-PartialSums(convert(V, list)); # Robert Israel, Dec 26 2024
CROSSREFS
KEYWORD
nonn
AUTHOR
EXTENSIONS
Corrected and extended by Ray Chandler, May 01 2007
STATUS
approved