HDU 5855 Less Time, More profit（最大权闭合子图）

Description
m家商店，n家工厂，第i家商店赢利vi元的前提是建立对应的ki家工厂，建立第i家工厂需要花费wi，用时ti，工厂的建立可以同时进行，问总赢利不小于L所需的最小时间
Input
第一行一整数T表示用例组数，每组用例首先输入三个整数n,m,L表示工厂数，商店数，以及赢利下限，之后m行每行两个整数wi和ti表示建立第i个工厂的花费和用时，最后n行每行首先输入该家商店的赢利wi，然后输入该商店赢利所需建立的工厂数量ki，之后ki个整数表示该商店赢利所需建立的工厂编号
(1<=T<=30,1<=n,m=200,1<=L,ti<=10^9,1<=wi,vi<=30000)
Output
对于每组用例，如果总赢利可以达到L则输出最少用时以及在最少用时下能够取得的最大赢利，否则输出impossible
Sample Input
2

1 1 2
1 5
3 1 1

1 1 3
1 5
3 1 1
Sample Output
Case #1: 5 2
Case #2: impossible
Solution
对于每家商店可以求出其赢利需要的最短时间，即其赢利所需建立工厂的用时最大值，之后枚举这些时间，对于每个时间，将可以赢利的商店看作一些点，将需要建立的工厂看作一些点，那么就有一个依赖关系，即选择商店点就要选择对应的工厂点，商店点有正点权，工厂点有负点权，进而问题转化为求一个最大权闭合子图问题，如果求出的最大权不小于L说明找到一个符合条件的t，更新一个最小值即可，如果对于最大的t所获得的总赢利还小于L则impossible
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
using namespace std;
typedef long long ll;
#define maxn 555 
#define maxm 111111
#define INF 0x3f3f3f3f
int head[maxn],cur[maxn],d[maxn],st[maxm],s,e,no;//s为源点,e为汇点,n为点数,no为边数 
struct point
{
    int u,v,flow,next;
    point(){};
    point(int x,int y,int z,int w):u(x),v(y),next(z),flow(w){};
}p[maxm];
void add(int x,int y,int z)//从x到y建容量为z的边 
{
    p[no]=point(x,y,head[x],z);//前向弧,标号为偶 
    head[x]=no++;
    p[no]=point(y,x,head[y],0);//后向弧,标号为奇 
    head[y]=no++;
}
void init()//初始化 
{
    memset(head,-1,sizeof(head));
    no=0;
}
bool bfs()
{
    int i,x,y;
    queue<int>q;
    memset(d,-1,sizeof(d));
    d[s]=0; 
    q.push(s);
    while(!q.empty())
    {
        x=q.front();    
        q.pop();
        for(i=head[x];i!=-1;i=p[i].next)
        {
            if(p[i].flow&& d[y=p[i].v]<0)
            {
                d[y]=d[x]+1;
                if(y==e)    
                    return true;
                q.push(y);
            }
        }
    }
    return false;
}
int dinic()//最大流 
{
    int i,loc,top,x=s,nowflow,maxflow=0;
    while(bfs())
    {
        memcpy(cur,head,sizeof(head));
        top=0;
        while(true)
        {
            if(x==e)
            {
                nowflow=INF;
                for(i=0;i<top;i++)
                {
                    if(nowflow>p[st[i]].flow)
                    {
                        nowflow=p[st[i]].flow;
                        loc=i;
                    }
                }
                for(i=0;i<top;i++)
                {
                    p[st[i]].flow-=nowflow;
                    p[st[i]^1].flow+=nowflow;
                }
                maxflow+=nowflow;
                top=loc;    
                x=p[st[top]].u;
            }
            for(i=cur[x];i!=-1;i=p[i].next)
                if(p[i].flow&&d[p[i].v]==d[x]+1) 
                    break;
            cur[x]=i;
            if(i!=-1)
            {
                st[top++]=i;
                x=p[i].v;
            }
            else 
            {
                if(!top)    
                    break;
                d[x]=-1;
                x=p[st[--top]].u;
            }
        }
    }
    return maxflow;
}
int T,n,m,w[maxn];
ll L,t[maxn];
struct node
{
    int v;
    ll t;
    vector<int>g;
}shop[maxn];
ll h[maxn];
int vis[maxn];
int Case=1;
bool cmp(ll a,ll b)
{
    return a>b;
}
int main()
{
    scanf("%d",&T);
    while(T--)
    {
        scanf("%d%d%I64d",&n,&m,&L);
        for(int i=1;i<=n;i++)scanf("%d%I64d",&w[i],&t[i]),h[i]=t[i];
        sort(h+1,h+n+1,cmp);
        for(int i=1;i<=m;i++)
        {
            int k,temp;
            scanf("%d%d",&shop[i].v,&k);
            shop[i].g.clear();
            shop[i].t=-1;
            while(k--)
            {
                scanf("%d",&temp);
                shop[i].g.push_back(temp);
                shop[i].t=max(shop[i].t,t[temp]);
            }
        }
        int anst=-1,ansv=0;
        for(int k=1;k<=n;k++)
        {
            memset(vis,0,sizeof(vis));
            init();
            s=0,e=n+m+1;
            int sum=0;
            for(int i=1;i<=m;i++)
                if(shop[i].t<=h[k])
                {
                    sum+=shop[i].v;
                    add(s,i,shop[i].v);
                    for(int j=0;j<shop[i].g.size();j++)
                    {
                        int temp=shop[i].g[j];
                        vis[temp]=1;
                        add(i,temp+m,INF);
                    }
                }
            for(int i=1;i<=n;i++)
                if(vis[i])
                    add(i+m,e,w[i]);
            int ans=sum-dinic();
            if(ans>=L)
                anst=h[k],ansv=ans;
            else break;
        }
        printf("Case #%d: ",Case++);
        if(anst==-1)printf("impossible\n");
        else printf("%d %d\n",anst,ansv);
    }
    return 0;
}