C++

[C] gestion des nombres aléatoires

Recherche :

Mot : Pseudo : Filtrer
Bas de page
Auteur	Sujet : [C] gestion des nombres aléatoires

Aricoh

gentil mais fo po pousser

Bonjour !

Comment faites-vous pour gérer les nombres aléatoires en C ???

Dans le cadre d'un projet perso, j'ai besoin de concevoir un tout petit code C générant une série de 10 nombres aléatoires compris entre 1 et 10. Un pote m'a filé un bout de code pour générer de l'aléatoire et j'ai fait ce petit exemple que voici :

Code :

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#define D10 ((double)rand() / ((double)RAND_MAX + 1) * 10) + (1)
int main(void) {
int i = 10;
srand (time(NULL));
while(i) {
printf("%d;", D10);
--i;
}
printf("\n" );
return 0;
}

Les chiffres générés sont bizarres, il me manque certainement un cast en int car je récupère des valeurs du genre 32765548, -284363, etc à chaque appel de D10.

Mais c'est pas ça qui m'inquiète, ce qui me gène bcp c'est que lorsque j'appelle le binaire du programme n fois, je retrouve les mêmes séries très régulièrement.

Or, je cherche surtout la manière de générer des nombres aléatoires qui soient VRAIMENT ALEATOIRES.

Comment faire ?

Aricoh.

PS : pros du C, soyez indulgent, je fais pas du C très souvent, mon domaine à moi ---> Perl

---------------
Samsung Galaxy S1 -> Samsung Galaxy S2 -> Samsung Note 2 -> Huawei Ascend Mate 7 -> ZTE Axon 7 -> OnePlus 6T -> Oppo Find X2 PRO -> Google Pixel 9 PRO XL

Publicité

Taz

bisounours-codeur

http://dejean.benoit.free.fr/rand.c.html

Aricoh

gentil mais fo po pousser

++Taz a écrit :

http://dejean.benoit.free.fr/rand.c.html

Merci Taz

Free merde apparemment, connexion impossible à l'url

---------------
Samsung Galaxy S1 -> Samsung Galaxy S2 -> Samsung Note 2 -> Huawei Ascend Mate 7 -> ZTE Axon 7 -> OnePlus 6T -> Oppo Find X2 PRO -> Google Pixel 9 PRO XL

Taz

bisounours-codeur

Code :

#include <stdio.h>
#include <stdlib.h> /* pour srand, rand et RAND_MAX */
#include <time.h> /* pour time */
#include <math.h>
unsigned my_bad_rand(unsigned N)
{
/**
* les bits de poids faibles ont une distribution tres peu aléatoire
* d'ou des séquences 0 1 0 1 0 1 0 1 0
* ce n'est donc pas une bonne méthode
* cette version ne peut renvoyée qu'un entier compris entre 0 et min(RAND_MAX, N)
* d'ou un usage restreint
*/
return (unsigned)rand()%N;
}
unsigned my_better_rand(unsigned N)
{
/**
* beaucoup plus aléatoire et renvoie un entier entre 0 et N
* (la je fais pas l'explication :oD )
*/
return (unsigned)((double)rand()/((double)RAND_MAX+1)*N);
}
double average(const unsigned *tab, unsigned size)
{
double sum=0.0;
unsigned i;
for(i=0U; i<size; ++i)
{
sum+=tab[i];
}
return sum/size;
}
int main()
{
/**
* initialisation du générateur de nombres pseudos-aleatoires
* on utilise généralement le temps courant (nombre de secondes ecoulees
* depuis 01/01/1970)
* srand => Seed == 'graine'
*/
srand(time(NULL));
/* petite routine de test */
const unsigned N= time(NULL) % RAND_MAX;
printf("comparaison entre my_bad_rand et my_better_rand\n"
"nombre d\'iterations (bigger is better) : " );
unsigned i, size;
scanf("%u", &size);
unsigned *tab;
if((tab=malloc(size*sizeof(unsigned))) == NULL)
{
perror("malloc" );
exit(EXIT_FAILURE);
}
for(i=0U; i<size; ++i)
{
tab[i]=my_bad_rand(N);
}
double bad_avg=average(tab, size) / N;
for(i=0U; i<size; ++i)
{
tab[i]=my_better_rand(N);
}
double better_avg=average(tab, size) / N;
free(tab);
printf("my_bad_rand moyenne : %f\n", bad_avg);
printf("my_better_rand moyenne : %f\n", better_avg);
printf("my_%s_rand gagne !!! (plus aleatoire)\n", (fabs(bad_avg-0.5) < fabs(better_avg-0.5) ? "bad" : "better" ));
return EXIT_SUCCESS;
}

Aricoh

gentil mais fo po pousser

:jap: ++Taz

---------------
Samsung Galaxy S1 -> Samsung Galaxy S2 -> Samsung Note 2 -> Huawei Ascend Mate 7 -> ZTE Axon 7 -> OnePlus 6T -> Oppo Find X2 PRO -> Google Pixel 9 PRO XL

Taz

bisounours-codeur

