login
A076178
a(n) = 2*n^2 - A077071(n).
3
0, 0, 0, 2, 2, 4, 6, 10, 10, 12, 14, 18, 20, 24, 28, 34, 34, 36, 38, 42, 44, 48, 52, 58, 60, 64, 68, 74, 78, 84, 90, 98, 98, 100, 102, 106, 108, 112, 116, 122, 124, 128, 132, 138, 142, 148, 154, 162, 164, 168, 172, 178, 182, 188, 194, 202, 206, 212, 218, 226, 232
OFFSET
0,4
FORMULA
a(n) = A001105(n) - A077071(n). - Omar E. Pol, Nov 13 2024
PROG
(PARI) a(n)=2*n^2-sum(k=0, n, -valuation(polcoeff(pollegendre(2*n), 2*k), 2))
(Python)
def A076178(n): return ((n+1)*n.bit_count()-n<<1)+sum((m:=1<<j)*((k:=n>>j)-(r if n<<1>=m*(r:=k<<1|1) else 0)) for j in range(1, n.bit_length()+1)) # Chai Wah Wu, Nov 12 2024
CROSSREFS
Equals 2 * A078903(n).
Sequence in context: A166289 A071058 A376005 * A347693 A308907 A285908
KEYWORD
nonn
AUTHOR
Benoit Cloitre, Nov 01 2002
STATUS
approved