前回からちょっとだけ進捗がありました。
上図のような半径をrとする半円があり、
∠αの弦をa、∠βの弦をb、直径をc、
mα＋nα＝180˚, m≧n＞0, m, n∈Nを満たす。
各(m, n)において、(a, b, c)が既約な正の
整数となるようなものを求めよ。
という自作問題をプログラミングしていまして、
| m | n | a | b | c |
| 1 | 1 | 3 | 4 | 5 |
| 1 | 1 | 4 | 3 | 5 |
| 1 | 1 | 5 | 12 | 13 |
| 1 | 1 | 12 | 5 | 13 |
| 1 | 1 | 8 | 15 | 17 |
| 1 | 1 | 15 | 8 | 17 |
| 1 | 1 | 7 | 24 | 25 |
| 1 | 1 | 24 | 7 | 25 |
| 1 | 1 | 20 | 21 | 29 |
| 1 | 1 | 21 | 20 | 29 |
| 2 | 1 | 2 | 7 | 8 |
| 2 | 1 | 3 | 7 | 9 |
| 2 | 1 | 6 | 1 | 9 |
| 2 | 1 | 3 | 17 | 18 |
| 2 | 1 | 5 | 23 | 25 |
| 2 | 1 | 10 | 17 | 25 |
| 2 | 1 | 15 | 7 | 25 |
| 2 | 1 | 4 | 31 | 32 |
| 2 | 1 | 12 | 23 | 32 |
| 2 | 1 | 20 | 7 | 32 |
| 3 | 1 | 845 | 828 | 2197 |
| 3 | 1 | 2312 | 495 | 4913 |
| 3 | 1 | 4375 | 10296 | 15625 |
| 3 | 1 | 16428 | 27755 | 50653 |
| 3 | 1 | 15129 | 54280 | 68921 |
| 3 | 1 | 40931 | 194220 | 226981 |
| 3 | 1 | 67600 | 201663 | 274625 |
| 3 | 1 | 93925 | 550116 | 614125 |
| 3 | 1 | 260100 | 157157 | 614125 |
| 3 | 1 | 308919 | 146960 | 704969 |
| 4 | 1 | 8 | 17 | 32 |
| 4 | 1 | 27 | 17 | 81 |
| 4 | 1 | 27 | 127 | 162 |
| 4 | 1 | 64 | 449 | 512 |
| 4 | 1 | 192 | 17 | 512 |
| 4 | 1 | 125 | 433 | 625 |
| 4 | 1 | 125 | 1151 | 1250 |
| 4 | 1 | 375 | 431 | 1250 |
| 4 | 1 | 343 | 2017 | 2401 |
| 4 | 1 | 216 | 2449 | 2592 |
| 5 | 1 | 2734375 | 1476984 | 9765625 |
| 5 | 1 | 25431849 | 51872200 | 115856201 |
| 5 | 1 | 152304251 | 520632300 | 844596301 |
| 5 | 1 | 285610000 | 372960063 | 1160290625 |
| 5 | 1 | 678608125 | 3192445284 | 4437053125 |
| 5 | 1 | 2081208020 | 5707904499 | 10510100501 |
| 5 | 1 | 2445710415 | 14490678832 | 18424351793 |
| 5 | 1 | 7514860625 | 53348005584 | 64097340625 |
| 5 | 1 | 10609215000 | 43190896463 | 64097340625 |
| 5 | 1 | 46578742132 | 6863126325 | 154963892093 |
| 6 | 1 | 32 | 7 | 128 |
| 6 | 1 | 243 | 782 | 1458 |
| 6 | 1 | 1024 | 5983 | 8192 |
| 6 | 1 | 3125 | 5543 | 15625 |
| 6 | 1 | 3125 | 25774 | 31250 |
| 6 | 1 | 7776 | 81863 | 93312 |
| 6 | 1 | 16807 | 76751 | 117649 |
| 6 | 1 | 16807 | 213982 | 235298 |
| 6 | 1 | 50421 | 63902 | 235298 |
| 6 | 1 | 32768 | 487807 | 524288 |
| 7 | 1 | 42750938169 | 4205436520 | 194754273881 |
| 7 | 1 | 566724117971 | 933420304380 | 3142742836021 |
| 7 | 1 | 4902943703125 | 15256342319316 | 32057708828125 |
| 7 | 1 | 21230403012020 | 18712692999099 | 107213535210701 |
| 7 | 1 | 31229276289135 | 140277922759408 | 235260548044817 |
| 7 | 1 | 157999944640625 | 916846383342096 | 1347646586640625 |
| 7 | 1 | 223058745375000 | 533250316479023 | 1347646586640625 |
| 7 | 1 | 668074738214539 | 4716260663909940 | 6364290927201661 |
| 7 | 1 | 1636648392655612 | 6237995487261915 | 11514990476898413 |
| 7 | 1 | 2446656141033741 | 20243147384053660 | 25748143198497941 |
| 8 | 1 | 2187 | 3007 | 13122 |
| 8 | 1 | 16384 | 70529 | 131072 |
| 8 | 1 | 78125 | 543551 | 781250 |
| 8 | 1 | 279936 | 2638369 | 3359232 |
| 8 | 1 | 823543 | 2371777 | 5764801 |
| 8 | 1 | 823543 | 9694847 | 11529602 |
| 8 | 1 | 2097152 | 29441537 | 33554432 |
| 8 | 1 | 6291456 | 2074529 | 33554432 |
| 8 | 1 | 4782969 | 27069761 | 43046721 |
| 8 | 1 | 4782969 | 77720959 | 86093442 |
| 9 | 1 | 35423768255078125 | 43493220003156084 | 231616946283203125 |
| 9 | 1 | 398766628935964815 | 1093500518290131952 | 3004041937984268273 |
| 9 | 1 | 3321948836069140625 | 13910996789480983824 | 28334269484119140625 |
| 9 | 1 | 4689810121509375000 | 2102004302995808783 | 28334269484119140625 |
| 9 | 1 | 21886796498646512179 | 121874910240833598420 | 208500535066053616021 |
| 9 | 1 | 63516687470571646108 | 126614861389951203555 | 446885265417950510117 |
| 9 | 1 | 119497132584228944181 | 823817094582142462780 | 1257565061957837936381 |
| 9 | 1 | 608996735987235225632 | 2115333500682367703295 | 4891005035897482905857 |
| 9 | 1 | 559364406039571484375 | 4572130042340783039304 | 6444850765238541015625 |
| 9 | 1 | 2303004088924671288025 | 21684566268331317474552 | 28833611193336884526073 |
| 3 | 2 | 2 | 9 | 16 |
| 3 | 2 | 21 | 10 | 54 |
| 3 | 2 | 35 | 44 | 125 |
| 3 | 2 | 7 | 235 | 343 |
| 3 | 2 | 161 | 18 | 343 |
| 3 | 2 | 161 | 286 | 686 |
| 3 | 2 | 153 | 329 | 729 |
| 3 | 2 | 328 | 299 | 1024 |
| 3 | 2 | 63 | 962 | 1458 |
| 3 | 2 | 890 | 187 | 2000 |
| 5 | 2 | 27 | 118 | 243 |
| 5 | 2 | 448 | 475 | 2048 |
| 5 | 2 | 875 | 237 | 3125 |
| 5 | 2 | 125 | 4193 | 6250 |
| 5 | 2 | 5831 | 11511 | 33614 |
| 5 | 2 | 3584 | 39589 | 65536 |
| 5 | 2 | 29889 | 17149 | 118098 |
| 5 | 2 | 30613 | 48713 | 161051 |
| 5 | 2 | 31000 | 77207 | 200000 |
| 5 | 2 | 22627 | 184425 | 322102 |
| 7 | 2 | 32 | 87 | 256 |
| 7 | 2 | 16807 | 539285 | 823543 |
| 7 | 2 | 1003833 | 216727 | 4782969 |
| 7 | 2 | 413343 | 5667493 | 9565938 |
| 7 | 2 | 1127357 | 10725128 | 19487171 |
| 7 | 2 | 7569397 | 3941297 | 38974342 |
| 7 | 2 | 11510083 | 8739370 | 62748517 |
| 7 | 2 | 248832 | 50054273 | 71663616 |
| 7 | 2 | 8539739 | 65279611 | 125497034 |
| 7 | 2 | 12909375 | 85000619 | 170859375 |
| 9 | 2 | 2187 | 5522 | 19683 |
| 9 | 2 | 78125 | 2502682 | 3906250 |
| 9 | 2 | 14000231 | 65934 | 80707214 |
| 9 | 2 | 14680064 | 137823851 | 268435456 |
| 9 | 2 | 310000000 | 169944073 | 2000000000 |
| 9 | 2 | 331281907 | 2131650810 | 4715895382 |
| 9 | 2 | 439239619 | 5978595177 | 10604499373 |
| 9 | 2 | 3074677333 | 2764642646 | 21208998746 |
| 9 | 2 | 3267818624 | 17194374007 | 41322093568 |
| 9 | 2 | 19058917376 | 21809310453 | 137438953472 |
| 4 | 3 | 28 | 223 | 512 |
| 4 | 3 | 270 | 289 | 1250 |
| 4 | 3 | 858 | 191 | 2592 |
| 4 | 3 | 936 | 2143 | 6561 |
| 4 | 3 | 4070 | 2159 | 14641 |
| 4 | 3 | 10764 | 239 | 28561 |
| 4 | 3 | 195 | 25087 | 50625 |
| 4 | 3 | 6578 | 20641 | 57122 |
| 4 | 3 | 12950 | 22543 | 76832 |
| 4 | 3 | 33872 | 22831 | 131072 |
| 8 | 3 | 2430 | 353 | 13122 |
| 8 | 3 | 302526 | 2157503 | 5764801 |
| 8 | 3 | 18700000 | 54097759 | 200000000 |
| 8 | 3 | 6119938 | 200073121 | 428717762 |
| 8 | 3 | 93937129 | 174695903 | 815730721 |
| 8 | 3 | 1104150528 | 1536790303 | 8589934592 |
| 8 | 3 | 1070572178 | 4372390753 | 13951514882 |
| 8 | 3 | 2337437456 | 2628945119 | 16983563041 |
| 8 | 3 | 644342688 | 9499728863 | 22039921152 |
| 8 | 3 | 15826803872 | 15065687551 | 109751747072 |
| 5 | 4 | 1673 | 10907 | 33614 |
| 5 | 4 | 9736 | 13425 | 65536 |
| 5 | 4 | 27063 | 12223 | 118098 |
| 5 | 4 | 59210 | 3401 | 200000 |
| 5 | 4 | 77295 | 199066 | 759375 |
| 5 | 4 | 270929 | 215984 | 1419857 |
| 5 | 4 | 653923 | 145062 | 2476099 |
| 5 | 4 | 75220 | 2362043 | 6400000 |
| 5 | 4 | 876854 | 2910139 | 10307264 |
| 5 | 4 | 1517977 | 3121757 | 12872686 |
| 7 | 4 | 5875 | 20158 | 78125 |
| 7 | 4 | 82728 | 74335 | 559872 |
| 7 | 4 | 164297 | 33933 | 823543 |
| 7 | 4 | 257337 | 3240713 | 9565938 |
| 7 | 4 | 4473491 | 7413831 | 38974342 |
| 7 | 4 | 525083 | 23161315 | 62748517 |
| 7 | 4 | 22042501 | 10494781 | 125497034 |
| 7 | 4 | 9304904 | 65405329 | 210827008 |
| 7 | 4 | 75225375 | 1461083 | 341718750 |
| 7 | 4 | 54980608 | 113578907 | 536870912 |
| 9 | 4 | 17408 | 163917 | 524288 |
| 9 | 4 | 8857289 | 11683330 | 80707214 |
| 9 | 4 | 280100000 | 152945093 | 2000000000 |
| 9 | 4 | 193422251 | 486749276 | 2357947691 |
| 9 | 4 | 3314532611 | 840140586 | 21208998746 |
| 9 | 4 | 5308790625 | 18068472334 | 76886718750 |
| 9 | 4 | 15131416049 | 12407846922 | 118587876497 |
| 9 | 4 | 22868393984 | 2277731797 | 137438953472 |
| 9 | 4 | 339345823 | 90057353906 | 237175752994 |
| 9 | 4 | 36825790752 | 72405415511 | 396718580736 |
| 6 | 5 | 203511 | 860078 | 3543122 |
| 6 | 5 | 728916 | 994607 | 5971968 |
| 6 | 5 | 1724125 | 951118 | 9653618 |
| 6 | 5 | 3440178 | 565823 | 15059072 |
| 6 | 5 | 1850100 | 24382921 | 85766121 |
| 7 | 5 | 5925 | 16124 | 78125 |
| 8 | 5 | 92232 | 896737 | 3359232 |
| 8 | 5 | 577612 | 881857 | 5764801 |
| 8 | 5 | 5165568 | 2241857 | 33554432 |
| 8 | 5 | 108916275 | 335727167 | 1631461442 |
m=5, n=3
がなかなか進まないので、他に手を付けました。
cに関しては、
初期値
c=xm
インクリメント
x++;
c=xm
と
c=2xm
だけを調査するという方針にしたので、cはxの増加で簡単に大きな値になる。
しかし、aに関しては、
初期値
a=xm-n
インクリメント
a+=xm-n
であるため、m-nが小さいと、cは極端に大きくなるが、aは小さいために、随分と時間がかかってしまう。
というわけで、とりあえずC言語で精度をちょっと緩めて、ある程度近いものを拾って、そこで抽出されたものを、多倍長電卓で精密検査するという方針で、解を求めた。
なので、m=8, n=3以降、解の個数が10個に達していないものもある。
とりあえず、解を190個ほど見つけたわけで、この見つかった解から、新しい法則などが見つかれば、それをつてに効率よく探すことが出来るのだろうか。
ではでは