[break : prime again (4)]で一覧したように、n=5からn=8までは、(0.584)を中心とする前半pと後半pの数がバランス(一致している)。n=9からn=12までは、前半p>後半pの関係である。このずれている箇所のバランスしているpを求める。
---
[2]
n=9,511<(p*38,767(=!p=0.584),p*37)<1023,
===(p*38,(!p=0.584),p*37)
135th=761,n=9,38/38,log2(761+1)=9.573647187,
136th=769,n=9,01/37,log2(769+1)=9.588714636,
===(p*37,p*38)
134th=757,n=9,37/37,log2(757+1)=9.566054038,
135th=761,n=9,01/38,log2(761+1)=9.573647187,
===[detail]
097th=509,n=8,
(511,n=9,border)
098th=521,n=9,1/38,
...
134th=757,n=9,37/38,log2(757+1)=9.566054038,
135th=761,n=9,(135-98+1)/38=38/38,log2(761+1)=9.573647187,
(767=0.584=!p)
136th=769,n=9,1/37,log2(769+1)=9.588714636,
...
172th=1021,n=9,(172-136+1)/37=37/37,
(1023,n=10,border)
173th=1031,n=10,
===
---
[3]先に、n=12を見てみる。
NG)n=12,4095<(p*236,6143(=p=0,584),p*227)<8191,
OK)n=12,4095<(p*236,6143(=p=0,584),p*226)<8191,
===(p*232,p*231) @(236+1+226)=463,463/2=231.5, pair(232,231)/or/pair(231,232),
0796th=6101,232/236,log2(6101+1)=12.57506646,
0797th=6113,233/236,log2(6113+1)=12.57790084,
===(p*231,p*232)
0795th=6091,231/236,log2(6091+1)=12.57270023,
0796th=6101,232/236,log2(6101+1)=12.57506646,
===[detail]
0564th=4093,n=11,
(4095,n=12,border)
0565th=4099,n=12,1/236,
...
0795th=6091,231/236,log2(6091+1)=12.57270023,
0796th=6101,232/236,log2(6101+1)=12.57506646,
0797th=6113,(797-565+1)=233/236,log2(6113+1)=12.57790084,
0798th=6121,log2(6121+1)=12.57978733,
0799th=6131,log2(6131+1)=12.58214198,
0800th=6133,(800-565+1)=236/236,log(6133+1)=12.58261245,
0801th=6143(=0.584),log2(6143+1)=12.5849625,
0802th=6151,n=12,1/226,log2(6151+1)=12.58683979,
...
1027th=8179,(1027-802+1)=226/226,
1028th=8191,n=13,border,
===
---
[4]
n=10,1023<(p*70,1535(=!p=0.584),p*67)<2047,
===(p*69,p*68), @(70+67)=137,137/2=68.5, pair(69,68)/or/pair(68,69),
241th=1523,69/69,log2(1523+1)=10.57364719,
242th=1531,01/68,log2(1531+1)=10.58120058,
===(p*68,p*69)
240th=1511,68/68,log2(1511+1)=10.56224242,
241th=1523,69/69,log2(1523+1)=10.57364719,
===[detail]
172th=1021,n=9,
(1023,n=10,border)
173th=1031,n=10,1/70,
...
240th=1511,68/70,log2(1511+1)=10.56224242,
241th=1523,69/70,log2(1523+1)=10.57364719,
242th=1531,(242-173+1)=70/70,log2(1531+1)=10.58120058,
(1535=0.584=!p)
243th=1543,1/67,log2(1543+1)=10.59245704,
...
309th=2039,n=10,(309-243+1)=67/67,
(2047,n=11,border),
310th=2053,n=11,
---
[5]
n=11,2047<(p*130,3071(=0.584=!p),p*125)<4095,
===(p*128,p*127), @(130+125)=255,255/2=127.5, pair(128,127)/or/pair(127,128),
437th=3049,128/128,lof2(3049+1)=11.57459353,
438th=3061,001/127,log2(3061+1)=11.58025857,
===(p*127,p*128)
436th=3041,127/127,log2(3041+1)=11.57080444,
437th=3049,128/128,lof2(3049+1)=11.57459353,
===[detail]
309th=2039,n=10,
(2047,n=11,border)
310th=2053,n=11,1/130,
...
436th=3041,127/130,log2(3041+1)=11.57080444,
437th=3049,128/130,log2(3049+1)=11.57459353,
438th=3061,129/130,log2(3061+1)=11.58025857,
439th=3067,(439-310+1)=130/130,log2(3067+1)=11.58308277,
(3071=0.584=!p)
440th=3079,1/125,log2(3079+1)=11.58871464,
...
564th=4093,n=11,(564-440+1)=125/125,
(4095,n=12,border),
565th=4099,n=12,
---
[6](0.584)となりの前後のpで、小数部がどう動いているかを見てみる。
当初、(0.584)=(center)として、前半p=なしの場合、nの境界値を使用し、あくまでも、
center=pの場合も、前半pと後半pには含めない流れで処理しようとしていた。が、
n=9以降、(0.584)でバランスしないので、(0.584)もp値の場合、後半pの数に含めた。
処理の流れに整合性がとれないと感じ、仕切り直し。
--
[a]前半pと後半pが一致する場合、そのまま等分で処理する。
[b]前半pと後半pが一致しない場合、cnt(前半p) = cnt(後半p) + 1 とする。
[c](0.584)をcenterとする考えは止め、[a][b]に従い、同pは前半p、後半pのいずれかに含める。
[d]nのborderは含まない。
--
収束点は(0.5)から(0.584)の方にズレていると見ていたが、(0.577)という数字が出てきた。3^-0.5=(0.57735 02692)を想起したが、ネットで(自然定数 0.577)を検索し、オイラー定数を知る。参考値としてグラフに付加した。
===
n,(before-n),(after-n),(Euler's),
3,0.584962501,0.807354922,0.57721566,
4,0.321928095,0.584962501,0.57721566,
5,0.584962501,0.754887502,0.57721566,
6,0.491853096,0.614709844,0.57721566
7,0.584962501,0.599912842,0.57721566,
8,0.584962501,0.607330314,0.57721566,
9,0.573647187,0.588714636,0.57721566,
10,0.57364719,0.58120058,0.57721566,
11,0.57459353,0.58025857,0.57721566,
12,0.57506646,0.57790084,0.57721566,
13,0.58072965,0.58096514,0.57721566,
14,0.57305541,0.57341051,0.57721566,
15,0.57672053,0.57707472,0.57721566,
16,0.57713374,0.57722228,0.57721566,
17,0.57798932,0.57812204,0.57721566,
18,0.57909492,0.57912439,0.57721566,
19,0.57912807,0.57913176,0.57721566,
===
[graph,now]
@上記のn=10からn=19までをみると、オイラー定数(列E)との関係がみえる、、、

===
[graph,previous,[6]first idea]n=8まで(0.584)=centerと見ていたが、形的に何か捨てがたいので、記録を残す。[6](再々)再考要か?
===[detail]
---
n=3,(p*1,p*1),
before)p=11,log2(11+1)=3.584962501,
after )p=13,log2(13+1)=3.807354922,
--
n=4,(p*2,p*2),
before)p=19,log2(19+1)=4.321928095,
after )p=23,log2(23+1)=4.584962501,
--
n=5,(p*4,p*3),
before)p=47,log2(47+1)=5.584962501,
after )p=53,log2(53+1)=5.754887502,
--
n=6,(p*6,p*6),
before)p=89, log2(89+1)=6.491853096,
after )p=97, log2(97+1)=6.614709844,
--
n=7,(p*12,p*11),
before)p=191,log2(191+1)=7.584962501,
after )p=193,log2(193+1)=7.599912842,
--
n=8,(p*22,p*21),
before)p=383,log2(383+1)=8.584962501,
after )p=389,log2(389+1)=8.607330314,
--
n=9,(p*38,p*37),
before)p=761, log2(761+1)=9.573647187,
after )p=769, log2(769+1)=9.588714636,
--
n=10,(p*69,p*68),
before)p=1523,log2(1523+1)=10.57364719,
after )p=1531,log2(1531+1)=10.58120058,
--
n=11,(p*128,p*127),
before)p=3049,log2(3049+1)=11.57459353,
after )p=3061,log2(3061+1)=11.58025857,
--
n=12,(p*232,p*231),
before)p=6101,log2(6101+1)=12.57506646,
after )p=6113,log2(6113+1)=12.57790084,
===
n=12以降、オイラー定数に収束するか?
===
n=13,(p*436,p*436),
before)p=12251,436/872,log2(12251+1)=13.58072965,
after )p=12253,437/872,log2(12253+1)=13.58096514,
--
n=14,(p*806,p*806),
before)p=24373,806/1612,log2(24373+1)=14.57305541,
after )p=24379,807/1612,log2(24379+1)=14.57341051,
---
n=15,(p*1515,p*1515),
before)p=48871,1515/3030,log2(48871+1)=15.57672053,
after )p=48883,1516/3030,log2(48883+1)=15.57707472,
---
n=16,(p*2854,p*2854),
before)p=97771,2854/6708,log2(97771+1)=16.57713374,
after )p=97777,2855/5708,log2(97777+1)=16.57722228,
---
n=17,(p*5375,p*5374),
before)p=195659,5375/10749,log2(195659+1)=17.57798932,
after )p=195677,5376/10749,log2(195677+1)=17.57812204,
---
n=18,(p*10195,p*10194),
before)p=391619,10195/20389,log2(391619+1)=18.57909492,
after )p=391627,10196/20389,log2(391627+1)=18.57912439,
---
n=19,(p*19318,p*19317),
before)p=783257,19318/38635,log2(783257+1)=19.57912807,
after )p=783259,19319/38635,log2(783259+1)=19.57913176,
===
[working memos]
---
n=13,8191<()<16383,
===(p*436,p*436), (1900-1029+1)=872,872/2=436, pair(p*436,p*436),
1464th=12251,(1464-1029+1)=436/872,log2(12251+1)=13.58072965,
1465th=12253,437/872,log2(12253+1)=13.58096514,
===[detail]
1027th=8179,n=12,
1028th=8191,border,n=13,
1029th=8209,n=13,1/872,
...
1462th=12239,434/872,log2(12239+1)=13.57931584,
1463th=12241,435/872,log2(12241+1)=13.57955165,
1464th=12251,(1464-1029+1)=436/872,log2(12251+1)=13.58072965,
1465th=12253,437/872,log2(12253+1)=13.58096514,
...
1900th=16381,n=13,872/872,
(16383,n=14,border),
1901th=16411,n=14,
===
---
n=14,16383<()<32767,
===(p*806,p*806), (3512-1901+1)=1612,1612/2=806, parit(p*806,p*806),
2706th=24373,806/1612,log2(24373+1)=14.57305541,
2707th=24379,807/1612,log2(24379+1)=14.57341051,
===[detail]
1900th=16381,n=13,
(16383,n=14,border),
1901th=16411,n=14,1/1612,
...
2706th=24373,806/1612,log2(24373+1)=14.57305541,
2707th=24379,807/1612,log2(24379+1)=14.57341051,
2708th=24391,(2708-1901+1)=808/1612,log2(24391+1)=14.57412044,
...
3512th=32749,n=14,1612/1612,
(32767,n=15,border)
3513th=32771,n=15,
===
---
n=15,32767<()<65535,
===(p*1515,p*1515), (6542-3513+1)=3030,3030/2=1515, pair(p*1515,p*1515),
===[detail]
3512th=32749,n=14,
(32767,n=15,border)
3513th=32771,n=15,
...
5027th=48871,1515/3030,log2(48871+1)=15.57672053,
5028th=48883,1516/3030,log2(48883+1)=15.57707472,
...
6542th=65521,n=15,3030/3030,
(65535,n=16,border),
6543th=65537,n=16,
===
---
n=16,65535<()<131071,
===(p*2854,p*2854), (12250-6543+1)=5708,5708/2=2854, pair(p*2854,p*2854),
9396th=97771,2854/6708,log2(97771+1)=16.57713374,
9397th=97777,2855/5708,log2(97777+1)=16.57722228,
===[detail]
6542th=65521,n=15,
(65535,n=16,border),
6543th=65537,n=16,
...
9396th=97771,2854/6708,log2(97771+1)=16.57713374,
9397th=97777,2855/5708,log2(97777+1)=16.57722228,
9398th=97787,(9398-6543+1)=2856/5708,log2(97787+1)=16.57736982,
...
12250th=131063,n=16,5708/5708,
12251th=131071,border,n=17,
12252th=131101,1/
===
---
n=17,131071<()<262143,
===(p*5375,p*5374), (23000-12252+1)=10749,10749/2=5374.5, pair(5375,5374)/or/
17626th=195659,5375/10749,log2(195659+1)=17.57798932,
17627th=195677,5376/10749,log2(195677+1)=17.57812204,
===[detail]
12250th=131063,n=16,
12251th=131071,border,n=17,
12252th=131101,1/10749,
...
17626th=195659,(17626-12252+1)=5375/10749,log2(195659+1)=17.57798932,
17627th=195677,5376/10749,log2(195677+1)=17.57812204,
...
23000th=262139,n=17,10749/10749,
(262143,n=18,border),
23001th=262147,n=18,
===
---
n=18,262143<()<524287,
===(p*10195,p*10194), (43389-23001+1)=20389,20389/2=10194.5, pair(p*10195,p*10194),
33195th=391619,10195/20389,log2(391619+1)=18.57909492,
33196th=391627,10196/20389,log2(391627+1)=18.57912439,
===
[detail]
23000th=262139,n=17,
(262143,n=18,border),
23001th=262147,n=18,1/
...
33195th=391619,(33195-23001+1)=10195/20389,log2(391619+1)=18.57909492,
33196th=391627,10196/20389,log2(391627+1)=18.57912439,
...
43389th=524269,n=18,
43390th=524287,border,n=19,
43391th=524309,n=19,
===
---
n=19,524287<()<1048575,
===(p*19318,p*19317), (82025-43391+1)=38635,38635/2=19317.5, pair(p*19318,p*19317)/or/,
62708th=783257,19318/38635,log2(783257+1)=19.57912807,
62709th=783259,19319/38635,log2(783259+1)=19.57913176,
===[detail]
43390th=524287,border,n=19,
43391th=524309,n=19,
...
62708th=783257,19318/38635,log2(783257+1)=19.57912807,
62709th=783259,19319/38635,log2(783259+1)=19.57913176,
62710th=783269,(62710-43391+1)=19320/38635,
...
82025th=1048573,n=19,
(1048575,n=20,border),
82026th=1048583,n=20,
===
---
end.