%I #19 May 04 2021 18:10:56
%S 64,72,81,100,108,125,128,216,225,288,324,500,576,864,972,1125,1152,
%T 1225,1800,2025,2700,3125,3200,3528,4500,7776,8100,10125,13068,13689,
%U 15488,17496,18496,21125,24500,28800,34848,42336,44100,48672,55225,69696,93636,95256
%N Numbers k that can be expressed as k = w+x = y*z with w*x = k*(y+z) where w, x, y, and z are all positive integers.
%C All terms are powerful (A001694).
%e 64 is a term because a solution exists at k=64, w=32, x=32, y=8, z=8:
%e k = w + x = y*z with w*x = k*(y+z)
%e becomes
%e 64 = 32 + 32 = 8*8 with 32*32 = 64*(8+8).
%o (PARI) is(k) = fordiv(k, y, if(issquare(k^2 - 4*k*(y+k/y)), return(1))); 0; \\ _Jinyuan Wang_, May 01 2021
%Y Cf. A001694, A057369, A057370, A057372, A057373, A057421.
%K nonn
%O 1,1
%A _Naohiro Nomoto_, Sep 23 2000
%E More terms from _Jinyuan Wang_, May 01 2021
|