QOJ.ac

QOJ

ID题目提交者结果用时内存语言文件大小提交时间测评时间
#282372#3249. 分组作业kilo_tobo_tarjen#TL 0ms0kbC++209.3kb2023-12-11 21:21:582023-12-11 21:21:58

Judging History

你现在查看的是最新测评结果

  • [2023-12-11 21:21:58]
  • 评测
  • 测评结果:TL
  • 用时:0ms
  • 内存:0kb
  • [2023-12-11 21:21:58]
  • 提交

answer

#include <bits/stdc++.h>
// #include <bits/extc++.h>
// using namespace __gnu_cxx;
// using namespace __gnu_pbds;
using i64 = long long;
using namespace std;
const int N = 1e4 + 5;
// const int B = 1250;
// const int M = 105;
// const int base = 17131;
// const int mod = 998244353;
// const int mod = 1e9 + 7;
// const long double pi = acosl(-1);

namespace maxflow
{

    typedef long long flow_t;
    const flow_t inf_flow = 1e9;
    const int inf_dep = 1e9;

    struct FlowEdge
    {
        int from, to;
        flow_t cap, low = 0, flow = 0;
    };

    int num_node(const std::vector<FlowEdge> &edges)
    {
        int n = 0;
        for (const auto &e : edges)
            n = std::max({n, e.from, e.to});
        return n;
    }

    flow_t get_flow(const std::vector<FlowEdge> &edges, int s)
    {
        flow_t flow = 0;
        for (const auto &e : edges)
        {
            if (e.from == s)
                flow += e.flow;
        }
        return flow;
    }

    struct MaxFlow
    {
        struct Edge
        {
            int from, to;
            flow_t cap;
        };
        int n;
        std::vector<std::vector<int>> eid;
        std::vector<Edge> edge;
        void build(const std::vector<FlowEdge> &edges)
        {
            n = num_node(edges);
            eid.assign(n + 1, {});
            edge.clear();
            int num_edges = 0;
            for (const auto &e : edges)
            {
                eid[e.from].push_back(num_edges++);
                edge.push_back({e.from, e.to, e.cap - e.flow});
                eid[e.to].push_back(num_edges++);
                edge.push_back({e.to, e.from, e.flow});
            }
        }

        std::vector<int> dis;
        std::vector<int> cur;
        bool bfs(int s, int t)
        {
            if (s > n || t > n)
                return false;
            dis.assign(n + 1, inf_dep);
            cur.assign(n + 1, 0);
            std::queue<int> que;
            dis[s] = 0;
            que.push(s);
            while (que.size())
            {
                int u = que.front();
                que.pop();
                for (auto i : eid[u])
                {
                    const auto &e = edge[i];
                    if (e.cap && dis[e.to] > dis[u] + 1)
                    {
                        dis[e.to] = dis[u] + 1;
                        que.push(e.to);
                    }
                }
            }
            return dis[t] < inf_dep;
        }

        flow_t dfs(int s, int t, flow_t flim)
        {
            if (s == t)
                return flim;
            flow_t flow = 0;
            for (int &i = cur[s]; i < eid[s].size() && flow < flim; i++)
            {
                auto &e = edge[eid[s][i]];
                if (dis[e.to] == dis[s] + 1 && e.cap)
                {
                    auto detf = dfs(e.to, t, std::min(flim - flow, e.cap));
                    flow += detf;
                    e.cap -= detf;
                    edge[eid[s][i] ^ 1].cap += detf;
                }
                if (flow == flim)
                    break;
            }
            return flow;
        }
        flow_t maxflow(int s, int t)
        {
            flow_t flow = 0;
            while (bfs(s, t))
            {
                flow += dfs(s, t, inf_flow);
            }
            return flow;
        }

        std::vector<FlowEdge> to_edge()
        {
            std::vector<FlowEdge> edges;
            for (int i = 0; i < edge.size(); i += 2)
                edges.push_back({
                    .from = edge[i].from,
                    .to = edge[i].to,
                    .cap = edge[i].cap + edge[i ^ 1].cap,
                    .low = 0,
                    .flow = edge[i ^ 1].cap,
                });
            return edges;
        }
    };

    struct Processor
    {
        std::vector<bool> neg;
        std::vector<flow_t> low;
        std::vector<flow_t> excess;
        void init(std::vector<FlowEdge> &edges)
        {
            int n = num_node(edges);
            neg.clear();
            neg.reserve(edges.size());
            low.clear();
            low.reserve(edges.size());
            excess.assign(n + 1, 0);
        }
        void rmv_low(std::vector<FlowEdge> &edges)
        {
            for (auto &e : edges)
            {
                low.push_back(e.low);
                if (e.flow >= e.low)
                {
                    e.flow -= e.low;
                }
                else
                {
                    excess[e.from] -= e.low - e.flow;
                    excess[e.to] += e.low - e.flow;
                    e.flow = 0;
                }
                e.cap -= e.low;
                e.low = 0;
            }
        }
        void add_low(std::vector<FlowEdge> &edges)
        {
            reverse(low.begin(), low.end());
            for (auto &e : edges)
            {
                e.low = low.back();
                e.flow += e.low;
                e.cap += e.low;
                low.pop_back();
            }
        }
    };

