Pagini recente » Cod sursa (job #1760824) | Cod sursa (job #1162661) | Cod sursa (job #1803446) | Cod sursa (job #811316) | Cod sursa (job #2522055)
#include<bits/stdc++.h>
#define R 65535//2^16-1
using namespace std;
FILE *f=fopen("radixsort.in","r");
FILE *g=fopen("radixsort.out","w");
int v[2][10000002],w[10000002],a[R+2];
inline void rez()
{
int n,i,A,B,C;
fscanf(f,"%d%d%d%d",&n,&A,&B,&C);
v[0][1]=B;
for(i=2;i<=n;i++)
v[0][i]=(1LL*v[0][i-1]*A+B)%C;
for(i=1;i<=n;i++)
{
w[i]=v[0][i]&R;
a[w[i]]++;
}
for(i=1;i<=R;i++)
a[i]+=a[i-1];
for(i=n;i>=1;i--)
{
v[1][a[w[i]]]=v[0][i];
a[w[i]]--;
}
for(i=0;i<=R;i++)
a[i]=0;
for(i=1;i<=n;i++)
{
w[i]=(v[1][i]>>16)&R;
a[w[i]]++;
}
for(i=1;i<=R;i++)
a[i]+=a[i-1];
for(i=n;i>=1;i--)
{
v[0][a[w[i]]]=v[1][i];
a[w[i]]--;
}
for(i=1;i<=n;i+=10)
fprintf(g,"%d ",v[0][i]);
}
int main()
{
rez();
return 0;
}