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
KEYWORD
nonn
AUTHOR
Max Alekseyev, Oct 19 2014
STATUS
approved