Pagini recente » Cod sursa (job #2633447) | Cod sursa (job #2429051) | Cod sursa (job #2259379) | Cod sursa (job #2429050) | Cod sursa (job #1249332)
#include <cstdio>
using namespace std;
const int NMAX = 1000000;
inline int get_digit(const int value, const int base, const int exp) {
return (value / exp) % base;
}
void radix_sort(int *A, const int len) {
const int BASE = 1024;
// max_A element
long long max_A = A[0];
int i;
for (i = 1; i < len; i++)
if (A[i] > max_A) max_A = A[i];
int B[NMAX];
for (long long exp = 1; exp < max_A; exp *= BASE) {
int bucket[BASE] = {0};
// Count elements by digit
for (i = 0; i < len; i++)
bucket[ get_digit(A[i], BASE, exp) ]++;
// Add the count from behind to use later in ordering
for (i = 1; i < BASE; i++)
bucket[i] += bucket[i-1];
// Put elements in order
for (i = len-1; i >= 0; i--)
B[ --bucket[ get_digit(A[i], BASE, exp) ] ] = A[i];
// Copy elements from B to A
for (i = 0; i < len; i++)
A[i] = B[i];
}
}
int main() {
int N, A, B, C;
FILE * in = fopen("radixsort.in", "r");
fscanf(in, "%d%d%d%d", &N, &A, &B, &C);
fclose(in);
// Generate elements of array
int v[NMAX] = {B};
for (int i = 1; i < N; i++)
v[i] = (1LL * A * v[i-1] + B) % C;
// Sort generated array
radix_sort(v, N);
// Output elements
FILE * out = fopen("radixsort.out", "w");
for (int i = 0; i < N; i += 10)
fprintf(out, "%d ", v[i]);
fclose(out);
return 0;
}