最大权闭合子图的解法

给出一些权值为正和一些权值为负的点，然后点与点之间有依赖关系，如选A点就必须选择B点之类，问最后选择一些点，使得权值之和最大。

利用最小割去求解此问题，将正权值的点与源点建一条流量为权值的边，负权值的点与汇点建一条流量为权值的绝对值的边。有依赖关系的点之间建立一条流量为无穷的边。

显然求出的最小割中，不会有无穷的这条割边，割只会是与源汇直接关联的边（两种边），最后答案就是正权值之和-最大流（最小割）。

我们可以这样思考，我们先假设把所有的正权值的点选上，不考虑所要支付的代价。要是最后权值之和最大，也就要删掉尽可能小的权值或者支付尽可能小的代价。假如最小割中有一条边是关联源点和某一个正权值的点，就相当于我不选择这个点，所以sum要减去这个权值，也就是这条割边的流量。最小割中的另一种边就是与汇点相连的负权值的点，减去这条割边的流量，就相当于我要支付这个权值的代价（由于点与点之间的依赖关系）。这样看来，求的sum-最小割就是我们得到最后答案。

下面附上代码。

#include <cstdio>
#include <cstring>
#include <vector>
#include <algorithm>
#define MAXN 205  
#define MAXM 20008
#define INF 0x3f3f3f3f 
#define max(a,b) (a>b?a:b)  
#define min(a,b) (a<b?a:b)  

using namespace std;  

struct edge  
{  
    int u,v,next,w;
    edge(){}
    edge(int _u, int _v, int _w, int _ne):u(_u), v(_v), w(_w), next(_ne){} 
}E[MAXM];
bool mark[MAXN];
int head[MAXN],size;  
int gap[MAXN],cur[MAXN],pre[MAXN],dis[MAXN];  
int N,M,scr,sink,vn,cas,cnt;  
int U[MAXM],V[MAXM],D[MAXM];  
int W[MAXM];
vector <int>ans;
void add(int u,int v,int w)  
{
     E[size] = edge(u, v, w, head[u]);
     head[u] = size++;
     E[size] = edge(v, u, 0, head[v]);
     head[v] = size++; 
}

void Init()  
{  
     memset (mark, false, sizeof (mark));
    memset(head,-1,sizeof(head));
    size = 0;    
}

int Sap(int s,int t,int n)//s表示起点，t表示终点，n表示点的个数 
{  
    int ans=0,aug=INF;//aug表示增广路的流量  
    int i,v,u=pre[s]=s;  
    for(i=0;i<=n;i++)  
    {  
        cur[i]=head[i];  
        dis[i]=gap[i]=0;  
    }  
    gap[s]=n;  
    bool flag;  
    while(dis[s]<n)  
    {  
        flag=false;  
        for(int &j=cur[u];j!=-1;j=E[j].next)
        {  
            v=E[j].v;  
            if(E[j].w>0&&dis[u]==dis[v]+1)  
            {  
                flag=true;//找到容许边  
                aug=min(aug,E[j].w);  
                pre[v]=u;  
                u=v;  
                if(u==t)  
                {  
                    ans+=aug;  
                    while(u!=s)  
                    {  
                        u=pre[u];  
                        E[cur[u]].w-=aug;  
                        E[cur[u]^1].w+=aug;//注意  
                    }  
                    aug=INF;  
                }  
                break;//找到一条就退出  
            }  
        }  
        if(flag) continue;  
        int mindis=n;  
        for(i=head[u];i!=-1;i=E[i].next)  
        {  
            v=E[i].v;  
            if(E[i].w>0&&dis[v]<mindis)  
            {  
                mindis=dis[v];  
                cur[u]=i;  
            }  
        }  
        if((--gap[dis[u]])==0) break;  
        gap[dis[u]=mindis+1]++;  
        u=pre[u];  
    }  
    return ans;  
}

void find(int u)
{
     mark[u] = true;
     for (int i = head[u];~i;i= E[i].next)
     {
         if (E[i].w > 0 && !mark[E[i].v])
            ans.push_back(E[i].v), find(E[i].v);
     }
}
int main()
{
    int i, n, end, st, cnt, v, d, sum;  
    while (scanf ("%d", &n), n)
    {
          sum = 0, st = 0, end = n+1;
          Init();
          for (i = 1;i <= n;i++)
          {
              scanf ("%d", &d);
              if (d < 0) add(i, end, -d);
              else sum += d, add(st, i, d);
              scanf ("%d", &cnt);
              while (cnt--)
              {
                    scanf ("%d", &v);
                    add(i, v, INF);
              }
          }
          int flow = sum - Sap(st, end, n+2);
          if (flow <= 0) puts("Refused");
          else
          {
              ans.clear();
              find(0);
              sort (ans.begin(), ans.end());
              printf ("%d %d\n", flow, ans.size());
              for (i = 0;i < ans.size();i++)
              {
                  if(i) printf (" ");
                  printf ("%d", ans[i]);
              }
              puts("");
          }
    }
}