login
Numbers of the form x + y + z with distinct positive integers x,y,z such that (x+y+z) | x*y*z.
2

%I #11 Apr 13 2021 20:42:01

%S 6,10,12,14,15,16,18,20,21,22,24,26,27,28,30,32,33,34,35,36,38,39,40,

%T 42,44,45,46,48,49,50,51,52,54,55,56,57,58,60,62,63,64,65,66,68,69,70,

%U 72,74,75,76,77,78,80,81,82,84,85,86,87,88,90,91,92,93,94,95,96,98,99,100

%N Numbers of the form x + y + z with distinct positive integers x,y,z such that (x+y+z) | x*y*z.

%C From _Robert Israel_, Apr 12 2021: (Start)

%C All terms are composite.

%C Conjecture: Consists of all composite numbers except 4, 8, 9, and 25. (End)

%H Robert Israel, <a href="/A343311/b343311.txt">Table of n, a(n) for n = 1..10000</a>

%e 10 is in the sequence since 10 = 1+4+5 = 2+3+5, (1+4+5) | 1*4*5 and (2+3+5) | 2*3*5.

%e 12 is in the sequence since 12 = 1+3+8 = 2+4+6 = 3+4+5, (1+3+8) | 1*3*8, (2+4+6) | 2*4*6 and (3+4+5) | 3*4*5.

%p filter:= proc(n) local x,y,z;

%p if isprime(n) then return false fi;

%p x:= min(numtheory:-factorset(n));

%p y:= n/x;

%p z:= n - x - y;

%p if z > 0 and nops({x,y,z}) = 3 then return true fi;

%p for x from 1 to n/3 do

%p for y from x+1 while x+2*y+1 <= n do

%p z:= n-x-y;

%p if x*y*z mod n = 0 then return true fi;

%p od od;

%p false

%p end proc:

%p select(filter, [$1..100]); # _Robert Israel_, Apr 12 2021

%t Table[If[Sum[Sum[(1 - KroneckerDelta[i, j]) (1 - KroneckerDelta[n - j, 2 i]) (1 - KroneckerDelta[n - i, 2 j]) (1 - Ceiling[i*j*(n - i - j)/n] + Floor[i*j*(n - i - j)/n]), {i, j, Floor[(n - j)/2]}], {j, Floor[n/3]}] > 0, n, {}], {n, 100}] // Flatten

%Y Cf. A343270.

%K nonn

%O 1,1

%A _Wesley Ivan Hurt_, Apr 11 2021