// #include // using namespace atcoder; // using mint = modint1000000007; // using mint = modint998244353; #include #define rep(i, n) for (int i = 0; i < (int)n; ++i) #define rrep(i, n) for (int i = (int)n-1; i >= 0; --i) using namespace std; using ll = long long; template inline bool chmax(T& a, const T& b) { if (a < b){ a = b; return true; } return false; } template inline bool chmin(T& a, const T& b) { if (b < a) { a = b; return true; } return false; } /** * @brief 多次元 vector の作成 * @author えびちゃん */ namespace detail { template auto make_vec(vector& sizes, T const& x) { if constexpr (N == 1) { return vector(sizes[0], x); } else { int size = sizes[N-1]; sizes.pop_back(); return vector(size, make_vec(sizes, x)); } } } template auto make_vec(int const(&sizes)[N], T const& x = T()) { vector s(N); for (int i = 0; i < N; ++i) s[i] = sizes[N-i-1]; return detail::make_vec(s, x); } __attribute__((constructor)) void fast_io() { ios::sync_with_stdio(false); cin.tie(nullptr); } int main() { int n, m; cin >> n >> m; auto a = make_vec({n, n}); rep(i, n) rep(j, n) cin >> a[i][j]; vector sum(n); rep(i, n) rep(j, n) sum[i] += a[i][j]; auto cost = make_vec({n, 1 << (n + 2)}, 0); rep(i, n) rep(bit, 1 << (n + 2)) { vector used(n); rep(j, n) if (bit >> j & 1) cost[i][bit] += a[i][j], used[j] = true; if ((bit >> n & 1) && !used[i]) cost[i][bit] += a[i][i], used[i] = true; if ((bit >> (n + 1) & 1) && !used[n - 1 - i]) cost[i][bit] += a[i][n - 1 - i], used[n - 1 - i] = true; } auto dp = make_vec({1 << (n + 2), n + 2}, LLONG_MAX); rep(bit, 1 << n) { ll now = 0; rep(j, n) { if (bit >> j & 1) now += a[0][j]; } chmin(dp[bit | ((bit & 1) << n) | (bit >> (n - 1) & 1) << (n + 1)][bit == (1 << n) - 1], now); } for (int i = 1; i < n; ++i) { auto ndp = make_vec({1 << (n + 2), n + 2}, LLONG_MAX); for (int bit = 0; bit < 1 << (n + 2); ++bit) { for (int x = 0; x <= n; ++x) { if (dp[bit][x] == LLONG_MAX) continue; for (int nbit = bit; nbit >= 0; --nbit) { nbit &= bit; int nx = x + (sum[i] == cost[i][nbit]); chmin(ndp[nbit][nx], dp[bit][x] + cost[i][nbit]); } chmin(ndp[bit][x + 1], dp[bit][x] + sum[i]); } } dp = move(ndp); } ll ans = LLONG_MAX; rep(bit, 1 << (n + 2)) rep(x, n + 1) { int y = x + __builtin_popcount(bit); if (y >= m) chmin(ans, dp[bit][x]); } cout << ans << '\n'; }