BZOJ1143 [CTSC2008] 祭祀river

#include<cstdio>
#include<cstring>
#include<algorithm>

using namespace std;

const int maxn=210;
const int maxm=1010;
const int INF=0x3f3f3f3f;

struct Node{
    int data,next,low;
}node[maxn*maxn*2+2*maxn];

#define now node[point].data
#define then node[point].next
#define www node[point].low

int n,m,cnt;
int s,t,ans;
int head[maxn],cur[maxn];
int dis[maxn],que[maxn];
bool f[maxn][maxn];

void add(int u,int v,int w){
    node[cnt].data=v;node[cnt].next=head[u];node[cnt].low=w;head[u]=cnt++;
    node[cnt].data=u;node[cnt].next=head[v];node[cnt].low=0;head[v]=cnt++;
}

bool BFS(){
    memset(dis,-1,sizeof(dis));
    int H=0,T=1;que[1]=0;dis[s]=0;
    while(H<T){
        H++;
        for(int point=head[que[H]];point!=-1;point=then)
            if(www && dis[now]<0){
                dis[now]=dis[que[H]]+1;
                que[++T]=now;
            }
    }
    return dis[t]>0;
}

int dfs(int x,int low){
    if(x==t) return low;
    int Low;
    for(int &point=cur[x];point!=-1;point=then)
        if(www && dis[now]==dis[x]+1){
            Low=dfs(now,min(low,www));
            if(Low){
                www-=Low,node[point^1].low+=Low;
                return Low;
            }
        }
    return 0;
}

int main(){
#ifndef ONLINE_JUDGE
    freopen("1143.in","r",stdin);
    freopen("1143.out","w",stdout);
#endif
    
    int u,v;
    
    scanf("%d%d",&n,&m);
    t=2*n+1;
    for(int i=s;i<=t;i++) head[i]=-1;
    for(int i=1;i<=n;i++) add(s,i,1),add(i+n,t,1);
    for(int i=1;i<=m;i++){
        scanf("%d%d",&u,&v);
        f[u][v]=true;
    }
    for(int k=1;k<=n;k++)
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                f[i][j]=f[i][j]|(f[i][k]&f[k][j]);
    
    for(int i=1;i<=n;i++)
        for(int j=1;j<=n;j++)
            if(f[i][j] && i!=j) add(i,j+n,1);
    
    int flag;
    
    while(BFS()){
        for(int i=s;i<=t;i++) cur[i]=head[i];
        while(flag=dfs(s,INF))
            ans+=flag;
    }
    
    printf("%d",n-ans);
    return 0;
}