Pagini recente » Cod sursa (job #1797462) | Cod sursa (job #1997893) | Cod sursa (job #2517618) | Cod sursa (job #2571821) | Cod sursa (job #1842877)
#include <cstdio>
#define NMAX 10000003
#define grupa(nr, step) ((nr >> step) & 15)
using namespace std;
int s[2][16][NMAX];
int ls[2][16];
int n, a, b, c;
int pc = 0, lc = 0, *vc, indc = -1;
int mc = 1, mp;
inline void next(int& val)
{
while(pc == lc)
{
vc = s[mp][++indc];
lc = ls[mp][indc];
pc = 0;
}
val = vc[pc++];
}
int main()
{
int i, gr = 0, aux, prev;
freopen("radixsort.in", "r", stdin);
freopen("radixsort.out", "w", stdout);
scanf("%d%d%d%d", &n, &a, &b, &c);
prev = b;
s[0][prev & 15][ls[0][prev & 15]++] = prev;
//printf("%d ", prev);
for(i = 1; i < n; i++)
{
aux = (a * prev + b) % c;
//printf("%d ", aux);
s[0][aux & 15][ls[0][aux & 15]++] = aux;
prev = aux;
}
//printf("\n");
for(gr = 4; gr != 32; gr += 4)
{
for(i = 0; i < 16; i++) ls[mc][i] = 0;
for(i = 0; i < n; i++)
{
next(aux);
s[mc][grupa(aux, gr)][ls[mc][grupa(aux, gr)]++] = aux;
}
pc = 0; lc = 0; indc = -1;
mc = mp;
mp = (mp == 1) ? 0 : 1;
}
mp = mc;
for(i = 0; i < n; i++)
{
next(aux);
if(i % 10 == 0) printf("%d ", aux);
}
return 0;
}