做法：
把图中所有的圈缩成一个点，那么就是求是否存在一个点，使得所有的点都能到达。
遍历所有入度为0的点，对所有遍历到的出度为0的点的标记+1；
若出度为0的点的数目大于两个，则输出0。否则若标记的值等于入度点的和，那么输出这个点缩点之前含的点。
#include<stdio.h>
#include<iostream>
#include<string.h>
#include<algorithm>
#include<queue>
#include<stack>
#include<map>
#include<vector>
#include<stdlib.h>
#define INF_MAX 0x7fffffff
#define INF 999999
#define max3(a,b,c) (max(a,b)>c?max(a,b):c)
#define min3(a,b,c) (min(a,b)<c?min(a,b):c)
#define mem(a,b) memset(a,b,sizeof(a))
using namespace std;
#define ll __int64
#define maxn 10001
#define maxm 50001
struct node
{int u;int v;int w;bool friend operator < (node a, node b){return a.w < b.w;}
}edge[maxn];
ll gcd(ll n,ll m){if(n<m) swap(n,m);return n%m==0?m:gcd(m,n%m);}
ll lcm(ll n,ll m){if(n<m) swap(n,m);return n/gcd(n,m)*m;}
vector<int>vec[maxn];
vector<int>vect[maxn];
stack<int>st;
int dnf[maxn],low[maxn],vis[maxn],sum[maxn],instack[maxn];
int du[maxn],du2[maxn];
int n,m;
int times,num;
void init()
{int i;for(i=0;i<=n;i++)dnf[i]=low[i]=instack[i]=du[i]=du2[i]=vis[i]=sum[i]=0;times=1;num=1;for(i=0;i<=n;i++)vec[i].clear();for(i=0;i<=n;i++)vect[i].clear();while(!st.empty())st.pop();
}
void tarjan(int x)
{int i;dnf[x]=low[x]=times++;instack[x]=1;st.push(x);int n=vec[x].size();for(i=0;i<n;i++){int y=vec[x][i];if(!dnf[y]){tarjan(y);low[x]=min(low[x],low[y]);}else if (instack[y]){low[x]=min(low[x],dnf[y]);}}if(low[x]==dnf[x]){int y=-1;while(y!=x){y=st.top();st.pop();instack[y]=0;sum[num]++;vis[y]=num;}num++;}
}
void jiantu()
{int i,j;for(i=1;i<=n;i++){int len=vec[i].size();for(j=0;j<len;j++){int y=vec[i][j];if(vis[i]==vis[y])continue;vect[vis[i]].push_back(vis[y]);du[vis[y]]++;du2[vis[i]]++;}}
}
void bfs(int x)
{int visit[maxn],i;for(i=0;i<=num;i++){visit[i]=0;}queue<int>q;q.push(x);visit[x]=1;while(!q.empty()){int y=q.front();q.pop();int len=vect[y].size();int leap=0;for(i=0;i<len;i++){leap=1;if(!visit[vect[y][i]]){q.push(vect[y][i]);visit[vect[y][i]]=1;}}if(leap==0){vis[y]++;}}
}
int main()
{int i;while(~scanf("%d%d",&n,&m)){init();int a,b;for(i=0;i<m;i++){scanf("%d%d",&a,&b);vec[a].push_back(b);}for(i=1;i<=n;i++)if(!dnf[i])tarjan(i);int ns=0;jiantu();mem(vis,0);int t1,t2,ip;t1=t2=0;for(i=1;i<num;i++){if(du[i]==0)t1++;if(du2[i]==0)t2=i,ip=i;}if(t2>1){cout<<"0"<<endl;continue;}for(i=1;i<num;i++){if(!du[i]){bfs(i),ns++;}}int sums=0;if(vis[ip]==t1)sums=sum[ip];cout<<sums<<endl;}return 0;
}