强连通+最小路径覆盖

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

强连通+最小路径覆盖
*/
#include<iostream>
#include<string>
#include<cstdio>
using namespace std;
int head[5005],n,shead[5005];
int DFN[5005],Low[5005],Stack[5005],Index,top,scc,Belong[5005],ecnt,secnt;
bool Instack[5005],vis[5005];
int link[50005];
struct edge
{
int u,v,next;
}E[1000000],E1[1000000];
void Insert(int u,int v)
{
E[ecnt].u=u;
E[ecnt].v=v;
E[ecnt].next=head[u];
head[u]=ecnt++;
}
void SInsert(int u,int v)//缩完点以后再次构图插边
{
E1[secnt].u=u;
E1[secnt].v=v;
E1[secnt].next=shead[u];
shead[u]=secnt++;
}
bool find(int x)
{
int i,v;
for(i=shead[x];i!=-1;i=E1[i].next)
{
v=E1[i].v;
if(!vis[v])
{
vis[v]=1;
if(link[v]==-1||find(link[v]))
{
link[v]=x;
return 1;
}
}
}
return 0;
}
int MaxMatch()
{
int i,ans=0;
memset(link,-1,sizeof(link));
for(i=1;i<=scc;i++)//注意点是重新构图以后1—scc进行找增广路
{
memset(vis,0,sizeof(vis));
if(find(i))ans++;
}
return ans;
}
void Tarjan(int u)
{
int i,v;
Low[u]=DFN[u]=++Index;
Stack[++top]=u;
Instack[u]=1;
for(i=head[u];i!=-1;i=E[i].next)
{
v=E[i].v;
if(!DFN[v])
{
Tarjan(v);
if(Low[u]>Low[v])
Low[u]=Low[v];
}
else if(Instack[v]&&Low[u]>DFN[v])
Low[u]=DFN[v];
}
if(Low[u]==DFN[u])
{
scc++;
do{
v=Stack[top--];
Instack[v]=0;
Belong[v]=scc;
}while(u!=v);
}
return;
}

void Init()
{
int i,u,v,m;
ecnt=0;
secnt=0;
memset(shead,-1,sizeof(shead));
memset(head,-1,sizeof(head));
scanf("%d%d",&n,&m);
for(i=1;i<=m;i++)
{
scanf("%d%d",&u,&v);
Insert(u,v);
}
}
void Solve()
{
int i,v,u;
scc=top=Index=0;
memset(DFN,0,sizeof(DFN));
memset(Low,0,sizeof(Low));
memset(Instack,0,sizeof(Instack));
for(i=1;i<=n;i++)
{
if(!DFN[i])Tarjan(i);
}
//for(i=1;i<=n;i++)
//printf("~%d ",Belong[i]);
//printf("\n");
for(i=0;i<ecnt;i++)
{
u=E[i].u;
v=E[i].v;
if(Belong[u]!=Belong[v])
{
SInsert(Belong[u],Belong[v]);
//printf("!! %d %d\n",Belong[u],Belong[v]);
}
}
int sum=MaxMatch();
printf("%d\n",scc-sum);

}
int main()
{
int cas;
scanf("%d",&cas);
while(cas--)
{
Init();
Solve();
}
return 0;
}
/*
6
6
1 2
2 3
3 4
4 1
3 6
3 5

8 9
1 2
2 3
3 1
1 4
4 6
1 5
5 6
6 7
5 8

*/