#include<cstdio>
using namespace std;
const int nmax = 10000005;
const int mask = (1<<8)-1;
int n,a,b,c,i,k,v[nmax],w[nmax],shift,buc[mask+5];
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;
	for(shift=0;shift<32;shift+=8)
	{
	    for(i=1;i<=n;i++)
	    {
	        k=(v[i]>>shift)&mask;
	        buc[k]++;
	    }
	    for(i=1;i<=mask;i++) buc[i]+=buc[i-1];
	    for(i=n;i;i--)
	    {
	        k=(v[i]>>shift)&mask;
	        w[buc[k]--]=v[i];
	    }
	    for(i=0;i<=mask;i++) buc[i]=0;
	    for(i=1;i<=n;i++) v[i]=w[i];
	}
	for(i=1;i<=n;i+=10) printf("%d ",v[i]);
	return 0;
}