|
|
A064463
|
|
Numbers k such that A001414(k) is a square and sets a new record for squares.
|
|
1
|
|
|
1, 4, 14, 39, 46, 94, 158, 291, 334, 446, 718, 878, 1454, 1678, 2174, 2446, 2734, 3694, 4414, 4798, 6046, 7438, 7934, 9518, 10078, 11246, 11854, 15838, 17294, 21214, 22894, 27374, 28318, 29278, 32254, 34318, 36446, 38638, 42046, 48046, 51838
(list;
graph;
refs;
listen;
history;
text;
internal format)
|
|
|
OFFSET
|
1,2
|
|
LINKS
|
|
|
PROG
|
(PARI) sopfr(n, s, fac, i) = fac=factor(n); for(i=1, matsize(fac)[1], s=s+fac[i, 1]*fac[i, 2]); return(s); aa=-1; j=[]; for(n=1, 150000, x=sopfr(n); if(issquare(x), b=x; if(b>aa, aa=b; j=concat(j, n)))); j
(PARI) sopfr(n)= { local(f, s=0); f=factor(n); for(i=1, matsize(f)[1], s+=f[i, 1]*f[i, 2]); return(s) } { n=0; a=-1; for (m=1, 10^9, if (issquare(b=sopfr(m)), if(b>a, a=b; write("b064463.txt", n++, " ", m); if (n==600, break))) ) } \\ Harry J. Smith, Sep 15 2009
|
|
CROSSREFS
|
|
|
KEYWORD
|
easy,nonn
|
|
AUTHOR
|
|
|
STATUS
|
approved
|
|
|
|