hdu4971 最大权闭合图

题意：有一些任务 然后每一个任务有一个收益，需要一些技术 每一种技术之间存在需要关系和每一种技术需要话费

然后问最大收益是多少

解法：最大圈闭合图的建图方法之外 ，需要先进行强联通，然后根据强联通分量之间进行建任务的建图

比赛的时候完全没有看出来

#include<queue>
#include<stdio.h>
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<string.h>
#include<vector>
#include<stack>
using namespace std;
#define maxn 11111
#define inf 11111111

struct EdgeNode
{
    int from;
    int to;
    int next;
}e[11111];


struct edge{int to,cap,rev;
    edge(int to=0,int cap=0,int rev=0):to(to),cap(cap),rev(rev){}
};
vector<edge>g[maxn];
int lv[maxn],iter[maxn];
void add(int u,int v,int w)
{
    edge a;a.to=v;a.cap=w;a.rev=(int)g[v].size();
    g[u].push_back(a);
    a.to=u;a.cap=0;a.rev=(int)g[u].size()-1;
    g[v].push_back(a);
}
void bfs(int s)
{
    memset(lv,-1,sizeof lv);
    queue<int>q;
    q.push(s);
    lv[s]=0;
    int v;
    while(!q.empty())
    {
        v=q.front();
        q.pop();
        for(int i=0;i<g[v].size();++i)
        {
            if(lv[g[v][i].to]<0&&g[v][i].cap>0)
            {
                lv[g[v][i].to]=lv[v]+1;
                q.push(g[v][i].to);
            }
        }
    }
}
int dfs(int u,int v,int f)
{
    if(u==v)return f;
    for(int &i=iter[u];i<g[u].size();++i)
    {
        edge &e=g[u][i];
        if(lv[u]<lv[e.to]&&e.cap>0)//
        {
            int d=dfs(e.to,v,min(f,e.cap));
            if(d>0)
            {
                e.cap-=d;
                g[e.to][e.rev].cap+=d;
                return d;
            }
        }
    }
    return 0;
}
int max_flow(int s,int t)
{
    int flow=0;
    for(;;)
    {
        memset(lv,-1,sizeof lv);
        bfs(s);
        if(lv[t]<0)return flow;
        memset(iter,0,sizeof iter);
        //        printf("%d\n",flow);
        int f;
        while((f=dfs(s,t,inf))>0){
            flow+=f;
        }
    }
}
int head[maxn],cnt;
void add1(int x,int y){
    e[cnt].from=x;e[cnt].to=y;e[cnt].next=head[x];head[x]=cnt++;
}

void init(){
    cnt=1;
    memset(head,-1,sizeof(head));
}

int dfn[maxn],low[maxn],ins[maxn];
int cixu,scc_count,n,ru[maxn],chu[maxn],sd[maxn];
stack<int>S;

void dfs(int u){
    int v,x;
    dfn[u]=low[u]=++cixu;
    S.push(u);
    ins[u]=1;
    for(int i=head[u];i!=-1;i=e[i].next)
    {
        v=e[i].to;
        if(!dfn[v]){
            dfs(v);
            low[u]=min(low[u],low[v]);
        }
        else if(ins[v])
            //如果已经搜过，并且在栈里面，如果不在栈里面，说明该点属于其他的强连通分量里面
            low[u]=min(low[u],dfn[v]);
        //此处low代替dfn也是可以的，但是在tarjan求无向图的时候不能
    }
    if(dfn[u]==low[u])
        //表示u为此强连通分量第一个被发现的点，属于该强连通分量的所有的其他点都在栈中，他的上面
    {
        scc_count++;
        do
        {
            x=S.top();
            S.pop();
            sd[x]=scc_count;
            ins[x]=0;
        }while(x!=u);
    }
}

void tarjan()
{
    for(int i=1;i<=n;i++)sd[i]=i;   //缩点标记数组
    memset(ins,0,sizeof(ins));
    memset(dfn,0,sizeof(dfn));
    memset(low,0,sizeof(low));
    cixu=scc_count=0;
    for(int i=1;i<=n;i++)     //可能所给的图为不连通的图，分块，因此需要遍历所有的点
    {
        if(!dfn[i])
            dfs(i);
    }
}
int x,y;
void SD()
{
    int xx,yy,len=cnt;
    init();
    memset(ru,0,sizeof(ru));
    memset(chu,0,sizeof(chu));
    for(int i=1;i<len;i++){
        xx=sd[e[i].from]; yy=sd[e[i].to];
        if(xx==yy)continue;
        add(xx+x,yy+x,inf);
//        add1(xx,yy);
        //        printf("xxx%d %d\n",x,y);
//        chu[xx]++; ru[yy]++;
    }
}

int va[maxn],value[maxn];
vector<int>mp[maxn];
int main()
{
    int t,m,_;
    int a,b,c;
    int ans;
//    init();
//    n=3;
//    add1(1,2);
//    add1(2,1);
//    tarjan();
//    SD();
//    printf("%d\n",scc_count);
    scanf("%d",&_);
    for(int z=1;z<=_;++z){

        scanf("%d%d",&x,&y);
        ans=0;

        for(int i=1;i<=x;++i){
            scanf("%d",&a);
            ans+=a;
            add(0,i,a);
//            while(b--){
//                scanf("%d",&c);
//                add(i,c+x,inf);
//            }
        }

        memset(va,0,sizeof va);
        memset(value,0,sizeof value);
        for(int i=1;i<=y;++i){
            scanf("%d",&va[i]);
        }
        for(int i=0;i<=x+1;++i)mp[i].clear();
        
        for(int i=1;i<=x;++i){
            scanf("%d",&b);
            while(b--){
                scanf("%d",&c);
                mp[i].push_back(c+1);
//                add(i,c+x+1,inf);
            }
        }
        init();
        
        for(int i=1;i<=y;++i)
            for(int j=1;j<=y;++j){
                scanf("%d",&c);
                if(c){
                    add1(i,j);
//                    printf("%d %d\n",i,j);
                }
            }
        n=y;
        tarjan();
        SD();
        n=x+scc_count+1;
//        printf("%d\n",scc_count);
        
        for(int i=1;i<=x;++i){
            for(int j=0;j<mp[i].size();++j){
                add(i,x+sd[mp[i][j]],inf);
            }
        }

        for(int i=1;i<=y;++i){
            value[sd[i]]+=va[i];
        }
        
        for(int i=1;i<=scc_count;++i){
            add(x+i,n,value[i]);
        }
        
        printf("Case #%d: %d\n",z,ans-max_flow(0,n));
        for(int i=0;i<=n;++i)g[i].clear();
    }
    return 0;
    
}