a48[n]：=a48[n]=除数和[NestWhile[#/2&，n，EvenQ]，MoebiusMu[#]2^（n/#）&]/（2n）；(*A000048号*)
a37[n]：=a37[n]=除数和[n，Moebius Mu[n/#]2^#&]/n；(*A001037号*)
CI0[{n_Integer}]：=CI0[}]=CI[Transpose[If[EvenQ[n]，p2=IntegerExponent[n，2]；sub=除数[n/2^p2]；{2^（p2+1）sub，a48/@（2^p2-sub）}，sub=除数[n]；{sub，a37/@sub}]]2^（n-1）；（*均匀烫发*）
CI1[{n_Integer}]：=CI1[}]=CI[sub=除数[n]；转座[If[EvenQ[n]，{sub，a37/@sub}，{2sub，（a37/@sub）/2}]]2^（n-1）；（*奇数perm.*）
compress[x:{{_，_}…}]：=（s=Sort[x]；对于[i=Length[s]，i>1，i-=1，如果[s[[i，1]]==s[[i-1，1]]，s[i-1，2]]+=s[i，2]]；s=Delete[s，i]，Null]]；s）
cix[{a_，b}，{c_，d_}]：={LCM[a，c]，（a c d）/LCM[a，c]}；
取消保护[次数]；次数[CI[a_List]，CI[b_List]]：=（*combine*）CI[compress[Flatten[Outer[cix，a，b，1]，1]]；保护[次数]；
CI0[p_List]：=CI0[p]=展开[CI0[Drop[p，-1]]CI0[{Last[p]}]+CI1[Drop[p，-1]]CI1[{Last[p]{]]
CI1[p_List]：=CI1[p]=展开[CI0[Drop[p，-1]]CI1[{Last[p]}]+CI1[Drop[p，-1]]CI0[{Last[p]{]]
pc[p_List]：=模块[{ci，mb}，mb=删除重复项[p]；ci=计数[p，#]&/@mb；n/（Times@@（ci！）Times@@（mb^ci））]（*分区计数*）
行[n_Integer]：=行[n]=系数[（总计[（（CI0[#]+CI1[#]）pc[#]]）&/@IntegerPartitions[n]]）/（n！2^n）]/。CI[l_List]：>j^（总计[l][[2]）
表[LinearSolve[表[二项式[i，j]，{i，1，2^n}，{j，1，2 ^n}]，表[数组[n，k]，{k，1，2A^n}]]，{n，1，6}]//展平
| 