Форум программистов, компьютерный форум, киберфорум
C# для начинающих
Форум программистов и сисадминов Киберфорум > Форум программистов > .NET > C# .NET > C# для начинающих
Войти
Регистрация
Восстановить пароль
Блоги Сообщество Поиск Заказать работу  
 
Рейтинг 4.80/5: Рейтинг темы: голосов - 5, средняя оценка - 4.80
1 / 1 / 0
Регистрация: 22.02.2016
Сообщений: 9

Долгое выполнение функции

05.09.2016, 12:12. Показов 1021. Ответов 13
Метки нет (Все метки)
Студворк — интернет-сервис помощи студентам
Доброго времени суток. Такой вопрос.
Есть две функции расчета:
C#
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24

            double p1;
            double p2;
            p1 = 0;//Convert.ToDouble(numericUpDown3.Value);
            p2 = 0;//Convert.ToDouble(numericUpDown4.Value);                                  
            double res=0;
            double t = PA(p1,p2);
            if (g1<=4&g2<=4) res = t*PBPsT(g1+1,g2)+(1-t)*PBPsT(g1,g2+1);
            if (g1==6&g2==6) res = PAPgt(p1,p2);
            if (g1==6&g2==5) res = t*1+(1-t)*PBPsT(6,6);
            if (g1==5&g2==6) res = t*PBPsT(6,6)+(1-t)*0;
            if (g1==5&g2==5) res = t*PBPsT(6,5)+(1-t)*PBPsT(5,6);
            if (g1==5&g2==4) res = t*1+(1-t)*PBPsT(5,5);
            if (g1==5&g2==3) res = t*1+(1-t)*PBPsT(5,4);
            if (g1==5&g2==2) res = t*1+(1-t)*PBPsT(5,3);
            if (g1==5&g2==1) res = t*1+(1-t)*PBPsT(5,2);
            if (g1==5&g2==0) res = t*1+(1-t)*PBPsT(5,1);        
            if (g1==4&g2==5) res = t*PBPsT(5,5)+(1-t)*0;            
            if (g1==3&g2==5) res = t*PBPsT(4,5)+(1-t)*0;
            if (g1==2&g2==5) res = t*PBPsT(3,5)+(1-t)*0;
            if (g1==1&g2==5) res = t*PBPsT(2,5)+(1-t)*0;
            if (g1==0&g2==5) res = t*PBPsT(1,5)+(1-t)*0;        
    
 
            return res;
C#
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27

        public double PBPsT(int g1, int g2)  
        {
            double p1;
            double p2;  
            p1 = 0;//Convert.ToDouble(numericUpDown3.Value);
            p2 = 0;//Convert.ToDouble(numericUpDown4.Value);                                  
            double res=0;
            double t = PB(p1,p2);
            if (g1<=4&g2<=4) res = t*PAPsT(g1+1,g2)+(1-t)*PAPsT(g1,g2+1);
            if (g1==6&g2==6) res = PAPgt(p1,p2);
            if (g1==6&g2==5) res = t*1+(1-t)*PAPsT(6,6);
            if (g1==5&g2==6) res = t*PAPsT(6,6)+(1-t)*0;
            if (g1==5&g2==5) res = t*PAPsT(6,5)+(1-t)*PAPsT(5,6);
            if (g1==5&g2==4) res = t*1+(1-t)*PAPsT(5,5);
            if (g1==5&g2==3) res = t*1+(1-t)*PAPsT(5,4);
            if (g1==5&g2==2) res = t*1+(1-t)*PAPsT(5,3);
            if (g1==5&g2==1) res = t*1+(1-t)*PAPsT(5,2);
            if (g1==5&g2==0) res = t*1+(1-t)*PAPsT(5,1);
            
            if (g1==4&g2==5) res = t*PAPsT(5,5)+(1-t)*0;            
            if (g1==3&g2==5) res = t*PAPsT(4,5)+(1-t)*0;
            if (g1==2&g2==5) res = t*PAPsT(3,5)+(1-t)*0;
            if (g1==1&g2==5) res = t*PAPsT(2,5)+(1-t)*0;
            if (g1==0&g2==5) res = t*PAPsT(1,5)+(1-t)*0;        
 
            return res;         
        }
