Pagini recente » Cod sursa (job #3161464) | Cod sursa (job #60521) | Cod sursa (job #1871015) | Cod sursa (job #301143) | Cod sursa (job #2945116)
///radixsort-LSD variant
#include <iostream>
#include <fstream>
#define Nmax 10000000
#define MaxBit 32
#define BitsPerStep 8
#define Mask ((1<<BitsPerStep)-1)
using namespace std;
ifstream fin ("radixsort.in");
ofstream fout ("radixsort.out");
int n,a,b,c;
int v1[Nmax], v2[Nmax], ind[Mask+1], fr[Mask+1];
void radixsort(int v[], int aux[], int bit){
if (bit==MaxBit)
return;
for (int i=0;i<=Mask;i++)
fr[i]=0;
for (int i=0;i<n;i++)
fr[(v[i]>>bit)&Mask]++;
ind[0]=0;
for (int i=1;i<=Mask;i++)
ind[i]=ind[i-1]+fr[i-1];
for (int i=0;i<n;i++)
aux[ind[(v[i]>>bit)&Mask]++]=v[i];
radixsort(aux, v, bit+BitsPerStep);
}
int main()
{
fin>>n>>a>>b>>c;
v1[0]=b;
for (int i=1;i<n;i++){
v1[i]=((long long)a*v1[i-1]+b)%c;
}
radixsort(v1, v2, 0);
for (int i=0;i<n;i+=10)
fout<<v1[i]<<' ';
return 0;
}