#include<iostream>
#include<vector>
using namespace std;
 
class ModInt{
private:
  long long val;
  static const int default_mod = 1e9 + 7;
  static int mod;
  static vector<ModInt> facts;
public:
  ModInt(long long val)
  {
    if(mod == 0){
      mod = default_mod;
    }
    this->val = ((val % mod) + mod) % mod;
  }
  ModInt()
    : ModInt(0)
  {
  }
  static void setMod(int mod){
    if(ModInt::mod != 0){
      throw "The modulo is already set.";
    }
    ModInt::mod = mod;
  }

  const ModInt operator* (const ModInt& another) const{
    return ModInt(val * another.val);
  }
  ModInt operator+ (const ModInt& another){
    return ModInt(val + another.val);
  }
  ModInt operator- (const ModInt & another){
    return ModInt(val - another.val + mod);
  }
  ModInt operator/(const ModInt& rhs){
    return (*this) * rhs.inv();
  }
  ModInt& operator*=(const ModInt& rhs){
    (*this) = (*this) * rhs;
    return *this;
  }
  ModInt& operator+=(const ModInt& rhs){
    (*this) = (*this) + rhs;
    return *this;
  }
  ModInt& operator-=(const ModInt& rhs){
    (*this) = (*this) - rhs;
    return *this;
  }
  ModInt& operator/=(const ModInt& rhs){
    (*this) = (*this) / rhs;
    return *this;
  }
  ModInt pow(long long exponent) const{
    exponent %= mod - 1;
    if(exponent == 0){
      return ModInt(1);
    }else{
      if(exponent % 2){
        return ((*this) * (*this)).pow(exponent / 2) * (*this);
      }else{
        return ((*this) * (*this)).pow(exponent / 2);
      }
    }
  }
  ModInt inv() const{
    return this->pow(mod - 2);
  }
  static ModInt fact(long long n){
    if(n >= mod){
      return 0;
    }
    if(facts.size() < static_cast<size_t>(n + 1)){
      if(facts.size() == 0){
        facts.push_back(1);
      }
      for(int i = facts.size(); i <= n; i++){
        facts.push_back(facts[i - 1] * i);
      }
    }
    return facts[n];
  }
  static ModInt comb(long long n, long long r){
    if(r > n){
      return 0;
    }
    return ModInt::fact(n) * ModInt::fact(r).inv() * ModInt::fact(n - r).inv();
  }
 
  friend ostream& operator<< (ostream& os, const ModInt& modint){
    os << modint.val;
    return os;
  }
  friend istream& operator>>(istream& is, ModInt& modint){
    long long val;
    is >> val;
    modint.val = val % modint.mod;
    return is;
  }
  friend ModInt operator*(const long long lhs, const ModInt& rhs){
    return ModInt(lhs) * rhs;
  }
};
int ModInt::mod = 0;
vector<ModInt> ModInt::facts;

ModInt calc_dp(int n, vector<ModInt>& dp, vector<bool>& checked, int N, vector<vector<int>>& a){
  if(n == 0){
    checked[n] = true;
    return dp[n] = 1;
  }
  if(checked[n]){
    return dp[n];
  }

  int n_bit_cnt = __builtin_popcount(n);
  
  ModInt ret = 0;
  for(int i = 0; i < N; i++){
    if((n >> i) & 1){
      ret += calc_dp(n & (~(1 << i)), dp, checked, N, a) * a[n_bit_cnt - 1][i];
    }
  }
  checked[n] = true;
  return dp[n] = ret;
}

int main(){
  int N;
  cin >> N;

  vector<vector<int>> a(N, vector<int>(N));
  for(auto& aa: a){
    for(auto& aaa: aa){
      cin >> aaa;
    }
  }

  // dp[00110] := #{(pair of woman 2, pair of woman 3)| paif of woman * in man {0, 1}}
  // dp[00110] = dp[00010] * a[1][2] + dp[00100] * a[1][3]
  vector<ModInt> dp(1 << N);
  vector<bool> checked(1 << N, false);
  cout << calc_dp((1 << N) - 1, dp, checked, N, a) << endl;
}


Measured memory_bytes: 16896