%I#42 2018年1月1日04:19:07

%S 1,4,4,4,16,4,8,4,16，16,4,4，16,8,16,4，16，16,16,16,16，64,4,8，4,16,18,8，

%电话：8,32,16,32,4,16,16,16，6,64,16,32，16,64,4,4,16，8,16,4,16，16,16，16，

%U 64,8,16,8,32,32,16,16,64,32,64,4,16,16,16，16,64,16,32,16，64

%N（1+x+x^2+x^3）^N中奇数系数的个数。

%H G.C.Greubel，n表，n=0..10000时的a（n）</a>

%H S.R.Finch、P.Sebah和Z.-Q.Bai，<a href=“http://arXiv.org/abs/0802.2654“>帕斯卡三角中的奇数项，arXiv:0802.2654[math.NT]，2008。

%F a（n）=2^A036555（n）。

%F a（n）=gcd（4^n，C（4*n，n））_Peter Luschny_，2011年11月8日

%e自2015年3月1日起生效：（开始）

%e写成一个不规则三角形，其中的行长度是A011782的项，序列开始于：

%e 1；

%e 4；

%e 4,4；

%e 4、16、4、8；

%e第4、16、16、4、4、16和8、16页；

%e第4,16,16,16、6,64,4,8,4,16,16,8,8,32,16,32页；

%e第4、16、16、66、64、16、64、62、64、44、16、8、16、14、16、18、16、16,64、8、8、32、16、24、32、64、64；

%e。。。

%e（结束）

%p序列（igcd（4^n，二项式（4*n，n）），n=0..77）；#_Peter Luschny_，2011年11月8日

%t多项式模型[（1+x+x^2+x^3）^n，2]/。x->1

%t A036555=总数/@整数位数[3范围[0，100]，2]；表[2^A036555[[n]]，{n，1，20}]（*或*）表[GCD[4^n，二项式[4*n，n]]、{n，0，50}]（*_G.C.Greubel_，2017年12月31日*）

%o（PARI）a（n）=｛my（pol=pol（[1,1,1,1]，xx）*Mod（1,2））；subst（升力（pol^n），xx，1）；｝\\_Michel Marcus_，2015年3月1日

%o（PARI）a（n）=2^重量（3*n）；\\_Joerg Arndt_，2015年3月10日

%Y参见A000120、A001316、A008287、A036555、A071053。

%K nonn公司

%0、2

%A _斯蒂芬芬奇，2008年1月25日