login
A181904
a(n) = 2*(4^n - 1) / A027760(n)
1
3, 5, 63, 17, 1023, 195, 16383, 4369, 262143, 31775, 4194303, 12291, 67108863, 89478485, 1073741823, 16843009, 17179869183, 172229265, 274877906943, 6663706835, 4398046511103, 254959218035, 70368744177663, 206208774147, 1125899906842623, 1501199875790165
OFFSET
1,1
COMMENTS
The Clausen numbers A027760 divide 2*(4^n - 1).
LINKS
Thomas Clausen, Lehrsatz aus einer Abhandlung ueber die Bernoulli'schen Zahlen, Astr. Nachr. 17 (1840), 351-352.
Ludwig Seidel, Über eine einfache Entstehungsweise der Bernoullischen Zahlen und einiger verwandten Reihen, Sitzungsberichte der mathematisch-physikalischen Classe der königlich bayerischen Akademie der Wissenschaften zu München, Vol. 7 (1877), 157-187.
MATHEMATICA
a[n_] := 2*(4^n-1)/Times @@ Select[Divisors[n]+1, PrimeQ]; Table[a[n], {n, 1, 26}] (* Jean-François Alcover, Jul 29 2013 *)
PROG
(Sage)
def A181904(n) :
def A027760(n) :
S = divisors(n)
S = map(lambda i: i+1, S)
S = filter(lambda i: is_prime(i), S)
return mul(i for i in S)
return 2*(4^n-1) / A027760(n)
[A181904(n) for n in (1..26)]
CROSSREFS
KEYWORD
nonn
AUTHOR
Peter Luschny, Apr 01 2012
STATUS
approved