#include <bits/stdc++.h>
using namespace std;
const int N = 5e2 + 5;
char s[N][N];
long long dp[2][N][N];
int x[N];
long long get(int x) { return 1ll * x * (x - 1) / 2; }
long long fac[N];
void init(int mod) {
  for (int i = 0; i < N; i++) fac[i] = get(i) % mod;
}
int main() {
  int n, m, mod;
  scanf("%d %d %d", &n, &m, &mod);
  for (int i = 1; i <= m; i++) {
    scanf(" %s", s[i] + 1);
    for (int j = 1; j <= n; j++) x[j] += s[i][j] - '0';
  }
  init(mod);
  int fi = 0, se = 0;
  for (int i = 1; i <= n; i++) {
    if (x[i] == 1)
      fi++;
    else if (x[i] == 2)
      se++;
  }
  dp[0][fi][se] = 1;
  int f = 0;
  for (int i = m + 1; i <= n; i++) {
    f ^= 1;
    memset(dp[f], 0, sizeof(dp[f]));
    for (int j = 0; j <= n; j++) {
      for (int k = 0, z; k + j <= n; k++) {
        if (!dp[f ^ 1][j][k]) continue;
        z = n - (j + k);
        if (z >= 2) dp[f][j + 2][k] += dp[f ^ 1][j][k] * fac[z] % mod;
        if (j >= 1 && z >= 1)
          dp[f][j][k + 1] += dp[f ^ 1][j][k] * j % mod * z % mod;
        if (j >= 2) dp[f][j - 2][k + 2] += dp[f ^ 1][j][k] * fac[j] % mod;
      }
    }
  }
  printf("%lld\n", dp[f][0][n] % mod);
  return 0;
}
