#include <stdio.h>

#define MAX 10000000
int v[MAX];

#define BITS 4
#define MASK 0xf
#define TINY 0xf

inline void insertion(int lef, int rig) {
  for (int i = lef; i < rig; ++i) {
    int t, j; for (t = v[i], j = i - 1; j >= lef; --j)
      if (t < v[j]) v[j + 1] = v[j]; else break;
    v[j + 1] = t;
  }
}

void radix(int lef, int rig, char bit) {
  int x[MASK + 1], y[MASK + 1];
  y[0] = lef; for (int i = 1; i <= MASK; ++i) y[i] = 0;
  for (int i = lef; i < rig; ++i) ++y[v[i] >> bit & MASK];
  for (int i = 1; i <= MASK; ++i) y[i] += y[i - 1];
  for (int i = 1; i <= MASK; ++i) x[i] = y[i - 1]; x[0] = lef;
  for (int i = 0; i <= MASK; ++i)
    for (; x[i] < y[i]; ++x[i])
      for (int j = v[x[i]] >> bit & MASK; j != i; j = v[x[i]] >> bit & MASK) {
        int tmp = v[x[i]]; v[x[i]] = v[x[j]]; v[x[j]++] = tmp; }
  if (bit <= 0) return;
  for (int i = 1; i <= MASK; ++i) x[i] = y[i - 1]; x[0] = lef;
  for (int i = 0; i <= MASK; ++i) {
    if (y[i] - x[i] < 2) continue;
    if (y[i] - x[i] < TINY) insertion(x[i], y[i]);
    else radix(x[i], y[i], bit - BITS);
  }
}

inline int read() {
  int n, a, b, c; scanf("%d%d%d%d", &n, &a, &b, &c);
  v[0] = b; a %= c, b %= c;
  for (int i = 1; i < n; ++i) v[i] = (1ULL * a * v[i - 1] + b) % c;
  return n;
}

int main() {
  freopen("radixsort.in", "r", stdin);
  freopen("radixsort.out", "w", stdout);
  int n = read();
  radix(0, n, 31 - 31 % BITS);
  for (int i = 0; i < n; i += 10) printf("%d ", v[i]);
}