#include <cstdio>
#define N 10000010
using namespace std;
const int P = 1<<16;
int n,a,b,c,x[2][N],B[P],i,mask=P-1,shift;
int main()
{
    freopen("radixsort.in","r",stdin);
    freopen("radixsort.out","w",stdout);
    scanf("%d%d%d%d",&n,&a,&b,&c);
    x[0][1]=b;
    for(i=2;i<=n;i++)x[0][i]=((1LL*a*x[0][i-1])%c+b)%c;
    a=0;c=1;b=0;
    for(shift=0;shift<=16;shift+=16)
    {
        for(i=1;i<=n;i++){b=(x[a][i]>>shift)&mask;B[b]++;}
        for(i=1;i<=mask;i++)B[i]+=B[i-1];
        for(i=n;i;i--)
        {
            b=(x[a][i]>>shift)&mask;
            x[c][B[b]]=x[a][i];
            B[b]--;
        }
        a=1-a;c=1-c;
        for(i=0;i<=mask;i++)B[i]=0;
    }
    for(i=1;i<=n;i+=10)
        printf("%d ",x[a][i]);
    return 0;
}

/*


#include <cstdio>

using namespace std;
const int N = 10000010;
int n,X[N],Y[N],p[260],*x,*y,buket;
int A,B,C;
int main()
{
    freopen("radixsort.out","w",stdout);
    freopen("radixsort.in","r",stdin);
    scanf("%d%d%d%d",&n,&A,&B,&C);
    x=X;y=Y;
    x[1]=B;
    for(int i=2;i<=n;i++)
        x[i]=(1LL*A*x[i-1]%C+B)%C;
    for(int step = 0;step<4;step++)
    {
        int shift=8*step;
        for(int i=0;i<256;i++)
            p[i]=0;
        for(int i=1;i<=n;i++)
        {
            buket = (x[i]>>shift)&255;
            p[buket]++;
        }
        for(int i=1;i<256;i++)
            p[i]+=p[i-1];
        for(int i=n;i>=1;i--)
        {
            buket=(x[i]>>shift)&255;
            y[p[buket]]=x[i];
            p[buket]--;
        }
        int *aux=x;x=y;y=aux;
    }
    for(int i=1;i<=n;i+=10)
        printf("%d ",x[i]);
    return 0;
}
*/