强连通+最小路径覆盖

/*

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

*/