login
The OEIS Foundation is supported by donations from users of the OEIS and by a grant from the Simons Foundation.

 

Logo


Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
A100307 Modulo 2 binomial transform of 3^n. 7

%I

%S 1,4,10,40,82,328,820,3280,6562,26248,65620,262480,538084,2152336,

%T 5380840,21523360,43046722,172186888,430467220,1721868880,3529831204,

%U 14119324816,35298312040,141193248160,282472589764,1129890359056

%N Modulo 2 binomial transform of 3^n.

%C 3^n may be retrieved through 3^n = Sum_{k=0..n}(-1)^A010060(n-k)*mod(binomial(n,k),2)*a(k).

%H Gheorghe Coserea, <a href="/A100307/b100307.txt">Table of n, a(n) for n = 0..200</a>

%H V. Shevelev, <a href="http://arxiv.org/abs/1011.6083">On Stephan's conjectures concerning Pascal triangle modulo 2 and their polynomial generalization</a>, J. of Algebra Number Theory: Advances and Appl., 7 (2012), no.1, 11-29.

%F a(n) = Sum_{k=0..n} mod(binomial(n, k), 2)*3^k.

%F From _Vladimir Shevelev_, Dec 26-27 2013: (Start)

%F Sum_{n>=0}1/a(n)^r = Product_{k>=0}(1 + 1/(3^(2^k)+1)^r),

%F Sum_{n>=0}(-1)^A000120(n)/a(n)^r = Product_{k>=0}(1 - 1/(3^(2^k)+1)^r), where r>0 is a real number.

%F In particular,

%F Sum_{n>=0}1/a(n) = Product_{k>=0}(1 + 1/(3^(2^k)+1)) = 1.391980...;

%F Sum_{n>=0}(-1)^A000120(n)/a(n) = 2/3.

%F a(2^n) = 3^(2^n)+1, n>=0.

%F Note that analogs of Stephan's limit formulas (see Shevelev link) reduce to the relations:

%F a(2^t*n+2^(t-1)) = 8*(3^(2^(t-1)+1))/(3^(2^(t-1))-1) * a(2^t*n+2^(t-1)-2), t>=2.

%F In particular, for t=2,3,4, we have the following formulas:

%F a(4*n+2) = 10 * a(4*n);

%F a(8*n+4) = 41/5 * a(8*n+2);

%F a(16*n+8)= 3281/410 * a(16*n+6), etc. (End)

%F From _Tom Edgar_, Oct 11 2015: (Start)

%F a(n) = Product_{b_j != 0} a(2^j) where n = Sum_{j>=0} b_j*2^j is the binary representation of n.

%F a(2*k+1) = 4*a(2*k).

%F (End)

%t Table[Sum[Mod[Binomial[n,k],2]3^k,{k,0,n}],{n,0,30}] (* _Harvey P. Dale_, Aug 28 2013 *)

%o (Sage) [sum((binomial(n,k)%2)*3^k for k in [0..n]) for n in [0..100]] # _Tom Edgar_, Oct 11 2015

%o (PARI) a(n) = subst(lift((Mod(1,2)+'x)^n), 'x, 3); \\ _Gheorghe Coserea_, Jun 11 2016

%Y Cf. A001316, A001317, A038183, A100308, A100309, A100310, A100311.

%K easy,nonn

%O 0,2

%A _Paul Barry_, Dec 06 2004

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recent
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 14 19:09 EDT 2020. Contains 336483 sequences. (Running on oeis4.)