結果

問題 No.957 植林
ユーザー minatominato
提出日時 2024-04-22 03:34:30
言語 C++23
(gcc 12.3.0 + boost 1.83.0)
結果
TLE  
実行時間 -
コード長 18,022 bytes
コンパイル時間 2,733 ms
コンパイル使用メモリ 209,868 KB
実行使用メモリ 29,784 KB
最終ジャッジ日時 2024-10-14 02:26:38
合計ジャッジ時間 11,438 ms
ジャッジサーバーID
（参考情報） 		judge5 / judge3
このコードへのチャレンジ
（要ログイン）

テストケース

テストケース表示
入力 結果 実行時間
実行使用メモリ
testcase_00 AC 2 ms
10,496 KB
testcase_01 AC 1 ms
5,248 KB
testcase_02 AC 2 ms
5,248 KB
testcase_03 AC 71 ms
23,212 KB
testcase_04 AC 68 ms
22,020 KB
testcase_05 AC 67 ms
24,004 KB
testcase_06 AC 82 ms
24,984 KB
testcase_07 AC 73 ms
22,556 KB
testcase_08 AC 48 ms
23,608 KB
testcase_09 AC 49 ms
23,616 KB
testcase_10 AC 51 ms
24,800 KB
testcase_11 AC 48 ms
23,872 KB
testcase_12 AC 48 ms
23,608 KB
testcase_13 AC 40 ms
20,740 KB
testcase_14 AC 54 ms
25,604 KB
testcase_15 AC 45 ms
23,260 KB
testcase_16 AC 38 ms
20,868 KB
testcase_17 AC 40 ms
22,292 KB
testcase_18 TLE -
testcase_19 TLE -
testcase_20 TLE -
testcase_21 -- -
testcase_22 -- -
testcase_23 -- -
testcase_24 -- -
testcase_25 -- -
testcase_26 -- -
testcase_27 -- -
testcase_28 -- -
testcase_29 -- -
testcase_30 -- -
testcase_31 -- -
testcase_32 -- -
testcase_33 -- -
testcase_34 -- -
testcase_35 -- -
testcase_36 -- -
testcase_37 -- -
testcase_38 -- -
testcase_39 -- -
testcase_40 -- -
testcase_41 -- -
testcase_42 -- -
testcase_43 -- -
testcase_44 -- -
testcase_45 -- -
testcase_46 -- -
testcase_47 -- -
権限があれば一括ダウンロードができます

ソースコード

diff # 