Так вот: если щелкнуть по кнопке:
C#
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21

        void Button3Click(object sender, EventArgs e)
        {   
            double win_set;
            double p1;
            double p2;
            int g1;
            int g2;
            p1 = Convert.ToDouble(numericUpDown3.Value);
            p2 = Convert.ToDouble(numericUpDown4.Value); 
            g1 = Convert.ToInt16(textBox85.Text);
            g2 = Convert.ToInt16(textBox86.Text);
            if (textBox85.Text=="6"&textBox86.Text=="6")
                win_set = PAPgt(p1,p2);
            else
                win_set = PA(p1,p2)*PBPsT(g1+1,g2)+(1-PA(p1,p2))*PBPsT(g1,g2+1);
            var watch = System.Diagnostics.Stopwatch.StartNew();
            MessageBox.Show(Convert.ToString(win_set));
            watch.Stop();
            var elapsedMs = watch.ElapsedMilliseconds;
            MessageBox.Show(Convert.ToString(elapsedMs/1000));
        }
при g1 и g2>=3 функция PBPst выполняется достаточно быстро, а вот при меньших значениях очень долго выполнение происходит.
С чем это связано, не подскажете?
PA и PAPqT выполняются быстро. Проверено.
0
cpp_developer
Эксперт
20123 / 5690 / 1417
Регистрация: 09.04.2010
Сообщений: 22,546
Блог
 05.09.2016, 12:12
Ответы с готовыми решениями:

Долгое выполнение запроса
Суть такая update и select делает перебор всей базы, хотя выборка по индексам стоит лог выдает это # Query_time: 0.220863 Lock_time:...

Долгое выполнение пролога
Добрый день! Есть ИМ, товара несколько десятков тысяч. Чтобы ускорить его, был написан скрипт, который рано утром проходится по всем...

Долгое выполнение exec()
Всем привет, есть скрипт файл a.php в нем делаем различные выборы , выбор записываем в сессию после нажатия на кнопку идет...

13
Администратор
Эксперт .NET
 Аватар для OwenGlendower
18298 / 14222 / 5368
Регистрация: 17.03.2014
Сообщений: 28,900
Записей в блоге: 1
05.09.2016, 21:45
RealFlush, из-за рекурсии наверное.
1
1 / 1 / 0
Регистрация: 22.02.2016
Сообщений: 9
06.09.2016, 10:16  [ТС]
Да, скорее всего. А есть какой-нибудь вариант сделать быстрей?

