login
A249016
Indices of records in A055460 (number of primes dividing the squarefree part of n!).
3
1, 2, 3, 5, 13, 17, 21, 23, 33, 42, 43, 56, 57, 75, 84, 99, 101, 105, 109, 119, 133, 139, 157, 162, 163, 182, 183, 207, 208, 219, 220, 255, 257, 263, 267, 303, 305, 307, 315, 340, 341, 343, 383, 385, 387, 397, 411, 420, 421, 423, 469, 483, 485, 489, 505, 519, 523, 547, 552, 553, 581, 602, 603, 609, 618
OFFSET
1,2
LINKS
PROG
(PARI) np=vector(10^6); b=-1; r=0; for(n=1, 10^6, f=factor(n); for(i=1, matsize(f)[1], if(f[i, 2]%2, r += (-1)^np[f[i, 1]]; np[f[i, 1]]=1-np[f[i, 1]]; )); if(r>b, b=r; print1(n, ", "))
CROSSREFS
Cf. A249017 (values of records), A055204
Sequence in context: A236183 A093077 A360418 * A241123 A038983 A235635
KEYWORD
nonn
AUTHOR
Max Alekseyev, Oct 19 2014
STATUS
approved