HDU-2017"百度之星"程序设计大赛-初赛(B)-1005-度度熊的交易计划

最新推荐文章于 2021-08-07 17:53:26 发布

f_zyj

最新推荐文章于 2021-08-07 17:53:26 发布

阅读量653

点赞数 2

分类专栏：网络流文章标签：网络流

本文链接：https://blog.youkuaiyun.com/f_zyj/article/details/77162755

版权

网络流专栏收录该内容

23 篇文章

订阅专栏

本文详细解析了一道典型的费用流问题，并分享了一个因代码细节疏忽而导致长时间调试的经历。通过对模版代码的深入分析，指出了一处容易忽视的bug，并最终解决了超时问题。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

ACM模版

描述

题解

裸的费用流问题，直接拷模版就行了……可是我 $6000ms$ 超时了，一直没有找到问题所在，今天又看了看代码，发现有一个地方少加了 $~$ 导致进入了死循环，也就是当数据读到文件尾后它并不会停止而是继续执行，从而挂了……真是一个打死都想不到的 $bug$ 所在，这种 $bug$ 是最容易被自己忽视的了！

不想说啥了，感觉自己好废啊……本来流就没有搞过几道题，这还出了这问题，太让人失望了我。

代码

#include <cstdio>
#include <iostream>
#include <cstring>
#include <queue>

using namespace std;

const int MAXN = 505;
const int MAXM = 400005;
const int INF = 0x3f3f3f3f;

struct Edge
{
    int v, c, w, next;
    Edge() {}
    Edge(int v, int c, int w, int next) : v(v), c(c), w(w), next(next) {}
} E[MAXM];

int flow, cost, s, t, T, cnt;
int hed[MAXN];
int cap[MAXN];
int vis[MAXN];
int dis[MAXN];
int cur[MAXN];
queue<int> q;

void addedges(int u, int v, int c, int w)
{
    E[cnt] = Edge(v, c, w, hed[u]);
    hed[u] = cnt++;
    E[cnt] = Edge(u, 0, -w, hed[v]);
    hed[v] = cnt++;
}

bool spfa()
{
    memset(dis, INF, sizeof dis);
    cur[s] = -1;
    vis[s] = ++T;
    cap[s] = INF;
    dis[s] = 0;
    q.push(s);

    while (!q.empty())
    {
        int u = q.front();
        q.pop();
        vis[u] = T - 1;
        for (int e = hed[u]; ~e; e = E[e].next)
        {
            int v = E[e].v, c = E[e].c, w = E[e].w;
            if (c && dis[v] > dis[u] + w)
            {
                dis[v] = dis[u] + w;
                cap[v] = min(cap[u], c);
                cur[v] = e;
                if (vis[v] != T)
                {
                    vis[v] = T;
                    q.push(v);
                }
            }
        }
    }
    if (dis[t] > 0)
    {
        return false;
    }

    cost += cap[t] * dis[t];
    flow += cap[t];
    for (int e = cur[t]; ~e; e = cur[E[e ^ 1].v])
    {
        E[e].c -= cap[t];
        E[e ^ 1].c += cap[t];
    }

    return true;
}

int mfmc()
{
    flow = cost = 0;
    while (spfa()) ;
    return cost;
}

void init()
{
    cnt = T = 0;
    memset(hed, -1, sizeof hed);
    memset(vis, 0, sizeof vis);
}

int n, m;
int a[MAXN], b[MAXN], c[MAXN], d[MAXN];
int g[MAXN][MAXN];

template <class T>
inline void scan_d(T &ret)
{
    char c;
    ret = 0;
    while ((c = getchar()) < '0' || c > '9');
    while (c >= '0' && c <= '9')
    {
        ret = ret * 10 + (c - '0'), c = getchar();
    }
}

void work()
{
    init();

    for (int i = 1; i <= n; i++)
    {
        scan_d(a[i]), scan_d(b[i]), scan_d(c[i]), scan_d(d[i]);
    }

    s = 0, t = n + 1;
    for (int i = 0; i <= t; i++)
    {
        for (int j = 0; j <= t; j++)
        {
            if (i == j)
            {
                g[i][j] = 0;
            }
            else
            {
                g[i][j] = INF;
            }
        }
    }
    int u, v, k;
    for (int i = 1; i <= m; i++)
    {
        scan_d(u), scan_d(v), scan_d(k);
        if (g[u][v] > k)
        {
            g[u][v] = g[v][u] = k;
        }
    }

    for (int i = 1; i <= n; i++)
    {
        addedges(i, t, b[i], a[i]);
        addedges(s, i, d[i], -c[i]);
    }

    for (int i = 1; i <= n; i++)
    {
        for (int j = 1; j <= n; j++)
        {
            if (i != j && g[i][j] != INF)
            {
                addedges(i, j, INF, g[i][j]);
            }
        }
    }
    int ans = -mfmc();

    printf("%d\n", ans);
}

int main()
{
    while (scanf("%d%d", &n, &m) != EOF)
    {
        work();
    }

    return 0;
}