// Slower C++ solution:
#include<bits/stdc++.h>

#define FIO ios_base::sync_with_stdio(false);cin.tie(NULL);cout.tie(NULL);

#define PI 3.141592653589793238462643383279502884L

#define make_unique(vec) vec.resize(distance(vec.begin(),unique(vec.begin(),vec.end())))

#define Sort(vec) sort(vec.begin(),vec.end())

#define RSort(vec) sort(vec.rbegin(),vec.rend())

#define endl "\n"

#define show(i); cout<<i<<" ";

#define tr1(i); cout<<i<<endl;

#define tr2(i,j); cout<<i<<" "<<j<<endl;

#define tr3(i,j,k); cout<<i<<" "<<j<<" "<<k<<endl;

#define vvi vector<vector<int> > 

#define vvl vector<vector<lli> >

#define all(st) st.begin(),st.end()

#define vpll vector<pair<lli,lli> >

#define vpii vector<pair<int,int> >

#define pi pair<int,int>

#define pl pair<lli,lli>

#define al vector<list<int> >

#define vs vector<string>

#define vi vector<int>

#define vb vector<bool>

#define vl vector<lli>

#define vc vector<char>

#define rep(i,l,r) for(int i=l;i<r;i++)

#define repit(st) for(auto it=st.begin();it!=st.end();it++)

#define ulli unsigned long long int

#define eb emplace_back

#define pb push_back

#define lli long long int

#define inf INT_MAX

#define minf INT_MIN

#define f first

#define se second

#define mod 1000000007

#define mx 500001

using namespace std;

void add_me(lli &a,lli b){

    a+=b;

    if(a>=mod) a-=mod;

}   

void sub_me(lli &a,lli b){

    a-=b;

    if(a<0) a+=mod;

}

bool sortbysec(const pair<int,int> &a,const pair<int,int> &b){

	if(a.se==b.se) return a.f<b.f;

    return (a.se < b.se); 

}

struct disjointset{

    int parent[mx];

   	disjointset(){

        for(int i=1;i<mx;i++) parent[i]=i;

    }

    int find(int x){

      	int par=parent[x];

        if(par==x) return x;

        return parent[x]=find(parent[x]);

    }

    void union_(int x,int y){

   	    int par1=find(x);

     	int par2=find(y);

       	if(par1==par2) return;

       	parent[par2]=par1;

    }

};

lli mi(lli n,lli m){

	lli pw=n;

	lli ex=m-2;

	lli ans=1;

	while(ex){

		if(ex&1) ans = ans*pw%m;

		pw = pw*pw%m;

		ex>>=1;

	}

	return ans%m;

}

// const int pm = 1000000;

// int spf[pm+1];

// bool is_prime[pm+1];

// int ETF[pm+1];

// void prime_solver(){

// 	memset(is_prime,true,sizeof(is_prime));

// 	for(int i=1;i<=pm;i++) spf[i]=i;

// 	for(int i=2;i*i<=pm;i++){

// 		if(is_prime[i]){

// 			for(int j=i*i;j<=pm;j+=i){

// 				is_prime[j]=false;

// 				if(spf[j]==j) spf[j]=i;

// 			}

// 		}

// 	}

// 	for(int i=1;i<=pm;i++) ETF[i]=i;

// 	for(int i=2;i*i<=pm;i++){

// 		if(ETF[i]==i){

// 			for(int j=i;j<=pm;j+=i){

// 				ETF[j]/=i;

// 				ETF[j]*=(i-1);

// 			}

// 		}

// 	}

// }

int target;

int n;

vector<vi> vec(21,vi(21));

int dp[1<<21];

int solve(int p,int mw){

	if(mw==target) return 1;

	if(dp[mw]!=-1) return dp[mw];

	int ans=0;

	int i,j;

	for(int j=0;j<n;j++){

		if((mw&(1<<j))==0 && vec[p][j]){

			ans+=solve(p+1,(mw|(1<<j)));

			ans%=mod;

		}

	}

	return dp[mw]=ans;

}

int main(){

	FIO;

	cin>>n;

	memset(dp,-1,sizeof(dp));

	rep(i,0,n){

		rep(j,0,n) cin>>vec[i][j];

	}

	target=((1<<n)-1);

	int ans=solve(0,0);

	cout<<ans<<endl;

}

// Faster C++ solution for the same problem (measured speedup: 9.708003063784234x):
#include<bits/stdc++.h>

#include <ext/pb_ds/assoc_container.hpp> 

#include <ext/pb_ds/tree_policy.hpp> 

#define FIO ios_base::sync_with_stdio(false);cin.tie(NULL);cout.tie(NULL);

