OFFSET
2,1
COMMENTS
EXAMPLE
Of the 30 members of the reduced residue system mod 31, exactly one-third, namely 2, 3, 5, 7, 11, 13, 17, 19, 23, 29, are prime. 31 is the least number with this property, so a(3) = 31.
MAPLE
f:= proc(n) uses numtheory;
phi(n)/(pi(n) - nops(factorset(n)));
end proc:
N:= 13: # to get a(2)..a(N)
R:= Array(2..N): count:= 0:
for k from 3 while count < N-1 do
v:= f(k);
if v::integer and v <= N and R[v] = 0 then
R[v]:= k;
count:= count+1;
fi
od:
convert(R, list);
MATHEMATICA
With[{s = Table[EulerPhi[n]/Count[Prime@ Range@ PrimePi@ n, _?(GCD[#, n] == 1 &)], {n, 3, 10^4}]}, Array[2 + FirstPosition[s, #][[1]] &, Max@ Select[s, IntegerQ] - 1, 2]] (* Michael De Vlieger, Apr 23 2019 *)
CROSSREFS
KEYWORD
nonn,more
AUTHOR
J. M. Bergot and Robert Israel, Apr 23 2019
STATUS
approved