f（n）={my（f=系数（n））；总和（k=1，#f~，my（i=primepi（f[k，1]）-1，b=Vecrev（二进制（f[k，2]））；总和（l=1，#b，Mod（b[l]，2）*'x^i*'y^（l-1））);}g（q）={prod（i=0，极度数（q，'x），my（p=素数（1+i），qx=polceof（q，i，'x））；触头（j=0，极性（qx），如果（polcoeffe（qx，j，'y），p^（2^j），1);););}T（n，k）=g（f（n）*f（k））t=0；对于（d=1100，对于（n=1，d，打印（t++“”t（n，d+1-n）））退出