#include <iostream>
#include <iomanip>
using namespace std;

void Display (int arg[], int length); //prototipo
void Display (double arg[], int length); //prototipo

  int main () {

    // algoritmo di generazione di numeri interi  pseudo-casuali tra 0 e (m-1) : 
    // r_i = mod(a*r_(i-1) +c, m);  
    // mod (x,y) = x - int(x/y)*y   ritorna il resto della divisione tra x e y 

    int seed ;
    cout  << "Inserisci un intero come seme: ";
    cin>> seed;
    int m = 9; // periodo della sequenza; deve essere il piu' grande possibile 
                // [2^(32) -1 in una macchina a 32 bit]
    int c = 1; // parametro da ottimizzare
    int a = 4; //  parametro da ottimizzare (una scelta piu' opportuna sarebbe  a = 16807)
    int N = 25;  // numero di eventi da simulare
    int numeriInteriGenerati[N]; 
    double numeriDoubleGenerati[N];  //tra 0 e 1 
    numeriInteriGenerati[0] = seed;
    for (int n=1; n<N; n++)
      {
	// x % y da il resto della divisione di x con y
	numeriInteriGenerati[n] = (a * numeriInteriGenerati[n-1] + c) %  m;
      }
    
    for (int n=0; n<N; n++)
      {
	numeriDoubleGenerati[n] = double(numeriInteriGenerati[n])/double(m-1);
      }
    


    cout << "Numeri interi pseudo-random generati = "<< N <<" con periodo = " << m << " e con seed = " << seed << endl; 
    Display (numeriInteriGenerati, N);
    
    cout << "Numeri double pseudo-random generati tra 0 e 1 = "<< N <<" con periodo = " << m << " e con seed = " << seed << endl; 
    Display (numeriDoubleGenerati, N);
    return (0);
}
void Display (int arg[], int length) {
  for (int n=0; n<length; n++) { cout << setw(5) << arg[n]; } 
  cout << endl;  
}

void Display (double arg[], int length) {
  for (int n=0; n<length; n++) { cout << setw (10) << setprecision (4) << arg[n]; } 
  cout << endl;  
}