c juste un petit programme sans doute buggé pour montrer l'interet

BifaceMcLeOD

The HighGlandeur

Aricoh a écrit :

Bonjour !

Comment faites-vous pour gérer les nombres aléatoires en C ???

(...)

Mais c'est pas ça qui m'inquiète, ce qui me gène bcp c'est que lorsque j'appelle le binaire du programme n fois, je retrouve les mêmes séries très régulièrement.

Or, je cherche surtout la manière de générer des nombres aléatoires qui soient VRAIMENT ALEATOIRES.

Comment faire ?

Vraiment aléatoires, c'est impossible avec la plupart des machines courantes. Comme leur documentation l'indique, les fonctions rand() et consort donnent accès à un générateur de nombres pseudo-aléatoires. Leur but est de donner des nombres apparemment non liés entre eux. En fait, ce sont les termes d'une suite mathématique toute bête.

Donc si tu pars toujours avec le même premier terme de la suite, tu obtiendras toujours les mêmes termes suivants d'une exécuation à l'autre. C'est assez pratique, ceci dit, lorsque tu as un bug dans un programe qui manipule des nombres pseudo-aléatoires, et que tu souhaites reproduire le bug...

Cependant, il est assez facile d'obtenir des séquences de nombres pseudo-aléatoires qui changent à chaque exécution. Il suffit de fixer le premier terme de la suite (i.e. initialiser le générateur de nombres pseudo-aléatoires) à une valeur différente à chaque exécution.

Pour initialiser le générateur, on utilise la fonction "srand(int seed)". Et pour trouver une valeur différente à chaque exécution, on utilise généralement la valeur que retourne la fonction "time()".

Le seul problème résiduel dépend du nombre de tirages que tu fais avec le générateur de nombres pseudo-aléatoires. Si ce nombre de tirages est inférieur à, disons, 5000, alors il n'y a aucun souci. Au-delà, tu risques de retrouver une certaine régularité dans les tirages, et tu devras alors utiliser un autre générateur de nombres pseudo-aléatoires, qui fonctionne avec une résolution supérieure. Le générateur de la librairie standard C est un générateur 16 bits. Mais si tu es sur UNIX, tu dois pouvoir également utiliser un générateur 48 bits (fais un "man rand48" pour plus d'infos).

Message édité par BifaceMcLeOD le 03-06-2003 à 12:23:21

Taz

bisounours-codeur

/dev/random

gatorette

Juste une petite précision (j'ai eu l'occasion de m'y intéresser récemment), tu peux également utiliser d'autres algorithmes de génération de nombre aléatoire. Je pense notamment au Mersenne Twister et tu trouveras pas mal d'infos (si ça t'intéresse) sur le site du pLab.

---------------
each day I don't die is cheating

Taz

bisounours-codeur

ou une lava lamp

Publicité

Aricoh

gentil mais fo po pousser

Je crois que je vais faire au plus rapide et plus sûr : un module Perl qui gère ces tirages aléatoires.

Après moults tests, ce qui pose problème, c'est l'initialisation de l'horloge système. Dès lors que j'insère un srand() dans mon code C, dès que je le lance n fois via un script Perl, je récupère toujours la même série de nombres. Idem si le source C ne contient aucun srand(). Idem si je place mon srand() dans le script Perl et pas dans le source C.

Le but de la manip' était de pouvoir exécuter un binaire faisant ces jets aléatoires, les jets étant retournés au programme appelant.

---------------
Samsung Galaxy S1 -> Samsung Galaxy S2 -> Samsung Note 2 -> Huawei Ascend Mate 7 -> ZTE Axon 7 -> OnePlus 6T -> Oppo Find X2 PRO -> Google Pixel 9 PRO XL

usa_satriani

S.P.Q.R.

Salut !
mon prof m'a filé ces deux truc : alea.h

Code :

#ifndef ALEA_H
#define ALEA_H
void init_alea ();
void init_alea_graine (long graine);
double alea_double();
int alea_int(int maximum);
#endif

et alea.c

Code :

#include <stdlib.h>
#include <sys/time.h>
#include "alea.h"
void init_alea () {
struct timeval nowtime;
gettimeofday(&nowtime, NULL);
srandom(nowtime.tv_sec^nowtime.tv_usec);
}
void init_alea_graine (long graine) {
srandom(graine);
}
double alea_double() {
double res;
res = (random() % 100000) / 100000.;
return res;
}
int alea_int(int maximum) {
return (alea_double() * maximum);
}

problèmes : à lécole on bosse sous Sun et moi chui sous windows alors je sais pas si c ca mais dev c ++ me renvoie cette erreur :

Citation :

In file included from f:\langc\projet\projet.c:2:
f:\langc\projet\alea.c: In function `init_alea':
f:\langc\projet\alea.c:6: storage size of `nowtime' isn't known

---------------
Ce monde n'est qu'une vaste entreprise à se foutre du monde. Céline

Angel_Dooglas

Le dernier des humains

Aricoh a écrit :

