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
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
|
#include <stdio.h>
#include <stdlib.h>
#include "config.h"
#include "BinHeap.h"
#define FATHER(i) (i >> 1)
#define LEFT(i) (i << 1)
#define RIGHT(i) ((i << 1) + 1)
/**********************************\
* *
* Tas binaires *
* *
\**********************************/
/*
* Implémentation quasi directe de l'algorithme Entasser(A, i) avec A = this.
* Seul le test de condition est inversé.
*/
void BinHeap::PackUp(int i)
{
int l = LEFT(i), r = RIGHT(i), min;
binheap_t *binheap = (binheap_t *) Datas;
min = ((l <= Key) && (binheap[l - 1].Key < binheap[i - 1].Key)) ? l : i;
if ((r <= Key) && (binheap[r - 1].Key < binheap[min - 1].Key))
min = r;
if (min != i) {
SWAP(binheap[i - 1], binheap[min - 1]);
PackUp(min);
}
}
/*
* Il n'y a aucune méthode "rapide" pour effectuer l'union de deux tas
* binaires. Nous utiliserons donc l'union générique.
*/
PriorityList *BinHeap::Union(PriorityList * P)
{
return GenericUnion(P);
}
/*
* Comme pour l'Union, nous allons devoir utiliser une méthode lente pour
* effectuer le diminuer clef (voir rapport)
*/
bool BinHeap::Lower_Key(Cell x, Key_t NKey)
{
int i = ((binheap_t *)x) - ((binheap_t*)Datas) + 1;
binheap_t * binheap = (binheap_t*) Datas;
if (((binheap_t*)x)->Key < NKey) return false;
while ((i > 1) && (binheap[FATHER(i) - 1].Key > NKey)) {
SWAP(binheap[i - 1], binheap[FATHER(i) - 1]);
i = FATHER(i);
}
binheap[i - 1].Key = NKey;
return true;
}
BinHeap::BinHeap(void)
{
Key = 0;
Datas = NULL;
}
BinHeap::~BinHeap(void)
{
if (Datas)
free(Datas);
}
Key_t BinHeap::ReadKey(Cell C)
{
return ((binheap_t *) C)->Key;
}
Datas_t BinHeap::ReadDatas(Cell C)
{
return ((binheap_t *) C)->Datas;
}
int BinHeap::rn(void)
{
return n();
}
void BinHeap::Dump(ostream & os)
{
int i;
binheap_t *binheap = (binheap_t *) Datas;
for (i = 0; i < Key; i++) {
os << binheap[i].Key << " ";
}
os << endl;
}
bool BinHeap::IsEmpty(void)
{
return (Key == 0);
}
Cell BinHeap::Min(void)
{
return Datas;
}
Cell BinHeap::Insert(Key_t IKey, Datas_t const &IDatas)
{
binheap_t newcell = { IKey, IDatas }, *binheap;
int i = Key++;
if (!Datas || ((((Key >> GRANUL) + 1) << GRANUL) != (((i >> GRANUL) + 1) << GRANUL))) {
if (!(Datas = realloc(Datas, (((Key >> GRANUL) + 1) << GRANUL) * sizeof(binheap_t)))) {
exception(5, _("Not enough memory"));
}
}
i = Key;
binheap = (binheap_t *) Datas;
while ((i > 1) && (binheap[FATHER(i) - 1].Key > IKey)) {
binheap[i - 1] = binheap[FATHER(i) - 1];
i = FATHER(i);
}
binheap[i - 1] = newcell;
return &(binheap[i - 1]);
}
/*
* Implémentation directe de l'algorithme EXTRAIRE-MIN-TAS(A) avec A = this.
*/
Key_t BinHeap::Extract_Min(Datas_t & Datas)
{
binheap_t *binheap = (binheap_t *) this->Datas;
Key_t ret;
if (Key < 1)
exception(5, _("negative overflow"));
ret = binheap[0].Key;
Datas = binheap[0].Datas;
binheap[0] = binheap[Key - 1];
int i = Key--;
if (!this->Datas || ((((Key >> GRANUL) + 1) << GRANUL) != (((i >> GRANUL) + 1) << GRANUL))) {
this->Datas = realloc(this->Datas, (((Key >> GRANUL) + 1) << GRANUL) * sizeof(binheap_t));
}
PackUp(1);
return ret;
}
/*
* Implémentation directe de l'algorithme TAS-BINOMIAL-SUPPRIMER(T,x) avec T = this.
* La valeur de la clef supprimée est renvoyée par la routine et les données satellites
* sont écrites dans la variable Datas.
*/
Key_t BinHeap::Delete(Datas_t & Datas, Cell x)
{
Key_t K;
K = ReadKey(x);
Lower_Key(x, M_INFINITY);
Extract_Min(Datas);
return K;
}
|
