Pagini recente » Cod sursa (job #1477862) | Cod sursa (job #1979022)
#include <iostream>
#include <fstream>
#include <queue>
#include <vector>
using namespace std;
ifstream in("radixsort.in");
ofstream out("radixsort.out");
#define pb push_back
#define ll long long
const int NMax = 1e7 + 5;
const int nrBytes = 4;
const int Radix = 0xff;
const int bits = 8;
ll N,A,B,C;
int v1[NMax],v2[NMax];
void countingSort(int*,int*,int);
int main() {
in>>N>>A>>B>>C;
v1[1] = B;
B %= C;
A %= C;
for (int i=2;i <= N;++i) {
v1[i] = (A * v1[i-1] + B) % C;
}
for (int i=0;i < nrBytes;++i) {
if (i & 1) {
countingSort(v2,v1,i);
}
else {
countingSort(v1,v2,i);
}
}
for (int i=1;i <= N;i += 10) {
out<<v1[i]<<' ';
}
out<<'\n';
in.close();
out.close();
return 0;
}
void countingSort(int* one,int* two,int byte) {
int occ[(1<<bits) + 1] = {};
for (int i=1;i <= N;++i) {
++occ[(one[i]>>(bits*byte)) & Radix];
}
int total = 0;
for (int i=0;i <= Radix;++i) {
int old = occ[i];
occ[i] = total + 1;
total += old;
}
for (int i=1;i <= N;++i) {
two[ occ[(one[i]>>(bits*byte)) & Radix]++ ] = one[i];
}
}