結果

問題 No.957 植林
ユーザー ooaiu
提出日時 2025-07-25 11:37:16
言語 C++23
(gcc 13.3.0 + boost 1.87.0)
結果
TLE  
実行時間 -
コード長 2,634 bytes
コンパイル時間 3,792 ms
コンパイル使用メモリ 291,364 KB
実行使用メモリ 30,744 KB
最終ジャッジ日時 2025-07-25 11:37:26
合計ジャッジ時間 10,297 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge1
このコードへのチャレンジ
（要ログイン）
ファイルパターン 結果
sample AC * 3
other AC * 15 TLE * 1 -- * 29
権限があれば一括ダウンロードができます

ソースコード

diff # 

#include <bits/stdc++.h>
using namespace std;
template <class F> class y_combinator {
    F f;

   public:
    y_combinator(F&& f) : f(std::forward<F>(f)) {}
    template <class... Args> auto operator()(Args&&... args) const { return f(*this, std::forward<Args>(args)...); }
};
using ll = long long;
using ld = long double;
template <class T, class U = std::less<T>> using prique = std::priority_queue<T, std::vector<T>, U>;
template <class T> T floor(T a, T b) noexcept { return a / b - (a % b && (a ^ b) < 0); }
template <class T> T ceil(T a, T b) noexcept { return floor(a + b - 1, b); }
template <class T> bool chmin(T& x, const T& y) noexcept { return (x > y ? x = y, true : false); }
template <class T> bool chmax(T& x, const T& y) noexcept { return (x < y ? x = y, true : false); }
#define overload4(a, b, c, d, e, ...) e
#define rep1(a) for (long long _i = 0; _i < (a); _i++)
#define rep2(i, a) for (long long i = 0; i < (a); i++)
#define rep3(i, a, b) for (long long i = (a); i < (b); i++)
#define rep4(i, a, b, c) for (long long i = (a); i < (b); i += (c))
#define rep(...) overload4(__VA_ARGS__, rep4, rep3, rep2, rep1)(__VA_ARGS__)
#define rrep(i, a, b, c) for (long long i = (a); i > (b); i += (c))
#define all(x) std::begin(x), std::end(x)
#define rall(x) std::rbegin(x), std::rend(x)
#define pb push_back
#ifndef LOCAL
#define debug(...)
#endif

#include <atcoder/maxflow>
void run_case() {
    int H, W;
    cin >> H >> W;
    vector<vector<int>> G(H, vector<int>(W));
    rep(i, H) rep(j, W) cin >> G[i][j];
    vector<int> R(H), C(W);
    rep(i, H) cin >> R[i];
    rep(i, W) cin >> C[i];
    const int N = H * W;
    const int SRC = N;
    const int SNK = N + 1;
    atcoder::mf_graph<ll> MFG(N + 2 + H + W);
    rep(i, H) rep(j, W) {
        MFG.add_edge(i * W + j, SNK, G[i][j]);
        MFG.add_edge(SRC, i * W + j, 0);
    }
    const ll INF = 1e18;
    ll tot = 0;
    {
        const int OFFSET = N + 2;
        rep(i, H) {
            tot += R[i];
            MFG.add_edge(SRC, OFFSET + i, R[i]);
            rep(j, W) {
                MFG.add_edge(OFFSET + i, i * W + j, INF);
            }
        }
    }
    {
        const int OFFSET = N + 2 + H;
        rep(j, W) {
            tot += C[j];
            MFG.add_edge(SRC, OFFSET + j, C[j]);
            rep(i, H) {
                MFG.add_edge(OFFSET + j, i * W + j, INF);
            }
        }
    }
    ll ans = tot - MFG.flow(SRC, SNK);
    cout << ans << "\n";
}
int main() {
    std::ios_base::sync_with_stdio(false);
    std::cin.tie(nullptr);
    std::fixed(std::cout).precision(16);
    int T = 1;
    while (T--) run_case();
    return 0;
}
0