A249017 Records in A055460 (number of primes dividing the squarefree part of n!).

0, 1, 2, 3, 4, 5, 6, 7, 9, 10, 11, 12, 14, 15, 16, 17, 18, 19, 20, 22, 23, 24, 25, 26, 27, 28, 30, 31, 32, 34, 36, 37, 38, 39, 41, 42, 43, 44, 46, 47, 49, 50, 51, 52, 53, 54, 55, 57, 58, 59, 60, 61, 63, 64, 65, 67, 68, 69, 70, 72, 74, 75, 76, 77, 78, 79, 80, 81, 82, 85, 86, 87, 88, 91, 93, 95, 98, 99, 100

Offset

1,3

Links

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

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(b, ", ")) )

Crossrefs

Cf. A055460, A055204, A249016 (indices of records).

Sequence in context: A283970 A376268 A121166 * A214879 A187226 A026470

Adjacent sequences: A249014 A249015 A249016 * A249018 A249019 A249020

Keyword

nonn

Author

Max Alekseyev, Oct 19 2014

Status

approved