{ nb = vector(1 000); for (n=1, #a=vector(10 000), if (n==1, v = 1, v > #nb, break, v = nb[v]; ); print (n " " a[n]=v); s = 1; forstep (k=n, 1, -1, s = lcm(s, a[k]); if (s>#nb, break, nb[s]++ ); ); ); } quit