BZOJ1143: [CTSC2008]祭祀river（最长反链）

传送门

题意：
给一个$DAG$,选出最大的点集，使得他们两两不在同一条链上。

题解：

最长反链裸题。

首先有一个结论是最长反链等于最小链覆盖。

这个理解起来也比较简单，首先最小链覆盖可以覆盖完所有点，所以最长反链$\le$最小链覆盖。
反过来可以用归纳法证明，可以参考vfk的博客。

然后这个可以用网络流优化。

首先考虑最小路径覆盖的建图方法，每个点拆成两个点$i_1,i_2$，然后如果有边$(i,j)$，那么连$(i_1,j_2,1)$。
最后跑最大匹配，再用总的点去减。
这样做的原理是每个点出发选择一个点走过去相当于是减少了一条路径。

同样，最小链覆盖只不过是加了可以经过重复点的性质，那么每个点连$(i_2,i_1,INF)$即可。原理同最大匹配，不同是匹配的条件换成了两个点连通。连的反边相当于是把条件换为连通。

#include <bits/stdc++.h>
using namespace std;
inline int rd() {
    char ch=getchar(); int i=0,f=1;
    while(!isdigit(ch)) {if(ch=='-')f=-1; ch=getchar();}
    while(isdigit(ch)) {i=(i<<1)+(i<<3)+ch-'0'; ch=getchar();}
    return i*f;
}
const int N=1e5+50;
const int M=1e6+50;
const int INF=0x3f3f3f3f;
int n,m;
const int src=2e5+40,des=2e5+41;
struct NetworkFlow {
    queue <int> q;
    int g[N*2],nt[M],vt[M],c[M],ec;
    int cur[N*2],lev[N*2];
    inline void init(){ec=1; memset(g,0,sizeof(g));}
    inline void add(int x,int y,int o){
        nt[++ec]=g[x]; g[x]=ec; vt[ec]=y; c[ec]=o;
        nt[++ec]=g[y]; g[y]=ec; vt[ec]=x; c[ec]=0;
    }
    inline bool bfs() {
        while(!q.empty()) q.pop();
        memset(lev,0,sizeof(lev));
        lev[src]=1; q.push(src);
        while(!q.empty()) {
            int u=q.front(); q.pop();
            for(int e=g[u];e;e=nt[e]) {
                if(!c[e]||lev[vt[e]]) continue;
                lev[vt[e]]=lev[u]+1; q.push(vt[e]);
                if(vt[e]==des) return true; 
            }
        } return false;
    }
    inline int dinic(int x,int f) {
        if(x==des) return f;
        int rs=0;
        for(int &e=cur[x];e;e=nt[e]) {
            if(!c[e]||lev[vt[e]]!=lev[x]+1) continue;
            int o=dinic(vt[e],min(f-rs,c[e]));
            c[e]-=o; c[e^1]+=o; rs+=o;
            if(rs==f) return rs;
        } return lev[x]=0,rs;
    }
    inline int maxflow() {
        int rs=0;
        while(bfs()){
            int tt; memcpy(cur,g,sizeof(cur));
            while((tt=dinic(src,INF))) rs+=tt,memcpy(cur,g,sizeof(cur));
        } return rs;
    }
}G;
int main() {
    n=rd(); m=rd(); G.init();
    for(int i=1;i<=n;i++) G.add(src,i,1),G.add(i+n,des,1),G.add(i+n,i,INF);
    for(int i=1;i<=m;i++) {int x=rd(),y=rd();G.add(x,y+n,INF);}
    printf("%d\n",n-G.maxflow());
}