#include <bits/stdc++.h>

using namespace std;

long long int A[16][16];
long long int dp[5][1<<16];
vector <int> state[17];
int main(void)
{
    cin.tie(0);
    ios::sync_with_stdio(false);

    int n = 16;

    for(int i=0;i<n;i++)
    {
        for(int j=0;j<n;j++)
        {
            cin >> A[i][j];
        }
    }

    for(int i=0;i<n;i++)
    {
        for(int j=i;j<n;j++)
        {
            A[j][i] = -A[i][j];
        }
    }

    for(int i=0;i<(1<<16);i++)
    {
        int bcnt = __builtin_popcount(i);
        state[bcnt].push_back(i);
    }

    for(int X=0;X<16;X++)
    {
        memset(dp,0,sizeof(dp));
        dp[0][(1<<X)] = 1;
        for(int i=0;i<4;i++)
        {
            int x = (1<<i);
            for(auto s : state[x])
            {
                if(dp[i][s]==0) continue;
                for(auto s2 : state[2*x])
                {
                    if((s&s2)==s)
                    {
                        int s3 = s^s2;
                        vector <int> a,b;
                        for(int j=0;j<16;j++)
                        {
                            if((s&(1<<j))) a.push_back(j);
                            if((s3&(1<<j))) b.push_back(j);
                        }
                        int N = a.size();
                        int M = b.size();
                        vector<vector<long long int>> dp2(N+1,vector<long long int>((1<<M), 0));
                        dp2[0][0] = 1;
                        for(int j=0;j<N;j++)
                        {
                            for(int k=0;k<(1<<M);k++)
                            {
                                if(dp2[j][k]==0) continue;
                                for(int l=0;l<M;l++)
                                {
                                    if((k&(1<<l))) continue;
                                    if(A[a[j]][b[l]]==1)
                                    {
                                        dp2[j+1][k|(1<<l)] += (2*dp2[j][k]);
                                    }
                                }
                            }
                        }
                        dp[i+1][s2] += (dp2[N][(1<<M)-1]*dp[i][s]);
                    }
                }
            }
        }

        cout << dp[4][(1<<16)-1] << '\n';
    }


    return 0;
}