(08-09-2018 10:30 AM)Dieter Wrote: What's this?
Oh, I see. It's just that: \(\frac{113}{355} < \frac{1}{\pi} < \frac{219}{688}\)
Thus we can't simply return \(e, f\):
Code:
if b+d > n:
return e, f
Instead we should check whether \(\frac{a}{b}\) or \(\frac{c}{d}\) is closer to \(\frac{1}{\pi}\).
Here's the output of:
a, b, c, d
Code:
>>> farey(1/pi, 1000)
0 1 1 2
0 1 1 3
1 4 1 3
2 7 1 3
3 10 1 3
4 13 1 3
5 16 1 3
6 19 1 3
7 22 1 3
7 22 8 25
7 22 15 47
7 22 22 69
7 22 29 91
7 22 36 113
7 22 43 135
7 22 50 157
7 22 57 179
7 22 64 201
7 22 71 223
7 22 78 245
7 22 85 267
7 22 92 289
7 22 99 311
7 22 106 333
113 355 106 333
113 355 219 688
(219, 688)
As we know \(\frac{113}{355}\) is closer for a long time. We have to wait until \(\frac{33215}{104348}\) to get closer:
Code:
>>> farey(1/pi, 200000)
0 1 1 2
0 1 1 3
1 4 1 3
2 7 1 3
3 10 1 3
4 13 1 3
5 16 1 3
6 19 1 3
7 22 1 3
7 22 8 25
7 22 15 47
7 22 22 69
7 22 29 91
7 22 36 113
7 22 43 135
7 22 50 157
7 22 57 179
7 22 64 201
7 22 71 223
7 22 78 245
7 22 85 267
7 22 92 289
7 22 99 311
7 22 106 333
113 355 106 333
113 355 219 688
113 355 332 1043
113 355 445 1398
113 355 558 1753
113 355 671 2108
113 355 784 2463
113 355 897 2818
113 355 1010 3173
113 355 1123 3528
113 355 1236 3883
113 355 1349 4238
113 355 1462 4593
113 355 1575 4948
113 355 1688 5303
113 355 1801 5658
113 355 1914 6013
113 355 2027 6368
113 355 2140 6723
113 355 2253 7078
113 355 2366 7433
113 355 2479 7788
113 355 2592 8143
113 355 2705 8498
113 355 2818 8853
113 355 2931 9208
113 355 3044 9563
113 355 3157 9918
113 355 3270 10273
113 355 3383 10628
113 355 3496 10983
113 355 3609 11338
113 355 3722 11693
113 355 3835 12048
113 355 3948 12403
113 355 4061 12758
113 355 4174 13113
113 355 4287 13468
113 355 4400 13823
113 355 4513 14178
113 355 4626 14533
113 355 4739 14888
113 355 4852 15243
113 355 4965 15598
113 355 5078 15953
113 355 5191 16308
113 355 5304 16663
113 355 5417 17018
113 355 5530 17373
113 355 5643 17728
113 355 5756 18083
113 355 5869 18438
113 355 5982 18793
113 355 6095 19148
113 355 6208 19503
113 355 6321 19858
113 355 6434 20213
113 355 6547 20568
113 355 6660 20923
113 355 6773 21278
113 355 6886 21633
113 355 6999 21988
113 355 7112 22343
113 355 7225 22698
113 355 7338 23053
113 355 7451 23408
113 355 7564 23763
113 355 7677 24118
113 355 7790 24473
113 355 7903 24828
113 355 8016 25183
113 355 8129 25538
113 355 8242 25893
113 355 8355 26248
113 355 8468 26603
113 355 8581 26958
113 355 8694 27313
113 355 8807 27668
113 355 8920 28023
113 355 9033 28378
113 355 9146 28733
113 355 9259 29088
113 355 9372 29443
113 355 9485 29798
113 355 9598 30153
113 355 9711 30508
113 355 9824 30863
113 355 9937 31218
113 355 10050 31573
113 355 10163 31928
113 355 10276 32283
113 355 10389 32638
113 355 10502 32993
113 355 10615 33348
113 355 10728 33703
113 355 10841 34058
113 355 10954 34413
113 355 11067 34768
113 355 11180 35123
113 355 11293 35478
113 355 11406 35833
113 355 11519 36188
113 355 11632 36543
113 355 11745 36898
113 355 11858 37253
113 355 11971 37608
113 355 12084 37963
113 355 12197 38318
113 355 12310 38673
113 355 12423 39028
113 355 12536 39383
113 355 12649 39738
113 355 12762 40093
113 355 12875 40448
113 355 12988 40803
113 355 13101 41158
113 355 13214 41513
113 355 13327 41868
113 355 13440 42223
113 355 13553 42578
113 355 13666 42933
113 355 13779 43288
113 355 13892 43643
113 355 14005 43998
113 355 14118 44353
113 355 14231 44708
113 355 14344 45063
113 355 14457 45418
113 355 14570 45773
113 355 14683 46128
113 355 14796 46483
113 355 14909 46838
113 355 15022 47193
113 355 15135 47548
113 355 15248 47903
113 355 15361 48258
113 355 15474 48613
113 355 15587 48968
113 355 15700 49323
113 355 15813 49678
113 355 15926 50033
113 355 16039 50388
113 355 16152 50743
113 355 16265 51098
113 355 16378 51453
113 355 16491 51808
113 355 16604 52163
113 355 16717 52518
113 355 16830 52873
113 355 16943 53228
113 355 17056 53583
113 355 17169 53938
113 355 17282 54293
113 355 17395 54648
113 355 17508 55003
113 355 17621 55358
113 355 17734 55713
113 355 17847 56068
113 355 17960 56423
113 355 18073 56778
113 355 18186 57133
113 355 18299 57488
113 355 18412 57843
113 355 18525 58198
113 355 18638 58553
113 355 18751 58908
113 355 18864 59263
113 355 18977 59618
113 355 19090 59973
113 355 19203 60328
113 355 19316 60683
113 355 19429 61038
113 355 19542 61393
113 355 19655 61748
113 355 19768 62103
113 355 19881 62458
113 355 19994 62813
113 355 20107 63168
113 355 20220 63523
113 355 20333 63878
113 355 20446 64233
113 355 20559 64588
113 355 20672 64943
113 355 20785 65298
113 355 20898 65653
113 355 21011 66008
113 355 21124 66363
113 355 21237 66718
113 355 21350 67073
113 355 21463 67428
113 355 21576 67783
113 355 21689 68138
113 355 21802 68493
113 355 21915 68848
113 355 22028 69203
113 355 22141 69558
113 355 22254 69913
113 355 22367 70268
113 355 22480 70623
113 355 22593 70978
113 355 22706 71333
113 355 22819 71688
113 355 22932 72043
113 355 23045 72398
113 355 23158 72753
113 355 23271 73108
113 355 23384 73463
113 355 23497 73818
113 355 23610 74173
113 355 23723 74528
113 355 23836 74883
113 355 23949 75238
113 355 24062 75593
113 355 24175 75948
113 355 24288 76303
113 355 24401 76658
113 355 24514 77013
113 355 24627 77368
113 355 24740 77723
113 355 24853 78078
113 355 24966 78433
113 355 25079 78788
113 355 25192 79143
113 355 25305 79498
113 355 25418 79853
113 355 25531 80208
113 355 25644 80563
113 355 25757 80918
113 355 25870 81273
113 355 25983 81628
113 355 26096 81983
113 355 26209 82338
113 355 26322 82693
113 355 26435 83048
113 355 26548 83403
113 355 26661 83758
113 355 26774 84113
113 355 26887 84468
113 355 27000 84823
113 355 27113 85178
113 355 27226 85533
113 355 27339 85888
113 355 27452 86243
113 355 27565 86598
113 355 27678 86953
113 355 27791 87308
113 355 27904 87663
113 355 28017 88018
113 355 28130 88373
113 355 28243 88728
113 355 28356 89083
113 355 28469 89438
113 355 28582 89793
113 355 28695 90148
113 355 28808 90503
113 355 28921 90858
113 355 29034 91213
113 355 29147 91568
113 355 29260 91923
113 355 29373 92278
113 355 29486 92633
113 355 29599 92988
113 355 29712 93343
113 355 29825 93698
113 355 29938 94053
113 355 30051 94408
113 355 30164 94763
113 355 30277 95118
113 355 30390 95473
113 355 30503 95828
113 355 30616 96183
113 355 30729 96538
113 355 30842 96893
113 355 30955 97248
113 355 31068 97603
113 355 31181 97958
113 355 31294 98313
113 355 31407 98668
113 355 31520 99023
113 355 31633 99378
113 355 31746 99733
113 355 31859 100088
113 355 31972 100443
113 355 32085 100798
113 355 32198 101153
113 355 32311 101508
113 355 32424 101863
113 355 32537 102218
113 355 32650 102573
113 355 32763 102928
113 355 32876 103283
113 355 32989 103638
113 355 33102 103993
33215 104348 33102 103993
(33215, 104348)
Sorry for the slip.
Thomas
Search
Member List
Calendar
Help
