OFFSET
1,1
PROG
(PARI) lista() = {my(nb = 10^7); my(ip = 1); while (factorback(primes(ip)) < nb, ip++); my(vp = vector(ip), list = List([2]), listi = List([1]), lastnb=1, lasti=1, nba = 3, nbb = nb); vp[lasti]++; for (i=1, 2, forsquarefree(isqf=nba, nbb, my(newi = omega(isqf[2])); vp[newi]++; if (newi != lasti, if (vp[newi] > vp[lasti], lasti = newi; listput(list, isqf[1]); listput(listi, vp[newi]); ); ); ); nba = nbb + 1; nbb = nbb + nb; ); Vec(list); } \\ Michel Marcus, Dec 22 2022
CROSSREFS
KEYWORD
nonn,more
AUTHOR
Michel Marcus, Dec 22 2022
EXTENSIONS
Name clarified by Peter Munn, Dec 31 2022
STATUS
approved