login
a(n) = (2^n-1-numerator(B(2^(n+1))))/2^n where numerator(B(k)) denotes the numerator of the k-th Bernoulli number.
2

%I #10 Apr 08 2025 19:01:01

%S 1,1,453,481832565077,3336994692120829058949545155598207748189,

%T 4183667308477313795108662587270037418426429056133620947343416969739965511482242658521291610836247844507748569589345

%N a(n) = (2^n-1-numerator(B(2^(n+1))))/2^n where numerator(B(k)) denotes the numerator of the k-th Bernoulli number.

%C Appears to always be an integer. Next term is too large to include.

%F a(n) = (2^n-1-A000367(2^n))/2^n

%t Table[(2^n-1-Numerator[BernoulliB[2^(n+1)]])/2^n,{n,6}] (* _Harvey P. Dale_, Nov 16 2024 *)

%o (PARI) a(n)=(2^n-1-numerator(bernfrac((2^(n+1)))))/2^n

%Y Cf. A000367.

%K nonn

%O 1,3

%A _Benoit Cloitre_, Oct 19 2002