在遙遠的東方，有一個神祕的民族，自稱Y族。他們世代居住在水面上，奉龍王為神。每逢重大慶典， Y族都
會在水面上舉辦盛大的祭祀活動。我們可以把Y族居住地水系看成一個由岔口和河道組成的網路。每條河道連線著
兩個岔口，並且水在河道內按照一個固定的方向流動。顯然，水系中不會有環流（下圖描述一個環流的例子）。

由於人數眾多的原因，Y族的祭祀活動會在多個岔口上同時舉行。出於對龍王的尊重，這些祭祀地點的選擇必
須非常慎重。準確地說，Y族人認為，如果水流可以從一個祭祀點流到另外一個祭祀點，那麼祭祀就會失去它神聖
的意義。族長希望在保持祭祀神聖性的基礎上，選擇儘可能多的祭祀的地點。
Input
第一行包含兩個用空格隔開的整數N、M，分別表示岔口和河道的數目，岔口從1到N編號。
接下來M行，每行包含兩個用空格隔開的整數u、v，
描述一條連線岔口u和岔口v的河道，水流方向為自u向v。
N≤100M≤1000
Output
第一行包含一個整數K，表示最多能選取的祭祀點的個數。

Sample Input
4 4

1 2

3 4

3 2

4 2
Sample Output
2

【樣例說明】

在樣例給出的水系中，不存在一種方法能夠選擇三個或者三個以上的祭祀點。包含兩個祭祀點的測試點的方案有兩種：

選擇岔口1與岔口3（如樣例輸出第二行），選擇岔口1與岔口4。

水流可以從任意岔口流至岔口2。如果在岔口2建立祭祀點，那麼任意其他岔口都不能建立祭祀點

但是在最優的一種祭祀點的選取方案中我們可以建立兩個祭祀點，所以岔口2不能建立祭祀點。對於其他岔口

至少存在一個最優方案選擇該岔口為祭祀點，所以輸出為1011。

其實這題是一個二分圖的最大獨立集問題

我們只要i，j不可以相互到達建立一個inf的邊
原點到所有的i建立一條1的邊
所有的j向t匯點建立一條1的邊

跑一下最大流即可

#include<bits/stdc++.h>
using namespace std;
const int maxn=120000;
struct node
{
	int v,nxt;
	int w;
}edge[maxn*4];
queue<int>q;
int mapp[1000][1000];
int s,t;
int cnt=0,head[maxn];
const int inf=0x3f3f3f3f;
void add_Edge(int u,int v,int w)
{
	edge[cnt].v=v;
	edge[cnt].w=w;
	edge[cnt].nxt=head[u];
	head[u]=cnt++;
}
int dis[maxn];
int bfs()
{
	int x;
	while(!q.empty())
		q.pop();
	memset(dis,0,sizeof(dis));
	dis[s]=1;
	q.push(s);
	while(!q.empty())
	{
		x=q.front();
		q.pop();
		for(int i=head[x];i!=-1;i=edge[i].nxt)
		{
			if(edge[i].w&&!dis[edge[i].v])
			{
				dis[edge[i].v]=dis[x]+1;
				if(edge[i].v==t)
					return 1;
				q.push(edge[i].v);
			}
		}
	}
	return 0;
}
int dfs(int x,int low)
{
	if(x==t)
		return low;
	int temp=low;
	int k;
	for(int i=head[x];i!=-1;i=edge[i].nxt)
	{
		if(edge[i].w&&dis[edge[i].v]==dis[x]+1)
		{
			k=dfs(edge[i].v,min(temp,edge[i].w));
			if(!k)
				dis[edge[i].v]=0;
			edge[i].w-=k;
			edge[i^1].w+=k;
			if(!(temp-=k))
				break;
		}
	}
	return low-temp;
}
int main ()
{
	int n,m;
	cnt=0;
	memset(head,-1,sizeof(head));
	scanf("%d%d",&n,&m);
	for(int i=1;i<=m;i++)
	{
		int xx,yy;
		scanf("%d%d",&xx,&yy);
		mapp[xx][yy]=1;
	}
		for(int k=1;k<=n;k++)
		{
			for(int i=1;i<=n;i++)
			{
				for(int j=1;j<=n;j++)
				{
					mapp[i][j]|=(mapp[i][k]&mapp[k][j]);
				}
			}
		}
		for(int i=1;i<=n;i++)
		{
			for(int j=1;j<=n;j++)
			{
				if(mapp[i][j])
				{
					add_Edge(i,j+n,inf),add_Edge(j+n,i,0);
				}
			}
		}
		s=2*n+1;
		t=2*n+2;
		for(int i=1;i<=n;i++)
		{
			add_Edge(s,i,1);
			add_Edge(i,s,0);
		}
		for(int i=n+1;i<=2*n;i++)
			add_Edge(i,t,1),add_Edge(t,i,0); 
		int flow=0;
		while(bfs())
		{
			flow+= dfs(s,inf);
		}
		printf("%d\n",n-flow);
}

/*
1 9
1 6
2 9
3 9
3 8
3 6
4 9
4 6

*/