login
A356655
Clausen numbers based on the strictly proper divisors of n, 1 < d < n.
0
1, 1, 1, 1, 3, 1, 3, 1, 15, 1, 3, 1, 105, 1, 3, 1, 15, 1, 21, 1, 165, 1, 3, 1, 1365, 1, 3, 1, 15, 1, 231, 1, 255, 1, 3, 1, 25935, 1, 3, 1, 165, 1, 21, 1, 345, 1, 3, 1, 23205, 1, 33, 1, 15, 1, 399, 1, 435, 1, 3, 1, 465465, 1, 3, 1, 255, 1, 483, 1, 15, 1, 33, 1
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
Sequence in context: A212183 A115716 A079412 * A306861 A262940 A278601
KEYWORD
nonn
AUTHOR
Peter Luschny, Aug 20 2022
STATUS
approved