#line 1 "library-cpp/other/template.hpp"
// clang-format off
#include <algorithm>
#include <array>
#include <bit>
#include <bitset>
#include <cassert>
#include <chrono>
#include <climits>
#include <cmath>
#include <complex>
#include <deque>
#include <forward_list>
#include <fstream>
#include <functional>
#include <iomanip>
#include <ios>
#include <iostream>
#include <limits>
#include <list>
#include <map>
#include <numeric>
#include <optional>
#include <queue>
#include <random>
#include <ranges>
#include <set>
#include <sstream>
#include <stack>
#include <string>
#include <tuple>
#include <type_traits>
#include <unordered_map>
#include <unordered_set>
#include <utility>
#include <vector>
using namespace std;
using uint = unsigned int;
using ll = long long;
using ull = unsigned long long;
using i128 = __int128_t;
using ld = long double;
using pii = pair<int, int>;
using pll = pair<long long, long long>;
template <class T> using maxheap = priority_queue<T>;
template <class T> using minheap = priority_queue<T, vector<T>, greater<T>>;
template <class T> using vec = vector<T>;
template <class T> using vvec = vector<vector<T>>;
#define OVERLOAD_REP(_1, _2, _3, name, ...) name
#define REP0(n) for (auto minato = decay_t<decltype(n)>{}; minato < (n); ++minato)
#define REP1(i, n) for (auto i = decay_t<decltype(n)>{}; (i) < (n); (i)++)
#define REP2(i, l, r) for (auto i = (l); (i) < (r); (i)++)
#define rep(...) OVERLOAD_REP(__VA_ARGS__, REP2, REP1, REP0)(__VA_ARGS__)
#define OVERLOAD_RREP(_1, _2, _3, name, ...) name
#define RREP1(i, n) for (auto i = (n) - 1; (i) >= decay_t<decltype(n)>{}; (i)--)
#define RREP2(i, l, r) for (auto i = (r) - 1; (i) >= (l); (i)--)
#define rrep(...) OVERLOAD_RREP(__VA_ARGS__, RREP2, RREP1)(__VA_ARGS__)
#define ALL(x) begin(x), end(x)
template <class Container> int SZ(const Container& v) { return int(v.size()); }
template <class T> void UNIQUE(vector<T>& v) { v.erase(unique(v.begin(), v.end()), v.end()); }
template <class T> T MAX(const vector<T>& v) { return *max_element(v.begin(), v.end()); }
template <class T> T MIN(const vector<T>& v) { return *min_element(v.begin(), v.end()); }
template <class T> T SUM(const vector<T>& v) { return accumulate(v.begin(), v.end(), T(0)); }
template <class T> T ABS(T x) { return max(x, -x); }
template <class T> constexpr T POW(T x, ull n) { T ret = 1; while (n > 0) { if (n & 1) ret *= x; x *= x; n >>= 1; } return ret; }
template <class T> constexpr T POW(T x, ull n, T mod) { T ret = 1; while (n > 0) { if (n & 1) ret = ret * x % mod; x = x * x % mod; n >>= 1; } return ret; }
constexpr long long floor_div(long long n, long long d) { return n / d - ((n ^ d) < 0 && n % d); }
template <class T1, class T2> bool chmax(T1& a, T2 b) { if (a < b) { a = b; return true; } return false; }
template <class T1, class T2> bool chmin(T1& a, T2 b) { if (a > b) { a = b; return true; } return false; }
constexpr int topbit(ull x) { return 63 - std::countl_zero(x); }
constexpr int botbit(ull x) { return std::countr_zero(x); }
constexpr int popcnt(ull x) { return std::popcount(x); }
constexpr int kthbit(ull x, int k) { return (x >> k) & 1; }
constexpr int popparity(ull x) { return __builtin_parityll(x); }
constexpr int parity_sign(ull x) { return (x & 1) ? -1 : 1; }
constexpr long long TEN(int x) { return x == 0 ? 1 : TEN(x - 1) * 10; }
template <typename S> void rearrange(const vector<S>& id) { (void)id; }
template <typename S, typename T> void rearrange_exec(const vector<S>& id, vector<T>& v) { vector<T> w(v.size()); for (size_t i = 0; i < id.size(); i++) { w[i] = v[id[i]]; } v.swap(w); }
template <typename S, typename Head, typename... Tail> void rearrange(const vector<S>& id, Head& a, Tail& ...tail) { rearrange_exec(id, a); rearrange(id, tail...); }
istream& operator>>(istream& is, __int128_t& x) {
    x = 0;
    string s;
    is >> s;
    int n = int(s.size()), it = 0;
    if (s[0] == '-') it++;
    for (; it < n; it++) x = (x * 10 + s[it] - '0');
    if (s[0] == '-') x = -x;
    return is;
}
ostream& operator<<(ostream& os, __int128_t x) {
    if (x == 0) return os << 0;
    if (x < 0) os << '-', x = -x;
    deque<int> deq;
    while (x) deq.emplace_front(x % 10), x /= 10;
    for (int e : deq) os << e;
    return os;
}
template <class T> vector<T> &operator++(vector<T>& v) { for (auto& e : v) { e++; } return v;} 
template <class T> vector<T> operator++(vector<T>& v, int) { auto res = v; for (auto& e : v) { e++; } return res; }
template <class T> vector<T> &operator--(vector<T>& v) { for (auto& e : v) { e--; } return v; }
template <class T> vector<T> operator--(vector<T>& v, int) { auto res = v; for (auto& e : v) { e--; } return res; }
template <class T1, class T2> pair<T1, T2> operator-(const pair<T1, T2>& x) { return pair<T1, T2>(-x.first, -x.second); }
template <class T1, class T2> pair<T1, T2> operator-(const pair<T1, T2>& x, const pair<T1, T2>& y) { return pair<T1, T2>(x.first - y.first, x.second - y.second); }
template <class T1, class T2> pair<T1, T2> operator+(const pair<T1, T2>& x, const pair<T1, T2>& y) { return pair<T1, T2>(x.first + y.first, x.second + y.second); }
template <class T1, class T2> pair<T1, T2> operator+=(pair<T1, T2>& l, const pair<T1, T2>& r) { return l = l + r; }
template <class T1, class T2> pair<T1, T2> operator-=(pair<T1, T2>& l, const pair<T1, T2>& r) { return l = l - r; }
constexpr char ln = '\n';
const string YESNO[2] = {"NO", "YES"};
const string YesNo[2] = {"No", "Yes"};
void YES(bool t = true) { cout << YESNO[t] << "\n"; }
void NO(bool t = 1) { YES(!t); }
void Yes(bool t = true) { cout << YesNo[t] << "\n"; }
void No(bool t = 1) { Yes(!t); }
template <class T> void drop(T x) { cout << x << "\n"; exit(0); }
#define INT(...)     \
    int __VA_ARGS__; \
    IN(__VA_ARGS__)
