Pagini recente » Cod sursa (job #730493) | Cod sursa (job #2808413) | Cod sursa (job #2593710) | Cod sursa (job #3308074) | Cod sursa (job #3327965)
#include <stdio.h>
#include <string.h>
#define NMAX 10000000
#define BITS_PER_STEP 4
#define BASE (1 << BITS_PER_STEP)
#define MASK (BASE - 1)
int v[NMAX], freq[BASE], aux[NMAX];
void radixSort( int v[], int begin, int end, int bit ) {
int i, cp;
if ( begin < end && bit >= 0 ) {
memset( freq, 0, sizeof( freq ) );
for ( i = begin; i <= end; i++ )
freq[(v[i] >> bit) & MASK]++;
for ( i = 1; i < BASE; i++ )
freq[i] += freq[i - 1];
for ( i = begin; i <= end; i++ ) {
--freq[cp = (v[i] >> bit) & MASK];
aux[freq[cp] + begin] = v[i];
}
for ( i = begin; i <= end; i++ )
v[i] = aux[i];
i = 0;
for ( i = 0; i < end; i++) {
cp = i;
while ( i < end && ((v[i] >> bit) & MASK) == ((v[i + 1] >> bit) & MASK) )
i++;
radixSort( v, cp, i, bit - BITS_PER_STEP );
}
}
}
int main() {
FILE *fin, *fout;
int n, a, b, c, i;
fin = fopen( "radixsort.in", "r" );
fscanf( fin, "%d%d%d%d", &n, &a, &b, &c );
fclose( fin );
v[0] = b;
for ( i = 1; i < n; i++ )
v[i] = ((long long)a * v[i - 1] + b) % c;
radixSort( v, 0, n - 1, 28 );
fout = fopen( "radixsort.out", "w" );
for ( i = 0; i < n; i += 10 )
fprintf( fout, "%d ", v[i] );
fputc( '\n', fout );
fclose( fout );
return 0;
}