Appel au générateur des nombres aleatoire dans un programme c++

Bonjour tout les monde

j'ai un programme qui consiste a générer 4000 nombre aleatoire entre 1 et 10 000 et pour chaque nombre aleatoire(n) il faut que j'exécute le P1, P2, P3. bien evidement P1 P2 et P3 vont me donner des résultats(x,y) mais pas tout ces résultats satisfaites ma condition qui exige que n(le nombre aleatoire) =x2-y2 et 0<=y<=x. et sur chaque nombre aléatoire bien sur. par la suite si la résultat fourni par chaque programme sur chaque nombre aléatoire satisfait ma condition je fais un compteur sur chaque programme qui s’incrémente a chaque fois que les sorties de programme satisfaites condition . le but de ce compteur c'est pour calculer la fiabilité a la fin qui est égale au nombre des résultats qui satisfait ma condition/ le nombre totale des nombres aleatoire(4000).

donc vue que je vais faire d'autres travaux pareil(avec d'autres programmes et d'autres condition ) donc l'idée été de faire quelque chose plus générique. mon problème ici c'est que le générateur des nombres aléatoires m'affiche toujours le même nombre aléatoire (ici c'est le n) pour chaque exécution de P1, P2, P3. ci-joint mon code et j’espère que quelqu'un peut m'aider a le corriger

#include <stdio.h>
#include <math.h>
#include <stdlib.h>
#include <time.h>
#include <iostream>
using namespace std;
//#include <iostream.h> //pour cout
typedef struct {
int n, x, y;
}statetype; // définir ma propre type de variable
bool domR (statetype s) // faire un test sur le nombre aléatoire (n) 
{
if (!((s.n/2) % 2) ==1)
return(true);
else
return(false);
}
bool R (statetype s, statetype sprime) // Le sprime est un nom de variable de type statetype . ici j'ai ecrit dans ce procédure la condition que les sorties(x,y) de chaque programme sur chaque n doit la verifier.
{
if ( s.n== pow(sprime.x,2)- pow(sprime.y,2) && 0<= sprime.y <= sprime.x)
{

return(true);
}
else
return(false);
}

bool oracle (statetype s, statetype sprime)
{ return ! domR(s) || R(s,sprime) ;
} 



void p1()
{unsigned int n, x, y; 
{unsigned int r;
x=0; y=0; r=0; 
while (r<n) {r=r+2*x+1; x=x+1;}}}


void p2() 
{unsigned int n, x, y; // input/output variables
{unsigned int r; x=0; r=0;
while (r<n) {r=r+2*x+1; x=x+1;}
if (r>n) {y=0; while (r>n) {r=r-2*y-1; y=y+1;}}}}



void p3() 
{unsigned int n, x, y; // input/output variables
{unsigned int r; x=0; r=0; 
while (r<n) {r=r+2*x+1; x=x+1;}
while (r>n) {int rsave; y=0; rsave=r;
while (r>n) {r=r-2*y-1; y=y+1;}
if (r<n) {r=rsave+2*x+1; x=x+1;}}}}


//------------------------------------------------------------------
int randomgeneration(int n)
{
int partSize = 1 + (n == RAND_MAX ? 0 : (RAND_MAX - n) / (n + 1));
int maxUsefull = partSize * n + (partSize - 1);
int draw;

do {
draw = rand();
} while (draw > maxUsefull);

return draw / partSize;
}
//---------------------------------------------------------------------
void tesdriver(int testdatasize)
{int c1=0; // compteur qui s'incrémente chaque fois que les résultats fourni par le programme satisfait la condition qui est exprimée dans le procedure bool R(N=x2-y2)et (0<=y<x)
int c2=0;
int c3=0;

int testindex=1; 

while (testindex <= testdatasize)
{

statetype inits,s, sprime; 

inits.n=randomgeneration(10000);

//P1 

s=inits;
p1();
if (oracle(inits,s)) {c1=c1+1;}
printf("n : %d \t" , inits);
printf("counter p1 : %d \t" , c1);


// P2

s=inits;
p2();
if (oracle(inits,s)) {c2=c2+1;}
printf("n : %d \t" , inits); 
printf("counter p2 : %d \t" , c2);


//P3
s=inits;
p3();
if (oracle(inits,s)) {c3=c3+1;}
printf("n : %d \t" , inits); 
printf("counter p3 : %d \t" , c3);

testindex ++;

}

//printf("\nc1=%d\t",c1);
//printf("c2=%d\t",c2);
//printf("c3=%d\t",c3);
//printf("reliability of p1 : %f \t" , float (c1)/float (testdatasize));
cout << "reliability of p1 : " << float (c1)/ float(testdatasize) << endl ; // calcul de realiabilité 
//printf("reliability of p2 : %f \t" , float (c2)/float (testdatasize));
cout << " reliability of p2 : " << float (c2)/ float(testdatasize) << endl ;
//printf("reliability of p3 : %f \t" , float (c3)/float (testdatasize));
cout << " reliability of p3 : " << float (c3)/ float(testdatasize) << endl ;
}
int main()
{
tesdriver(4000);
}