这题。。。。。。。做了半天,倒是让我更加了解了强连通分量,是我的第二道强连通。
石建让我做一下这道题,我用模板再变形,改了很久之后,能给出的样例都过了。
discussion上的样例都过了,倒是上面N多人说数据没过AC了。。。。。
你让我们这数据过了没AC的生活如何自理。。。
先放着,哪天有思路再来搞。。
/重新做了一遍,发现之前的思路太乱了,整理之后水过了、、、、、、、、、、、、、
#include <iostream>
#include <cstring>
#include <cstdio>
#include <cstdlib>
using namespace std;
#define MAXN 50010
#define MAXM 100010
struct Edge
{
int v, next;
}edge[MAXM]; //边结点数组
int num[MAXN];
int first[MAXN], stack[MAXN], DFN[MAXN], Low[MAXN], Belong[MAXM];
// first[]头结点数组,stack[]为栈,DFN[]为深搜次序数组,Belong[]为每个结点所对应的强连通分量标号数组
// Low[u]为u结点或者u的子树结点所能追溯到的最早栈中结点的次序号
int instack[10010]; // instack[]为是否在栈中的标记数组
int n, m, cnt, scnt, top, tot;
void init()
{
cnt = 0;
scnt = top = tot = 0; //初始化连通分量标号,次序计数器,栈顶指针为0
memset(first, -1, sizeof(first));
memset(DFN, 0, sizeof(DFN)); //结点搜索的次序编号数组为0,同时可以当是否访问的数组使用
}
void read_graph(int u, int v) //构建邻接表
{
edge[tot].v = v;
edge[tot].next = first[u];
first[u] = tot++;
}
void Tarjan(int v) //Tarjan算法求有向图的强连通分量
{
int min, t;
DFN[v] = Low[v] = ++tot; //cnt为时间戳
instack[v] = 1; //标记在栈中
stack[top++] = v; //入栈
for(int e = first[v]; e != -1; e = edge[e].next)
{ //枚举v的每一条边
int j = edge[e].v; //v所邻接的边
if(!DFN[j])
{ //未被访问
Tarjan(j); //继续向下找
if(Low[v] > Low[j]) Low[v] = Low[j]; // 更新结点v所能到达的最小次数层
}
else if(instack[j] && DFN[j] < Low[v])
{ //如果j结点在栈内,
Low[v] = DFN[j];
}
}
if(DFN[v] == Low[v])
{ //如果节点v是强连通分量的根
scnt++; //连通分量标号加1
do
{
t = stack[--top]; //退栈
instack[t] = 0; //标记不在栈中
Belong[t] = scnt; //出栈结点t属于cnt标号的强连通分量
num[scnt]++;
}while(t != v); //直到将v从栈中退出
}
}
void solve()
{
for(int i = 1; i <= n; i++) //枚举每个结点,搜索连通分量
if(!DFN[i]) //未被访问
Tarjan(i); //则找i结点的连通分量
}
int b1[MAXN],b2[MAXN];
int outdegree[MAXN];
int main()
{
scanf("%d%d",&n,&m);
init();
int temp=1;
int m1=m;
memset(num,0,sizeof(num));
memset(outdegree,0,sizeof(outdegree));
while(m--)
{
int u, v;
scanf("%d%d", &u, &v);
b1[temp]=u;b2[temp++]=v;
read_graph(u, v);
}
solve(); //求强连通分量
for(int i=1;i<=m1;i++)
{
if(Belong[b1[i]]!=Belong[b2[i]])
outdegree[Belong[b1[i]]]++;
}
int tmp=0;int k=0;
for(int i=1;i<=scnt;i++)
{
//cout<<outdegree[i]<<endl;
if(outdegree[i]==0)
{tmp++;k=i;}
}
if(tmp>1)
{
cout<<0<<endl;
}
else if(tmp==1)
{
cout<<num[k]<<endl;
}
return 0;
}