Pagini recente » Arbori de intervale si aplicatii in geometria computationala | Monitorul de evaluare | Cod sursa (job #456056) | Cod sursa (job #2473840) | Cod sursa (job #1251856)
#include <fstream>
#include <cstring>
#define Nmax 10100100
#define Bmax (1 << 8)
#define Byte(x) ((x >> Step) & (Bmax - 1))
using namespace std;
int N, A, B, C, V[Nmax], Tmp[Nmax];
void radixSort() {
int i, Step, Size[Bmax], Index[Bmax];
for(Step = 0; Step <= 24; Step += 8) {
memset(Size, 0, sizeof(Size));
for(i = 1; i <= N; i++) {
Tmp[i] = V[i];
Size[Byte(V[i])]++;
}
Index[0] = 1;
for(i = 1; i < Bmax; i++)
Index[i] = Index[i - 1] + Size[i - 1];
for(i = 1; i <= N; i++)
V[Index[Byte(Tmp[i])]++] = Tmp[i];
}
}
void Solve() {
int i;
V[1] = B;
for(i = 2; i <= N; i++)
V[i] = (1LL * A * V[i-1] + B) % C;
radixSort();
}
void Read() {
ifstream in("radixsort.in");
in >> N >> A >> B >> C;
in.close();
}
void Write() {
ofstream out("radixsort.out");
for(int i = 1; i <= N; i += 10)
out << V[i] << ' ';
out << '\n';
out.close();
}
int main() {
Read();
Solve();
Write();
return 0;
}