POJ-2987-Firing

最新推荐文章于 2020-05-19 08:38:46 发布

南宮逸辰

最新推荐文章于 2020-05-19 08:38:46 发布

阅读量478

点赞数

CC 4.0 BY-SA版权

分类专栏： POJ ACM 文章标签：最大权闭合图

本文链接：https://blog.youkuaiyun.com/z309241990/article/details/38662705

ACM 同时被 2 个专栏收录

716 篇文章

订阅专栏

POJ

378 篇文章

订阅专栏

本文深入探讨了最大闭权闭合图的求解方法，并通过实例展示了其在实际问题中的应用，包括使用Dinic算法进行优化。

这个题依然是求最大闭权闭合图，在胡博涛论文中有相关介绍，我用Dinic做的~

代码：

#include<cstdio>
#include<cstring>
#include<iostream>
#include<queue>
using namespace std;
const long long inf=1LL<<40;
const int maxn=5010;
const int maxm=2e5;
int e,st,des,head[maxn],pnt[maxm],nxt[maxm],level[maxn];
int n,m;
long long flow[maxm];
bool vis[maxn];
queue<int> q;
void AddEdge(int u,int v,long long f)
{
    pnt[e]=v;nxt[e]=head[u];flow[e]=f;head[u]=e++;
    pnt[e]=u;nxt[e]=head[v];flow[e]=0;head[v]=e++;
}
bool BFS()
{
    memset(level,0,sizeof(level));
    q.push(st);
    level[st]=1;
    while(!q.empty())
    {
        int u=q.front();
        q.pop();
        for(int i=head[u];i!=-1;i=nxt[i])
            if(flow[i]&&!level[pnt[i]])
            {
                level[pnt[i]]=level[u]+1;
                q.push(pnt[i]);
            }
    }
    return level[des];
}
long long DFS(int u,long long maxf)
{
    if(u==des||maxf<=0)
        return max(maxf,0LL);
    long long f=0;
    for(int i=head[u];i!=-1;i=nxt[i])
        if(flow[i]&&level[pnt[i]]==level[u]+1)
        {
            long long t=DFS(pnt[i],min(flow[i],maxf-f));
            if(t>0)
            {
                flow[i]-=t;
                flow[i^1]+=t;
                f+=t;
                if(f==maxf)
                    break;
            }
            else
                level[pnt[i]]=0;
        }
    return f;
}
long long maxflow()
{
    long long ans=0;
    while(BFS())
        ans+=DFS(st,inf);
    return ans;
}
void dfs(int u)
{
    for(int i=head[u];i!=-1;i=nxt[i])
        if(flow[i]&&!vis[pnt[i]])
        {
            vis[pnt[i]]=1;
            dfs(pnt[i]);
        }

}
void solve(long long profit)
{
    long long f=maxflow();
    memset(vis,0,sizeof(vis));
    vis[st]=1;
    dfs(st);
    int ans=0;
    for(int i=1;i<=n;i++)
        if(vis[i])
            ans++;
    printf("%d %I64d\n",ans,profit-f);
}
int main()
{
    while(scanf("%d%d",&n,&m)!=EOF)
    {
        e=st=0;des=n+1;
        memset(head,-1,sizeof(head));
        long long profit=0;
        for(int i=1;i<=n;i++)
        {
            int val;
            scanf("%d",&val);
            if(val>0)
            {
                AddEdge(st,i,val);
                profit+=val;
            }
            else
                AddEdge(i,des,-val);
        }
        for(int i=0;i<m;i++)
        {
            int u,v;
            scanf("%d%d",&u,&v);
            AddEdge(u,v,inf);
        }
        solve(profit);
    } 
    return 0;
}