#include using ll = long long; constexpr int N = 16; constexpr int MASK = 1 << N; ll memo[N][MASK]; int win[N][N]; int main() { for (int i = 0; i < N; i++) { std::fill(memo[i], memo[i] + MASK, -1); } auto next_combination = [](const int mask) { const int x = mask & (-mask), y = mask + x; return (((mask & ~y) / x) >> 1) | y; }; for (int i = 0; i < N; i++) { for (int j = 0, c; j < N; j++) { std::cin >> c; if (j >= i) { win[i][j] = (c == 1), win[j][i] = (c == -1); } } } auto dp = [&](auto&& self, const int p, const int mask, const int b) -> ll { if (b == 1) { return mask == (1 << p) ? 1 : 0; } if (memo[p][mask] != -1) { return memo[p][mask]; } std::vector bits; for (int i = 0; i < N; i++) { if (mask & (1 << i)) { bits.push_back(i); } } const int half = b / 2; int subi = (1 << half) - 1; ll ans = 0; for (; subi < (1 << b); subi = next_combination(subi)) { int sub = 0; for (int i = 0; i < b; i++) { if (subi & (1 << i)) { sub |= (1 << bits[i]); } } int other = mask - sub; if (other & (1 << p)) { std::swap(sub, other); } for (const int q : bits) { if ((sub & (1 << q)) or (win[q][p])) { continue; } ans += self(self, p, sub, b / 2) * self(self, q, other, b / 2); } } return memo[p][mask] = ans; }; for (int i = 0; i < N; i++) { std::cout << dp(dp, i, MASK - 1, N) << std::endl; } return 0; }