OFFSET
2,1
LINKS
Max Alekseyev, Table of n, a(n) for n = 2..47
Yasutoshi Kohmoto, UnitaryPhi, SeqFan list, Aug 31 2018.
Michel Marcus, solve_uphi pari code
FORMULA
a(pp) = pp iff pp is a prime power (A000961) > 1.
a(n) = A319313(2*n). - Max Alekseyev, Oct 20 2025
PROG
(PARI) \\ uses the "solve_uphi pari code", see links
a(n) = {my(lim = 1, v); while (1, v = solve_uphi(n-1, n, lim); if (#v, return (v[1])); lim *= 10; ); }
CROSSREFS
KEYWORD
nonn
AUTHOR
Michel Marcus, Sep 04 2018
EXTENSIONS
Terms a(30) onward from Max Alekseyev, Oct 20 2025
STATUS
approved
