OFFSET
0,2
COMMENTS
a(n) = n+1 for n+1 in A192453. - Robert Israel, May 24 2019
LINKS
Robert Israel, Table of n, a(n) for n = 0..10000
MAPLE
N:= 100: # for a(0)..a(N)
V:= Array(0..N):
count:= 1: V[0]:= 1:
for k from 2 while count < N+1 do
R:= select(t -> t <= N and V[t] = 0, {seq(x^4 mod k, x=1..k/2)});
V[convert(R, list)]:= k;
count:= count + nops(R);
od:
convert(V, list); # Robert Israel, May 24 2019
PROG
(PARI) a(n) = {my(k=n+1); while (!ispower(Mod(n, k), 4), k++); k; } \\ Michel Marcus, May 25 2019
CROSSREFS
KEYWORD
easy,nonn
AUTHOR
Henry Bottomley, Jun 08 2000
STATUS
approved