BZOJ 1143: [CTSC2008]祭祀river

显然这是一个最小路径覆盖
转化为二分图做
先拆点，把每个点拆成$x_1$和$x_2$这样
如果$x$到$y$有一条边，就连$x_1,y_2$这样一条边
然后跑一遍二分图最大匹配
匈牙利或网络流都行
我匈牙利忘记了只好dinic
最后答案就是原图的点减二分图最大匹配

---

#include<bits/stdc++.h>
#define g getchar()
#define T n*2+1
#define ll long long
#define inf 0x3f3f3f3f
using namespace std;
inline ll read(){
    ll x=0,f=1;char ch=g;
    for(;ch<'0'||ch>'9';ch=g)if(ch=='-')f=-1;
    for(;ch>='0'&&ch<='9';ch=g)x=x*10+ch-'0';
    return x*f;
}
inline void out(ll x){
    int a[25],t=0;
    if(x<0)putchar('-'),x=-x;
    for(;x;x/=10)a[++t]=x%10;
    for(int i=t;i;--i)putchar('0'+a[i]);
    if(t==0)putchar('0');
    putchar('\n');
}
struct re{int v,fl,next;}ed[100001];
int head[505],dep[505],dui[505],mp[201][201];
int e=1,n,ans,m;
inline void ins(int x,int y){
    ed[++e]=(re){y,1,head[x]};head[x]=e;
    ed[++e]=(re){x,0,head[y]};head[y]=e;
}
bool bfs(){
    int tou=1,wei=1;
    memset(dep,-1,sizeof(dep));
    dui[1]=0,dep[0]=0;
    for(;tou<=wei;++tou){
        int u=dui[tou];
        for(int i=head[u];i;i=ed[i].next)
        if(ed[i].fl&&dep[ed[i].v]==-1)
        dep[ed[i].v]=dep[u]+1,dui[++wei]=ed[i].v;
    }
    return dep[T]!=-1;
}
int dfs(int x,int fl){
    if(x==T)return fl;
    int used=0;
    for(int i=head[x];i;i=ed[i].next)
    if(ed[i].fl&&dep[ed[i].v]==dep[x]+1){
        int w=dfs(ed[i].v,min(fl-used,ed[i].fl));
        used+=w;ed[i].fl-=w;ed[i^1].fl+=w;
        if(used==fl)return fl;
    }
    if(!used)dep[x]=-1;
    return used;
}
void dinic(){for(;bfs();)ans+=dfs(0,inf);}
int main(){
    n=read();m=read();
    for(int i=1;i<=m;++i){
        int x=read(),y=read();
        mp[x][y]=1;
    }
    for(int i=1;i<=n;++i)mp[i][i]=1;
    for(int i=1;i<=n;++i)
    for(int j=1;j<=n;++j){
        for(int k=1;k<=n;++k)
        if(mp[i][j])continue;else mp[i][j]=mp[i][k]&&mp[k][j];
    }
    for(int i=1;i<=n;++i){
        for(int j=1;j<=n;++j)
        if(i!=j&&mp[i][j]){
            ins(i,j+n);
        }
        ins(0,i);
        ins(i+n,T);
    }
    dinic();
    out(n-ans);
    return 0;
}