Реализация программы на С++ (Найти вероятности следующих событий:)

#include <cstdlib>
#include <ctime>
#include <cmath>
#include <iostream>
#include <iomanip>
#include <vector>
#include <algorithm>
#include <iterator>
 
double negateProbability(double probability) {
  return 1 - probability;
}
 
class Shooter {
 public:
  Shooter(int ammo, double hitProbability) : ammo(ammo),
    hitProbability(hitProbability) {}
  int getAmmo() const {
    return ammo;
  }
  // шанс попасть
  double getHitProbability() const {
    return hitProbability;
  }
  // шанс промахнуться
  double getMissProbability() const {
    return negateProbability(getHitProbability());
  }
  // шанс, что как минимум все, кроме последнего
  // выстрелы промахнутся
  double getEmptyAmmoProbability() const {
    return pow(getMissProbability(), getAmmo() - 1);
  }
  // шанс, что останется хотя бы один патрон
  double getNonEmptyAmmoProbability() const {
    return negateProbability(getEmptyAmmoProbability());
  }
 private:
  int ammo;
  double hitProbability;
};
 
struct RandomShooter {
  int ammo;
  double minHitProbability, maxHitProbability;
  RandomShooter(int ammo, double minHitProbability, double maxHitProbability)
    : ammo(ammo), minHitProbability(minHitProbability),
    maxHitProbability(maxHitProbability) {}
  Shooter operator()() const {
    return Shooter(ammo, minHitProbability +
      (rand() % (int)((maxHitProbability - minHitProbability) * 1000) / 1000.));
  }
};
 
std::ostream &operator<<(std::ostream &stream, const Shooter &s) {
  return stream << "Shooter{" <<
    "ammo=" << s.getAmmo() << 
    std::fixed << std::setprecision(5) <<
    ",hitProbability=" << s.getHitProbability() << "}";
}
 
double multiplyByEmptyAmmoPorbability(double x, const Shooter &s) {
  return x * s.getEmptyAmmoProbability();
}
 
double multiplyByNonEmptyAmmoPorbability(double x, const Shooter &s) {
  return x * s.getNonEmptyAmmoProbability();
}
 
double accumulateMissChance(double x, const Shooter &s, double Pb) {
  return x + (s.getHitProbability() / s.getMissProbability()) * Pb; 
}
 
int main(int argc, char *argv[]) {
  srand(time(0));
  
  int ammo = 4;
  int numberOfShooters = 10;
  std::vector<Shooter> shooters;
  
  std::generate_n(std::back_inserter(shooters), numberOfShooters,
    RandomShooter(ammo, 0.2, 0.3));
   
  std::copy(shooters.begin(), shooters.end(),
    std::ostream_iterator<Shooter>(std::cout, "\n"));
 
  double Pa = 1.0;
  for (size_t i = 0; i < shooters.size(); ++i) {
    Pa = multiplyByEmptyAmmoPorbability(Pa, shooters[i]);
  }
  Pa = negateProbability(Pa);
  std::cout << "Probability that at least one shooter would have\n"
    "at least one shot left: " << Pa << std::endl;
  
  double Pb = 1.0;
  for (size_t i = 0; i < shooters.size(); ++i) {
    Pb = multiplyByNonEmptyAmmoPorbability(Pb, shooters[i]);
  }
  std::cout << "Probability that all shooters would have\n"
    "at least one shot left: " << Pb << std::endl;
    
  double Pc = 0.0;
  for (size_t i = 0; i < shooters.size(); ++i) {
    Pc = accumulateMissChance(Pc, shooters[i], Pb);
  }
  std::cout << "Probability that one and only one shooter would not have \n"
    "at least one shot left: " << Pc << std::endl;
}