Dmitriy40, да не убрал. У меня в 2-х местах именно разрешённые остатки, так проще формулы.
Я имел в виду, что к формуле (2) поправочных коэффициентов не надо.
Возможно, но я не склонен в таких вопросах верить на слово. Покажите, пожалуйста, как те формулы, что у меня в программе, следуют из формулы (2).
Код:
v=[0,6,12,30,42,72,90,96,120,126,132,156,162,180,210,222,240,246,252];
MC = 2^(#v-1);
forprime(p=3,#v, m=p-#Set(-v%p);
MC *= m/p/(1 - 1/p)^#v);
BC = prodeulerrat(( p^#v - #v*p^(#v-1) )/(p-1)^#v, 1, nextprime(#v+1));
PM = 1;
forprime(p=#v+1,v[#v]/2, m=p-#Set(-v%p);
PM *= m/(p - #v));
C = MC * BC * PM;
for(po = 1, 100,
Li2 = intnum(t=2, 10^po, C/log(t)^#v);
Li2 это и есть количество всех кортежей. Оно во 2-м столбце:
Код:
19 [0, 6, 12, 30, 42, 72, 90, 96, 120, 126, 132, 156, 162, 180, 210, 222, 240, 246, 252]
10^ HL-1 24/11 Posl/Pred Fact Pogresh
37 7.482339 e10 0.365
38 2.366044 e11 3.16
39 2.340443 e11 0.989
40 1.174290 e12 5.02
41 6.047189 e12 5.15
42 3.778631 e13 6.25
43 2.388353 e14 6.32
44 1.534438 e15 6.42
45 9.957602 e15 6.49
46 6.525521 e16 6.55
47 4.315983 e17 6.61
48 2.879911 e18 6.67
49 1.937988 e19 6.73
50 1.314755 e20 6.78
51 8.989148 e20 6.84
52 6.192116 e21 6.89
53 4.296164 e22 6.94
54 3.001411 e23 6.99
55 2.110869 e24 7.03
56 1.494110 e25 7.08
57 1.064117 e26 7.12
58 7.624074 e26 7.16
59 5.493968 e27 7.21
60 3.981084 e28 7.25