A091509 Numbers (in ascending order, without multiplicity) which occur in triples a,b,c such that a < b < c, (b * c) / (b + c) leaves a, (a * c) / (a + c) leaves b and (a * b) / (a + b) leaves c as a remainder.

5, 7, 11, 69, 87, 99, 111, 125, 153, 184, 207, 315, 324, 335, 365, 384, 391, 425, 455, 475, 494, 525, 575, 589, 612, 625, 663, 671, 675, 679, 704, 707, 711, 713, 720, 725, 731, 735, 736, 747, 749, 756, 763, 776, 783, 816, 847, 851, 864, 896, 901, 913, 923

Offset

1,1

Comments

First, second and third component of the triples are resp. in A092817, A092818, A092819. 5, 7 and 11 (forming the first triple) seem to be the only primes in the sequence. The greatest common divisor of the triples is in A092864. Some numbers occur in two (111, 475, 864, 931), some in three (575, 3015, 3200, 3311) triples. It is easy to prove that a triple has either only odd or only even components.

By means of the PARI program 337 triples (sum of first and second component < 20000) were computed, comprising 925 different terms. The 86th triple is the last one of these which has a first component < 1000, so it is almost certain that no terms are missing in the list given above. But to prove this one has to find a lower bound for the first component of the k-th triple.

Links

Table of n, a(n) for n=1..53.

Example

5, 7 and 11 form a triple because (7 * 11) / (7 + 11) = 77/18 leaves 5, (5 * 11) / (5 + 11) = 55/16 leaves 7 and (5 * 7) / (5 + 7) = 35/12 leaves 11 as a remainder, hence 11 (as well as 5 and 7) is a term.

Prog

(QBasic) INPUT "Type in the maximum value for c: "; maxc FOR a = 1 TO maxc - 2 FOR b = a + 1 TO maxc - 1 FOR c = b + 1 TO maxc IF (b * c) MOD (b + c) = a AND (a * c) MOD (a + c) = b AND (a * b) MOD (a + b) = c THEN PRINT a; b; c NEXT NEXT NEXT
(PARI) {m=12500; w=[]; for(n=3, m, for(a=1, (n-1)\2, b=n-a; c=a*b%(a+b); if(b<c, if((b*c)%(b+c)==a, if((a*c)%(a+c)==b, w=concat(w, [a, b, c])))))); w=vecsort(w); print1(w[1], ", "); for(j=2, length(w), if(w[j]>w[j-1], print1(w[j], ", "))); }

Crossrefs

Sequence in context: A098036 A127269 A071781 * A027728 A218275 A086140

Adjacent sequences: A091506 A091507 A091508 * A091510 A091511 A091512

Keyword

nonn

Author

Tamas Sandor Nagy, Mar 03 2004

Extensions

Edited and corrected by Klaus Brockhaus, Mar 07 2004

Status

approved