#define LL(...)     \
    ll __VA_ARGS__; \
    IN(__VA_ARGS__)
#define STR(...)        \
    string __VA_ARGS__; \
    IN(__VA_ARGS__)
#define CHR(...)      \
    char __VA_ARGS__; \
    IN(__VA_ARGS__)
#define LDB(...)             \
    long double __VA_ARGS__; \
    IN(__VA_ARGS__)
#define VEC(type, name, size) \
    vector<type> name(size);  \
    IN(name)
#define VEC2(type, name1, name2, size)     \
    vector<type> name1(size), name2(size); \
    for (int i = 0; i < size; i++) IN(name1[i], name2[i])
#define VEC3(type, name1, name2, name3, size)           \
    vector<type> name1(size), name2(size), name3(size); \
    for (int i = 0; i < size; i++) IN(name1[i], name2[i], name3[i])
#define VEC4(type, name1, name2, name3, name4, size)                 \
    vector<type> name1(size), name2(size), name3(size), name4(size); \
    for (int i = 0; i < size; i++) IN(name1[i], name2[i], name3[i], name4[i]);
#define VV(type, name, N, M)                       \
    vector<vector<type>> name(N, vector<type>(M)); \
    IN(name)
template <class T> void scan(T& a) { cin >> a; }
template <class T> void scan(vector<T>& a) { for (auto& i : a) scan(i); }
void IN() {}
template <class Head, class... Tail> void IN(Head& head, Tail&... tail) { scan(head); IN(tail...); }
std::vector<std::vector<int>> read_graph(int N, int M) {
    std::vector<std::vector<int>> g(N);
    for (int i = 0; i < M; i++) {
        int a, b;
        cin >> a >> b;
        a--, b--;
        g[a].emplace_back(b);
        g[b].emplace_back(a);
    }
    return g;
}
std::vector<std::vector<int>> read_tree(int N) { return read_graph(N, N - 1); }
void print() { cout << "\n"; }
template <class T> void print(const vector<T>& v) { for (auto it = v.begin(); it != v.end(); ++it) { if (it != v.begin()) { cout << " "; } cout << *it; } print(); }
template <class T, class... Args> void print(const T& x, const Args& ... args) { cout << x; if (sizeof...(Args)) cout << " "; print(args...); }
#ifdef MINATO_LOCAL
template <class T1, class T2> ostream& operator<<(ostream& os, pair<T1, T2> p);
template <class ...Args> ostream& operator<<(ostream& os, tuple<Args...> t);
template <class T> ostream& operator<<(ostream& os, vector<T> v);
template <class T, size_t N> ostream& operator<<(ostream& os, array<T, N> a);
template <class T, size_t N> enable_if_t<!is_same_v<char, remove_cv_t<T>>, ostream>& operator<<(ostream& os, T(&a)[N]);
template <class Key> ostream& operator<<(ostream& os, set<Key> s);
template <class Key, class T> ostream& operator<<(ostream& os, map<Key, T> mp);
template <class Key> ostream& operator<<(ostream& os, multiset<Key> s);
template <class T> ostream& operator<<(ostream& os, queue<T> q);
template <class T> ostream& operator<<(ostream& os, deque<T> q);
template <class T> ostream& operator<<(ostream& os, priority_queue<T> q);
template <class T> ostream& operator<<(ostream& os, priority_queue<T, vector<T>, greater<T>> q);
template <class T1, class T2> ostream& operator<<(ostream& os, pair<T1, T2> p) { return os << "(" << p.first << ", " << p.second << ")"; }
template <size_t N, class TUPLE> void debug_tuple(ostream& os, TUPLE _) { (void)os; (void)_; }
template <size_t N, class TUPLE, class T, class ...Args> void debug_tuple(ostream &os, TUPLE t) { os << (N == 0 ? "" : ", ") << get<N>(t); debug_tuple<N + 1, TUPLE, Args...>(os, t); }
template <class ...Args> ostream& operator<<(ostream& os, tuple<Args...> t) { os << "("; debug_tuple<0, tuple<Args...>, Args...>(os, t); return os << ")"; }
string debug_delim(int& i) { return i++ == 0 ? "" : ", "; }
#define debug_embrace(x) { int i = 0; os << "{";  { x } return os << "}"; }
template <class T> ostream& operator<<(ostream& os, vector<T> v) { debug_embrace( for (T e : v) { os << debug_delim(i) << e; } ) }
template <class T, size_t N> ostream& operator<<(ostream& os, array<T, N> a) { debug_embrace( for (T e : a) { os << debug_delim(i) << e; } ) }
template <class T, size_t N> enable_if_t<!is_same_v<char, remove_cv_t<T>>, ostream>& operator<<(ostream& os, T(&a)[N]) { debug_embrace( for (T e : a) { os << debug_delim(i) << e; } ) }
template <class Key> ostream& operator<<(ostream& os, set<Key> s) { debug_embrace( for (Key e : s) { os << debug_delim(i) << e; }) }
template <class Key, class T> ostream& operator<<(ostream& os, map<Key, T> mp) { debug_embrace( for (auto e : mp) { os << debug_delim(i) << e; }) }
template <class Key> ostream& operator<<(ostream& os, multiset<Key> s) { debug_embrace( for (Key e : s) { os << debug_delim(i) << e; }) }
template <class T> ostream& operator<<(ostream& os, queue<T> q) { debug_embrace( for (; !q.empty(); q.pop()) { os << debug_delim(i) << q.front(); } ) }
template <class T> ostream& operator<<(ostream& os, deque<T> q) { debug_embrace( for (T e : q) { os << debug_delim(i) << e; } ) }
template <class T> ostream& operator<<(ostream& os, priority_queue<T> q) { debug_embrace( for (; !q.empty(); q.pop()) { os << debug_delim(i) << q.top(); } ) }
template <class T> ostream& operator<<(ostream& os, priority_queue<T, vector<T>, greater<T>> q) { debug_embrace( for (; !q.empty(); q.pop()) { os << debug_delim(i) << q.top(); } ) }
void debug_out() { cerr << endl; }
template <class T, class... Args> void debug_out(const T& x, const Args& ... args) { cerr << " " << x; debug_out(args...); }
#define debug(...) cerr << __LINE__ << " : [" << #__VA_ARGS__ << "] =", debug_out(__VA_ARGS__)
#else
#define debug(...) (void(0))
#endif
///////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
// clang-format on
#line 2 "library-cpp/flow/ProjectSelectionProblem.hpp"

