OFFSET
1,3
COMMENTS
As in A166062, the offset is rather arbitrary.
The sequence contains numbers like 210 which are not in A006954.
One could also consider dividing by the largest prime divisor of A027642 instead of A089026, which yields 1, 1, 2, 1, 6, 1, 6, 1, 6, 1, 6, 1, 210, 1, 2, 1, 30, 1, 42, 1, 30, ... as an alternative version.
These are the Clausen numbers based on the proper divisors of n whereas the classical Clausen numbers A160014 are based on all divisors of n. (The proper divisors are the divisors of n that are less than n.) - Peter Luschny, Aug 20 2022
MAPLE
A027642 := proc(n) denom(bernoulli(n)) ; end:
A089026 := proc(n) if isprime(n) then n; else 1; end if; end proc:
A166120 := proc(n) A027642(n-1)/A089026(n) ; end proc: seq(A166120(n), n=1..80) ; # R. J. Mathar, Mar 25 2010
# Second program, assuming offset 0:
clausen := proc(n) if irem(n, 2)=1 then 1 else numtheory[divisors](n) minus {n};
map(i -> i+1, %); select(isprime, %); mul(i, i=%) fi end:
seq(clausen(n), n = 0..79); # Peter Luschny, Aug 20 2022
CROSSREFS
KEYWORD
nonn
AUTHOR
Paul Curtz, Oct 07 2009
EXTENSIONS
Extended by R. J. Mathar, Mar 25 2010
STATUS
approved