A372351 Odd bisection of A371094.
21, 21, 341, 45, 117, 69, 341, 93, 213, 117, 5461, 141, 309, 165, 725, 189, 405, 213, 1877, 237, 501, 261, 1109, 285, 597, 309, 5461, 333, 693, 357, 1493, 381, 789, 405, 3413, 429, 885, 453, 1877, 477, 981, 501, 87381, 525, 1077, 549, 2261, 573, 1173, 597, 4949, 621, 1269, 645, 2645, 669, 1365, 693, 11605, 717
Offset: 1
Keywords
Crossrefs
Programs
-
Mathematica
Table[With[{e = IntegerExponent[6*n - 2, 2]}, (6*n - 2)*2^e + (4^e - 1)/3], {n, 100}] (* Paolo Xausa, Apr 29 2024 *)
-
PARI
A371094(n) = { my(m=1+3*n, e=valuation(m,2)); ((m*(2^e)) + (((4^e)-1)/3)); }; A372351(n) = A371094(n+n-1);
-
Python
def A372351(n): return ((m:=6*n-2)<<(e:=(~m & m-1).bit_length()))+((1<<(e<<1))-1)//3 # Chai Wah Wu, Apr 28 2024
Formula
a(n) = A371094(2*n-1).