    bool excess_flow(std::vector<FlowEdge> &edges,
                     const std::vector<flow_t> &excess)
    {
        int n = num_node(edges), m = edges.size();
        for (int i = 1; i <= n; i++)
        {
            if (excess[i] > 0)
                edges.push_back({.from = n + 1, .to = i, .cap = excess[i]});
            if (excess[i] < 0)
                edges.push_back({.from = i, .to = n + 2, .cap = -excess[i]});
        }
        MaxFlow g;
        g.build(edges);
        g.maxflow(n + 1, n + 2);
        edges = g.to_edge();
        for (int i = m; i < edges.size(); i++)
            if (edges[i].flow != edges[i].cap)
                return false;
        edges.resize(m);
        return true;
    }

    std::optional<flow_t> feasible_flow(std::vector<FlowEdge> &edges, int s = 0,
                                        int t = 0)
    {
        if (s && t)
            edges.push_back({.from = t, .to = s, .cap = inf_flow});
        Processor p;
        p.init(edges);
        p.rmv_low(edges);
        if (!excess_flow(edges, p.excess))
            return std::nullopt;
        if (s && t)
            edges.pop_back();
        p.add_low(edges);
        return get_flow(edges, s);
    }

    std::optional<flow_t> maximum_flow(std::vector<FlowEdge> &edges, int s, int t)
    {
        edges.push_back({.from = t, .to = s, .cap = inf_flow});
        Processor p;
        p.init(edges);
        p.rmv_low(edges);
        if (!excess_flow(edges, p.excess))
            return std::nullopt;
        edges.pop_back();
        MaxFlow g;
        g.build(edges);
        g.maxflow(s, t);
        edges = g.to_edge();
        p.add_low(edges);
        return get_flow(edges, s);
    }

    std::optional<flow_t> minimum_flow(std::vector<FlowEdge> &edges, int s, int t)
    {
        edges.push_back({.from = t, .to = s, .cap = inf_flow});
        Processor p;
        p.init(edges);
        p.rmv_low(edges);
        if (!excess_flow(edges, p.excess))
            return std::nullopt;
        edges.pop_back();
        MaxFlow g;
        Processor q;
        excess_flow(edges, q.excess);
        g.build(edges);
        g.maxflow(t, s);
        edges = g.to_edge();
        p.add_low(edges);
        return get_flow(edges, s);
    }

} // namespace maxflow
  // ed.push_back({.from=x, .to=y, .cap=z});
using namespace maxflow;

int n, m, c[N], d[N], e[N];
void solve()
{
    cin >> n >> m;
    for (int i = 1; i <= n * 2; i++)
        cin >> c[i] >> d[i] >> e[i];
    vector<FlowEdge> E;
    int S = 6 * n + 2, T = 6 * n + 1;
    for (int i = 1; i <= n; i++)
    {
        E.push_back({.from = S, .to = 6 * i - 5, .cap = inf_flow + d[2 * i - 1]});
        E.push_back({.from = S, .to = 6 * i - 4, .cap = inf_flow});
        E.push_back({.from = S, .to = 6 * i - 3, .cap = inf_flow + d[2 * i]});

        E.push_back({.from = 6 * i - 2, .to = T, .cap = inf_flow + c[2 * i - 1]});
        E.push_back({.from = 6 * i - 1, .to = T, .cap = inf_flow});
        E.push_back({.from = 6 * i, .to = T, .cap = inf_flow + c[2 * i]});

        E.push_back({.from = 6 * i - 5, .to = 6 * i - 2, .cap = inf_flow * inf_flow});
        E.push_back({.from = 6 * i - 3, .to = 6 * i, .cap = inf_flow * inf_flow});
        E.push_back({.from = 6 * i - 4, .to = 6 * i - 2, .cap = inf_flow * inf_flow});
        E.push_back({.from = 6 * i - 4, .to = 6 * i - 1, .cap = inf_flow * inf_flow});
        E.push_back({.from = 6 * i - 4, .to = 6 * i, .cap = inf_flow * inf_flow});

        E.push_back({.from = 6 * i - 5, .to = 6 * i, .cap = inf_flow + e[2 * i - 1]});
        E.push_back({.from = 6 * i - 3, .to = 6 * i - 2, .cap = inf_flow + e[2 * i]});
    }
    for (int i = 1, a, b, x, y; i <= m; i++)
    {
        cin >> a >> b >> x >> y;
        a = (a + 1) / 2, b = (b + 1) / 2;
        E.push_back({.from = 6 * b - 4, .to = 6 * a - 2, .cap = x});
        E.push_back({.from = 6 * a - 5, .to = 6 * b, .cap = y});
    }
    i64 ans = *maximum_flow(E, S, T) - inf_flow * n * 3;
    cout << ans << '\n';
}

signed main()
{
    ios_base::sync_with_stdio(false);
    cin.tie(nullptr);
    int t = 1;
    // cin >> t;
    cout << fixed << setprecision(10);
    while (t--)
        solve();
}

詳細信息

Test #1:

score: 0
Time Limit Exceeded

input:

5000 10000
23060775 12 2
255978380 28 517
5 6624 26
151149 45131806 23849036
489 484971 24970
162846 1993316 188305
56311199 2003 211
1 50534913 517527
364913 882765 298
71 26 122914059
13 65459 18150033
20 607 8
380059068 3873712 228
9813 5449 6370
3309369 37410691 8
181 1 62340851
1705 4 107
8 209...

output:


result: