29/04/13 8071 Богородский
|
Последний раз редактировалось Yadryara 22.07.2024, 13:22, всего редактировалось 1 раз.
И ещё есть эффект чистоплотности. Его тоже интересно проверить. Проверил. Опять подтвердился. Код: 10^ OP 0 1 2 3 4 5 6
5 0.120 0.343 0.346 0.155 0.033 0.003 0.000 5 0.109 0.174 0.478 0.130 0.087 0.000 0.022
6 0.185 0.393 0.299 0.104 0.017 0.001 0.000 6 0.209 0.343 0.313 0.097 0.030 0.000 0.007
7 0.247 0.415 0.255 0.072 0.010 0.001 0.000 7 0.208 0.435 0.270 0.062 0.022 0.002 0.002
8 0.303 0.422 0.217 0.052 0.006 0.000 0.000 8 0.300 0.429 0.213 0.048 0.009 0.001 0.000
9 0.352 0.420 0.185 0.038 0.004 0.000 0.000 9 0.360 0.426 0.175 0.035 0.004 0.000 0.000
10 0.396 0.412 0.160 0.029 0.003 0.000 0.000 10 0.417 0.406 0.150 0.025 0.002 0.000 0.000
11 0.435 0.402 0.139 0.023 0.002 0.000 0.000 11 0.460 0.393 0.126 0.019 0.002 0.000 0.000
В верхней строчке каждой пары теоретические доли, в нижней — фактические. Для 5-36 эффект малозаметен: реальная доля чистых больше всего лишь в 1.06 раза. А вот для 7-108-1 уже в 1.24 раза: Код: 0 1 2 3 4 5 6 7 8
8 0.004 0.030 0.095 0.180 0.230 0.210 0.141 0.071 0.027 8 0.006 0.025 0.067 0.172 0.209 0.239 0.129 0.055 0.055
9 0.009 0.053 0.138 0.220 0.235 0.180 0.102 0.044 0.014 9 0.006 0.045 0.159 0.238 0.229 0.168 0.078 0.044 0.020
10 0.016 0.079 0.179 0.243 0.225 0.149 0.073 0.027 0.008 10 0.014 0.087 0.186 0.255 0.217 0.149 0.057 0.023 0.009
11 0.025 0.108 0.213 0.254 0.206 0.120 0.052 0.017 0.004 11 0.030 0.121 0.227 0.257 0.197 0.109 0.041 0.013 0.004
12 0.036 0.136 0.240 0.256 0.185 0.096 0.037 0.011 0.002 12 0.043 0.157 0.257 0.251 0.169 0.083 0.029 0.008 0.002
13 0.048 0.165 0.261 0.251 0.164 0.077 0.027 0.007 0.001 13 0.060 0.189 0.273 0.241 0.146 0.064 0.020 0.005 0.001
14 0.061 0.191 0.276 0.242 0.144 0.062 0.020 0.005 0.001 14 0.076 0.217 0.285 0.229 0.124 0.049 0.015 0.003 0.001
15 0.075 0.216 0.286 0.230 0.126 0.049 0.014 0.003 0.001 15 0.093 0.243 0.292 0.214 0.107 0.039 0.010 0.002 0.000
Напомню формулировку: сначала фактические доли всегда выигрывают, затем вблизи самого популярного варианта происходит перелом и уже только проигрывают. Три нижних левых выигрывают у верхних, остальные проигрывают. И именно для чистых выигрыш наибольший: . Программа: (PARI)
Код: {print();
cm = 7; k = 8;
f=vector(15);fo=vector(15);
f[8] = [ 1, 4, 11, 28, 34, 39, 21, 9, 9, 2 ]; f[9] = [ 4, 29, 102, 153, 147, 108, 50, 28, 13, 3 ]; f[10] = [ 41, 250, 536, 736, 627, 431, 166, 66, 26, 3 ]; f[11] = [ 419, 1671, 3130, 3546, 2722, 1501, 565, 186, 56, 12 ]; f[12] = [ 3087, 11269, 18443, 17949, 12140, 5932, 2078, 593, 125, 21 ]; f[13] = [ 24030, 75583, 108959, 96375, 58305, 25391, 8181, 2014, 382, 48 ]; f[14] = [ 177056, 502986, 661258, 530930, 288540, 114539, 33688, 7580, 1296, 170 ]; f[15] = [ 1299308, 3400588, 4082114, 2996106, 1497968, 540945, 145706, 29900, 4631, 549 ];
fo[1] = 0; fo[2] = 2; fo[3] = 4; fo[4] = 6; fo[5] = 10; fo[6] = 24; fo[7] = 57; fo[8] = 163; fo[9] = 642; fo[10] = 2887; fo[11] = 13813; fo[12] = 71642; fo[13] = 399280; fo[14] = 2318066; fo[15] = 13997875;
mor=[ 6, 210, 200560490130, 2.305567963946 E36, 6.107692946593 E127, 1.959034064500 E415, 1.952288231513 E1329, 5.949067958000 E4297, 4.305235385595 E13620, 6.523576146660 E43292, 4.458746156594 E136987, 1.470555276084 E433636, 1.753664729516 E1372340, 1.672547671403 E4340851, 4.012967119696 E13731288, 2.540148384222 E43424119, 2.048235031377 E137328734 ];
predok = 53;
vc=[554408, 115745660, 7407543216, 216318061332, 3429058539864, 32521729709928, 196088606111368, 781359053318046, 2106299320426628, 3885606017883462, 4914156988792810, 4230571480579892, 2435923197894882, 907258494170660, 205839739312028, 25609408974336, 1489970557864, 36686447616];
c=vc*1.0;
print1(" ");
for(i=1,9, printf(" %d", i-1)); print();print();
forprime(p=predok+1,sqrt(1e16),
while(p^2>10^k,
if(k<10,print1(k," "),print1(k," ")); for(i=1,9, printf(" %0.3f ", c[i]/vecsum(c)));
print();
if(k<10,print1(k," "),print1(k," ")); for(i=1,9, printf(" %0.3f ", f[k][i]/fo[k]));
print(); print();
k++);
if(c[1]>1e160000000,
kpon++; print(); print(kpon); print();
for(i=1,#c, c[i]=c[i]/(1e160000000);
print(c[i]);
); print(); );
for(i=1,#c-1, c[i]=c[i]*(p-i-cm+1)+c[i+1]*i);
c[#c]*=p-#c-cm+1);
print();
}quit; Об остальном позже.
|
|