(PARI)
F=矢量(30,n,n^3);/*修改以获得“贪婪表示”*/last_leq(v,F)的其他序列=
{/*返回排序数组F[]中的最后一个元素<=v*/
局部(j=1);
而(F[j]<=v,j+=1);
返回(F[j-1]);
}
贪婪(n,F)=
{
局部(v=n,ct=0);
而(v,v-=last_leq(v,F);ct+=1;);
返回(ct);
}
向量(min(100,F[#F-1]),n,贪婪(n,F))/*显示项*/
/*Joerg Arndt,2011年4月8日*/
(哈斯克尔)
a055401 n=s n$reverse$takeWhile(<=n)$tail a000578_list,其中
s _[]=0
s m(x:xs)|x>m=s m xs
|否则=m'+s r xs,其中(m',r)=divMod m x
(方案,带有备忘录-宏定义)
|