SQ[n_]:=SQ[n]=整数Q[Sqrt[n]];
p[n_]:=p[n]=分区p[n];
笔[n_]:=笔[n]=SQ[24n+1]&&模型[Sqrt[24n+1],6]==1;
Do[m=1;标签[bb];如果[p[m]>n,转到[cc]];执行[If[Pen[n-p[m]-x(3x-1)/2],打印[n,“”,m];转到[aa]],{x,0,(平方[24(n-p[m])+1]+1)/6}];m=m+1;转到[bb];标签[cc];打印[n,“”,0];标签[aa];继续,{n,1,80}]
|