Добавлено через 8 минут
И еще вдогонку вопрос.
Тоже функция ест, тоже рекурсивно выполняется, но расчет там за 5 секунд идет:
C#
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58

        public double PBPgt(double p1=0,double p2=0)
        {
            double res=0;
            p11 = CalcPoint1();
            p12 = 1-p11;
            p21 = CalcPoint2();
            p22=1-p21;
            if (p1==p2&p1>=6) res = (p11*p21)/(p11*p21+p12*p22);  
            else if (p1==p2+1&p1>=6) res = p21*1+p22*PBPgt(6,6);
            else if (p2==p1+1&p2>=6) res = p21*PBPgt(6,6)+p22*0;
            else if (p1==5&p2==5) res = p21*PAPgt(6,5)+p22*PAPgt(5,6);
            else if (p1==6&p2==4) res = p21*1+p22*PAPgt(6,5);
            else if (p1==6&p2==3) res = p21*1+p22*PBPgt(6,4);
            else if (p1==6&p2==2) res = p21*1+p22*PAPgt(6,3);
            else if (p1==6&p2==1) res = p21*1+p22*PBPgt(6,2);
            else if (p1==6&p2==0) res = p21*1+p22*PAPgt(6,1);
            else if (p1==5&p2==4) res = p21*PBPgt(6,4)+p22*PBPgt(5,5);
            else if (p1==5&p2==3) res = p21*PAPgt(6,3)+p22*PAPgt(5,4);
            else if (p1==5&p2==2) res = p21*PBPgt(6,2)+p22*PBPgt(5,3);
            else if (p1==5&p2==1) res = p21*PAPgt(6,1)+p22*PAPgt(5,2);
            else if (p1==5&p2==0) res = p21*PBPgt(6,0)+p22*PBPgt(5,1);
            else if (p1==4&p2==6) res = p21*PAPgt(5,6)+p22*0;
            else if (p1==4&p2==5) res = p21*PBPgt(5,5)+p22*PBPgt(4,6);
            else if (p1==4&p2==4) res = p21*PAPgt(5,4)+p22*PAPgt(4,5);
            else if (p1==4&p2==3) res = p21*PBPgt(5,3)+p22*PBPgt(4,4);
            else if (p1==4&p2==2) res = p21*PAPgt(5,2)+p22*PAPgt(4,3);
            else if (p1==4&p2==1) res = p21*PBPgt(5,1)+p22*PBPgt(4,2);
            else if (p1==4&p2==0) res = p21*PAPgt(5,0)+p22*PAPgt(4,1);
            else if (p1==3&p2==6) res = p21*PBPgt(4,6)+p22*0;
            else if (p1==3&p2==5) res = p21*PAPgt(4,5)+p22*PAPgt(3,6);
            else if (p1==3&p2==4) res = p21*PBPgt(4,4)+p22*PBPgt(3,5);
            else if (p1==3&p2==3) res = p21*PAPgt(4,3)+p22*PAPgt(3,4);
            else if (p1==3&p2==2) res = p21*PBPgt(4,2)+p22*PBPgt(3,3);
            else if (p1==3&p2==1) res = p21*PAPgt(4,1)+p22*PAPgt(3,2);
            else if (p1==3&p2==0) res = p21*PBPgt(4,0)+p22*PBPgt(3,1);
            else if (p1==2&p2==6) res = p21*PAPgt(3,6)+p22*0;
            else if (p1==2&p2==5) res = p21*PBPgt(3,5)+p22*PBPgt(2,6);
            else if (p1==2&p2==4) res = p21*PAPgt(3,4)+p22*PAPgt(2,5);
            else if (p1==2&p2==3) res = p21*PBPgt(3,3)+p22*PBPgt(2,4);
            else if (p1==2&p2==2) res = p21*PAPgt(3,2)+p22*PAPgt(2,3);
            else if (p1==2&p2==1) res = p21*PBPgt(3,1)+p22*PBPgt(2,2);
            else if (p1==2&p2==0) res = p21*PAPgt(3,0)+p22*PAPgt(2,1); 
            else if (p1==1&p2==6) res = p21*PBPgt(2,6)+p22*0;
            else if (p1==1&p2==5) res = p21*PAPgt(2,5)+p22*PAPgt(1,6);
            else if (p1==1&p2==4) res = p21*PBPgt(2,4)+p22*PBPgt(1,5);
            else if (p1==1&p2==3) res = p21*PAPgt(2,3)+p22*PAPgt(1,4);
            else if (p1==1&p2==2) res = p21*PBPgt(2,2)+p22*PBPgt(1,3);
            else if (p1==1&p2==1) res = p21*PAPgt(2,1)+p22*PAPgt(1,2);
            else if (p1==1&p2==0) res = p21*PBPgt(2,0)+p22*PBPgt(1,1);  
            else if (p1==0&p2==6) res = p21*PAPgt(1,6)+p22*0;
            else if (p1==0&p2==5) res = p21*PBPgt(1,5)+p22*PBPgt(0,6);
            else if (p1==0&p2==4) res = p21*PAPgt(1,4)+p22*PAPgt(0,5);
            else if (p1==0&p2==3) res = p21*PBPgt(1,3)+p22*PBPgt(0,4);
            else if (p1==0&p2==2) res = p21*PAPgt(1,2)+p22*PAPgt(0,3);
            else if (p1==0&p2==1) res = p21*PBPgt(1,1)+p22*PBPgt(0,2);
            else  res = p21*PAPgt(1,0)+p22*PAPgt(0,1);          
            return res;             
        }
А предыдущая функция аналогична по типу, но выполняется раз в 10 медленнее.
0
 Аватар для Kill100
434 / 299 / 82
Регистрация: 11.12.2010
Сообщений: 1,209
06.09.2016, 12:47
зачем вам в
C#
1

g1<=4&g2<=4
одинарное & ??? Ведь в этом случае даже если первое условие не верное будет проверено и второе. оно вам надо?
+ даётся мне что этот ужс можно записать циклом
2) Зачем писать так
C#
1

  if (g1==5&g2==2) res = t*1+(1-t)*PBPsT(5,3);
если можно так
C#
1

   if (g1==5&g2==2) return t*1+(1-t)*PBPsT(5,3);
тогда всё ниже не будет проверяться
0
 Аватар для Cupko
658 / 595 / 171
Регистрация: 17.07.2012
Сообщений: 1,682
Записей в блоге: 1
06.09.2016, 13:03
Kill100, ага, зачем нужно
C#
1
2
3
4
5
6

else if (p1==5&p2==5) res = p21*PAPgt(6,5)+p22*PAPgt(5,6);
else if (p1==5&p2==4) res = p21*PBPgt(6,4)+p22*PBPgt(5,5);
else if (p1==5&p2==3) res = p21*PAPgt(6,3)+p22*PAPgt(5,4);
else if (p1==5&p2==2) res = p21*PBPgt(6,2)+p22*PBPgt(5,3);
else if (p1==5&p2==1) res = p21*PAPgt(6,1)+p22*PAPgt(5,2);
else if (p1==5&p2==0) res = p21*PBPgt(6,0)+p22*PBPgt(5,1);
если можно сделать так:
C#
1
2
3
4

else if (p1 == 5)
{
    return p21*PAPgt(6,p2)+p22*PAPgt(5,p2 + 1);
}
С кодом тут беда, да.
0
1 / 1 / 0
Регистрация: 22.02.2016
Сообщений: 9
06.09.2016, 13:14  [ТС]
Я знаю, можно короче. Но скорости к расчету это не прибавит.
0
 Аватар для Kill100
