b: =proc(n,i)选项记忆`如果`(n=0,[[]],
`如果`(i<1,[],[b(n,i-1)[],seq(映射(p->[p[],[i,j]]),
b(n-i*j,i-1)[],j=1…n/i)]))
结束时间:
a: =程序(n)不适用(添加(添加(x^添加(添加)i[2]*j[2]*
igcd(i[1],j[1]),j=t),i=s)/mul(i[1]^i[2]!,i=s)
/mul(i[1]^i[2]*i[2]!,i=t),t=b(n$2)),s=b(n$2),x)
结束(8):
seq(a(n),n=0..10);