17年ICPC新疆網路賽F Islands(強連通 縮點)
阿新 • • 發佈:2019-02-16
題目大意
給你一個有向無環圖(DAG),問你最少加入幾條邊可以將這個圖變成強連通分析
裸的一道強連通題,縮點後統計入度為0和出度為0的個數,輸出這兩者間的最大值即可程式碼
#include<cstdio>
#include<iostream>
#include<cmath>
#include<cstring>
#include<cstdlib>
#include<queue>
#include<map>
#include<algorithm>
#include<set>
#include<stack>
using namespace std;
const int MAXN=40005;
const int MAXM=400005;
int n;
int m;
int DFN[MAXN];
int LOW[MAXN];
int vis[MAXN];
int belong[MAXN];//belong[i]表示i屬於縮點後的哪個節點
int cnt;
int out[MAXN];
int in[MAXN];
int tot;
struct Edge
{
int v;
int next;
}edge[MAXM];
int edgecount;
int head[MAXN];
void Init()
{
edgecount=0 ;
memset(head,-1,sizeof(head));
}
void Add_edge(int u,int v)
{
edge[++edgecount].v=v;
edge[edgecount].next=head[u];
head[u]=edgecount;
}
stack<int > St;
void Tarjan(int u)//從節點x開始搜尋
{
DFN[u]=LOW[u]=++tot;
vis[u]=1;//為1表示在佇列裡面
St.push(u);
for(int k=head[u];k!=-1 ;k=edge[k].next)
{
int v=edge[k].v;
if(!DFN[v])//還未訪問過
{
Tarjan(v);
LOW[u]=min(LOW[u],LOW[v]);
}
else if(vis[v])//被訪問過,還在佇列裡
{
LOW[u]=min(LOW[u],DFN[v]);
}
}
if(LOW[u]==DFN[u])
{
int x;
++cnt;
while(1)
{
x=St.top();
St.pop();
vis[x]=0;
belong[x]=cnt;
if(x==u)break;
}
}
}
void Solve()
{
tot=0;
cnt=0;//縮點後的點數
memset(DFN,0,sizeof(DFN));
memset(LOW,0,sizeof(LOW));
memset(vis,0,sizeof(vis));
memset(in,0,sizeof(in));
memset(out,0,sizeof(out));
while(!St.empty()) St.pop();
for(int i=1;i<=n;i++)
{
if(DFN[i]==0)Tarjan(i);
}
for(int u=1;u<=n;u++)
{
for(int k=head[u];k!=-1;k=edge[k].next)
{
int v=edge[k].v;
if(belong[u]!=belong[v])
{
out[belong[u]]++;
in[belong[v]]++;
}
}
}
int ans1=0;
int ans2=0;
for(int i=1;i<=cnt;i++)
{
if(in[i]==0)ans1++;
if(out[i]==0)ans2++;
}
ans2=max(ans1,ans2);
if(cnt==1)printf("0\n");
else printf("%d\n",ans2);
}
void In()
{
int a,b;
scanf("%d%d",&n,&m);
for(int i=1;i<=m;i++)
{
scanf("%d%d",&a,&b);
Add_edge(a,b);
}
}
int main()
{
int T;
scanf("%d",&T);
while(T--)
{
Init();
In();
Solve();
}
return 0;
}
/*
3
5 5
1 2
2 3
3 1
4 5
5 4
5 4
1 2
2 3
3 1
4 5
5 4
1 2
2 3
3 4
4 5
*/