/* #include #include typedef boost::multiprecision::cpp_int mlint; */ #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include #include using namespace std; typedef int64_t i64; typedef double db; typedef vector> mat; typedef pair pii; typedef pair pli; typedef pair pdi; #define esc(x) cout << (x) << endl,exit(0) #define inf(T) (numeric_limits::max() / 2 - 1) #define each(i,v) for(auto i = begin(v); i != end(v); ++i) #define reach(i,v) for(auto i = rbegin(v); i != rend(v); ++i) #define urep(i,s,t) for(int64_t i = (s); i <= (t); ++i) #define drep(i,s,t) for(int64_t i = (s); i >= (t); --i) #define rep(i,n) urep(i,0,(n) - 1) #define rep1(i,n) urep(i,1,(n)) #define rrep(i,n) drep(i,(n) - 1,0) #define all(v) begin(v),end(v) #define rall(v) rbegin(v),rend(v) #define vct vector #define prque priority_queue #define u_map unordered_map #define u_set unordered_set #define l_bnd lower_bound #define u_bnd upper_bound #define rsz resize #define ers erase #define emp emplace #define emf emplace_front #define emb emplace_back #define pof pop_front #define pob pop_back #define mkp make_pair #define mkt make_tuple #define fir first #define sec second struct setupper { constexpr setupper(uint8_t prec) { cin.tie(0); ios::sync_with_stdio(false); cout << fixed << setprecision(prec); #ifdef Local #define debug 1 #define print(x) cout << #x << " = " << (x) << " (L" << __LINE__ << ")" << endl cout << "-- Execution At Local ---\n" << "\n\n\n"; auto print_atexit = []() { time_t end_time = time(NULL); struct tm *ret = localtime(&end_time); cout << "\n----------------\n"; cout << "\nSuccessfully Executed At " << (ret -> tm_hour) << ":" << (ret -> tm_min) << ":" << (ret -> tm_sec) << "\n\n"; }; atexit(print_atexit); #endif } } setupper_obj(10); template ostream& operator << (ostream &s, pair p) { return s << p.fir << " " << p.sec; } template ostream& operator << (ostream &s, vector &v) { for(auto i = v.begin(); i != v.end(); ++i) { if(i != begin(v)) s << " "; s << *i; } return s; } template void hash_combine(size_t &seed, T const &key) { hash hasher; seed ^= hasher(key) + 0x9e3779b9 + (seed << 6) + (seed >> 2); } namespace std { template struct hash> { size_t operator()(pair const &p) const { size_t seed(0); hash_combine(seed,p.first); hash_combine(seed,p.second); return seed; } }; template struct hash> { size_t operator()(tuple const &t) const { size_t seed(0); hash_combine(seed,get<0>(t)); hash_combine(seed,get<1>(t)); hash_combine(seed,get<2>(t)); return seed; } }; } const auto add = [](auto &x, auto y) { x += y; }; const auto mul = [](auto &x, auto y) { x *= y; }; const auto lam_min = [](auto x, auto y) { return min(x,y); }; const auto lam_max = [](auto x, auto y) { return max(x,y); }; const auto chmax = [](auto &m, auto x) { if(m < x){ m = x; return true; } return false; }; const auto chmin = [](auto &m, auto x) { if(m > x){ m = x; return true; } return false; }; bool odd(i64 x) { return x & 1; } bool even(i64 x) { return !odd(x); } constexpr bool bit(i64 n, uint8_t e) { return n >> e & 1; } i64 mask(i64 n, uint8_t e) { return n & (1 << e) - 1; } int ilog(i64 x, uint64_t b = 2) { if(x) return ilog(x / b, b) + 1; return -1; } i64 qceil(i64 x, i64 y) { return x > 0 ? (x - 1) / y + 1 : x / y; } i64 gcd(i64 a, i64 b) { if(!b) return abs(a); return gcd(b, a % b); } i64 lcm(i64 a, i64 b) { if(a | b) return abs(a / gcd(a, b) * b); return 0; } i64 extgcd(i64 a, i64 b, i64 &x, i64 &y) { i64 d = a; if(b) { d = extgcd(b,a % b,y,x); y -= (a / b) * x; } else { x = 1,y = 0; } return d; } i64 bins(i64 ok, i64 ng, const auto &judge) { while(abs(ok - ng) > 1) { i64 mid = ok + ng >> 1; (judge(mid) ? ok : ng) = mid; } return ok; } template T cmprs(T &v) { T tmp = v,ret = v; sort(all(tmp)); tmp.erase(unique(all(tmp)),end(tmp)); each(i,ret) *i = l_bnd(all(tmp),*i) - begin(tmp) + 1; return ret; } constexpr int dx[9] = {1,0,-1,0,1,-1,-1,1,0}; constexpr int dy[9] = {0,1,0,-1,1,1,-1,-1,0}; constexpr i64 mod = 1e9 + 7; constexpr db gold = sqrt(1.25) + 0.5; constexpr db eps = 1e-15; struct nxtset { int ary[4][1<<16]; int sz[4]; constexpr nxtset() : ary(),sz() { rep(i,4) { sz[i]=0; for(int s=(1<<(1<>1); } } } } nx; int pw2[17]; int win[17][17]; i64 dp[5][17][1<<16]; int main() { pw2[0]=1; rep(i,16) pw2[i+1]=pw2[i]<<1; rep(i,16)rep(j,16) cin>>win[i][j]; rep(i,16)rep(j,i) win[i][j]=-win[j][i]; rep(i,16) dp[0][i][pw2[i]]=1; rep(i,4) { for(int g=0;g0?j:k][s|t]+=dp[i][j][s]*dp[i][k][t]; } } } } } rep(i,16) cout<