#include<cstdio>
#include<algorithm>
#include<iostream>
using namespace std;
const int nmax=10000000;
long long ve[nmax+5],v2[nmax+5],p2[35];
char si[nmax+5];
int main()
{
    freopen("radixsort.in","r",stdin);
    freopen("radixsort.out","w",stdout);
    long long int a,b,c,d;
    int i,n,st,dr,l,j,af,x;
    af=0;
    dr=n;
    cin>>n>>a>>b>>c;
        ve[1]=b;
    d=0;
    for(i=2;i<=n;i++)
    {
        d=(d*a+b);
      if(d>=c)
        d=d%c;
        ve[i]=d;
    }
    p2[0]=1;
    for(i=1;i<=29;i++)
        p2[i]=p2[i-1]*2;
        l=1;
    for(i=0;i<=29;i++)
    {
        st=dr=l;
        for(j=l;j<=n;j++)
            if(ve[j]<p2[i])
        {
            v2[st]=ve[j];
            st++;
        }
        dr=st;
         for(j=l;j<=n;j++)
            if((ve[j]>=p2[i]) && !(ve[j] & p2[i]))
        {
            v2[dr]=ve[j];
            dr++;
        }
            for(j=l;j<=n;j++)
            {
            if((ve[j]>=p2[i]) && (ve[j] & p2[i]))
        {
            v2[dr]=ve[j];
            dr++;
        }
            }
            for(j=l;j<=n;j++)
                ve[j]=v2[j];
            l=st;
    }
    x=n/10*10+1;
    if(x>n)
        x=x-10;
for(i=x;i>=1;i=i-10)
{
    j=ve[i];
    do
    {
        af++;
        si[af]=j%10+'0';
        j=j/10;
    }while(j!=0);
    af++;
    si[af]=' ';
}
af--;
printf("%s\n",si+1);
    return 0;
}