m=1；（*颜色元素的尺寸，此处为边*）
Fi1[p1_]：=模[{g，h}，系数[积[g=GCD[k1，p1]；h=GCD[2 k1，p1]；（1+2x^（k1/g））^（r1[[k1]]g）如果[可除[k1，h]，1，（1+2x^（2 k1/h））^（r2[[k1_]h/2）]，{k1，展平[位置[cs，n1_/；n1>0]}]，x，n-m]]；
FiSum[]：=（Do[Fi2[k2]=Fi1[k2]，{k2，Divisors[per]}]；DivisorSum[per，Divisor Sum[d1=#，MoebiusMu[d1/#]Fi2[#]&]/#&]）；
CCPol[r_List]：=（r1=r；r2=cs-r1；per=LCM@@表[If[cs[j2]==r1[j2]]，If[0==cs[[j2]，1，j2]、2j2]和{j2，n}]；If[EvenQ[Sum[If[CevenQ[j3]，r1[[j3]]，r2[[j3]]，{j3，n}]]，1，-1]倍@@二项式[cs，r1]2^（n-总计[cs]）b^FiSum[]）；
PartPol[p_List]：=（cs=计数[p，#]&/@Range[n]；总计[CPol[#]&@Tuples[Range[0，cs]]]）；
pc[p_List]：=模块[{ci，mb}，mb=删除重复项[p]；ci=计数[p，#]&&@mb；不/（Times@@（ci！）Times@@（mb^ci））]（*分区计数*）
行[n_Integer]：=行[n]=系数[（总计[（PartPol[#]pc[#]）和/@IntegerPartitions[n]]）/（n！2^n）]
数组[n，k_]：=行[n]/。b->k
表[数组[n，d+m-n]，{d，7}，{n，m，d+m-1}]//扁平