Etrange qu'il te rebalance la meme suite meme sans aucune graine (je ne me rappelle plus de l'usage exacte de rand, vu que rand c'est le mal, un tres mauvais generateur). Remarque qu'a priori il repart depuis le debut de la suite a chaque appel (me semblait pas que c'etait le cas mais bon...).
Si tu fais des appels en cascade a ton programme C (et que celui ci est rapide a executer en plus), il faut nourrir la graine avec autre chose que time (par exemple si tu disposes de lecture de "parasites" sur un signal, ca peut le faire) ou alors faire des sleep (le plus facile et le moins elegant mais si le temps machine n'est pas essentiel pourquoi pas) ou encore gerer le pointeur de la suite du rand de maniere a partir de la ou tu t'etais arrete (c'est un peu plus complique).

Edit: utilise mersenne twister pour des applis "serieuses".

Message édité par Angel_Dooglas le 06-06-2003 à 02:17:53

TheTooN

moi g fé un script pour ca pour le taf, ca exporte ds un fichier rapport.dat :

#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <conio.h>
main(){
printf(" -= p@ssWoRd g3n3r@t0r PoWeReD by t00n 2003 =-\n" );
srand(time(NULL));
int taille, gener, i, j;
FILE *rapport;
rapport = fopen("rapport.dat", "r+" );
taille = 5;
for(j=0; j<10; j++){
for(i=0; i<taille; i++){
gener = rand();
if(gener<(RAND_MAX / 2)){
while(gener>122){ gener = gener - 26; }
} else {
while(gener>57){ gener = gener - 9; }
}
printf("%c", gener);
fprintf(rapport, "%c", gener);
}
printf("\n" );
fprintf(rapport, "\n" );
}
fclose(rapport);

fclose(passgener);
puts("\nPress a key to finish." ); getch();
}

et il fo ke tu ais un fichier rapport.dat vide qui existe paske je me suis po pris la tete à faire des tests et tt et tt vu ke il y a ke moi ki utilise ca.

Message édité par TheTooN le 06-06-2003 à 10:54:36

Taz

bisounours-codeur

c'est loin d'etre un bon générateur dis-donc

TheTooN

ben té dur toi, ca genere des chaines alphanumériques de 5 caractères et ca marche plutot bien - pkoi cé po top ?

Taz

bisounours-codeur

niveau aléa c'est pas le top

essaye ça

Code :

unsigned my_better_rand(unsigned N)
{
/**
* beaucoup plus aléatoire et renvoie un entier entre 0 et N
* (la je fais pas l'explication :oD )
*/
return (unsigned)((double)rand()/((double)RAND_MAX+1)*N);
}

tu peux faire par exemple

Code :

c= 'a'+ my_better_rand(26);
d= '0'+ my_better_rand(10);
// ....
if(my_better_rand(2))
{
c= 'a'+ my_better_rand(26);
}
else
{
c= '0'+ my_better_rand(10);
}
// ou bien
c=my_better_rand(10+26);
c+= (c < 10 ? '0' : 'a'-10);

edit: et ça bouffera moins de CPU surtout

Message édité par Taz le 06-06-2003 à 11:25:26

LetoII

Le dormeur doit se réveiller

++Taz a écrit :

niveau aléa c'est pas le top

essaye ça

Code :

unsigned my_better_rand(unsigned N)
{
/**
* beaucoup plus aléatoire et renvoie un entier entre 0 et N
* (la je fais pas l'explication :oD )
*/
return (unsigned)((double)rand()/((double)RAND_MAX+1)*N);
}

tu peux faire par exemple

Code :

c= 'a'+ my_better_rand(26);
d= '0'+ my_better_rand(10);
// ....
if(my_better_rand(2))
{
c= 'a'+ my_better_rand(26);
}
else
{
c= '0'+ my_better_rand(10);
}
// ou bien
c=my_better_rand(10+26);
c+= (c < 10 ? '0' : 'a'-10);

edit: et ça bouffera moins de CPU surtout

De toute façon rand est pas terrible comem générateur alros bon...

---------------
Le Tyran

Taz

bisounours-codeur

alors autant optimiser son utilisation n'est ce pas

LetoII

Le dormeur doit se réveiller

++Taz a écrit :

alors autant optimiser son utilisation n'est ce pas

Je parlais au niveau vitesse, juste d'un point de vue qualité des nombres pseudo aléatoire fournis. Après j'ai pas regardé les codes en détails

---------------
Le Tyran

FORUM HardWare.fr

Programmation

C++

[C] gestion des nombres aléatoires

Sujets relatifs
[Java]Gestion de sources...	2 questions : gestion des exceptions et paramètres des fonctions
probleme avec la gestion d'evenement dans une balise div	Gestion Clients
Gestion de logs, timestamp, TStringList et memo	gestion bdd
[ASP] Gestion de forumulaire avec nombreux champs	PGCD de plus de 2 nombres
J'aimerai bien faire un truc de ce style..(Gestion graphique d'objets)	vous connaissez un outil de "gestion de version" a verrouillage
Plus de sujets relatifs à : [C] gestion des nombres aléatoires

Page générée en 0.139 secondes