#include <iostream>
#include <ctime>
#include <cstdlib>

using namespace std;

const long long N = 1000000000;

double montecarlo()
{
    srand(static_cast<unsigned>(time(nullptr)));

    long long licznik = 0;
    double x,y;

    for(long long i = 0; i < N; ++i)
    {
        x = rand() / 32767.0;       // losowanie liczby rzeczywistej od 0 do 1 (32767 to maksymalna losowana)
        y = rand() / 32767.0;
        if(x*x + y*y <= 1)
            ++licznik;
    }

    return 4.0 * licznik / N;       // UWAGA: zapis licznik / N * 4.0 NIE byłby poprawny. Dlaczego?
}

int main()
{
    cout << "PI oszaczowano na " << montecarlo() << endl;
    return 0;
}