solusi untuk $\frac{1}a + \frac{1}b + \frac{1}c = \frac{1}{2018}$

Aug 21 2020

Temukan, dengan bukti, semua triplet terurut dari bilangan bulat positif $(a,b,c)$ maka $\dfrac{1}a + \dfrac{1}b + \dfrac{1}c = \dfrac{1}{2018}.$

Secara umum, jika $d$ adalah bilangan bulat positif, lalu $(a,b,c) = (3d,3d,3d),(d,2d,6d),(d,6d,2d), (2d,d,6d),(2d,6d,d),(6d,d,2d),$ dan $(6d,2d,d)$ adalah semua triplet terurut dari bilangan bulat positif sehingga $\dfrac{1}a + \dfrac{1}b + \dfrac{1}c = \dfrac{1}d.$Namun, saya tidak yakin bagaimana menemukan semua kembar tiga seperti itu. Saya mencoba memanipulasi persamaan sebagai berikut:\begin{align} &1 + \dfrac{a}b + \dfrac{a}c = \dfrac{a}{2018}\\ &\dfrac{a}b + \dfrac{a}c = \dfrac{a-2018}{2018}\\ &\dfrac{ab}{c } = \dfrac{b(a-2018)-2018a}{2018}\\ &c[b(a-2018)-2018a]-2018ab = 0\\ &c[(a-2018)(b-2018)-2018^2]-2018ab = 0\\ &c[(a-2018)(b-2018)]-2018(ab-2018(a+b)+2018^2+2018(a+b)-2018^2)-2018^2 c =0 \\ &(c-2018)(a-2018)(b-2018)-2018^2(a+b+c)+2018^3 = 0\\ &(c-2018)(b-2018)(a-2018) = 2018^2(a+b+c) - 2018^3, \end{align}

tapi saya tidak tahu apakah ini berguna.

Sumber (dari komentar): Saya mendasarkan ini pada masalah kontes. Saya datang sendiri. Pertanyaan yang saya dasarkan adalah pertanyaan Putnam 2018 A1 .

Jawaban

3 RossMillikan Aug 21 2020 at 01:06

Saya menulis program Python kecil yang melakukan kekerasan. Itu ditemukan$670$ solusi dengan $a \le b \le c$ Yang pertama adalah $$a=2019, b= 4074343, c= 16600266807306$$ Ada $40$ dengan $a=2019$ dan yang terakhir adalah $$a=2019, b=c=8148684$$

Selama tidak ada yang tertawa terlalu keras, berikut kodenya. Saya biarkan saja$a$ berkisar dari $n+1$ untuk $3n$, hitung kisaran $b$ untuk menjadi seperti itu $\frac 1b \le \min(\frac 1a, \frac 1n-\frac 1a)$, hitung $c$menggunakan pembagian integer, lalu lihat apakah hasilnya genap. sukses menghitung keberhasilan. Ini adalah Python 2.

def prog(n, plev=0):
    succ=0
    astart=n+1
    aend=3*n
    if (plev > 19): print 'astart, aend',astart, aend
    for a in xrange(astart,aend+1):
        bstart=max(n*a/(a-n)+1,a)
        bend=2*n*a/(a-n)
        if (plev > 19):  print 'bstart, bend', bstart, bend
        for b in xrange(bstart, bend+1):
            c=n*a*b/(a*b-n*(a+b))
            if (n==a*b*c/(a*b+a*c+b*c)):
                print 'success',a,b,c
                succ+=1
    print 'successes',succ
3 Dr.Mathva Aug 21 2020 at 00:03

Jawaban parsial .$\,$ Membiarkan $d=\text{gcd}(a,b,c)$ dan dengan demikian, $a=dx, b=dy, c=dz$ - dimana $\text{gcd}(x,y,z)=1$. Jika kita berasumsi demikian$\text{gcd}(x,y)=\text{gcd}(y,z)=\text{gcd}(z,x)=1$ kemudian $\text{gcd}(xyz, xy+yz+zx)=1$ $$\frac1a+\frac1b+\frac1c=\frac1{2018}\iff \frac{abc}{ab+bc+ca}=\frac{dxyz}{xy+yz+zx}=2018$$ Ini menyiratkan itu $xy+yz+zx\mid d\iff d=(xy+yz+zx)\cdot k, k\in\mathbb Z$. Oleh karena itu, persamaannya menjadi$$k\cdot xyz=2018$$ Setelah Anda mendapatkan solusi bilangan bulat untuk persamaan ini - dan ini tidak akan memakan waktu lama $2018=2\cdot 1009$ -, gunakan nilai untuk mendapatkan kembali solusi $$(a,b,c)\equiv (kx\cdot (xy+yz+zx),ky\cdot (xy+yz+zx), kz\cdot (xy+yz+zx) )$$

Pengamatan . Karena simetri, Anda "hanya" harus mempertimbangkan kasusnya$\text{gcd}(x,y)>1$.

DmitryEzhov Aug 20 2020 at 23:46

$\dfrac{1}a + \dfrac{1}b + \dfrac{1}c = \dfrac{1}{2018}\implies \Big(a b - 2018 (a + b)\Big) \Big(a c - 2018 (a + c)\Big)=\Big(2018a\Big)^2$

Membiarkan $a<b<c$ dan $\Big(2018a\Big)^2=mn$, di tempat yang alami $m<n$.

Kemudian $b=\dfrac{2018 a + m}{a-2018}$ dan $c=\dfrac{2018 a + n}{a-2018}$

kode pari / gp:

abc2018()=
{
 s= 0;
 for(a=2019, 10^4,
  d= (2018*a)^2;
  D= divisors(d);
  for(i=1, #D\2,
   m= D[i];
   b= (2018*a+m)/(a-2018);
   if(a<b, if(b==floor(b),
    n= d/m;
    c= (2018*a+n)/(a-2018);
    if(b<c, if(c==floor(c),
     if(1/a + 1/b + 1/c == 1/2018,
      s++;
      print("("a", "b", "c")")
     )
    ))
   ))
  )
 );
 print("\nNumber of solutions: "s)
};

Keluaran:

(2019, 4074343, 16600266807306)
(2019, 4074344, 8300135440824)
(2019, 4074345, 5533424985330)
(2019, 4074346, 4150069757583)
(2019, 4074348, 2766714529836)
(2019, 4074351, 1844477711338)
(2019, 4074354, 1383359302089)
(2019, 4074360, 922240892840)
(2019, 4074378, 461122483591)
(2019, 4075015, 24670140810)
(2019, 4075351, 16456267338)
(2019, 4075688, 12337107576)
(2019, 4076360, 8230170840)
(2019, 4076361, 8226096498)
(2019, 4077034, 6170590959)
(2019, 4077369, 5488138674)
. . .
(3636, 4545, 2038180)
(3687, 4458, 921365323)
(3700, 4440, 22399800)
(3828, 4268, 187328922)
(3885, 4200, 22399800)
(3940, 4137, 83484660)
(3960, 4115, 657682344)
(4002, 4071, 39707177)
(4036, 4037, 16293332)
(4036, 4038, 8148684)
(4036, 4040, 4076360)
(4036, 4044, 2040198)
(4036, 4052, 1022117)
(4036, 5045, 20180)
(4036, 6054, 12108)

Number of solutions: 658