BZOJ 2718: [Violet 4]毕业旅行 二分图最大点独立集 传递闭包

2718: [Violet 4]毕业旅行

Time Limit: 10 Sec  Memory Limit: 128 MB
Submit: 541  Solved: 321
[Submit][Status][Discuss]

Description

Input

Output

最多可选多少景点

Sample Input

7 6
1 2
2 3
5 4
4 3
3 6
6 7

Sample Output

2

HINT

这个题一看就应该知道是二分图最大点独立集

用传递闭包维护一下连通性 再跑匈牙利算法就行了

用bitset优化更快呦 :-)

#include<cmath>
#include<ctime>
#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<complex>
#include<iostream>
#include<algorithm>
#include<iomanip>
#include<vector>
#include<bitset>
#include<string>
#include<queue>
#include<set>
#include<map>
using namespace std;
inline int read()
{
	int x=0,f=1;char ch=getchar();
	while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
	while(ch<='9'&&ch>='0'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
	return x*f;
}
const int N=210;
int n,m,ecnt,last[N],vis[N],match[N],clr;
struct EDGE{int to,nt;}e[(N*N)<<1];
inline void add(int u,int v)
{e[++ecnt]=(EDGE){v,last[u]};last[u]=ecnt;}
bitset<N>b[N];
bool hungary(int u)
{
	for(int i=last[u];i;i=e[i].nt)
	{
		if(clr==vis[e[i].to])continue;
		vis[e[i].to]=clr;
		if(!match[e[i].to]||hungary(match[e[i].to]))
		{
			match[e[i].to]=u;
			return 1;
		}
	}
	return 0;
}
int main()
{
	n=read();m=read();
	for(int i=1;i<=m;i++)
	{b[read()][read()]=1;}
	for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)
	if(b[i][j])b[i]|=b[j];
	for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)
	if(b[i][j])add(i,j);
	int ans=n;
	for(int i=1;i<=n;i++){clr++;if(hungary(i))ans--;}
	printf("%d\n",ans);
	return 0;
}