Prova questo codice, per i numeri a ed m ho messo i primi che mi venivano in mente. In teoria andrebbero scelti con cura per far si che il generatore abbia determinate proprieta' .
Ho fatto pure una verifica che il numero sia compreso tra 0 e 1, infatti la media è circa 0.5. Puoi verificare che si avvicina a 0.5 aumentando il numero di iterazioni . Puoi giocare anche sul seme. NB Il cuore del codice l'ho inplementato nelle sole 2 righe della funzione: generaRandom().
Carino no?
#include <cstdlib>
#include <iostream>
#include <cmath>
using namespace std;
int generaRandom(int a, int m, int n)
{
int n_piu_uno;
n_piu_uno=(a*n) - m*(a*n/m); // funzione resto della divisione tra ;
// (a*n) ed m ; Z_nres(a*Z_(n-1), m);
return n_piu_uno;
}
int main()
{
int a=2345;
int m=253256;
int nIterazioni=100;
double x;
int prossimo;
int n;
int seme=87432;
n=seme;
double somma=0;
for (int i=0; i<nIterazioni; i++)
{
prossimo=generaRandom(a, m, n);
double reale=double(prossimo)/double(m); //è compreso tra 0 e 1;
cout<<"iterazione: "<<i+1<<" intero estratto= "<<prossimo<<
" reale = "<<reale<<endl;
n=prossimo;
somma=somma+reale;
}
double media=double(somma)/double(nIterazioni);
cout<<"media= "<<media;
system("PAUSE");
return EXIT_SUCCESS;
}