434 / 299 / 82
Регистрация: 11.12.2010
Сообщений: 1,209
08.09.2016, 07:38
Цитата Сообщение от RealFlush Посмотреть сообщение
Я знаю, можно короче. Но скорости к расчету это не прибавит.
Прибавит Вы тут зря так думаете.
0
1 / 1 / 0
Регистрация: 22.02.2016
Сообщений: 9
19.09.2016, 16:14  [ТС]
Переделал я функции.
Всё равно долго
C#
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39

        public double PAPsT(double g1=0, double g2=0)
        {
            double p1;
            double p2;
            p1 = 0;//Convert.ToDouble(numericUpDown3.Value);
            p2 = 0;//Convert.ToDouble(numericUpDown4.Value);                                  
            double t = PA(p1,p2);
            if (g1<=4&g2<=4) return t*PBPsT(g1+1,g2)+(1-t)*PBPsT(g1,g2+1);
            else if (g1==6&g2==6) return PAPgt(p1,p2);
            else if (g1==6&g2==5) return t*1+(1-t)*PBPsT(6,6);
            else if (g1==5&g2==6) return t*PBPsT(6,6)+(1-t)*0;
            else if (g1==5&g2==5) return t*PBPsT(6,5)+(1-t)*PBPsT(5,6);
            else if (g1==5&g2<=4) return t*1+(1-t)*PBPsT(5,g2+1);       
            else  return t*PBPsT(g1+1,5)+(1-t)*0;       
        }
        public double PBPsT(double g1=0, double g2=0)  
        {
            double p1;
            double p2;  
            p1 = 0;//Convert.ToDouble(numericUpDown3.Value);
            p2 = 0;//Convert.ToDouble(numericUpDown4.Value);                                  
            double t = PB(p1,p2);
            if (g1<=4&g2<=4) return t*PAPsT(g1+1,g2)+(1-t)*PAPsT(g1,g2+1);
            else if (g1==6&g2==6) return PAPgt(p1,p2);
            else if (g1==6&g2==5) return t*1+(1-t)*PAPsT(6,6);
            else if (g1==5&g2==6) return t*PAPsT(6,6)+(1-t)*0;
            else if (g1==5&g2==5) return t*PAPsT(6,5)+(1-t)*PAPsT(5,6);
            else if (g1==5&g2==4) return t*1+(1-t)*PAPsT(5,5);
            else if (g1==5&g2==3) return t*1+(1-t)*PAPsT(5,4);
            else if (g1==5&g2==2) return t*1+(1-t)*PAPsT(5,3);
            else if (g1==5&g2==1) return t*1+(1-t)*PAPsT(5,2);
            else if (g1==5&g2==0) return t*1+(1-t)*PAPsT(5,1);
            
            else if (g1==4&g2==5) return t*PAPsT(5,5)+(1-t)*0;          
            else if (g1==3&g2==5) return t*PAPsT(4,5)+(1-t)*0;
            else if (g1==2&g2==5) return t*PAPsT(3,5)+(1-t)*0;
            else if (g1==1&g2==5) return t*PAPsT(2,5)+(1-t)*0;
            else  return t*PAPsT(1,5)+(1-t)*0;          
        }
При выполнении:
C#
1
2
3
4
5

System.Diagnostics.Stopwatch sw = new System.Diagnostics.Stopwatch();
sw.Start();
MessageBox.Show(Convert.ToString(PAPsT(4,4)));
sw.Stop();
MessageBox.Show((sw.ElapsedMilliseconds/6000.0).ToString());
10 - секунд

C#
1
2
3
4
5

System.Diagnostics.Stopwatch sw = new System.Diagnostics.Stopwatch();
sw.Start();
MessageBox.Show(Convert.ToString(PAPsT(0,0)));
sw.Stop();
MessageBox.Show((sw.ElapsedMilliseconds/6000.0).ToString());
- 1000 секунд.
0
 Аватар для Cupko
