%I#36 2024年4月15日06:45:56

%S 0,1,1,1,2,3,6,7,16,21,43,631292046851343238546258492，

%电话：16409307355929011253021718241562080307615454632990968，

%电话：5778267112014722170268642140890818307441591394983095908830293360293571311747793332290915478

%N恰好在N层之后重复的Barlow填料数量。

%H N.J.A.Sloane，N表，N=1..200的A（N）</a>

%H Dennis S.Bernstein和Omran Kouba，<a href=“https://arxiv.org/abs/1901.10703“>计算三种颜色的彩色项链和手镯</a>，arXiv:1901.10703[math.CO]，2019年。

%H E.Estevez-Rams、C.Azanza-Ricardo、J.Martinez-Garcia和B.Argon-Fernandez，<a href=“https://doi.org/10.107/S0108767304034294“>关于表示封闭堆积码序列的二进制码代数，Acta Cryst.A61（2005），201-208。

%H E.Esteves-Rams，C.L.Azana Ricardo，B.Aragon Fernandez，<a href=“https://doi.org/10.1524/zkri.220.7.592.67101“>计算封闭多型数的另一种表达式，Z.Krist.220（2005）592-595，表1

%H T.J.McLarnan，<a href=“http://dx.doi.org/10.1524/zkri.1981.155.3-4.269“>封闭包装和相关结构中的多型数</a>，Zeits.Krist.155，269-291（1981）。

%F a（n）=A011946（n/4）+A011947_安德烈·扎博洛茨基，2024年2月14日

%p与（数字理论）；读取转换；M： =200；

%pA:=proc（N，d），如果d mod 3=0，则2^（N/d）else（1/3）*（2^，N/d）+2*cos（Pi*N/d））；fi；结束；

%p E：=proc（N）如果N mod 2=0，则N*2^（N/2）+add（did（N/2，d）*phi（2*d）*2^（N/（2*d）），d=1.N/2）else（N/3）*（2^（（N+1）/2）+2*cos（Pi*（N+1）/2））；fi；结束；

%p PP:=程序（N）（1/（4*N））*（添加（did（N，d）*phi（d）*A（N，d），d=1..N）+E（N））；结束；

%对于从1到M的N，p做t1[N]：=PP（N）；日期：

%p p：=proc（N）局部s，d；s： =0；对于从1到N的d，如果N mod d=0，则s：=s+mobius（N/d）*t1[d]；fi；od：秒；结束；对于从1到M的N，进行lprint（N，P（N））；日期：#_N.J.A.Sloane，2006年8月10日

%t M=40；

%t did[m，n_]：=如果[Mod[m，n]==0，1，0]；

%tA[n_，d_]：=如果[Mod[d，3]==0，2^（n/d），（1/3）（2^；

%t EE[n_]：=如果[Mod[n，2]==0，n2^（n/2）+总和[did[n/2，d]EulerPhi[2d]*2^；

%t PP[n]：=PP[n]=（1/（4n））（Sum[did[n，d]EulerPhi[d]A[n，d]，{d，1，n}]+EE[n]）；

%tP[n_]：=模[{s=0，d}，对于[d=1，d<=n，d++，如果[Mod[n，d]==0，s+=MoebiusMu[n/d]PP[d]]；s] ；

%t数组[P，M]（*_Jean-François Alcover_，2020年4月21日，摘自Maple*）

%Y参考A114438。

%Y参见A011946、A011947、A011948、A011949、A011950、A011951、A011952、A011953、A011954、A011955、A011956。

%K nonn，简单

%O 1,6型

%A _N.J.A.Sloane和Michael OKeeffe（莫基夫（AT）asu.edu）

%E更多条款，来自N.J.A.Sloane，2006年8月10日