1 #include <algorithm>
  2 #include <iostream>
  3 #include <fstream>
  4 #include <cstring>
  5 #include <cstdlib>
  6 #include <cmath>
  7 using namespace std;
  8 ifstream fin("popular.in");
  9 ofstream fout("popular.out");
 10 struct ls
 11 {
 12  int nw;
 13  int nxt;
 14 };
 15 ls qxx[50001];//链式前向星
 16 int cows=0,gx=0,bian=0,qlt=0,cs=0;
 17 int qlts[10001]={0};//每个强连通分量包含的点数
 18 int DFN[10001]={0},LOW[10001]={0},tou[10001]={0};
 19 int outs[10001]={0};//每个强连通分量的出度
 20 int from[10001]={0};//每个点所属的强连通分量
 21 bool rz[10001]={0};//每个点是否入栈
 22 bool pd[10001][10001]={0};//两个强连通分量之间是否相通
 23 int zhan[10001]={0},top=0;
 24 void add(int fr,int to);
 25 void add2(int fr,int to);
 26 void tarjan(int nw);
 27 int main(void)
 28 {
 29  fin>>cows>>gx;
 30  int a=0,b=0;
 31  memset(tou,-1,sizeof(tou));
 32  for(int i=1;i<=gx;i++)
 33     {
 34      fin>>a>>b;
 35      add(a,b);
 36     }
 37  for(int i=1;i<=cows;i++)
 38     {
 39      if(!DFN[i])tarjan(i);
 40     }
 41  int v=0;
 42  for(int i=1;i<=cows;i++)
 43     {
 44         for(int j=tou[i];j>0;j=qxx[j].nxt)
 45         {
 46          v=qxx[j].nw;
 47          if(from[i]!=from[v]&&!pd[from[i]][from[v]])add2(from[i],from[v]);
 48         }
 49     }
 50  int total=0,ans=0;
 51  for(int i=1;i<=qlt;i++)
 52     {
 53      if(outs[i]==0)
 54        {
 55         total++;
 56         ans+=qlts[i];
 57        }
 58     }
 59  if(total==1)fout<<ans;
 60  else fout<<"0";
 61  return 0;
 62 }
 63
 64 void add(int fr,int to)
 65 {
 66  bian++;
 67  qxx[bian].nw=to;
 68  qxx[bian].nxt=tou[fr];
 69  tou[fr]=bian;
 70 }
 71
 72 void tarjan(int nw)
 73 {
 74  top++;
 75  DFN[nw]=LOW[nw]=++cs;
 76  zhan[top]=nw;
 77  rz[nw]=true;
 78  int v=0;
 79  for(int i=tou[nw];i>0;i=qxx[i].nxt)
 80     {
 81      v=qxx[i].nw;
 82      if(!DFN[v])
 83        {
 84         tarjan(v);
 85         LOW[nw]=min(LOW[v],LOW[nw]);
 86        }
 87      else
 88        {
 89         if(rz[v]==true)
 90           {
 91            LOW[nw]=min(LOW[nw],DFN[v]);
 92           }
 93        }
 94     }
 95  if(DFN[nw]==LOW[nw])
 96    {
 97     qlt++;
 98     do
 99     {
100      v=zhan[top];
101      from[v]=qlt;
102      qlts[qlt]++;
103      rz[v]=false;
104      top--;
105     }while(nw!=v);
106    }
107  return;
108 }
109
110 void add2(int fr,int to)
111 {
112  pd[fr][to]=true;
113  outs[fr]++;
114  return;
115 }