OFFSET
0,3
COMMENTS
The Stern diatomic sequence A002487 could be called the Stern-Fibonacci sequence, since it is given by A002487(n) = Sum_{k=0..floor(n/2)} (binomial(n-k,k) mod 2), where F(n+1) = Sum_{k=0..floor(n/2)} binomial(n-k,k). Now a(n) = Sum_{k=0..floor(n/2)} (binomial(n-k,k) mod 2)*2^k, where J(n+1) = Sum_{k=0..floor(n/2)} binomial(n-k,k)*2^k, with J(n) = A001045(n), the Jacobsthal numbers. - Paul Barry, Sep 16 2015
These numbers seem to encode Stern (0, 1)-polynomials in their binary expansion. See Dilcher & Ericksen paper, especially Table 1 on page 79, page 5 in PDF. See A125184 (A260443) for another kind of Stern-polynomials, and also A177219 for a reference to maybe a third kind. - Antti Karttunen, Nov 01 2016
LINKS
Ivan Panchenko, Table of n, a(n) for n = 0..1000
K. Dilcher and L. Ericksen, Reducibility and irreducibility of Stern (0, 1)-polynomials, Communications in Mathematics, Volume 22/2014 , pp. 77-102.
FORMULA
a(n) = Sum_{k=0..floor(n/2)} (binomial(n-k, k) mod 2)*2^k.
a(n) = Sum_{k=0..n} (binomial(k, n-k) mod 2)*2^(n-k). - Paul Barry, May 10 2005
a(n) = Sum_{k=0..n} A106344(n,k)*2^(n-k). - Philippe Deléham, Dec 18 2008
a(0)=1, a(1)=1, a(n) = a(n-1) XOR (a(n-2)*2), where XOR is the bitwise exclusive-OR operator. - Alex Ratushnyak, Apr 14 2012
PROG
(Python)
prpr = 1
prev = 1
print("1, 1", end=", ")
for i in range(99):
current = (prev)^(prpr*2)
print(current, end=", ")
prpr = prev
prev = current
# Alex Ratushnyak, Apr 14 2012
(Python)
def A101624(n): return sum(int(not k & ~(n-k))*2**k for k in range(n//2+1)) # Chai Wah Wu, Jun 20 2022
(Haskell)
a101624 = sum . zipWith (*) a000079_list . map (flip mod 2) . a011973_row
-- Reinhard Zumkeller, Jul 14 2015
KEYWORD
easy,nonn
AUTHOR
Paul Barry, Dec 10 2004
STATUS
approved