658 / 595 / 171
Регистрация: 17.07.2012
Сообщений: 1,682
Записей в блоге: 1
19.09.2016, 19:14
RealFlush, давайте весь код, с PA() PB() и PBgt().
Как бы, что вы там делаете с рекурсивными вычислениями мне не понять, но попрофилировать может удасться.
Есть еще такая штука как мемоизация - может помочь в вашем случае. Хотя я грешу на неоптимальнось вычислений или же кривой код.
0
1 / 1 / 0
Регистрация: 22.02.2016
Сообщений: 9
19.09.2016, 22:10  [ТС]
C#
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218

        public double CalcPoint1()
        {
            double g11;
            double g12;
            double g21;
            double g22;
            double pnt;
            double res=0;
            g11 = 1/Convert.ToDouble(textBox20.Text);
            g12 = 1/Convert.ToDouble(textBox21.Text);
            g21 = 1/Convert.ToDouble(textBox96.Text);
            g22 = 1/Convert.ToDouble(textBox95.Text);
            double margin1g=1-1/(g11+g12);
            double p11 = (1-margin1g)/Convert.ToDouble(textBox20.Text);
            //MessageBox.Show(Convert.ToString(p11));
            for (double i=0;i<=1;i+=0.001)
            {
                pnt = Math.Pow(i,4)*(15-34*i+28*Math.Pow(i,2)-8*Math.Pow(i,3))/(1-2*i+2*i*i)-p11;
            //  MessageBox.Show(Convert.ToString(pnt));
                if (pnt<=0.001& pnt>=-0.001) res=i;
                //else res= 0;
            }
       // MessageBox.Show(Convert.ToString(res));
        return res;
        
        }
        public double CalcPoint2()
        {
            double g11;
            double g12;
            double g21;
            double g22;
            double pnt;
            double res=0;
            g11 = 1/Convert.ToDouble(textBox20.Text);
            g12 = 1/Convert.ToDouble(textBox21.Text);
            g21 = 1/Convert.ToDouble(textBox96.Text);
            g22 = 1/Convert.ToDouble(textBox95.Text);
            double margin1g=1-1/(g11+g12);
            double p21 = (1-margin1g)/Convert.ToDouble(textBox96.Text);
            //MessageBox.Show(Convert.ToString(p11));
            for (double i=0;i<=1;i+=0.001)
            {
                pnt = Math.Pow(i,4)*(15-34*i+28*Math.Pow(i,2)-8*Math.Pow(i,3))/(1-2*i+2*i*i)-p21;
            //  MessageBox.Show(Convert.ToString(pnt));
                if (pnt<=0.001& pnt>=-0.001) res=i;
                //else res= 0;
            }
       // MessageBox.Show(Convert.ToString(res));
        return res;
        
        } 
        public double PA(double p1=0,double p2=0)
        {   
//          double res=0;
            p11 = CalcPoint1();
            p12 = 1-p11;
//          MessageBox.Show(Convert.ToString(p11));
            
//          if (p1==0&p2==0)
//          {
//              return p11*PA(1,0)+p12*PA(0,1);
//          }  
            if (p1==p2&p1>=3)
            {
                return Math.Pow(p11,2)/(Math.Pow(p11,2)+Math.Pow(p12,2));
            }
            else if (p1==p2+1&p1>=3)
            {
                return p11*1+p12*PA(3,3);
            }
            else if (p2==p1+1&p2>=3)
            {
                return p11*PA(3,3)+p12*0;
            }    
            else if (p1==3&p2==1)
            {
                return p11*1+p12*PA(3,2);
            }
            else if (p1==1&p2==3)
            {
                return p11*PA(2,3)+p12*0;
            }
            else if (p1==3&p2==0)
            {
                return p11*1+p12*PA(3,1);
            }
            else if (p1==0&p2==3)
            {
                return p11*PA(1,3)+p12*0;
            }
            else
            {
                return p11*PA(p1+1,p2)+p12*PA(p1,p2+1);
            }
            
        }
        public double PB(double p1=0,double p2=0)
        {   
//          double res=0;
            p21 = CalcPoint2();
            p22 = 1-p21;
 
            if (p1==p2&p1>=3)
            {
                return Math.Pow(p21,2)/(Math.Pow(p21,2)+Math.Pow(p22,2));
            }
            else if (p1==p2+1&p1>=3)
            {
                return p21*1+p22*PB(3,3);
            }
            else if (p2==p1+1&p2>=3)
            {
                return p21*PB(3,3)+p22*0;
            }    
            else if (p1==3&p2==1)
            {
                return p21*1+p22*PB(3,2);
            }
            else if (p1==1&p2==3)
            {
                return p21*PB(2,3)+p22*0;
            }
            else if (p1==3&p2==0)
            {
                return p21*1+p22*PB(3,1);
            }
            else if (p1==0&p2==3)
            {
                return p21*PB(1,3)+p22*0;
            }
            else
            {
                return p21*PB(p1+1,p2)+p22*PB(p1,p2+1);
            }           
        }   
        public double PAPgt(double p1=0,double p2=0)   
        {
//          double res=0;
            p11 = CalcPoint1();
            p12 = 1-p11;
            p21 = CalcPoint2();
            p22=1-p21;
            if (p1==0&p2==0) return p11*PBPgt(1,0)+p12*PBPgt(0,1);
            else if (p1==p2&p1>=6) return (p11*p21)/(p11*p21+p12*p22);
            else if (p1==p2+1&p1>=6) return p11*1+p12*PAPgt(6,6);
            else if (p2==p1+1&p2>=6) return p11*PAPgt(6,6)+p12*0;
            else if (p1==5&p2==5) return p11*PBPgt(6,5)+p12*PBPgt(5,6);
            else if (p1==6&p2<=4&p2%2==0) return p11*1+p12*PBPgt(6,p2+1);
            else if (p1==6&p2<=4&p2%2>0) return p11*1+p12*PAPgt(6,p2+1);
//          else if (p1==6&p2==3) return p11*1+p12*PAPgt(6,4);
//          else if (p1==6&p2==2) return p11*1+p12*PBPgt(6,3);
//          else if (p1==6&p2==1) return p11*1+p12*PAPgt(6,2);
//          else if (p1==6&p2==0) return p11*1+p12*PBPgt(6,1);
            else if (p1<=5&p2<=5&p2%2==0) return p11*PAPgt(p1+1,p2)+p12*PAPgt(p1,p2+1);
            else if (p1<=5&p2<=5&p2%2>0) return p11*PBPgt(p1+1,p2)+p12*PBPgt(p1,p2+1);          
            else if (p1<=4&p2==6&p2%2==0) return p11*PBPgt(p1+1,p2)+p12*0;
            else return p11*PAPgt(p1+1,p2)+p12*0;
 
        }
        public double PBPgt(double p1=0,double p2=0)
        {
//          double res=0;
            p11 = CalcPoint1();
            p12 = 1-p11;
            p21 = CalcPoint2();
            p22=1-p21;
            if (p1==0&p2==0) return p21*PAPgt(1,0)+p22*PAPgt(0,1);
            else if (p1==p2&p1>=6) return  (p11*p21)/(p11*p21+p12*p22);  
            else if (p1==p2+1&p1>=6) return  p21*1+p22*PBPgt(6,6);
            else if (p2==p1+1&p2>=6) return p21*PBPgt(6,6)+p22*0;
            else if (p1==6&p2<=4&p2%2==0) return p21*1+p22*PAPgt(p1,p1+1);
            else if (p1==6&p2<=4&p2%2>0) return p21*1+p22*PBPgt(p1,p1+1);           
            else if (p11<=5&p2<=5&p2%2==0) return p21*PBPgt(p1+1,p2)+p22*PBPgt(p1,p2+1);
            else if (p11<=5&p2<=5&p2%2>0) return p21*PAPgt(p1+1,p2)+p22*PAPgt(p1,p2+1);
            else if (p1<=4&p2==6&p1%2==0) return p21*PAPgt(p1+1,p2)+p22*0;
            else return p21*PBPgt(p1+1,p2)+p22*0;
            
        }
        public double PAPsT(double g1=0, double g2=0)
        {
            double p1;
            double p2;
            p1 = 0;//Convert.ToDouble(numericUpDown3.Value);
            p2 = 0;//Convert.ToDouble(numericUpDown4.Value);                                  
            double t = PA(p1,p2);
            if (g1<=4&g2<=4) return t*PBPsT(g1+1,g2)+(1-t)*PBPsT(g1,g2+1);
            else if (g1==6&g2==6) return PAPgt(p1,p2);
            else if (g1==6&g2==5) return t*1+(1-t)*PBPsT(6,6);
            else if (g1==5&g2==6) return t*PBPsT(6,6)+(1-t)*0;
            else if (g1==5&g2==5) return t*PBPsT(6,5)+(1-t)*PBPsT(5,6);
            else if (g1==5&g2<=4) return t*1+(1-t)*PBPsT(5,g2+1);       
            else  return t*PBPsT(g1+1,5)+(1-t)*0;       
        }
        public double PBPsT(double g1=0, double g2=0)  
        {
            double p1;
            double p2;  
            p1 = 0;//Convert.ToDouble(numericUpDown3.Value);
            p2 = 0;//Convert.ToDouble(numericUpDown4.Value);                                  
            double t = PB(p1,p2);
            if (g1<=4&g2<=4) return t*PAPsT(g1+1,g2)+(1-t)*PAPsT(g1,g2+1);
            else if (g1==6&g2==6) return PAPgt(p1,p2);
            else if (g1==6&g2==5) return t*1+(1-t)*PAPsT(6,6);
            else if (g1==5&g2==6) return t*PAPsT(6,6)+(1-t)*0;
            else if (g1==5&g2==5) return t*PAPsT(6,5)+(1-t)*PAPsT(5,6);
            else if (g1==5&g2==4) return t*1+(1-t)*PAPsT(5,5);
            else if (g1==5&g2==3) return t*1+(1-t)*PAPsT(5,4);
            else if (g1==5&g2==2) return t*1+(1-t)*PAPsT(5,3);
            else if (g1==5&g2==1) return t*1+(1-t)*PAPsT(5,2);
            else if (g1==5&g2==0) return t*1+(1-t)*PAPsT(5,1);
            
            else if (g1==4&g2==5) return t*PAPsT(5,5)+(1-t)*0;          
            else if (g1==3&g2==5) return t*PAPsT(4,5)+(1-t)*0;
            else if (g1==2&g2==5) return t*PAPsT(3,5)+(1-t)*0;
            else if (g1==1&g2==5) return t*PAPsT(2,5)+(1-t)*0;
            else  return t*PAPsT(1,5)+(1-t)*0;          
        }
