Форум программистов, компьютерный форум CyberForum.ru

метод Монте-Карло - C++

Войти
Регистрация
Восстановить пароль
Другие темы раздела
C++ Окна обозреветель решений ? http://www.cyberforum.ru/cpp-beginners/thread144232.html
Скажите как поменять чтобы окно обозреватель решений было с лева а не с права ? фаил приложил !
C++ поправить чуток Добрый день! кто-нибудь может подправить графику чтоб не все в перемешку двигалось, а в виде снежинок.... #include<graphics.h> #include<conio.h> #include<stdlib.h> int random(int N){return rand() %N;} http://www.cyberforum.ru/cpp-beginners/thread144225.html
Поменять слова местами C++
Поменять местами слова с чет и неч номерами не изменяя порядка следования внутри группы. Аа ее сс tt на ее аа tt сс. На си. Помогите плиз на экзамене.
C++ матрица
подскажите пожалуйста как вывести такую матрицу? 1 4 16 64 4 0 0 16 16 0 0 4 64 16 4 1
C++ Гири http://www.cyberforum.ru/cpp-beginners/thread144193.html
На Си. Самородок весит н грамм. Есть гири по 100 50 3 1 гр. Взвесить, используя мин колич гирь.
C++ Строки Помогите щя экзамен незнаю как делать Вводяться строка символов и слово.Определять сколько раз встречаеться в строке символ слова,находящийся в центре слова. подробнее

Показать сообщение отдельно
Mr.X
Эксперт С++
3040 / 1685 / 265
Регистрация: 03.05.2010
Сообщений: 3,867
12.06.2010, 18:33     метод Монте-Карло
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
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
//Трехмерное тело образовано объединением нескольких сфер произвольного размера 
//и взаимного расположения. Найти объем этого тела, используя метод Монте-Карло.
#include <iostream>
#include <cmath>
#include <vector>
#include <algorithm>
#include <functional>
#include <limits>
#include <numeric>
#include <ctime>
#include <cstdlib>
 
enum  T_coord_ind  
{
    X = 0,
    Y,
    Z,
    TOTAL_COORDS
};
 
T_coord_ind& operator++(T_coord_ind&  coord_ind)
{
    return coord_ind = T_coord_ind(coord_ind + 1);
}
 
typedef std::vector<double> T_coords;
 
class T_3_point
{
    T_coords  coords_;
public:
    typedef T_coords::iterator        iterator;
    typedef T_coords::const_iterator  const_iterator;
    typedef T_coords::value_type      value_type;   
 
    T_3_point
        (            
            value_type  x = 0,
            value_type  y = 0,
            value_type  z = 0
        ) : coords_(TOTAL_COORDS)
    {
        coords_[X] = x;
        coords_[Y] = y;
        coords_[Z] = z;
    }    
 
    iterator begin()
    {
        return coords_.begin();
    }
 
    iterator end()
    {
        return coords_.end();
    }
    
    T_3_point operator-(T_3_point  p) 
    {
        T_3_point  res_p;
        std::transform
            (
                begin(), 
                end(), 
                p.begin(),                
                res_p.begin(),                                
                std::minus<value_type>()
            );
        return  res_p;     
    }
 
    value_type& operator[](T_coord_ind  coord_ind)
    {
        return coords_[coord_ind];
    }
 
    value_type operator[](T_coord_ind  coord_ind) const
    {
        return coords_[coord_ind];
    }
 
    T_3_point::value_type  distance(T_3_point  p) const
    { 
        T_3_point  vector_from_this_in_p = p - *this;
 
        return sqrt
            (
                std::inner_product
                    (
                        vector_from_this_in_p.begin(), 
                        vector_from_this_in_p.end(),
                        vector_from_this_in_p.begin(),                        
                        0.0
                    )
            );
    }
 
    void input()
    {
        std::cout << std::endl
                  << "x = ";
        std::cin >> coords_[X];
        std::cout << "y = ";
        std::cin >> coords_[Y];
        std::cout << "z = ";
        std::cin >> coords_[Z];
    }
};//class T_3_point
 
typedef std::vector<T_3_point>  T_3_points;
 
class T_sfera
{
    T_3_point              center_;
    T_3_point::value_type  radius_;
public:
    T_sfera
        (
            T_3_point              center = T_3_point(), 
            T_3_point::value_type  radius = 0
        ) : center_(center), radius_(radius) 
    {}
 
    bool  is_vnutri_sfery(T_3_point  p) const
    {
        return  center_.distance(p) <= radius_;
    }
 
    void input()
    {
        std::cout << std::endl
                  << "Введите данные сферы:"
                  << std::endl
                  << "координаты центра:";
        center_.input();
        std::cout << "радиус: ";
        std::cin >> radius_;
    }
 
    T_3_point::value_type  get_min_coord(T_coord_ind  coord_ind)
    {       
        return center_[coord_ind] - radius_;
    }
 
    T_3_point::value_type  get_max_coord(T_coord_ind  coord_ind)
    {       
        return center_[coord_ind] + radius_;
    }
};//class T_sfera
 
typedef std::vector<T_sfera>  T_sfery;
 
//////////////////////////////////////////////////////////////////////////
class T_get_min_coord
{
    T_coord_ind            coord_ind_;
    T_3_point::value_type  min_coord_;
public:
    T_get_min_coord
        (
            T_coord_ind  coord_ind
        ) 
        : coord_ind_(coord_ind), min_coord_(std::numeric_limits<T_3_point::value_type>::max())
    {}
 
    void operator()(T_sfera sfera)
    {        
        if(sfera.get_min_coord(coord_ind_) < min_coord_)
        {
            min_coord_ = sfera.get_min_coord(coord_ind_);
        }
    }
 
