OFFSET
0,5
FORMULA
a(n) = Product_{d | n} (d + 1), where d + 1 is prime and 1 < d < n.
MAPLE
clausen := proc(n) numtheory[divisors](n) minus {1, n};
map(i -> i+1, %); select(isprime, %); mul(i, i=%) end:
seq(clausen(n), n = 0..80);
MATHEMATICA
a[n_] := Product[If[1 < d < n && PrimeQ[d + 1], d + 1, 1], {d, Divisors[n]}]; Array[a, 100, 0] (* Amiram Eldar, Aug 20 2022 *)
PROG
(PARI) a(n) = if (n, vecprod(select(isprime, apply(x->x+1, setminus(divisors(n), [1, n])))), 1); \\ Michel Marcus, Aug 21 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Aug 20 2022
STATUS
approved