結果

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

ソースコード

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;
    }

    auto dp = make_vec<ll>({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<ll>({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';
}
0