#include<cstdio>
using namespace std;
const int NMAX = 10000005;
const int MASK = (1<<8)-1;
int N,A,B,C,v[NMAX],w[NMAX],i,j,k,shift,Buckets[MASK+5];
int *V,*W,*Aux;
int main()
{
    freopen("radixsort.in","r",stdin);
    freopen("radixsort.out","w",stdout);
    scanf("%d%d%d%d",&N,&A,&B,&C);
    for(v[1]=B,i=2;i<=N;i++) v[i]=(1LL*A*v[i-1]+B)%C;
    V=v;
    W=w;
    for(i=0,shift=0;i<4;i++,shift+=8)
    {
        for(j=1;j<=N;j++)
        {
            k=(V[j]>>shift)&MASK;
            Buckets[k]++;
        }
        for(j=1;j<=MASK;j++)
            Buckets[j]+=Buckets[j-1];
        for(j=N;j;j--)
        {
            k=(V[j]>>shift)&MASK;
            W[Buckets[k]]=V[j];
            Buckets[k]--;
        }
        Aux=V;
        V=W;
        W=Aux;
        for(j=0;j<=MASK;j++) Buckets[j]=0;
    }
    for(i=1;i<=N;i+=10) printf("%d ",V[i]);
    return 0;
}