|
%I #11 Nov 14 2021 10:15:33
%S 11,14,17,19,20,27,32,34,36,43,46,49,52,64,67,73,82,97,100,142,148,
%T 163,193
%N Numbers having a unique representation as ab+ac+bc, with 0 < a < b < c.
%C Are there more terms?
%C No more terms < 10^6. - _Joerg Arndt_, Oct 01 2017
%D See A025052.
%e 11 is on the list because 11 = 1*2+1*3+2*3.
%t oneSol={}; Do[lim=Ceiling[(n-2)/3]; cnt=0; Do[If[n>a*b && Mod[n-a*b, a+b]==0 && Quotient[n-a*b, a+b]>b, cnt++; If[cnt>1, Break[]]], {a, 1, lim-1}, {b, a+1, lim}]; If[cnt==1, AppendTo[oneSol, n]], {n, 10000}]; oneSol
%o (Python)
%o from collections import Counter
%o def aupto(N):
%o acount = Counter()
%o for i in range(1, N-1):
%o for j in range(i+1, N//i + 1):
%o p, s = i*j, i+j
%o for k in range(j+1, (N-p)//s + 1):
%o acount.update([p + s*k])
%o return sorted([k for k in acount if acount[k] == 1])
%o print(aupto(10**5)) # _Michael S. Branicky_, Nov 14 2021
%Y Cf. A025052, A025058, A025060.
%Y Cf. A000926 (numbers not of the form ab+ac+bc, 0<a<b<c).
%K nonn,more
%O 1,1
%A _T. D. Noe_, Apr 08 2004
|