#line 6 "library-cpp/flow/ProjectSelectionProblem.hpp"

#include <atcoder/maxflow>

template <class Cap> struct ProjectSelectionProblem {
    struct Edge {
        int from, to;
        Cap cost;

        Edge(int from, int to, Cap cost) : from(from), to(to), cost(cost) {
        }
    };
    int n, mf_n, s, t;
    Cap base;
    std::vector<Edge> edges;
    atcoder::mf_graph<Cap> g;

    ProjectSelectionProblem() {
    }
    ProjectSelectionProblem(int n_)
        : n(n_), mf_n(n_ + 2), s(n_), t(n_ + 1), base(0) {
    }

    /**
     * 頂点 u が赤色かつ頂点 v が青色だとペナルティを受ける
     */
    void if_u_red_and_v_blue_then_penalty(int u, int v, Cap penalty) {
        assert(0 <= u && u < n && 0 <= v && v < n && u != v && penalty >= 0);
        edges.emplace_back(u, v, penalty);
    }

    /**
     * 頂点 v が赤色になるとペナルティを受ける
     */
    void red_penalty(int v, Cap penalty) {
        assert(0 <= v && v < n && penalty >= 0);
        edges.emplace_back(v, t, penalty);
    }

    /**
     * 頂点 v が青色になるとペナルティを受ける
     */
    void blue_penalty(int v, Cap penalty) {
        assert(0 <= v && v < n && penalty >= 0);
        edges.emplace_back(s, v, penalty);
    }

    /**
     * 頂点 u と v の色が異なるとペナルティを受ける
     */
    void difference_penalty(int u, int v, Cap penalty) {
        assert(0 <= u && u < n && 0 <= v && v < n && u != v && penalty >= 0);
        edges.emplace_back(u, v, penalty);
        edges.emplace_back(v, u, penalty);
    }

    /**
     * 頂点 u が赤色になると頂点 v も赤色になる
     */
    void if_u_red_then_v_red(int u, int v) {
        assert(0 <= u && u < n && 0 <= v && v < n && u != v);
        edges.emplace_back(u, v, std::numeric_limits<Cap>::max());
    }

    /**
     * 頂点 u が青色になると頂点 v も青色になる
     */
    void if_u_blue_then_v_blue(int u, int v) {
        assert(0 <= u && u < n && 0 <= v && v < n && u != v);
        edges.emplace_back(v, u, std::numeric_limits<Cap>::max());
    }

    /**
     * 頂点 v が赤色になる
     */
    void force_red(int v) {
        assert(0 <= v && v < n);
        edges.emplace_back(s, v, std::numeric_limits<Cap>::max());
    }

    /**
     * 頂点 v が青色になる
     */
    void force_blue(int v) {
        assert(0 <= v && v < n);
        edges.emplace_back(v, t, std::numeric_limits<Cap>::max());
    }

    /**
     * 頂点 u と v が同じ色になる
     */
    void force_same(int u, int v) {
        assert(0 <= u && u < n && 0 <= v && v < n && u != v);
        edges.emplace_back(u, v, std::numeric_limits<Cap>::max());
        edges.emplace_back(v, u, std::numeric_limits<Cap>::max());
    }

    /**
     * 頂点 v が赤色になると報酬を得る
     */
    void red_reward(int v, Cap reward) {
        assert(0 <= v && v < n && reward >= 0);
        base += reward;
        edges.emplace_back(s, v, reward);
    }

    /**
     * 頂点 v が青色になると報酬を得る
     */
    void blue_reward(int v, Cap reward) {
        assert(0 <= v && v < n && reward >= 0);
        base += reward;
        edges.emplace_back(v, t, reward);
    }

    /**
     * 頂点 u と v が同じ色になると報酬を得る
     */
    void same_reward(int u, int v, Cap reward) {
        assert(0 <= u && u < n && 0 <= v && v < n && u != v && reward >= 0);
        base += reward;
        edges.emplace_back(u, v, reward);
        edges.emplace_back(v, u, reward);
    }

    /**
     * 赤色の頂点が存在するとペナルティを受ける
     */
    void any_red_peanlty(const std::vector<int>& vs, Cap penalty) {
        int u = mf_n++;
        for (int v : vs) {
            assert(0 <= v && v < n);
            edges.emplace_back(v, u, std::numeric_limits<Cap>::max());
        }
        edges.emplace_back(u, t, penalty);
    }

    /**
     * 青色の頂点が存在するとペナルティを受ける
     */
    void any_blue_penalty(const std::vector<int>& vs, Cap penalty) {
        int u = mf_n++;
        for (int v : vs) {
            assert(0 <= v && v < n);
            edges.emplace_back(u, v, std::numeric_limits<Cap>::max());
        }
        edges.emplace_back(s, u, penalty);
    }

    /**
     * すべての頂点が赤になると報酬を得る
     */
    void all_red_reward(const std::vector<int>& vs, Cap reward) {
        base += reward;
        int u = mf_n++;
        for (int v : vs) {
            assert(0 <= v && v < n);
            edges.emplace_back(u, v, std::numeric_limits<Cap>::max());
        }
        edges.emplace_back(s, u, reward);
    }

    /**
     * すべての頂点が青になると報酬を得る
     */
    void all_blue_reward(const std::vector<int>& vs, Cap reward) {
        base += reward;
        int u = mf_n++;
        for (int v : vs) {
            assert(0 <= v && v < n);
            edges.emplace_back(v, u, std::numeric_limits<Cap>::max());
        }
        edges.emplace_back(u, t, reward);
    }

    /**
     * 報酬の最大値を求める
     */
    Cap solve() {
        g = atcoder::mf_graph<Cap>(mf_n);
        for (const auto& [u, v, cost] : edges) {
            g.add_edge(u, v, cost);
        }
        return base - g.flow(s, t);
    }

    std::vector<bool> red_selection() {
        auto selection = g.min_cut(s);
        selection.resize(n);
        return selection;
    }
};
#line 3 "A.cpp"
void solve() {
    INT(H, W);
    VV(ll, G, H, W);
    VEC(ll, R, H);
    VEC(ll, C, W);

    ProjectSelectionProblem<ll> psp(H + W);
    rep(i, H) {
        ll sum = 0;
        rep(j, W) {
            sum += G[i][j];
        }
        ll val = R[i] - sum;
        if (val > 0) {
            psp.red_reward(i, val);
        } else if (val < 0) {
            psp.red_penalty(i, -val);
        }
    }
    rep(j, W) {
        ll sum = 0;
        rep(i, H) {
            sum += G[i][j];
        }
        ll val = C[j] - sum;
        if (val > 0) {
            psp.red_reward(H + j, val);
        } else if (val < 0) {
            psp.red_penalty(H + j, -val);
        }
    }
    rep(i, H) rep(j, W) {
        psp.all_red_reward({i, H + j}, G[i][j]);
    }
    ll ans = psp.solve();
    print(ans);
}

int main() {
    cin.tie(nullptr);
    ios::sync_with_stdio(false);
    cout << fixed << setprecision(20);
    cerr << fixed << setprecision(7);

    int T = 1;
    // cin >> T;
    for (int test_case = 1; test_case <= T; test_case++) {
        // debug(test_case);
        solve();
    }
}
0