OFFSET
1,3
LINKS
Alois P. Heinz, Table of n, a(n) for n = 1..10000
N. J. A. Sloane, Transforms
MAPLE
a:= proc(n) option remember;
ithprime(n-1)-add(a(d), d=numtheory[divisors](n) minus {n})
end: a(1):=1:
seq(a(n), n=1..64); # Alois P. Heinz, Oct 04 2025
PROG
(PARI) a(n) = sumdiv(n, d, if (d==1, 1, prime(d-1))*moebius(n/d)); \\ Michel Marcus, Nov 04 2018
CROSSREFS
KEYWORD
nonn,look
AUTHOR
EXTENSIONS
Offset 1 and more terms from Michel Marcus, Nov 04 2018
STATUS
approved
