m=20；
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]；
边[v_]：=总和[Sum[GCD[v[i]]，v[[j]]]，{j，1，i-1}]，{i，2，Length[v]}]+总和[Quotient[v[i]]，2]，{i；
oddp[v_]：=（对于[i=1，i<=长度[v]，i++，如果[BitAnd[v[i]]，1]==0，返回[0]]]；1）；
b[n_]：=b[n]=（s=0；Do[If[oddp[p]==1，s+=permcount[p]*2^edges[p]]，{p，IntegerPartitions[n]}]；s/n！）；
B[x_]=总和[B[k]x^k，{k，0，m}]；
A[x_]=2-1/B[x]；
