Cod sursa(job #2092764)

Utilizator ArctopusKacso Peter-Gabor Arctopus Data 22 decembrie 2017 12:00:32
Problema Radix Sort Scor 0
Compilator cpp Status done
Runda Arhiva educationala Marime 1.32 kb
#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