|
（PARI）
Q（N）={
my（q=矢量（N））；q[1]=[1，0，0，0]；
对于（n=2，n，
我的（m=（平方（8*n+1）-1）\2）；
q[n]=矢量（（1+（m>>2））<<2）；q[n][1]=1；
对于（k=2，m，q[n][k]=q[n-k][k]+q[n-k][k-1]））；
收益（q）；
};
序列（N）={
my（a=矢量（N），q=q（N）；vmx=应用（vecmax，q））；
对于（n=1，n，a[n]=vecmin（选择（v->v==vmx[n]，q[n]，1））；
a；
};
| 