0
Эксперт .NET
 Аватар для Wolfdp
3790 / 1767 / 371
Регистрация: 15.06.2012
Сообщений: 6,543
Записей в блоге: 3
20.09.2016, 02:21
ить, осилил только один метод...
C#
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24

public static double PBPsT(int g1, int g2)
        {
            double p1;
            double p2;
            p1 = 0;//Convert.ToDouble(numericUpDown3.Value);
            p2 = 0;//Convert.ToDouble(numericUpDown4.Value);
            var t = PA(p1, p2);
            var u = (1 - t);
            if (g1 <= 4 & g2 <= 4) return t * (t + u * t * (t + 2 * u * PAPgt(p1, p2))) + u * Math.Pow(t, 2) * (t + 2 * u * PAPgt(p1, p2));
            if (g1 == 6 & g2 == 6) return PAPgt(p1, p2);
            if (g1 == 6 & g2 == 5) return t + u * PAPgt(p1, p2);
            if (g1 == 5 & g2 == 6) return t * PAPgt(p1, p2);
            if (g1 == 5 & g2 == 5) return t * (t + 2 * u * PAPgt(p1, p2));
            if (g1 == 5 & g2 == 4) return t + u * t * (t + 2 * u * PAPgt(p1, p2));
            if (g1 == 5 & g2 == 3) return t + u * (t + u * t * (t + 2 * u * PAPgt(p1, p2)));
            if (g1 == 5 & g2 == 2) return t + u * (t + u * (t + u * t * (t + 2 * u * PAPgt(p1, p2))));
            if (g1 == 5 & g2 == 1) return t + u * (t + u * (t + u * (t + u * t * (t + 2 * u * PAPgt(p1, p2)))));
            if (g1 == 5 & g2 == 0) return t + u * (t + u * (t + u * (t + u * (t + u * t * (t + 2 * u * PAPgt(p1, p2))))));
 
            if (g2 == 5 && (g1 >= 0 || g1 <= 4))
                return Math.Pow(t, 6 - g1) * (t + 2 * u * PAPgt(p1, p2));
 
            return 0;
        }
