|
%I #14 Jun 24 2018 09:01:16
%S 1,4,14,39,46,94,158,291,334,446,718,878,1454,1678,2174,2446,2734,
%T 3694,4414,4798,6046,7438,7934,9518,10078,11246,11854,15838,17294,
%U 21214,22894,27374,28318,29278,32254,34318,36446,38638,42046,48046,51838
%N Numbers k such that A001414(k) is a square and sets a new record for squares.
%H Harry J. Smith, <a href="/A064463/b064463.txt">Table of n, a(n) for n = 1..600</a>
%o (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
%o (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
%Y Cf. A001414.
%K easy,nonn
%O 1,2
%A _Jason Earls_, Oct 02 2001
|
