15 lines
282 B
Text
15 lines
282 B
Text
Hupto(n)={
|
|
my(r=Vec([1],n),v=primes(3),[v1,v2,v3]=v,i=1,j=1,k=1,t);
|
|
for(m=2,n,
|
|
r[m]=t=min(v1,min(v2,v3));
|
|
if(v1 == t, v1 = v[1] * r[i++]);
|
|
if(v2 == t, v2 = v[2] * r[j++]);
|
|
if(v3 == t, v3 = v[3] * r[k++]);
|
|
);
|
|
r
|
|
};
|
|
H(n)=Hupto(n)[n];
|
|
|
|
Hupto(20)
|
|
H(1691)
|
|
H(10^6)
|