题意就是有n个工厂，m个商店 每个工厂有建造时间ti，花费payi 每个商店和k个工厂有关，如果这k个工厂都建造了，那么能获利proi 问你求收益（∑pro−∑pay）≥L时，首先满足时间t最小，其次是收益p最大
首先二分时间的答案，然后看那些工厂能建造，然后工厂是花费，商店是收益，并且要与商店有关的工厂都建造了，才能获利，所以这是一个最大权闭合图的模型，把能建造的工厂和汇点t相连，容量为花费，然后如果这个商店有关的工厂都能建造，那么就把商店和有关的工厂连边，容量是INF，然后把源点和能获利的商店连边，容量是收益，然后跑最大流，然后把能获利的商店的利润都加起来得到sum，sum−flow就是最大的收益
代码：
#include <map> #include <set> #include <ctime> #include <stack> #include <queue> #include <cmath> #include <string> #include <vector> #include <cstdio> #include <cctype> #include <cstring> #include <sstream> #include <cstdlib> #include <iostream> #include <algorithm> #pragma comment(linker,"/STACK:102400000,102400000") using namespace std; #define MAX 1000005 #define MAXN 1000005 #define maxnode 205 #define sigma_size 26 #define lson l,m,rt<<1 #define rson m+1,r,rt<<1|1 #define lrt rt<<1 #define rrt rt<<1|1 #define middle int m=(r+l)>>1 #define LL long long #define ull unsigned long long #define mem(x,v) memset(x,v,sizeof(x)) #define lowbit(x) (x&-x) #define pii pair<int,int> #define bits(a) __builtin_popcount(a) #define mk make_pair #define limit 10000 //const int prime = 999983; const int INF = 0x3f3f3f3f; const LL INFF = 0x3f3f; const double pi = acos(-1.0); const double inf = 1e18; const double eps = 1e-4; const LL mod = 772002; const ull mx = 133333331; /*****************************************************/ inline void RI(int &x) { char c; while((c=getchar())<'0' || c>'9'); x=c-'0'; while((c=getchar())>='0' && c<='9') x=(x<<3)+(x<<1)+c-'0'; } /*****************************************************/ const int N = 505; const int M = 100005; struct Isap{ int tot; int d[N],pre[N],cur[N],gap[N]; int head[N]; struct Edge{ int v,next; int cap,flow; }edge[M*2]; void init(){ mem(head,-1); tot=0; } void add_edge(int a, int b, int c){ edge[tot]=(Edge){b,head[a],c,0}; head[a]=tot++; } void add(int a,int b, int c){ add_edge(a,b,c); add_edge(b,a,0); } void set_d(int t){ queue<int> q; mem(d,-1); mem(gap,0); d[t]=0; q.push(t); while(!q.empty()) { int u=q.front();q.pop(); ++gap[d[u]]; for(int i=head[u];i!=-1;i=edge[i].next) { int v=edge[i].v; if(d[v]==-1) { d[v]=d[u]+1; q.push(v); } } } } int sap(int s,int t,int num) { set_d(t); int ans=0,u=s; int flow=INF; memcpy(cur,head,sizeof(head)); while(d[s]<num){ int &i=cur[u]; for(;i!=-1;i=edge[i].next) { int v=edge[i].v; if(edge[i].cap>edge[i].flow&&d[u]==d[v]+1) { u=v; pre[v]=i; flow=min(flow,edge[i].cap-edge[i].flow); if(u==t){ while(u!=s){ int j=pre[u]; edge[j].flow+=flow; edge[j^1].flow-=flow; u=edge[j^1].v; } ans+=flow; flow=INF; } break; } } if(i==-1) { if(--gap[d[u]]==0) break; int dmin=num-1; cur[u]=head[u]; for(int j=head[u];j!=-1;j=edge[j].next) if(edge[j].cap>edge[j].flow) dmin=min(dmin,d[edge[j].v]); d[u]=dmin+1; ++gap[d[u]]; if(u!=s) u=edge[pre[u]^1].v; } } return ans; } }Sap; /*调用方式： Sap.init(); //建边前调用 Sap.add(u, v, c); //在u->v之间建一条容量为c的边 Sap.sap(s, t, n); //s为源点，t为汇点，n为点的数量*/ int p1[205],t1[205]; int p2[205]; vector<int> v[205]; int n,m,L; int check(int x){ Sap.init(); int sum=0; for(int i=1;i<=m;i++){ Sap.add(0,i,p2[i]); int flag=0; for(int j=0;j<v[i].size();j++){ if(t1[v[i][j]]>x) flag=1; } if(!flag){ for(int j=0;j<v[i].size();j++){ Sap.add(i,m+v[i][j],INF); } sum+=p2[i]; } } for(int i=1;i<=n;i++){ if(t1[i]<=x) Sap.add(m+i,n+m+1,p1[i]); } int ans=Sap.sap(0,n+m+1,n+m+2); if(sum-ans>=L) return sum-ans; else return -1; } int main(){ //freopen("froggy.in","r",stdin); //freopen("froggy.out","w",stdout); int t,kase=0; cin>>t; while(t--){ cin>>n>>m>>L; kase++; for(int i=1;i<=n;i++) scanf("%d%d",&p1[i],&t1[i]); for(int i=1;i<=m;i++){ int k; scanf("%d%d",&p2[i],&k); v[i].clear(); while(k--){ int a; scanf("%d",&a); v[i].push_back(a); } } int l=0,r=1e9; int ans=-1; while(l<=r){ int mid=(l+r)/2; int cnt=check(mid); if(cnt!=-1){ ans=cnt; r=mid-1; } else l=mid+1; } printf("Case #%d: ",kase); if(ans==-1) printf("impossible\n"); else printf("%d %d\n",l,ans); } return 0; }