woj1124(网络流好题)

来源:互联网 发布:ajax股票实时刷新数据 编辑:程序博客网 时间:2024/06/10 19:19
【题目大意】
有 N 支球队,互相之间已经进行了一些比赛,还剩下 M 场没有比。现在给出各支球队目前的总分以及还剩下哪 M 场没有比,问能否合理安排这 M 场比赛的结果,使得第 N 支球队最后的总分大于其他任何一支球队的总分。已知每场比赛胜者得 2 分,败者 0 分,平局则各得 1 分。(1 <= N <= 100, 0 <= M <= 1000)
【建模方法】

首先,所有跟球队 N 相关的比赛都要让球队 N 赢。如果此时仍有某支球队的总分大于等于球队 N 的总分,则已经不可能满足要求;否则按如下方法建图:每场比赛 i(不包括与球队 N 相关的比赛)作为一个点并加边(s, i, 2),每支球队 j(不包括球队 N)作为一个点并加边(j, t, score[N]-score[j]-1),每场比赛向与其关联的两支球队 u, v 连边(i, u, 2), (i, v, 2)。若最大流等于 2*比赛场数(不包括与球队 N 相关的比赛)则可以满足要求

注意数组要开够,m是1000!!!

#include<cstdio>#include<cstring>#include<algorithm>#include<queue>#include<vector>using namespace std;////////////////////////最大流开始//////////////////////////////////////typedef int cap_type;#define MAX_V 1200 + 30 + 16// 用于表示边的结构体(终点、容量、反向边)struct edge{int to, rev;cap_type cap;edge(int to, cap_type cap, int rev) : to(to), cap(cap), rev(rev){}};vector <edge> G[MAX_V];   // 图的邻接表表示int level[MAX_V];      // 顶点到源点的距离标号int iter[MAX_V];       // 当前弧,在其之前的边已经没有用了// 向图中加入一条从from到to的容量为cap的边void add_edge(int from, int to, int cap){G[from].push_back(edge(to, cap, G[to].size()));G[to].push_back(edge(from, 0, G[from].size() - 1));}// 通过BFS计算从源点出发的距离标号void bfs(int s){memset(level, -1, sizeof(level));queue<int> que;level[s] = 0;que.push(s);while (!que.empty()){int v = que.front();que.pop();for (int i = 0; i < G[v].size(); ++i){edge &e = G[v][i];if (e.cap > 0 && level[e.to] < 0){level[e.to] = level[v] + 1;que.push(e.to);}}}}// 通过DFS寻找增广路cap_type dfs(int v, int t, cap_type f){if (v == t){return f;}for (int &i = iter[v]; i < G[v].size(); ++i){edge &e = G[v][i];if (e.cap > 0 && level[v] < level[e.to]){cap_type d = dfs(e.to, t, min(f, e.cap));if (d > 0){e.cap -= d;G[e.to][e.rev].cap += d;return d;}}}return 0;}// 求解从s到t的最大流cap_type max_flow(int s, int t){cap_type flow = 0;for (;;){bfs(s);if (level[t] < 0){return flow;}memset(iter, 0, sizeof(iter));cap_type f;while ((f = dfs(s, t, 0x3f3f3f3f3f3f3f3f)) > 0){flow += f;}}}vector<int>g[105];int grade[MAX_V];///////////////////////////////最大流结束/////////////////////////////////////int main(){    int n,m,u,v;    while(scanf("%d%d",&n,&m)!=EOF)    {        for(int i=0;i<=n+m+50;i++) G[i].clear();        for(int i=0;i<=n;i++) g[i].clear();        int ma = 0,tot = 1,k;        for(int i=1;i<=n;i++) {            scanf("%d",&grade[i]);            if(ma<grade[i]) ma = grade[i],k=i;        }        for(int i=0;i<m;i++)        {            scanf("%d%d",&u,&v);            if(u==n||v==n) grade[n]+=2;            else            {                g[u].push_back(tot);                g[v].push_back(tot);                tot++;            }        }        if(k!=n&&grade[k]>=grade[n])        {            printf("NO\n");            continue;        }        int s = m+n+2; int t = s+1;        for(int i=1;i<n;i++) add_edge(s,i,grade[n]-grade[i]-1);        for(int i=1;i<tot;i++) add_edge(i+n,t,2);        for(int i=1;i<n;i++)        {            int len = g[i].size();            for(int j=0;j<len;j++)            {                int k = g[i][j];                add_edge(i,k+n,2);            }        }        if( max_flow(s,t)==2*(tot-1)) puts("YES");        else puts("NO");    }    return 0;}


原创粉丝点击