OFFSET
1,2
FORMULA
PROG
(PARI) pbc(n) = if (n==1, 1, prime(n-1));
a(n) = sumdiv(n, d, pbc(d)*eulerphi(n/d)); \\ Michel Marcus, Feb 28 2019
CROSSREFS
KEYWORD
nonn
AUTHOR
STATUS
approved