2015年8月23日日曜日

150823(6)

Ruby


どの2つの和も立方数(9)

どの2つの和も立方数となるような異なる4個の自然数の組を、
OEIS A018787(https://oeis.org/A018787/b018787.txt)のデータを使って、
以下のコードで作ってみた。

# OEIS A018787(https://oeis.org/A018787/b018787.txt)のデータの3000行分使って計算
3000.times{
  m0, m = gets.chomp.split(' ').map{|i| i.to_i}
  i = []
  j, s = 1, 0
  while s < 3
    k = m - j * j * j
    l = (k ** (1.0 / 3)).round
    if k == l * l * l
      i.push([j, l])
      s += 1
    end
    j += 1
  end
  i0 = i[1][0] ** 3 + i[1][1] ** 3
  j0 = i[0][0] ** 3 + i[1][0] ** 3 + i[2][0] ** 3 - i0
  if j0 > 0
    if j0 % 2 == 0
      a0 = j0 / 2
      ary = [a0, i[0][0] ** 3 - a0, i[1][0] ** 3 - a0, i[2][0] ** 3 - a0]
      p [1, ary, ary.inject(:+)]
      puts "#{ary[0]} + #{ary[1]} = #{i[0][0]}^3"
      puts "#{ary[0]} + #{ary[2]} = #{i[1][0]}^3"
      puts "#{ary[0]} + #{ary[3]} = #{i[2][0]}^3"
      puts "#{ary[1]} + #{ary[2]} = #{i[2][1]}^3"
      puts "#{ary[1]} + #{ary[3]} = #{i[1][1]}^3"
      puts "#{ary[2]} + #{ary[3]} = #{i[0][1]}^3"
    end
  end
  j1 = i[0][0] ** 3 + i[1][0] ** 3 + i[2][1] ** 3 - i0
  if j1 > 0
    if j1 % 2 == 0
      a1 = j1 / 2
      ary = [a1, i[0][0] ** 3 - a1, i[1][0] ** 3 - a1, i[2][1] ** 3 - a1]
      p [2, ary, ary.inject(:+)]
      puts "#{ary[0]} + #{ary[1]} = #{i[0][0]}^3"
      puts "#{ary[0]} + #{ary[2]} = #{i[1][0]}^3"
      puts "#{ary[0]} + #{ary[3]} = #{i[2][1]}^3"
      puts "#{ary[1]} + #{ary[2]} = #{i[2][0]}^3"
      puts "#{ary[1]} + #{ary[3]} = #{i[1][1]}^3"
      puts "#{ary[2]} + #{ary[3]} = #{i[0][1]}^3"
    end
  end
}

出力結果
[2, [1493580092, 2510949380, 7071213244, 51114333631], 62190076347]
1493580092 + 2510949380 = 1588^3
1493580092 + 7071213244 = 2046^3
1493580092 + 51114333631 = 3747^3
2510949380 + 7071213244 = 2124^3
2510949380 + 51114333631 = 3771^3
7071213244 + 51114333631 = 3875^3
[2, [11948640736, 20087595040, 56569705952, 408914669048], 497520610776]
11948640736 + 20087595040 = 3176^3
11948640736 + 56569705952 = 4092^3
11948640736 + 408914669048 = 7494^3
20087595040 + 56569705952 = 4248^3
20087595040 + 408914669048 = 7542^3
56569705952 + 408914669048 = 7750^3
[2, [40326662484, 67795633260, 190922757588, 1380087008037], 1679132061369]
40326662484 + 67795633260 = 4764^3
40326662484 + 190922757588 = 6138^3
40326662484 + 1380087008037 = 11241^3
67795633260 + 190922757588 = 6372^3
67795633260 + 1380087008037 = 11313^3
190922757588 + 1380087008037 = 11625^3
[2, [95589125888, 160700760320, 452557647616, 3271317352384], 3980164886208]
95589125888 + 160700760320 = 6352^3
95589125888 + 452557647616 = 8184^3
95589125888 + 3271317352384 = 14988^3
160700760320 + 452557647616 = 8496^3
160700760320 + 3271317352384 = 15084^3
452557647616 + 3271317352384 = 15500^3
[2, [186697511500, 313868672500, 883901655500, 6389291703875], 7773759543375]
186697511500 + 313868672500 = 7940^3
186697511500 + 883901655500 = 10230^3
186697511500 + 6389291703875 = 18735^3
313868672500 + 883901655500 = 10620^3
313868672500 + 6389291703875 = 18855^3
883901655500 + 6389291703875 = 19375^3
[2, [730327591966, 1085966964635, 2233317468578, 5518494262093], 9568106287272]
730327591966 + 1085966964635 = 12201^3
730327591966 + 2233317468578 = 14364^3
730327591966 + 5518494262093 = 18419^3
1085966964635 + 2233317468578 = 14917^3
1085966964635 + 5518494262093 = 18762^3
2233317468578 + 5518494262093 = 19791^3
[2, [22459099604, 937293264940, 1726359954988, 10593628833748], 13279741153280]
22459099604 + 937293264940 = 9864^3
22459099604 + 1726359954988 = 12048^3
22459099604 + 10593628833748 = 21978^3
937293264940 + 1726359954988 = 13862^3
937293264940 + 10593628833748 = 22592^3
1726359954988 + 10593628833748 = 23096^3
[2, [322613299872, 542365066080, 1527382060704, 11040696064296], 13433056490952]
322613299872 + 542365066080 = 9528^3
322613299872 + 1527382060704 = 12276^3
322613299872 + 11040696064296 = 22482^3
542365066080 + 1527382060704 = 12744^3
542365066080 + 11040696064296 = 22626^3
1527382060704 + 11040696064296 = 23250^3
[2, [323294653868, 950466050132, 1095982079404, 12371699929684], 14741442713088]
323294653868 + 950466050132 = 10840^3
323294653868 + 1095982079404 = 11238^3
323294653868 + 12371699929684 = 23328^3
950466050132 + 1095982079404 = 12696^3
950466050132 + 12371699929684 = 23706^3
1095982079404 + 12371699929684 = 23792^3
[2, [512297971556, 861255637340, 2425426142692, 17532216435433], 21331196187021]
512297971556 + 861255637340 = 11116^3
512297971556 + 2425426142692 = 14322^3
512297971556 + 17532216435433 = 26229^3
861255637340 + 2425426142692 = 14868^3
861255637340 + 17532216435433 = 26397^3
2425426142692 + 17532216435433 = 27125^3
[2, [764713007104, 1285606082560, 3620461180928, 26170538819072], 31841319089664]
764713007104 + 1285606082560 = 12704^3
764713007104 + 3620461180928 = 16368^3
764713007104 + 26170538819072 = 29976^3
1285606082560 + 3620461180928 = 16992^3
1285606082560 + 26170538819072 = 30168^3
3620461180928 + 26170538819072 = 31000^3
[2, [293917302828, 482234256172, 4855120970004, 30329251118828], 35960523647832]
293917302828 + 482234256172 = 9190^3
293917302828 + 4855120970004 = 17268^3
293917302828 + 30329251118828 = 31286^3
482234256172 + 4855120970004 = 17476^3
482234256172 + 30329251118828 = 31350^3
4855120970004 + 30329251118828 = 32768^3
[2, [1088819887068, 1830482098020, 5154914454876, 37262349216999], 45336565656963]
1088819887068 + 1830482098020 = 14292^3
1088819887068 + 5154914454876 = 18414^3
1088819887068 + 37262349216999 = 33723^3
1830482098020 + 5154914454876 = 19116^3
1830482098020 + 37262349216999 = 33939^3
5154914454876 + 37262349216999 = 34875^3
[2, [380440602076, 658453015076, 6618119383076, 47649632272924], 55306645273152]
380440602076 + 658453015076 = 10128^3
380440602076 + 6618119383076 = 19128^3
380440602076 + 47649632272924 = 36350^3
658453015076 + 6618119383076 = 19378^3
658453015076 + 47649632272924 = 36420^3
6618119383076 + 47649632272924 = 37860^3
[2, [1493580092000, 2510949380000, 7071213244000, 51114333631000], 62190076347000]
1493580092000 + 2510949380000 = 15880^3
1493580092000 + 7071213244000 = 20460^3
1493580092000 + 51114333631000 = 37470^3
2510949380000 + 7071213244000 = 21240^3
2510949380000 + 51114333631000 = 37710^3
7071213244000 + 51114333631000 = 38750^3
[2, [5842620735728, 8687735717080, 17866539748624, 44147954096744], 76544850298176]
5842620735728 + 8687735717080 = 24402^3
5842620735728 + 17866539748624 = 28728^3
5842620735728 + 44147954096744 = 36838^3
8687735717080 + 17866539748624 = 29834^3
8687735717080 + 44147954096744 = 37524^3
17866539748624 + 44147954096744 = 39582^3
[2, [1987955102452, 3342073624780, 9411784827764, 68033178062861], 82774991617857]
1987955102452 + 3342073624780 = 17468^3
1987955102452 + 9411784827764 = 22506^3
1987955102452 + 68033178062861 = 41217^3
3342073624780 + 9411784827764 = 23364^3
3342073624780 + 68033178062861 = 41481^3
9411784827764 + 68033178062861 = 42625^3
[2, [1380092762632, 9296973645368, 20258245135264, 64838664561361], 95773976104625]
1380092762632 + 9296973645368 = 22020^3
1380092762632 + 20258245135264 = 27866^3
1380092762632 + 64838664561361 = 40457^3
9296973645368 + 20258245135264 = 30918^3
9296973645368 + 64838664561361 = 42009^3
20258245135264 + 64838664561361 = 43985^3
[2, [179672796832, 7498346119520, 13810879639904, 84749030669984], 106237929226240]
179672796832 + 7498346119520 = 19728^3
179672796832 + 13810879639904 = 24096^3
179672796832 + 84749030669984 = 43956^3
7498346119520 + 13810879639904 = 27724^3
7498346119520 + 84749030669984 = 45184^3
13810879639904 + 84749030669984 = 46192^3
[2, [2580906398976, 4338920528640, 12219056485632, 88325568514368], 107464451927616]
2580906398976 + 4338920528640 = 19056^3
2580906398976 + 12219056485632 = 24552^3
2580906398976 + 88325568514368 = 44964^3
4338920528640 + 12219056485632 = 25488^3
4338920528640 + 88325568514368 = 45252^3
12219056485632 + 88325568514368 = 46500^3
[2, [2586357230944, 7603728401056, 8767856635232, 98973599437472], 117931541704704]
2586357230944 + 7603728401056 = 21680^3
2586357230944 + 8767856635232 = 22476^3
2586357230944 + 98973599437472 = 46656^3
7603728401056 + 8767856635232 = 25392^3
7603728401056 + 98973599437472 = 47412^3
8767856635232 + 98973599437472 = 47584^3
[2, [3281395462124, 5516555787860, 15535455497068, 112298190987307], 136631597734359]
3281395462124 + 5516555787860 = 20644^3
3281395462124 + 15535455497068 = 26598^3
3281395462124 + 112298190987307 = 48711^3
5516555787860 + 15535455497068 = 27612^3
5516555787860 + 112298190987307 = 49023^3
15535455497068 + 112298190987307 = 50375^3
[2, [9152310064644, 23911487556092, 33995215320580, 69801278935356], 136860291876672]
9152310064644 + 23911487556092 = 32096^3
9152310064644 + 33995215320580 = 35074^3
9152310064644 + 69801278935356 = 42900^3
23911487556092 + 33995215320580 = 38688^3
23911487556092 + 69801278935356 = 45422^3
33995215320580 + 69801278935356 = 46996^3
[2, [10475359466035, 16113643762388, 32986624058644, 94672005155564], 154247632442631]
10475359466035 + 16113643762388 = 29847^3
10475359466035 + 32986624058644 = 35159^3
10475359466035 + 94672005155564 = 47199^3
16113643762388 + 32986624058644 = 36618^3
16113643762388 + 94672005155564 = 48028^3
32986624058644 + 94672005155564 = 50352^3
[2, [4098383772448, 6890045098720, 19403409141536, 140257731483464], 170649569496168]
4098383772448 + 6890045098720 = 22232^3
4098383772448 + 19403409141536 = 28644^3
4098383772448 + 140257731483464 = 52458^3
6890045098720 + 19403409141536 = 29736^3
6890045098720 + 140257731483464 = 52794^3
19403409141536 + 140257731483464 = 54250^3
[2, [5040832810500, 8474454157500, 23865344698500, 172510876004625], 209891507671125]
5040832810500 + 8474454157500 = 23820^3
5040832810500 + 23865344698500 = 30690^3
5040832810500 + 172510876004625 = 56205^3
8474454157500 + 23865344698500 = 31860^3
8474454157500 + 172510876004625 = 56565^3
23865344698500 + 172510876004625 = 58125^3
[2, [19453816042236, 24533160957764, 37448437276611, 152611348285764], 234046762562375]
19453816042236 + 24533160957764 = 35300^3
19453816042236 + 37448437276611 = 38463^3
19453816042236 + 152611348285764 = 55620^3
24533160957764 + 37448437276611 = 39575^3
24533160957764 + 152611348285764 = 56162^3
37448437276611 + 152611348285764 = 57495^3
[2, [6117704056832, 10284848660480, 28963689447424, 209364310552576], 254730552717312]
6117704056832 + 10284848660480 = 25408^3
6117704056832 + 28963689447424 = 32736^3
6117704056832 + 209364310552576 = 59952^3
10284848660480 + 28963689447424 = 33984^3
10284848660480 + 209364310552576 = 60336^3
28963689447424 + 209364310552576 = 62000^3
[2, [19718844983082, 29321108045145, 60299571651606, 148999345076511], 258338869756344]
19718844983082 + 29321108045145 = 36603^3
19718844983082 + 60299571651606 = 43092^3
19718844983082 + 148999345076511 = 55257^3
29321108045145 + 60299571651606 = 44751^3
29321108045145 + 148999345076511 = 56286^3
60299571651606 + 148999345076511 = 59373^3
[2, [2351338422624, 3857874049376, 38840967760032, 242634008950624], 287684189182656]
2351338422624 + 3857874049376 = 18380^3
2351338422624 + 38840967760032 = 34536^3
2351338422624 + 242634008950624 = 62572^3
3857874049376 + 38840967760032 = 34952^3
3857874049376 + 242634008950624 = 62700^3
38840967760032 + 242634008950624 = 65536^3
[2, [7337958991996, 12336294303940, 34740870667772, 251124721129103], 305539845092811]
7337958991996 + 12336294303940 = 26996^3
7337958991996 + 34740870667772 = 34782^3
7337958991996 + 251124721129103 = 63699^3
12336294303940 + 34740870667772 = 36108^3
12336294303940 + 251124721129103 = 64107^3
34740870667772 + 251124721129103 = 65875^3
[2, [606395689308, 25306918153380, 46611718784676, 286027978511196], 358553011138560]
606395689308 + 25306918153380 = 29592^3
606395689308 + 46611718784676 = 36144^3
606395689308 + 286027978511196 = 65934^3
25306918153380 + 46611718784676 = 41586^3
25306918153380 + 286027978511196 = 67776^3
46611718784676 + 286027978511196 = 69288^3
[2, [8710559096544, 14643856784160, 41239315639008, 298098793735992], 362692525255704]
8710559096544 + 14643856784160 = 28584^3
8710559096544 + 41239315639008 = 36828^3
8710559096544 + 298098793735992 = 67446^3
14643856784160 + 41239315639008 = 38232^3
14643856784160 + 298098793735992 = 67878^3
41239315639008 + 298098793735992 = 69750^3
[2, [8728955654436, 25662583353564, 29591516143908, 334035898101468], 398018953253376]
8728955654436 + 25662583353564 = 32520^3
8728955654436 + 29591516143908 = 33714^3
8728955654436 + 334035898101468 = 69984^3
25662583353564 + 29591516143908 = 38088^3
25662583353564 + 334035898101468 = 71118^3
29591516143908 + 334035898101468 = 71376^3
[2, [10244465851028, 17222601797420, 48501451640596, 350593214375029], 426561733664073]
10244465851028 + 17222601797420 = 30172^3
10244465851028 + 48501451640596 = 38874^3
10244465851028 + 350593214375029 = 71193^3
17222601797420 + 48501451640596 = 40356^3
17222601797420 + 350593214375029 = 71649^3
48501451640596 + 350593214375029 = 73625^3
[1, [16955834122364, 117750072394364, 128741810605636, 163827996114820], 427275713237184]
16955834122364 + 117750072394364 = 51262^3
16955834122364 + 128741810605636 = 52620^3
16955834122364 + 163827996114820 = 56544^3
117750072394364 + 128741810605636 = 62700^3
117750072394364 + 163827996114820 = 65544^3
128741810605636 + 163827996114820 = 66386^3
[2, [49809860503772, 84896046012956, 95887784224228, 196682022496228], 427275713237184]
49809860503772 + 84896046012956 = 51262^3
49809860503772 + 95887784224228 = 52620^3
49809860503772 + 196682022496228 = 62700^3
84896046012956 + 95887784224228 = 56544^3
84896046012956 + 196682022496228 = 65544^3
95887784224228 + 196682022496228 = 66386^3
[2, [3043524816608, 5267624120608, 52944955064608, 381197058183392], 442453162185216]
3043524816608 + 5267624120608 = 20256^3
3043524816608 + 52944955064608 = 38256^3
3043524816608 + 381197058183392 = 72700^3
5267624120608 + 52944955064608 = 38756^3
5267624120608 + 381197058183392 = 72840^3
52944955064608 + 381197058183392 = 75720^3
[2, [11948640736000, 20087595040000, 56569705952000, 408914669048000], 497520610776000]
11948640736000 + 20087595040000 = 31760^3
11948640736000 + 56569705952000 = 40920^3
11948640736000 + 408914669048000 = 74940^3
20087595040000 + 56569705952000 = 42480^3
20087595040000 + 408914669048000 = 75420^3
56569705952000 + 408914669048000 = 77500^3
[2, [9170818419052, 39057725580948, 116836877387052, 349224335725716], 514289757112768]
9170818419052 + 39057725580948 = 36400^3
9170818419052 + 116836877387052 = 50134^3
9170818419052 + 349224335725716 = 71032^3
39057725580948 + 116836877387052 = 53820^3
39057725580948 + 349224335725716 = 72954^3
116836877387052 + 349224335725716 = 77532^3

0 件のコメント:

コメントを投稿

注: コメントを投稿できるのは、このブログのメンバーだけです。