[bzoj1143][CTSC2008]祭祀river

Description

给出一个n个点，m条边的DAG，你可以从中选择一些点，被选择的点所能到达的点都无法被选择。
求最多能选出多少点。
n<=100,m<=1000

Solution

回归bzoj O(∩_∩)O~~
啥时候CTSC的题那么水了~~
首先预处理处每两个点之间能否到达，能到达的点连一条边。
就相当于求最大独立集了。
最大独立集=n-最大匹配，网上一堆资料蛤
然而本蒟蒻不会证
练手写匈牙利~

Code

#include 
#include 
#include 
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define rep(i,a) for(int i=last[a];i;i=next[i])
using namespace std;
const int N=105;
int n,m,x,y,l,ans,id,to[N],bz[N];
int t[N*N],next[N*N],last[N];
bool f[N][N];
void add(int x,int y) {
    t[++l]=y;next[l]=last[x];last[x]=l;
}
bool find(int x) {
    rep(i,x) if (bz[t[i]]!=id) {
        bz[t[i]]=id;
        if (!to[t[i]]||find(to[t[i]])) {
            to[t[i]]=x;
            return 1;
        }
    }
    return 0;
}
int main() {
    scanf("%d%d",&n,&m);ans=n;
    fo(i,1,m) scanf("%d%d",&x,&y),f[x][y]=1;
    fo(k,1,n) fo(i,1,n) fo(j,1,n) f[i][j]=f[i][j]||f[i][k]&&f[k][j];
    fo(i,1,n) fo(j,1,n) if (f[i][j]) add(i,j);
    fo(i,1,n) id++,ans-=find(i);
    printf("%d\n",ans);
}