結果

問題 No.1345 Beautiful BINGO
コンテスト
ユーザー lorent_kyopro
提出日時 2021-01-16 14:28:34
言語 C++17
(gcc 13.3.0 + boost 1.87.0)
結果
WA  
実行時間 -
コード長 3,229 bytes
コンパイル時間 2,174 ms
コンパイル使用メモリ 210,020 KB
最終ジャッジ日時 2025-01-17 22:52:50
ジャッジサーバーID
（参考情報） 		judge2 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 2
other AC * 39 WA * 3 TLE * 19
権限があれば一括ダウンロードができます

ソースコード

diff # 

// #include <atcoder/all>
// using namespace atcoder;
// using mint = modint1000000007;
// using mint = modint998244353;
#include <bits/stdc++.h>
#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<typename T>
inline bool chmax(T& a, const T& b) {
    if (a < b){
        a = b;
        return true;
    }
    return false;
}
template<typename T>
inline bool chmin(T& a, const T& b) {
    if (b < a) {
        a = b;
        return true;
    }
    return false;
}
/**
 * @brief 多次元 vector の作成
 * @author えびちゃん
 */
namespace detail {
    template<typename T, int N>
    auto make_vec(vector<int>& 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<T, N-1>(sizes, x));
        }
    }
}
template<typename T, int N>
auto make_vec(int const(&sizes)[N], T const& x = T()) {
    vector<int> s(N);
    for (int i = 0; i < N; ++i) s[i] = sizes[N-i-1];
    return detail::make_vec<T, N>(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<int>({n, n});
    rep(i, n) rep(j, n) cin >> a[i][j];

    vector<ll> sum(n);
    rep(i, n) rep(j, n) sum[i] += a[i][j];

    auto cost = make_vec<ll>({n, 1 << (n + 2)}, 0);
    rep(i, n) rep(bit, 1 << (n + 2)) {
        vector<bool> 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;
    }

    vector<ll> dp(1 << (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)], now);
    }

    for (int i = 1; i < n; ++i) {
        vector<ll> ndp(1 << (n + 2), LLONG_MAX);
        for (int bit = 0; bit < 1 << (n + 2); ++bit) {
            if (dp[bit] == LLONG_MAX) continue;
            for (int nbit = bit; nbit >= 0; --nbit) {
                nbit &= bit;
                chmin(ndp[nbit], dp[bit] + cost[i][nbit]);
            }
        }
        dp = move(ndp);
    }

    ll ans = LLONG_MAX;
    rep(bit, 1 << (n + 2)) {
        int x = __builtin_popcount(bit);
        int y = max(0, m - x);
        if (y == 0) {
            chmin(ans, dp[bit]);
        } else {
            auto b = a;
            rep(j, n) if (bit >> j & 1) {
                rep(i, n) b[i][j] = 0;
            }
            if (bit >> n & 1) rep(i, n) b[i][i] = 0;
            if (bit >> (n + 1) & 1) rep(i, n) b[i][n-1-i] = 0;
            rep(i, n) for (int j = 1; j < n; ++j) b[i][0] += b[i][j];
            vector<int> c(n);
            rep(i, n) c[i] = b[i][0];
            sort(c.begin(), c.end());
            rep(i, y) dp[bit] += c[i];
            chmin(ans, dp[bit]);
        }
    }
    cout << ans << '\n';
}
0