    operator T_3_point::value_type()
    {
        return min_coord_;
    }
};//class T_get_min_coord
//////////////////////////////////////////////////////////////////////////
class T_get_max_coord
{
    T_coord_ind            coord_ind_;
    T_3_point::value_type  max_coord_;
public:
    T_get_max_coord
        (
            T_coord_ind  coord_ind
        ) 
        : coord_ind_(coord_ind), max_coord_(std::numeric_limits<T_3_point::value_type>::min())
    {}
 
    void operator()(T_sfera sfera)
    {        
        if(sfera.get_max_coord(coord_ind_) > max_coord_)
        {
            max_coord_ = sfera.get_max_coord(coord_ind_);
        }
    }
 
    operator T_3_point::value_type()
    {
        return max_coord_;
    }
};
//////////////////////////////////////////////////////////////////////////
class T_get_random_point
{
    T_3_point  min_XYZ_point_;    
    T_3_point  max_XYZ_point_;
public:
    T_get_random_point
        (
            T_3_point  min_XYZ_point,
            T_3_point  max_XYZ_point
        ) : min_XYZ_point_(min_XYZ_point), max_XYZ_point_(max_XYZ_point) 
    {}
    
    double  rand_between_0_and_1()
    {
        return static_cast<double>(rand()) / static_cast<double>(RAND_MAX);
    }
 
    T_3_point operator() ()
    {
        T_3_point::value_type  x = min_XYZ_point_[X] 
                       + rand_between_0_and_1() * (max_XYZ_point_[X] - min_XYZ_point_[X]);
 
        T_3_point::value_type  y = min_XYZ_point_[Y] 
                       + rand_between_0_and_1() * (max_XYZ_point_[Y] - min_XYZ_point_[Y]);
 
        T_3_point::value_type  z = min_XYZ_point_[Z] 
                       + rand_between_0_and_1() * (max_XYZ_point_[Z] - min_XYZ_point_[Z]);
 
        return T_3_point(x, y, z);
    }
};
//////////////////////////////////////////////////////////////////////////
class  T_is_v_objedinenii_sfer
{
    T_sfery  sfery_;    
public:
    T_is_v_objedinenii_sfer(T_sfery  sfery) : sfery_(sfery)
    {}
 
    bool operator() (T_3_point  point)
    {           
        for(T_sfery::const_iterator  cur_sfera_it = sfery_.begin();
            cur_sfera_it != sfery_.end(); ++cur_sfera_it)
        {
            if(cur_sfera_it->is_vnutri_sfery(point)) return true;
        }
        return false;
    }
};
//////////////////////////////////////////////////////////////////////////           
T_3_point::value_type get_objem_objedin_sfer
    (
        T_sfery  sfery,
        int      n
    )
{    
    //Находим минимальный параллелепипед, содержащий все заданные сферы.
    T_3_point  min_XYZ_point;
 
    for(T_coord_ind  cur_coord_ind = X; cur_coord_ind < TOTAL_COORDS; ++cur_coord_ind)
    {
        min_XYZ_point[cur_coord_ind] 
            = std::for_each(sfery.begin(), sfery.end(), T_get_min_coord(cur_coord_ind));    
    }
 
    T_3_point  max_XYZ_point;
 
    for(T_coord_ind  cur_coord_ind = X; cur_coord_ind < TOTAL_COORDS; ++cur_coord_ind)
    {
        max_XYZ_point[cur_coord_ind] 
            = std::for_each(sfery.begin(), sfery.end(), T_get_max_coord(cur_coord_ind));    
    }
 
    //Вектор разностей максимальных и минимальных координат.
    T_3_point  dliny_storon_parallelepipeda = max_XYZ_point - min_XYZ_point;
    
    T_3_point::value_type  objem_parallelepipeda = std::accumulate
                                        (
                                            dliny_storon_parallelepipeda.begin(),
                                            dliny_storon_parallelepipeda.end(),
                                            1.0,                                            
                                            std::multiplies<T_3_point::value_type>()
                                        );    
 
    //Генерируем вектор случайных точек.
    T_3_points  random_points(n);
 
    std::generate(random_points.begin(), random_points.end(), 
                  T_get_random_point(min_XYZ_point, max_XYZ_point));
 
    T_3_point::value_type  points_v_objedin_sfer_count 
        = std::count_if(random_points.begin(), random_points.end(), 
                        T_is_v_objedinenii_sfer(sfery));
 
    return objem_parallelepipeda * points_v_objedin_sfer_count / n;
}
//////////////////////////////////////////////////////////////////////////           
 
int main()
{
    std::locale::global(std::locale(""));
    srand(static_cast<unsigned>(time(0)));
 
    for(;;)
    {
        int sfer_total;
 
        std::cout << std::endl
                  << std::endl
                  << "Введите число сфер: ";
        std::cin >> sfer_total;
        if(sfer_total < 1) break;
        
        T_sfery  sfery(sfer_total);
        std::for_each(sfery.begin(), sfery.end(), 
                      std::mem_fun_ref(&T_sfera::input));
 
        for(;;)
        {
            std::cout << std::endl
                      << "Введите число испытаний : ";
            int n;
            std::cin >> n;
            if(n < 1) break;
            std::cout << "Ждите..."
                      << std::endl;
            std::cout << "Объем объединения "
                      << sfer_total
                      << " заданных сфер равен: "
                      << get_objem_objedin_sfer(sfery, n)
                      << std::endl;  
        }
    }  
    return 0;
}
 
КиберФорум - форум программистов, компьютерный форум, программирование
Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2017, vBulletin Solutions, Inc.
Рейтинг@Mail.ru