#define PI 3.141592653589793238462643383279502884L

#define make_unique(vec) vec.resize(distance(vec.begin(),unique(vec.begin(),vec.end())))

#define Sort(vec) sort(vec.begin(),vec.end())

#define RSort(vec) sort(vec.rbegin(),vec.rend())

#define endl "\n"

#define tr1(i) cout<<i<<endl;

#define tr2(i,j) cout<<i<<" "<<j<<endl;

#define tr3(i,j,k) cout<<i<<" "<<j<<" "<<k<<endl;

#define vvi vector<vector<int> > 

#define vvl vector<vector<ll> >

#define all(st) st.begin(),st.end()

#define vpl vector<pair<ll,ll> >

#define vpi vector<pair<int,int> >

#define pi pair<int,int>

#define pl pair<ll,ll>

#define al vector<list<int> >

#define vs vector<string>

#define vb vector<bool>

#define vi vector<int>

#define vl vector<ll>

#define vc vector<char>

#define rep(i,l,r) for(int i=l;i<r;i++)

#define repit(st) for(auto it=st.begin();it!=st.end();it++)

#define cs(ans); cout<<setprecision(6)<<fixed<<ans<<endl;

#define ull unsigned long long int

#define eb emplace_back

#define pb push_back

#define ll long long int

#define minf -(1e16)

//#define inf 1e16

#define f first

#define se second

const ll mod=1e9+7;

using namespace std;

using namespace __gnu_pbds; 

ll mi(ll n,ll m){ll pw=n%m;ll ex=m-2;ll ans=1;while(ex){if(ex&1) ans = ans*pw%m;pw = pw*pw%m;ex>>=1;}return ans%mod;}

ll pw(ll a,ll n){ll pw=a;ll ex=n;ll ans=1;while(ex){if(ex&1) ans = ans*pw;pw = pw*pw;ex>>=1;}return ans;}

ll pwm(ll a,ll n){ll pw=a%mod;ll ex=n;ll ans=1;while(ex){if(ex&1) ans = ans*pw%mod;pw = pw*pw%mod;ex>>=1;}return ans%mod;}

template<typename T>istream &operator>>(istream &is, vector<T> &v) { for (T &t : v) is >> t; return is; }

template<typename T>ostream &operator<<(ostream &os, const vector<T> &v) { for (const T &t : v) {os << t <<" ";} os << '\n'; return os; }

template<typename T> using ordered_set=tree<T,null_type,less<T>,rb_tree_tag,tree_order_statistics_node_update>;

//st.order_of_key(num)->number of elements < num

//st.find_by_order(i)->iterator of ith element in sorted order starting from zero

ll gcd(ll a,ll b){return b==0?a:gcd(b,a%b);}

ll lcm(ll a,ll b){return a/gcd(a,b)*b;}

ll max(ll a,ll b){return a>b?a:b;}

ll min(ll a,ll b){return a<b?a:b;}

ll Abs(ll x){return x>0?x:(-x);}

ll Ceil(ll a,ll b){return a/b+(a%b!=0);}

#define bs binary_search

#define lb lower_bound

#define ub upper_bound

#define mkp make_pair

//#define int long long int

//Binary_Search->dp->greedy

//auto check = [&](ll x){};

int dx[]={1,0,-1,0};

int dy[]={0,1,0,-1};

char move[]={'D','R','U','L'};

int X[]={2,1,-1,-2,-2,-1,1,2}; 

int Y[]={1,2,2,1,-1,-2,-2,-1}; 

struct comp{

	bool operator()(const pair<int,int>&a,const pair<int,int>&b) const {

		if(a.first>b.first) return true;

		if(a.first<b.first) return false; 

		if(a.second<b.second) return true;

		return false;

	}

};	

void solve() {

	//ifstream fin("input.txt");

	//ofstream fout("output.txt");

	int n;

	cin>>n;

	vvi m(n,vi(n));

	rep(i,0,n){

		rep(j,0,n) cin>>m[i][j];

	}

	vector<ll> dp(1<<n);

	dp[0]=1;

	for(int mask=0;mask<(1<<n)-1;mask++){

		int matched=__builtin_popcount(mask);

		for(int j=0;j<n;j++){

			if(m[matched][j] && ((mask&(1<<j))==0)){

				dp[mask|(1<<j)]+=dp[mask];

				dp[mask|(1<<j)]%=mod;

			}

		}

	}

	cout<<dp[(1<<n)-1]<<endl;

	

	

	

	

}	



int32_t main(){

int t=1;

//scanf("%d",&t);

while(t--){

	solve();     

}

return 0;

}