題目傳送門
題面有點坑（對於我這種不讀題面的人來說）
題面上有個大大的圖。我看了一下好像是環誒。
我就以為這題有環。
一開始
我想：
既然有環，那麼我用強聯通縮點之後每個環最多隻能選出一個點咯。
然後我再在剩下的點裡面求點集，使兩兩不能互相到達。

然後我看錯了：
題目說顯然，水系中不會有環流（下圖描述一個環流的例子）
微笑。

然而這道題是裸的最大獨立集。
學（mo）了一發。
師兄blog
最大獨立集=點數-最大匹配。

程式碼實現：（有強聯通版本。。）

#include<cstdio>
#include<cstring>
#include<cstdlib>
 

#include<iostream>
#include<algorithm>
using namespace std;
struct node {
    int x,y,next;
}a[1100];int len,last[110];
void ins(int x,int y) {
    len++;
    a[len].x=x;a[len].y=y;
    a[len].next=last[x];last[x]=len;
}
int belong[110],sta[1100],tp,cnt;
int low[110],dfn[110],id;
bool v[110];
void
 
 dfs(int x) {
    low[x]=dfn[x]=++id;
    sta[++tp]=x;v[x]=true;
    for(int k=last[x];k;k=a[k].next) {
        int y=a[k].y;
        if(dfn[y]==-1) {
            dfs(y);
            low[x]=min(low[x],low[y]);
        }
        else 
            if(v[y]==true)
                low[x]=min(low[x],dfn[y]);
    }
    if
 
(dfn[x]==low[x]) {
        int i;cnt++;
        do {
            i=sta[tp--];
            belong[i]=cnt;
            v[i]=false;
        }while(i!=x);
    }
}
int map[110][110],match[110];
bool chw[110];
bool get_match(int x) {
    for(int i=1;i<=cnt;i++)
        if(map[x][i]==true)
            if(chw[i]==false) {
                chw[i]=true;
                if(match[i]==0||get_match(match[i])==true) {
                    match[i]=x;return true;
                }
            }
    return false;
}
int main() {
    int n,m;scanf("%d%d",&n,&m);
    len=0;memset(last,0,sizeof(last));
    for(int i=1;i<=m;i++) {
        int x,y;scanf("%d%d",&x,&y);
        ins(x,y);
    }
    memset(dfn,-1,sizeof(dfn));
    memset(sta,0,sizeof(sta));
    memset(v,false,sizeof(v));
    memset(low,0,sizeof(low));
    tp=id=cnt=0;
    for(int i=1;i<=n;i++)
        if(dfn[i]==-1)
            dfs(i);
    memset(map,false,sizeof(map));
    for(int i=1;i<=len;i++) 
        if(belong[a[i].x]!=belong[a[i].y])
            map[belong[a[i].x]][belong[a[i].y]]=true;
    for(int k=1;k<=cnt;k++)
        for(int i=1;i<=cnt;i++)
            for(int j=1;j<=cnt;j++)
                if(map[i][k]==true&&map[k][j]==true)
                    map[i][j]=true; //要加floyd使得間接連通也算聯通
    int ans=0;memset(match,0,sizeof(match));
    for(int i=1;i<=cnt;i++) {
        memset(chw,false,sizeof(chw));
        if(get_match(i)==true)
            ans++;
    }
    printf("%d\n",cnt-ans);
    return 0;
}

不加強聯通：

#include<cstdio>
#include<cstring>
#include<cstdlib>
#include<iostream>
#include<algorithm>
using namespace std;
int map[110][110],match[110];
bool chw[110];int n;
bool get_match(int x) {
    for(int i=1;i<=n;i++)
        if(map[x][i]==true)
            if(chw[i]==false) {
                chw[i]=true;
                if(match[i]==0||get_match(match[i])==true) {
                    match[i]=x;return true;
                }
            }
    return false;
}
int main() {
    int m;scanf("%d%d",&n,&m);
    memset(map,false,sizeof(map));
    for(int i=1;i<=m;i++) {
        int x,y;scanf("%d%d",&x,&y);
        map[x][y]=true;
    }
    for(int k=1;k<=n;k++)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                if(map[i][k]==true&&map[k][j]==true)
                    map[i][j]=true;
    int ans=0;
    memset(match,0,sizeof(match));
    for(int i=1;i<=n;i++) {
        memset(chw,false,sizeof(chw));
        if(get_match(i)==true)
            ans++;
    }
    printf("%d\n",n-ans);
    return 0;
}