Pagini recente » Cod sursa (job #2403887) | Cod sursa (job #1215115) | Cod sursa (job #3202502) | Cod sursa (job #3040390) | Cod sursa (job #2092764)
#include <fstream>
#include <iostream>
#include <queue>
#define ll long long
using namespace std;
const int NLIM = 1e7 + 10;
const int BUCKNR = 70000;
ll N, A, B, C;
int *v;
int *v2;
int cnt[BUCKNR];
int cnt2[BUCKNR];
ll x;
ifstream fcin("radixsort.in");
ofstream fcout("radixsort.out");
int main()
{
ios::sync_with_stdio( false );
v = new int[NLIM];
v2 = new int[NLIM];
fcin >> N >> A >> B >> C;
// generate input
v[0] = B % C;
for( int i = 1; i < N; ++i )
v[i] = ( (ll)v[i - 1] * A + B ) % C;
int NR = 2;
int BYTES = 16;
ll radix = 1<<BYTES - 1;
for( int i = 0; i < NR; ++i )
{
//cout << i << "\n";
for( int j = 0; j < BUCKNR; ++j ) cnt[j] = cnt2[j] = 0;
// counting sort:
for( int j = 0; j < N; ++j )
{
x = ( v[j] >> ( i * BYTES ) ) & radix;
++cnt[x];
}
int sum = 0; cnt2[0] = 0;
for( int j = 1; j < BUCKNR; ++j )
{
sum += cnt[j - 1];
cnt2[j] = sum;
}
int vi = 0;
// get the numbers from the buckets:
for( int j = 0; j < N; ++j )
{
x = ( v[j] >> ( i * BYTES ) ) & radix;
v2[cnt2[x]] = v[j];
++cnt2[x];
}
int* emp = v;
v = v2;
v2 = emp;
}
for( int i = 0; i < N; ++i )
if( i % 10 == 0 )
fcout << v[i] << " ";
fcout << "\n";
}
// 100 12 38 123
// 10 1 1 100