Pagini recente » Cod sursa (job #2446975) | Cod sursa (job #2881825) | Cod sursa (job #26253) | Cod sursa (job #245483) | Cod sursa (job #2522067)
#include<bits/stdc++.h>
#include<cstdlib>
#define R 65535//2^16-1
using namespace std;
ifstream f("radixsort.in");
ofstream g("radixsort.out");
int v[2][10000002],w[10000002],a[R+2];
string sol;
inline void rez(int x)
{
if(x!=0)
{
rez(x/10);
sol+=x%10+'0';
}
}
inline void rez()
{
int n,i,A,B,C;
f>>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)
{
if(v[0][i]!=0)
rez(v[0][i]);
else
sol+='0';
sol+=' ';
}
g<<sol;
}
int main()
{
rez();
return 0;
}