#include <cstdio>
#include <cstring>
#include <iostream>
#include <string>
#include <cmath>
#include <bitset>
#include <vector>
#include <map>
#include <set>
#include <queue>
#include <deque>
#include <algorithm>
#include <complex>
#include <unordered_map>
#include <random>
using namespace std;
typedef long long int ll;
typedef pair<int, int> P;
int main()
{
	int a[16][16];
	for(int i=0; i<16; i++){
		for(int j=0; j<16; j++){
			cin>>a[i][j];
		}
	}
	ll dp[16][1<<16]={};
	for(int i=0; i<16; i++){
		for(int j=0; j<i; j++){
			if(a[j][i]>0) dp[j][(1<<i)+(1<<j)]=2;
			else dp[i][(1<<i)+(1<<j)]=2;
		}
	}
	for(int k=2; k<=4; k++){
		for(int i=0; i<(1<<16); i++){
			int c=0;
			for(int j=0; j<16; j++){
				if(i&(1<<j)) c++;
			}
			if(c!=(1<<k)) continue;
			for(int j=i; j>0; j=(j-1)&i){
				int c=0;
				for(int l=0; l<16; l++){
					if(j&(1<<l)) c++;
				}
				if(c!=(1<<(k-1))) continue;
				vector<int> v1, v2;
				for(int l=0; l<16; l++){
					if(j&(1<<l)) v1.push_back(l);
					else v2.push_back(l);
				}
				for(auto x:v1){
					for(auto y:v2){
						if(x<y){
							if(a[x][y]>0) dp[x][i]+=(dp[x][j]*dp[y][i-j]);
							else dp[y][i]+=(dp[x][j]*dp[y][i-j]);
						}else{
							if(a[y][x]>0) dp[y][i]+=(dp[x][j]*dp[y][i-j]);
							else dp[x][i]+=(dp[x][j]*dp[y][i-j]);
						}
					}
				}
			}
		}
	}
	for(int i=0; i<16; i++){
		cout<<dp[i][(1<<16)-1]<<endl;
	}
	return 0;
}