Метод Монте-Карло (вероятность попадания при стрельбе)

#include <iostream>
#include <vector>
#include <algorithm>
#include <numeric>
 
int main()
{
using namespace std;
int num_of_tanks;
cout « "Enter number of tanks:\n";
cin » num_of_tanks;
vector<double> tanks(num_of_tanks);
cout « "enter posibilities for each tank:\n";
for (int i = 0; i < tanks.size(); i++)
{
cout « "posibility for " « i + 1 « ") tank: ";
cin » tanks[i];
}
if (accumulate(tanks.begin(), tanks.end(), 0) > 1)
cout « "To correct ansver sum of pisibilities must be no more, than 1\n";
 
//расчет вероятности двойного попадания
double double_hit = 0;
for (int i = 0; i < tanks.size(); i++)
double_hit += (tanks[i] * tanks[i]);
 
//вероятность попадания в соседние баки
double neibour_hir = 0;
for (int i = 0; i < tanks.size() - 1; i++)
neibour_hir += tanks[i] * tanks[i + 1];
for (int i = tanks.size()-1; i > 0; i--)
neibour_hir += tanks[i] * tanks[i - 1];
 
cout « "inflammation posibility = "
« neibour_hir + double_hit « "%" « endl;
 
cin.get();
cin.get();
return 0;
}

double_hit += (tanks[i] * tanks[i]);

#include <vector>
#include <ctime>
#include <cstdlib>  
 
int main(){   
    int attempts, num_of_tanks;     
    std::cout<<"Input number of iterations (please more): ";
    std::cin>>attempts; 
    std::cout<<"Input number of tanks: ";
    std::cin>>num_of_tanks;
    std::vector<double> probs(num_of_tanks);
    std::vector<int>    Fbounds(num_of_tanks);
    std::cout << "enter hit probability for each tank:\n";
    for (int i = 0; i < probs.size(); ++i)
    {
        std::cout << "hit probability for (" << i + 1 << ") tank: ";
        std::cin >> probs[i];
        if (i)
            Fbounds[i]=Fbounds[i-1]+ probs[i]*RAND_MAX;
        else
            Fbounds[i]=probs[i]*RAND_MAX;
    }
    Fbounds[num_of_tanks-1]=RAND_MAX;
    srand(time(NULL));
    int triggered_counter=0;
    for (int i=0; i<attempts; ++i)
    {
        int random_key, random_hit1, random_hit2;
 
        random_key=rand();
        for(random_hit1=0; random_hit1<num_of_tanks && random_key>Fbounds[random_hit1]; ++random_hit1);  
        random_key=rand();
        for(random_hit2=0; random_hit2<num_of_tanks && random_key>Fbounds[random_hit2]; ++random_hit2);  
        if (abs(random_hit1-random_hit2)<=1)//если баки соседние
            ++triggered_counter;
    }
    double explosion= triggered_counter*1.0/attempts;
    std::cout<<"explosion probability "<< explosion*100<<"%"<<std::endl;
 
    std::cin.get();
    return 0;
}