A249016 Indices of records in A055460 (number of primes dividing the squarefree part of n!).

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

Max Alekseyev, Table of n, a(n) for n = 1..38282

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

Adjacent sequences: A249013 A249014 A249015 * A249017 A249018 A249019

Keyword

nonn,changed

Author

Max Alekseyev, Oct 19 2014

Status

approved