Добавлено через 2 минуты
в догонку: вообще вьіше написаное крайне смахивает на писанину китайца, но заморачиватся ради четьірех строчек не хочется.
1
1168 / 886 / 517
Регистрация: 09.04.2014
Сообщений: 2,101
20.09.2016, 10:46
странно вы время меряете - вы включаете всемя показа 1 окна в длительность расчета, попробуйте так
C#
1
2
3
4
5
6

System.Diagnostics.Stopwatch sw = new System.Diagnostics.Stopwatch();
sw.Start();
string result=Convert.ToString(PAPsT(4,4));
sw.Stop();
MessageBox.Show(result);
MessageBox.Show((sw.ElapsedMilliseconds/6000.0).ToString());
Добавлено через 23 минуты
Интересно, что вы ститаете... Для параметров 4, 4 количество вызовов методов поражает:
PAPsT: 7
PA: 483
CalcPoint1: 13911
PBPsT: 6
PB: 414
CalcPoint2: 13842
PAPgt: 6624
PBPgt: 6804
Добавлено через 10 минут
В CalcPoint1 можно заметно ускорить если не использовать Math.Pow в цикле и миримизировать число операций умножения (у меня время выполнения с параметрами 4,4 уменьшилось с 4700 мс до 350 мс):
C#
1
2
3
4
5
6
7
8

...
 for (double i = 0; i <= 1; i += 0.001)
            {
                double i2 = i * i;
                pnt = i2*i2 * (15 - 34 * i + 28 * i2 - 8 * i2*i) / (1 - 2 * i + 2 * i2) - p11;
                if (pnt <= 0.001 && pnt >= -0.001) res = i;
            }
