%I M4010 N1662#51 2024年7月5日16:12:42

%S 1,5,5215221459844846465614145444229148550030864，

%电话：33333 107860769639668174695272746749919580928，

%电话：333017109925390903792693181442327872412934945330155633255521860084295633556503802059558812644803074048576025077880237087769467308716151571214336

%N N个未标记点上二进制关系数的一半。

%D M.D.McIlroy，有限集上关系结构数的计算，麻省理工学院电子研究实验室，季度进展报告，第17期，1955年9月15日，第14-22页。

%D N.J.A.斯隆，《整数序列手册》，学术出版社，1973年（包括该序列）。

%D N.J.A.斯隆和Simon Plouffe，《整数序列百科全书》，学术出版社，1995年（包括该序列）。

%H Chai Wah Wu，n的表格，n=1..59的a（n）</a>

%H R.L.Davis，<a href=“https://dx.doi.org/10.1090/S002-9939-1953-005294-2“>有限关系的结构数，Proc.Amer.Math.Soc.4（1953），486-495。

%H M.D.McIlroy，有限集上关系结构数的计算，麻省理工学院电子研究实验室，季报，第17期，1955年9月15日，第14-22页。[带注释的扫描副本]

%H W.Oberschelp，<a href=“http://gdz.sub.uni-goettingen.de/dms/resolveppn/？PPN=GDZPPN002298732“>Kombinatorische Anzahlbestimmungen in Relationen，《数学年鉴》，174（1967），53-78。

%F a（n）=A000595（n）/2.-_Sean A.Irvine_，2012年3月16日

%t permcount[v_]：=模[{m=1，s=0，k=0，t}，对于[i=1，i<=长度[v]，i++，t=v[i]]；k=如果[i>1&&t==v[[i-1]]，k+1，1]；m*=t*k；s+=t]；s/m] ；

%t边[v_]：=总和[2 GCD[v[i]]，v[[j]]]，{i，2，长度[v]}，{j，1，i-1}]+总和[v]；

%t a[n_]：=模块[{s=0}，Do[s+=permcount[p]*2^edges[p]，{p，IntegerPartitions[n]}]；s/（2 n！）]；

%t数组[a，12]（*_Jean-François Alcover_，2019年8月1日，在A000595*中的_Andrew Howroyd_之后）

%o（Python）

%o来自itertools导入产品

%o来自数学导入prod、阶乘、gcd

%o从分数导入分数

%o从sympy.utilities.iterables导入分区

%o def A001173（n）：对于分区（n）中的p，返回int（sum（分数（1<<sum（p[r]*p[s]*gcd（r，s）表示乘积中的r，s（p.keys（），repeat=2）），prod（q**p[q]*阶乘（p[q]）表示q in p））>>1#_Chai Wah Wu_，2024年7月2日

%Y参考A000595、A001174。

%K nonn很好

%O 1,2号机组

%A _N.J.A.斯隆_

%E更多条款摘自_Vladeta Jovovic_，2000年4月18日

%E a（13）-a（14）（基于A000595），来自Pontus von Brömssen，2022年8月4日