OFFSET
2,1
COMMENTS
If p is a prime, the p-th term is p.
LINKS
Ivan Neretin, Table of n, a(n) for n = 2..2000
MAPLE
for b from 2 do if isprime(b) then print(b): next: end if: L := ifactors(b)[2]: last := 0: for n from 2 to 10000 do C := [seq(0, i=1..nops(L))]: for i from 1 while 2^i <= n do for j from 1 to nops(L) do C := subsop(j=C[j]+floor(n/(L[j][1])^i), C): end do: end do: for j from 1 to nops(L) do C := subsop(j=floor(C[j]/(L[j][2])), C): end do: new := min(op(C)): if last + 1 < new then print(last+1): break: end if: last := new: end do: end do:
CROSSREFS
KEYWORD
nonn
AUTHOR
Asger Grunnet, Jul 13 2002
STATUS
approved