HDU 3861 The King’s Problem

http://acm.hdu.edu.cn/showproblem.php?pid=3861

#include<iostream>
#include<vector>
using namespace std;
#define N 5005
#define M 100005

int top,bcnt;
int stack[N],index;
bool instack[N];
int belong[N];
int dfn[N],low[N];
int out[N],in[N];
vector<int> V[N];
struct node{
	int next,v;
	node(){};
	node(int a,int b){
		next=a,v=b;
	}
}E[M];
int head[N],NE;
int pre[2*N];
bool h[2*N];
void init(){
	NE=0;bcnt=0;top=0;index=0;
	memset(head,-1,sizeof(head));
	memset(instack,0,sizeof(instack));
	memset(dfn,0,sizeof(dfn));
	memset(low,0,sizeof(low));
}
void insert(int u,int v){
	E[NE]=node(head[u],v);
	head[u]=NE++;
}

void tarjan(int u){
	dfn[u]=low[u]=++index;
	instack[u]=true;
	stack[top++]=u;
	for(int i=head[u];i!=-1;i=E[i].next){
		int v=E[i].v;
		if(!dfn[v]){
			tarjan(v);
			if(low[v]<low[u])
				low[u]=low[v];
		}
		else if(instack[v]&&dfn[v]<low[u])
			low[u]=dfn[v];
	}
	if(low[u]==dfn[u]){
		bcnt++;
		int v;
		do{
			v=stack[--top];
			instack[v]=false;
			belong[v]=bcnt;
		}while(u!=v);
	}
}

bool dfs(int u){
	for(int i=0;i<V[u].size();i++){
		int v=V[u][i];
		if(h[v]) continue;
		h[v]=1;
		if(pre[v]==-1||dfs(pre[v])){
			pre[v]=u;
			return true;
		}
	}
	return false;
}
int u[M],v[M];
int main(void){
	int t,n,m;
	scanf("%d",&t);
	while(t--){
		scanf("%d%d",&n,&m);
		init();
		for(int i=1;i<=m;i++){
			scanf("%d%d",&u[i],&v[i]);
			insert(u[i],v[i]);
		}
		for(int i=1;i<=n;i++)
			if(!dfn[i])
				tarjan(i);
		if(bcnt==1){
			printf("1\n");
			continue;
		}
		for(int i=1;i<=bcnt;i++)
			V[i].clear();
		for(int i=1;i<=m;i++){
			int x=belong[u[i]];
			int y=belong[v[i]];
			if(x!=y){
				V[x].push_back(y+bcnt);
			}
		}
		memset(pre,-1,sizeof(pre));
		int cn=0;
		for(int i=1;i<=bcnt;i++){
			memset(h,0,sizeof(h));
			if(dfs(i))
				cn++;
		}
		printf("%d\n",bcnt-cn);
	}
}

强连通缩点之后进行二分匹配。。。最小路径覆盖点。。。好题啊！比赛的时候想不到要用二分匹配，水啊。。。