1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
|
/**************************************************************************
randomlist.cpp - Some "random functions" :-)
Copyright (C) 1997,98 Antonio Larrosa Jimenez
This program is free software; you can redistribute it and/or modify
it under the terms of the GNU General Public License as published by
the Free Software Foundation; either version 2 of the License, or
(at your option) any later version.
This program is distributed in the hope that it will be useful,
but WITHOUT ANY WARRANTY; without even the implied warranty of
MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
GNU General Public License for more details.
You should have received a copy of the GNU General Public License
along with this program; if not, write to the Free Software
Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
Send comments and bug fixes to larrosa@kde.org
or to Antonio Larrosa, Rio Arnoya, 10 5B, 29006 Malaga, Spain
***************************************************************************/
#include "randomlist.h"
#include <stdlib.h>
#include <stdio.h>
#define RAND_UNIFORM (double)rand()/(double)RAND_MAX
int random_discrete(double *distrib,int n)
{
int i=0;
double g=0.0;
double z=0.0;
while ((z==0.0)||(z==1.0)) z=RAND_UNIFORM;
while ((g<z)&&(i<n))
{
g+=distrib[i];
i++;
};
return i-1;
}
double *generate_discrete_uniform_distrib(int n)
{
double *distrib=new double[n];
for (int i=0;i<n;i++)
distrib[i]=1.0/(double)n;
return distrib;
}
void remove_lmn_from_discrete_distrib(int k,double *distrib,int n,int used=0)
{
int i;
if (used==0) //guess how many values of the v.a. are strictly positive
{
for (i=0;i<n;i++) if (distrib[i]>0) used++;
};
used--; // k will no longer be used :-)
if (used==0) return;
double piece=distrib[k]/(double)used;
distrib[k]=0.0;
for (i=0;i<n;i++) if (distrib[i]>0) distrib[i]+=piece;
}
void show_distrib(double *distrib,int n)
{
printf("(");
for (int j=0;j<n;j++) printf("%g,",distrib[j]);
printf(")\n");
}
int *generate_random_list(int n)
{
if (n==0) return NULL;
int *list=new int[n];
double *distrib=generate_discrete_uniform_distrib(n);
int used=n;
int x;
int i=1;
while (used>0)
{
x=random_discrete(distrib,n);
list[x]=i;
i++;
remove_lmn_from_discrete_distrib(x,distrib,n,used);
used--;
};
delete distrib;
return list;
}
int *generate_list(int n)
{
int *list=new int[n];
for (int i=0;i<n;i++)
list[i]=i+1;
return list;
}
|
