-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathMain.java
More file actions
160 lines (134 loc) · 4.09 KB
/
Copy pathMain.java
File metadata and controls
160 lines (134 loc) · 4.09 KB
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
public class Main {
public static void main(String[] args) {
InputGetter inputHandler = new InputGetter(null);
execute(inputHandler.getInput());
}
public static void execute(double[] input) {
// Calcolo somma array / 2
double middleArraySum = arraySum(input, 0, input.length - 1) / 2;
// Calcololo mediana inferiore
double result = calcInferiorMedian(input, 0, input.length - 1, middleArraySum);
System.out.println(result);
}
public static double calcInferiorMedian(double[] array, int p, int q, double middleArraySum) {
if (p > q) {
return -1;
}
int m = (p + q) / 2;
// cerco la mediana candidata
select(array, m, p, q);
// calcolo somme per effettuare verifica della mediana
double leftSum = arraySum(array, 0, m - 1);
double rightSum = arraySum(array, 0, m);
if (leftSum < middleArraySum && middleArraySum <= rightSum) {
return array[m];
} else {
if (leftSum > middleArraySum) {
return calcInferiorMedian(array, p, m - 1, middleArraySum);
} else {
return calcInferiorMedian(array, m + 1, q, middleArraySum);
}
}
}
/**
* Algoritmo partition
*/
public static int partition(double[] array, int p, int q) {
double pivot = choosePivot(array, p, q);
int i = p - 1;
int pivotIndex = seek(array, pivot, p, q);
swapCells(array, pivotIndex, q);
for (int j = p; j <= q; j++) {
if (array[j] <= pivot) {
i++;
swapCells(array, i, j);
}
}
return i;
}
/**
* algoritmo select
*/
public static double select(double[] array, int i, int p, int q) {
int x = partition(array, p, q);
if (i == x) {
return array[x];
} else if (i < x) {
return select(array, i, p, x - 1);
} else {
return select(array, i, x + 1, q);
}
}
/**
* Metodo che ordina una piccola porzione di array
*/
public static double[] sortAPart(double[] array, int begin, int endOfArray, int k) {
for (int i = begin; i < begin + (k - 1) && i < endOfArray; i++) {
for (int j = i + 1; j < begin + k && j <= endOfArray; j++) {
if (array[i] > array[j]) {
swapCells(array, i, j);
}
}
}
return array;
}
/**
* Metodo per la scelta di un pivot (mediano dei mediani)
*/
public static double choosePivot(double[] array, int p, int q) {
int n = q - p + 1;
int dim = n / 5;
if (n % 5 != 0) {
dim = dim + 1;
}
double[] B = new double[dim];
int i = p;
for (int j = 0; j < dim; j++) {
sortAPart(array, i, q, 5);
double median;
if (q - i < 5) {
median = array[i + (int) Math.floor((q - i) / 2)];
} else {
median = array[i + 2];
}
B[j] = median;
i = i + 5;
}
if (dim == 1) {
return B[0];
} else {
return select(B, (int) Math.floor(dim / 2), 0, dim - 1);
}
}
/*
* metodo per lo scambio di due celle
*/
public static double[] swapCells(double[] array, int i, int j) {
double temp;
temp = array[i];
array[i] = array[j];
array[j] = temp;
return array;
}
/**
* Metodo che ricerca un elemento in un array.
*/
public static int seek(double[] array, double x, int p, int q) {
for (int i = p; i <= q; i++) {
if (array[i] == x) {
return i;
}
}
return -1;
}
/**
* Metodo che calcola la somma degli elementi di un array
*/
public static double arraySum(double[] array, int pos1, int pos2) {
double sum = 0;
for (int i = pos1; i <= pos2; i++) {
sum = sum + array[i];
}
return sum;
}
}