OFFSET
1,1
COMMENTS
The sequence of primes indexed by this sequence is (5, 89, 113, 887, 1327, 5591, 16141, 15683, 31907, 43331, 34061, 31397, 188029, 461717, 155921, ...). See A390785 for a guide to related sequences.
FORMULA
floor(prime(a(n)+1)/6) - floor(prime(a(n))/6) = n.
MAPLE
N:= 40: # for a(1) .. a(N)
V:= Vector(N): p:= 2: m:= 0: count:= 0:
for k from 1 while count < N do
p:= nextprime(p); mp:= floor(p/6);
v:= mp - m; m:= mp;
if v > 0 and v <= N and V[v] = 0 then
V[v]:= k; count:= count+1;
fi;
od:
convert(V, list); # Robert Israel, Nov 25 2025
MATHEMATICA
p[n_] := p[n] = Prime[n]; m = 6; z = 50000;
t = Table[Floor[p[n + 1]/m] - Floor[p[n]/m], {n, 1, z}];
Flatten[Table[First[Position[t, k]], {k, 1, 40}]]
CROSSREFS
KEYWORD
nonn
AUTHOR
Clark Kimberling, Nov 24 2025
STATUS
approved