...
анлогично в CalcPoint2
C#
1
2
3
4
5
6

 for (double i = 0; i <= 1; i += 0.001)
            {
                double i2 = i * i;
                pnt = i2*i2 * (15 - 34 * i + 28 * i2 - 8 * i2*i) / (1 - 2 * i + 2 * i2) - p21;
                if (pnt <= 0.001 && pnt >= -0.001) res = i;
            }
1
1 / 1 / 0
Регистрация: 22.02.2016
Сообщений: 9
21.09.2016, 13:12  [ТС]
Wolfdp, Сесе.
nedel, Спасибо большое, помогли.
0
1168 / 886 / 517
Регистрация: 09.04.2014
Сообщений: 2,101
21.09.2016, 15:00
и еще, если в CalcPoint1 и CalcPoint2 в цикле вы ищете первый подходящий результат, то стоит прерывать цикл после того как вы его нашлиб чтобы не считать лишнего
C#
1
2
3
4
5

 if (pnt <= 0.001 && pnt >= -0.001)
{
     res = i;
     break;
}
0
Надоела реклама? Зарегистрируйтесь и она исчезнет полностью.
raxper
Эксперт
30234 / 6612 / 1498
Регистрация: 28.12.2010
Сообщений: 21,154
Блог
 21.09.2016, 15:00
Помогаю со студенческими работами здесь

Заливка затравкой. Долгое выполнение
Здравствуйте, алгоритм вроде как работает, но делает это овер медлено, возможно проблема со списками, ну или с тем что к канве обращаюсь,...

ОЧЕНЬ долгое выполнение запроса
Очень долгое выполнение запроса, в каждой таблице (кроме regions), примерно 60 000 строк. Такой запрос, который выбирает все товары у...

Долгое выполнение вставки в цикле
Все привет, подскажите что не так сделал, уж слишком долго он вставляет данные в таблицу DECLARE @CountGroup int = 500 DECLARE...

Долгое выполнение ajax запросов
Всем привет. У меня проблема с долгим выполнением ajax запросов. Они обрабатываются по 2-3 секунды каждый запрос. = ( Сперва...

Долгое выполнение вывода новостей
Добрый вечер! Есть 2 скрипта: popular.php &lt;?php $sql_result = $db-&gt;query( &quot;SELECT id, title, date, alt_name, category,...


Искать еще темы с ответами

Или воспользуйтесь поиском по форуму:
14
Ответ Создать тему
Новые блоги и статьи
1С: Контроль уникальности заводского номера
Maks 23.03.2026
Алгоритм контроля уникальности заводского (или серийного) номера на примере документа выдачи шин для спецтехники с табличной частью. Данные берутся из регистра сведений, по которому настроено. . .
Хочу заставить корпорации вкладываться в здоровье сотрудников: делаю мат модель здравосохранения
anaschu 22.03.2026
e7EYtONaj8Y Z4Tv2zpXVVo https:/ / github. com/ shumilovas/ med2. git
 1С: Программный отбор элементов справочника по группе
Maks 22.03.2026
Установка программного отбора элементов справочника "Номенклатура" из модуля формы документа. В качестве фильтра для отбора справочника служит группа номенклатуры. Отбор по наименованию группы. . .
Как я обхитрил таблицу Word
Alexander-7 21.03.2026
Когда мигает курсор у внешнего края таблицы, и нам надо перейти на новую строку, а при нажатии Enter создается новый ряд таблицы с ячейками, то мы вместо нервных нажатий Энтеров мы пишем любые буквы. . .
Krabik - рыболовный бот для WoW 3.3.5a
AmbA 21.03.2026
без регистрации и смс. Это не торговля, приложение не содержит рекламы. Выполняет свою непосредственную задачу - автоматизацию рыбалки в WoW - и ничего более. Однако если админы будут против -. . .
1С: Программный отбор элементов справочника по значению перечисления
Maks 21.03.2026
Установка программного отбора элементов справочника "Сотрудники" из модуля формы документа. В качестве фильтра для отбора служит значение перечислений. / / Событие "НачалоВыбора" реквизита на форме. . .
Переходник USB-CAN-GPIO
Eddy_Em 20.03.2026
Достаточно давно на работе возникла необходимость в переходнике CAN-USB с гальваноразвязкой, оный и был разработан. Однако, все меня терзала совесть, что аж 48-ногий МК используется так тупо: просто. . .
Оттенки серого
Argus19 18.03.2026
Оттенки серого Нашёл в интернете 3 прекрасных модуля: Модуль класса открытия диалога открытия/ сохранения файла на Win32 API; Модуль класса быстрого перекодирования цветного изображения в оттенки. . .
КиберФорум - форум программистов, компьютерный форум, программирование
Поддержать - Реклама - Условия использования - Обратная связь
Powered by vBulletin
Copyright ©2000 - 2026, CyberForum.ru