OFFSET
1,1
COMMENTS
The unitary version of A100410.
No more terms below 10^7.
EXAMPLE
MATHEMATICA
uphi[n_] :=(Times @@ (Table[#[[1]]^#[[2]] - 1, {1}] & /@ FactorInteger[n] ))[[1]] ; seqQ[n_] := (n - uphi[n] - 1 == PrimePi[n]); Select[Range[2, 10^7], seqQ]
PROG
(PARI) uphi(n) = my(f=factor(n)); prod(i=1, #f~, f[i, 1]^f[1, 2]-1);
isok(k) = k - primepi(k) - 1 == uphi(k); \\ Michel Marcus, Mar 07 2018
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Amiram Eldar, Mar 06 2018
EXTENSIONS
a(18)-a(22) from Robert G. Wilson v, Mar 07 2018
STATUS
approved