%I #16 Aug 08 2023 12:19:34
%S 1,3,5,7,9,11,13,17,19,21,23,25,27,29,31,33,37,39,41,43,47,49,51,53,
%T 57,59,61,67,69,71,73,79,81,83,87,89,93,95,97,101,103,107,109,111,113,
%U 115,121,123,125,127,129,131,133,137,139,141,145,147,149,151,155,157,159
%N Numbers k with the property that (k-x)*(k-y)*(k-z) = x*y*z has no integer solutions 0 < x,y,z < k.
%C All odd primes are terms of this sequence.
%e 15 is not a term of this sequence because (15-x)*(15-y)*(15-z) = x*y*z has the solution (5,5,12).
%o (Python)
%o def exis(n):
%o for x in range(1,n):
%o for y in range(x+1):
%o for z in range(y+1):
%o if x*y*z==(k-x)*(k-y)*(k-z):
%o return True
%o return False
%o for k in range(1, 200, 2):
%o if not exis(k):
%o print(str(k), end=',')
%o (PARI) is(n)=for(x=1,n-1,for(y=1,x, my(t=(n-x)*(n-y),z=t*n/(x*y+t)); if(denominator(z)==1 && 0 < z && z < n, return(0)))); 1 \\ _Charles R Greathouse IV_, Dec 09 2014
%Y Cf. A065091 (odd primes).
%K nonn
%O 1,2
%A Nurdin Takenov (greanvert(AT)gmail.com), Sep 14 2009
%E Edited by _Charles R Greathouse IV_, Dec 09 2014