(PARI)\\要生成图形:
对于(n=1,尾,u1=1;u2=n;而(u2!=u1)||(u2%2)==1),u3=u2;u2=地板(u2/2)+地板(u1/2);u1=u3;);b[n]=u2;)fr(m,k)=(3*总和(i=1,k,b[i]))-k^2+m*k;bound(m)=if((m%2)==1,p=(m-1)/2;5/3*(4^p-1),2/3*(4 ^(m/2)-1));m=5;分形=矢量(边界(m));对于(i=1,界(m),分形[i]=fr(m,i););Mm=vecmax(分形)指数=矢量(界限(m));对于(i=1,界(m),指数[i]=i);psplotraw(指数,分形,1);
(PARI)
up_to=10000;
A073504list(up_to)={my(v=向量(up_tto));对于(n=1,up_to,u1=1;u2=n;当((u2!=u1)||((u2%2)==1),u3=u2;u2=(u2\2)+(u1\2);u1=u3);v[n]=u2);(v);};
v073504=A073504list(up_to);
|