FZU - 2295 - 最大权闭合子图

题目链接:https://vjudge.net/problem/FZU-2295

 

解题思路:

如果不考虑最后K对是相互对立的话，那么就是模板最大权闭合子图。因为K最大是5，所以最多是有2^5种情况，这个数可以看做常数，那么就是做32遍最大权闭合子图就是了。

 

#include <iostream>
#include <cstdio>
#include <algorithm>
#include <cstring> 
#include <queue>
#define inf 0x3f3f3f3f
using namespace std;
const int mx = 2e2 + 10;
typedef long long ll;
int n,m,K,a[mx][mx],b[mx][mx],ret;
int vi[mx],wi[mx],job[10][2],cur[mx];
int head[mx],tot,c[mx],d[mx],S,T;
bool vis[100];
struct node
{
	int y,v;
	int nxt;
	node(){}
	node(int yy,int vv,int h){ y=yy,v=vv,nxt=h;}
}edge[mx*mx];
void init()
{
	ret = tot = 0;
	memset(head,-1,sizeof(head));
	memset(vis,0,sizeof(vis));
}
void AddEdge(int x,int y,int v)
{
	edge[tot] = node(y,v,head[x]);
	head[x] = tot++;
	edge[tot] = node(x,0,head[y]);
	head[y] = tot++;
}
int dep[mx];
bool bfs()
{
	memset(dep,0,sizeof(dep));
	queue<int> que;
	que.push(S);
	dep[S] = 1;
	while(!que.empty())
	{
		int no = que.front();
		que.pop();
		for(int i=head[no];~i;i=edge[i].nxt)
		{
			int y = edge[i].y;
			if(edge[i].v&&!dep[y]){
				dep[y] = dep[no] + 1;
				que.push(y);
			}
		}
	}
	return dep[T];
}
int dfs(int x,int flow)
{
	if(x==T || !flow) return flow;
	int used = 0;
	for(int& i=cur[x];~i;i=edge[i].nxt)
	{
		int y = edge[i].y;
		if(dep[x]+1==dep[y]){
			int w = dfs(y,min(flow-used,edge[i].v));
			edge[i].v -= w;
			edge[i^1].v += w;
			used += w;
			if(used==flow) return flow;
		}
	}
	//if(!used) dep[x] = 0;
	return used;
}
int maxflow()
{
	int ans = 0;
	while(bfs()){
		copy(head,head+mx,cur);
		ans += dfs(S,inf);
	}
	return ans;
}
void add(int x)
{
	ret += wi[x];
	AddEdge(S,x+n,wi[x]);
	for(int i=1;i<=d[x];i++) AddEdge(x+n,b[x][i],inf);
}
int main()
{
	int t;
	scanf("%d",&t);
	while(t--){
		scanf("%d%d%d",&n,&m,&K);
		S = 0,T = n+m+1;
		for(int i=1;i<=n;i++){
			scanf("%d%d",vi+i,c+i);
			for(int j=1;j<=c[i];j++){
				scanf("%d",a[i]+j);
			}
		}
		for(int i=1;i<=m;i++){
			scanf("%d%d",wi+i,d+i);
			for(int j=1;j<=d[i];j++){
				scanf("%d",b[i]+j);
			}
		}
		for(int i=0;i<K;i++) scanf("%d%d",job[i]+0,job[i]+1);
		int ans = 0;
		for(int i=0;i<(1<<K);i++){
			init();
			for(int j=1;j<=n;j++){
				AddEdge(j,T,vi[j]);
				for(int k=1;k<=c[j];k++){
					AddEdge(j,a[j][k],inf);
				}
			}
			for(int j=0;j<K;j++){
				if((1<<j)&i) vis[job[j][0]] = 1;	
				else vis[job[j][1]] = 1;	
			} 
			for(int j=1;j<=m;j++) if(!vis[j]) add(j);
			ans = max(ans,ret-maxflow());
		}
		printf("%d\n",ans);
	}
	return 0;
}