#include <iostream>
#include <string>
#include <algorithm>
#include <vector>
#include <array>
#include <queue>
#include <deque>
#include <map>
#include <set>
#include <sstream>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <list>
#include <numeric>
#include <stack>
#include <iomanip>
#include <random>
#include <complex>
#include <functional>
#include <tuple>

using namespace std;

#define Rep(i,a,b) for(ll i = a; i < b; ++i)
#define rep(i,b) Rep(i,0,b)
#define allof(a) (a).begin(), (a).end()

#define Yes(q) ((q) ? "Yes" : "No")
#define YES(q) ((q) ? "YES" : "NO")
#define Possible(q) ((q) ? "Possible" : "Impossible")
#define POSSIBLE(q) ((q) ? "POSSIBLE" : "IMPOSSIBLE")

using ll = long long;

constexpr int inf = 1e9 + 7;
constexpr ll infll = 1ll << 60ll;
constexpr ll mod = 1e9 + 7;
// 0~3までは右左下上 4~7までは斜め
constexpr int dx[] = { 1, 0, -1, 0, 1, 1, -1, -1 };
constexpr int dy[] = { 0, -1, 0, 1, 1, -1, -1, 1 };

namespace {
    template<typename T> void chmax(T& a, T b) { a = std::max(a, b); }
    template<typename T> void chmin(T& a, T b) { a = std::min(a, b); }
    template<typename T> void chadd(T& a, T b) { a = a + b; }

    void Cout(long long x, const char* end = "\n") { std::cout << x << end; }
    template <typename T> void Cout(const T& x, const char* end = "\n") { std::cout << x << end; }
    template <typename T> void Cout(const std::vector<T>& x, const char* sep = " ", const char* end = "\n") { for (std::size_t i = 0, sz = x.size(); i < sz; i++) { std::cout << x[i] << (i == sz - 1 ? end : sep); } }

    // 標準入出力
    struct inp {
        std::size_t sz;
        inp(std::size_t _sz = 1) : sz(_sz) {}
        template <typename T> operator T () const { T a; std::cin >> a; return a; }
        template <typename T> operator std::vector<T>() const { vector<T> a(sz); for (std::size_t i = 0; i < sz; i++) std::cin >> a[i]; return a; }
        template <typename T, typename U> operator std::pair<T, U>() const { T f; U s; std::cin >> f >> s; return std::pair<T, U>(f, s); }
    };

    inp inp1; // input one

    // Bit演算を楽にするライブラリ ２進数
    // @author tada
    class Bit {
    public:
        // コンストラクタでサイズ数
        Bit(int length) {
            length_ = length;
            bit_ = -1; // while(bit.next()) をされるのを想定
        }

        // 一つ先の値に移動する もう次に進めないならばfalseを返す
        bool Next() {
            if (++bit_ >= (1 << length_)) return false;
            return true;
        }

        // Bit操作 指定したインデックスを変更する
        void Change(int i, int value = 1) {
            if (value == 1) bit_ |= (1 << i);
            else bit_ &= ~(1 << i);
        }

        // 指定したインデックスが1かどうか
        bool IsSelected(int i) {
            return bit_ & (1 << i);
        }

        // 1の個数を返す 引数で範囲 [l, r) 準備中
        int GetSelectedNum(){ //int l, int r) {
            return __builtin_popcount(bit_);
        }

        // ビット情報を初期状態に戻す
        void Reset() {
            bit_ = -1;
        }

        // ビット列をint型で返す
        int GetBit() {
            return bit_;
        }

        // ビット長を返す
        int GetLength() {
            return length_;
        }

    private:
        int bit_;
        int length_;
    };
}

int main() {

    int n = inp1;
    vector<vector<int>> a(n, vector<int>(n));
    rep(i, n)rep(j, n) cin >> a[i][j];

    // bitDP

    // dp[i][s] := i番目までの男を見たときに，女性集合sとマッチするパターン数
    // 女性集合sとは 0110 の場合は，2,3番目の女性がマッチングしている
    vector<vector<ll>> dp(n + 1, vector<ll>((1 << n) + 1, 0LL));

    dp[0][0] = 1LL;

    Bit bit(n);

    rep(i, n) {
        bit.Reset();
        while (bit.Next()) {
            int b = bit.GetBit();
            // ペアが決まっている男性の数と女性の数が一致しないと意味がない
            if (i + 1 != bit.GetSelectedNum()) continue;
            // ひとり少ないやつから遷移
            rep(j, n) {
                dp[i + 1][b] += dp[i][b & ~(1 << j)] * a[i][j];
                dp[i + 1][b] %= mod;
            }
        }
    }

    //Cout(dp[n]);

    Cout(dp[n][(1 << n) - 1]);

    return 0;
}


Measured memory_bytes: 377216