

A286516


a(n) = b(2*n1)/b(2*n) where b(n) = A195441(n1) = denominator(Bernoulli_{n}(x)  Bernoulli_{n}).


9



1, 2, 3, 2, 5, 3, 7, 2, 3, 5, 11, 1, 13, 7, 15, 2, 17, 3, 19, 5, 7, 11, 23, 1, 5, 13, 3, 7, 29, 5, 31, 2, 11, 17, 7, 1, 37, 19, 13, 5, 41, 21, 43, 11, 3, 23, 47, 1, 7, 5, 17, 13, 53, 3, 11, 7, 19, 29, 59, 1, 61, 31, 7, 2, 65, 11, 67, 17, 23, 5, 71, 1, 73, 37
(list;
graph;
refs;
listen;
history;
text;
internal format)



OFFSET

1,2


COMMENTS

a(n) is an integer for all n, a(n) is odd if n is not a power of 2, a(2^k)=2 for all k>=1, a(n)=1 infinitely often, and a(n)=p infinitely often for every prime p. See Cor. 2 and Cor. 3 in "The denominators of power sums of arithmetic progressions". See also "Powersum denominators".


LINKS



FORMULA



MATHEMATICA

b[n_] := Denominator[ Together[ BernoulliB[n, x]  BernoulliB[n]]]; Table[
b[2 n  1]/b[2 n], {n, 1, 74}]


CROSSREFS



KEYWORD

nonn


AUTHOR



STATUS

approved



