hdu_5585_Less Time, More profit(二分+最大权闭合图)

题目链接：hdu_5855_Less Time, More profit

题意：

有n个工厂，每建一个工厂要花费v_i,需要时间t_i,然后有m个商店，每个商店需要在指定的k个工厂中进货，才能盈利，如果其中一个不在，都不能盈利，问在满足利润大于等于k的情况下的最少时间t，和在t时刻的最大利润

题解：

二分时间+网络流的最大权闭合图。比赛的时候搞了半天的题意才看懂，4点过才开始写网络流，然而没时间调试了。然后就GG了

 #include<bits/stdc++.h>
 #define mst(a,b) memset(a,b,sizeof(a));
 #define F(i,a,b) for(int i=a;i<=b;++i)
 using namespace std;

 const int N=,inf=~0U>>,M=;
 struct edge{int t,f;edge*nxt,*pair;}*g[N],*d[N],pool[M],*cur=pool;
 struct ISAP{
     int n,m,i,S,T,h[N],gap[N],maxflow;
     void init(int ss,int tt){for(S=ss,T=tt,cur=pool,i=;i<=T;i++)g[i]=d[i]=NULL,h[i]=gap[i]=;}
     void add(int s,int t,int f){
         edge*p=cur++;p->t=t,p->f=f,p->nxt=g[s],g[s]=p;
         p=cur++,p->t=s,p->f=,p->nxt=g[t],g[t]=p;
         g[s]->pair=g[t],g[t]->pair=g[s];
     }
     int sap(int v,int flow){
         if(v==T)return flow;
         int rec=;
         for(edge*p=d[v];p;p=p->nxt)if(h[v]==h[p->t]+&&p->f){
         int ret=sap(p->t,min(flow-rec,p->f));
         p->f-=ret;p->pair->f+=ret;d[v]=p;
         if((rec+=ret)==flow)return flow;
         }
         if(!(--gap[h[v]]))h[S]=T;
         gap[++h[v]]++;d[v]=g[v];
         return rec;
     }
     int get_ans(){
         for(gap[maxflow=]=T,i=;i<=T;i++)d[i]=g[i];
         while(h[S]<T)maxflow+=sap(S,inf);
         return maxflow;
     }
 }G;

 int t,n,m,l,L,R,ansf,ic=,all;
 struct pat{
     int v,t;
 }pt[];
 struct shop{
     int v,k,vis[];
 }sp[];

 bool check(int x)
 {
     all=;
     G.init(,);
     F(i,,n)if(pt[i].t<=x)G.add(i+,,pt[i].v);
     int vv[];
     mst(vv,);
     F(i,,m)
     {
         int fg=;
         F(j,,sp[i].k)
         if(pt[sp[i].vis[j]].t>x)fg=;
         if(fg)vv[i]=;
     }
     F(i,,m)if(vv[i])G.add(,i+,sp[i].v);
     F(i,,m)if(vv[i])
     {
         all+=sp[i].v;
         F(j,,sp[i].k)G.add(i+,sp[i].vis[j]+,inf);
     }
     ansf=all-G.get_ans();
     if(ansf>=l)return ;
     else return ;
 }

 int main(){
     scanf("%d",&t);
     while(t--)
     {
         ansf=,all=;
         L=1e9,R=;
         scanf("%d%d%d",&n,&m,&l);
         F(i,,n)scanf("%d%d",&pt[i].v,&pt[i].t),L=min(L,pt[i].t),R=max(R,pt[i].t);
         F(i,,m)
         {
             scanf("%d%d",&sp[i].v,&sp[i].k);
             F(j,,sp[i].k)scanf("%d",&sp[i].vis[j]);
         }
         int mid,ans=;
         while(L<=R)
         {
             mid=(L+R)>>;
             if(check(mid))R=(ans=mid)-;
             else L=mid+;
         }
         check(ans);
         if(ansf>=l)printf("Case #%d: %d %d\n",ic++,ans,ansf);
         else printf("Case #%d: impossible\n",ic++);
     }
 }