#ifdef LOCAL
#include "stdafx.h"
#else
#include <bits/stdc++.h>
#define IL inline
#define LL long long
#define eb emplace_back
#define sz(v) static_cast<int>((v).size())
#define L(i, j, k) for (int i = (j); i <= (k); ++i)
#define R(i, j, k) for (int i = (j); i >= (k); --i)
#define FIO(FILE) freopen(FILE".in", "r", stdin), freopen(FILE".out", "w", stdout)
using namespace std;

using vi = vector<int>;
#endif

constexpr int N = 507;
int n, P;

IL constexpr int qpow (int b, int k = P - 2) {
  int r = 1;
  for (; k; k >>= 1, b = (LL)b * b % P) if (k & 1) r = (LL)r * b % P;
  return r;
}

int fc[N], ifc[N], C[N][N];

void init (int Z) {
  fc[0] = 1;
  L (i, 1, Z) {
    fc[i] = (LL)fc[i - 1] * i % P;
  }
  ifc[Z] = qpow(fc[Z]);
  R (i, Z - 1, 0) {
    ifc[i] = (LL)ifc[i + 1] * (i + 1) % P;
  }
  L (i, 0, Z) {
    C[i][0] = 1;
    L (j, 1, i) {
      C[i][j] = (C[i - 1][j] + C[i - 1][j - 1]) % P;
    }
  }
}

int f[N][N], g[N][N], h[N][N];

IL int com (int x) {
  int ret = ifc[n];
  L (i, 1, n) {
    ret = (LL)ret * (x + n - i) % P;
  }
  return ret;
}

IL int sgn (int x) {
  return x & 1 ? -1 : 1;
}

int main () {
  ios::sync_with_stdio(0), cin.tie(0);
  cin >> n >> P;
  init(n);
  L (i, 0, n) {
    L (j, 0, n) {
      f[i][j] = com((LL)i * j % P);
    }
  }
  L (i, 0, n) {
    L (j, 0, n) {
      int cur = 0;
      L (k, 0, j) {
        cur = (cur + (LL)sgn(k) * f[i][j - k] * C[j][k]) % P;
      }
      g[i][j] = (cur % P + P) % P;
    }
  }
  L (j, 0, n) {
    L (i, 0, n) {
      int cur = 0;
      L (k, 0, i) {
        cur = (cur + (LL)sgn(k) * g[i - k][j] * C[i][k]) % P;
      }
      h[i][j] = (cur % P + P) % P;
    }
  }
  L (i, 0, n) {
    L (j, 0, n) {
      f[i][j] = h[i][j];
    }
  }
  L (i, 0, n) {
    L (j, 0, n) {
      int cur = 0;
      L (k, 0, j) {
        cur = (cur + (LL)sgn(j - k) * C[n - k][n - j] * f[i][k]) % P;
      }
      g[i][j] = (cur % P + P) % P;
    }
  }
  L (j, 0, n) {
    L (i, 0, n) {
      int cur = 0;
      L (k, 0, i) {
        cur = (cur + (LL)sgn(i - k) * C[n - k][n - i] * g[k][j]) % P;
      }
      h[i][j] = (cur % P + P) % P;
    }
  }
  L (i, 0, n - 1) {
    L (j, 0, n - 1) {
      cout << h[n - i][n - j] << ' ';
    }
    cout << '\n';
  }
}
// I love WHQ!