#include <iostream>
#include <fstream>

using namespace std;

ifstream in("radixsort.in");
ofstream out("radixsort.out");

typedef long long ll;

const int NB = 16, N = 10000005, NC = 2, B = 1 << NB;
int poz[B], nr[B], cif, v[2][N];

int main()
{
    int n;
    ll a, b, c;
    in >> n >> a >> v[0][0] >> c;
    int p = 1;
    for(int i = 1; i < n; i++){
        v[0][i] = (a * v[0][i - 1] + v[0][0]) % c;
    }
    for(int k = 0; k < NC; k++){
        int nb = k * NB;
        for(int j = 0; j < B; j++){
            nr[j] = 0;
        }
        for(int i = 0; i < n; i++){
            int cif = (v[k & 1][i] >> nb) & (B - 1);
            nr[cif]++;
        }
        poz[0] = 0;
        for(int j = 1; j < B; j++){
            poz[j] = poz[j - 1] + nr[j - 1];
        }
        for(int i= 0; i < n; i++){
            int cif = (v[k & 1][i] >> nb) & (B - 1);
            v[1][poz[cif]++] = v[0][i];
        }
        for(int i = 0; i < n; i++){
            v[0][i] = v[1][i];
        }
        p *= B;
    }
    for(int i = 0; i < n; i += 10){
        out << v[0][i] << " ";
    }

    return 0;
}