Ambush UVALive - 8267

凯子哥走了，把燚哥拉进占星者号，大概就是以后燚哥负者数学类题目我改行做图论数据结构类的题目。。。

一道网络流题目，读完题以后，我擦好熟悉的题目前几天见过没补，三个人没人看出来是网络流，想着是图论，大鹏偷听到隔壁队伍说是最小割，其实我觉得隔壁队伍没啥参考建议，都一样水，然后大鹏就去试一下WA了，他继续改没改出来（当然改不出来。。。）今天补题问了大佬，大佬说是最大流，仔细想想确实如此最大流没毛病。

可以把ceil（敌军总数量 * 1.0 / 3)看成是每一条边的容量，这样求出来的最大流就是我们需要的士兵数量也就是最终答案

只可意会不可言传自己理解一下就懂了

1.前几遍T了是因为有一行memset的代码我看着他长就把他和一行短的代码交换了位置，交换了位置。。。这个可不能交换。。。

2.有一遍WA了是因为这是一幅无向图一条边的正反容量都是ceil（敌军总数量 * 1.0 / 3)

#include <bits/stdc++.h>

using namespace std;
typedef long long ll;
const int maxn = 1000 + 10;
const int INF = 0x3f3f3f3f;
bool vis[maxn];
int d[maxn];
int cur[maxn];
int n,m,s,t,c,a,b;
struct edge
{
    int from,to,cap,flow;
    edge(int a,int b,int c,int d):from(a),to(b),cap(c),flow(d) {}
};
vector<int> G[maxn];
vector<edge> edges;
void init()
{
    for(int i=1; i<=n; i++) G[i].clear();
    edges.clear();
}
void addedge(int u,int v)
{
    edges.push_back(edge(u,v,c,0));
    edges.push_back(edge(v,u,c,0));
    int t = edges.size();
    G[u].push_back(t-2);
    G[v].push_back(t-1);
}

bool bfs()
{
    queue<int> q;
    q.push(s);
    memset(vis,0,sizeof(vis));
    d[s] = 0,vis[s] = 1;
    while(!q.empty())
    {
        int x = q.front();
        q.pop();
        for(int i=0; i<(int)G[x].size(); i++)
        {
            edge e = edges[G[x][i]];
            if(!vis[e.to] && e.cap > e.flow)
            {
                vis[e.to] = 1;
                d[e.to] = d[x] + 1;
                q.push(e.to);
            }
        }
    }
    return vis[t];
}
int dfs(int x,int a)
{
    if(x == t || a == 0) return a;
    int flow = 0, f;
    for(int &i = cur[x]; i<(int)G[x].size(); i++)
    {
        edge &e = edges[G[x][i]];
        if(d[x] + 1 == d[e.to] && (f = dfs(e.to,min(a,e.cap - e.flow))) > 0)
        {
            e.flow += f;
            edges[G[x][i]^1].flow -= f;
            a -= f;
            flow += f;
            if(a == 0) break;
        }
    }
    return flow;
}
int max_flow()
{
    int flow = 0;
    while(bfs())
    {
        memset(cur,0,sizeof(cur));
        flow += dfs(s,INF);
    }
    return flow;
}
int main()
{
    int k;
    scanf("%d",&k);
    while(k--)
    {
        scanf("%d%d%d%d%d",&n,&m,&s,&t,&c);
        init();
        c = ceil(c * 1.0 / 3);
        for(int i=1; i<=m; i++)
        {
            scanf("%d%d",&a,&b);
            addedge(a,b);
        }
        int ans = max_flow();
        printf("%d\n",ans);
    }
    return 0;
}