16 lines
273 B
Text
16 lines
273 B
Text
Hupto(n)={
|
|
my(v=vector(n),x2=2,x3=3,x5=5,i=1,j=1,k=1,t);
|
|
v[1]=1;
|
|
for(m=2,n,
|
|
v[m]=t=min(x2,min(x3,x5));
|
|
if(x2 == t, x2 = v[i++] << 1);
|
|
if(x3 == t, x3 = 3 * v[j++]);
|
|
if(x5 == t, x5 = 5 * v[k++]);
|
|
);
|
|
v
|
|
};
|
|
H(n)=Hupto(n)[n];
|
|
|
|
Hupto(20)
|
|
H(1691)
|
|
H(10^6)
|