BZOJ1143 [CTSC2008] 祭祀river

AC通道：http://www.lydsy.com/JudgeOnline/problem.php?id=1143

题目大意：

　　给你n个点，点与点之间由有向边相连。如果u能到达v的话，那么他们就不能同时选。问最多选多少个点。

[原题很强，可惜这里只有第一问，就变成大水题了...]

首先当然先跑一遍Floyd跑个传递闭包。

然后就是最大独立集了，最大独立集

最长反链
=最小路径覆盖数
=n-二分图最大匹配

然后就可以跑最大匹配了！

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

using namespace std;

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

struct Node{
    int data,next,low;
}node[maxn*maxn*+*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=;head[v]=cnt++;
}

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

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

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=*n+;
    for(int i=s;i<=t;i++) head[i]=-;
    for(int i=;i<=n;i++) add(s,i,),add(i+n,t,);
    for(int i=;i<=m;i++){
        scanf("%d%d",&u,&v);
        f[u][v]=true;
    }
    for(int k=;k<=n;k++)
        for(int i=;i<=n;i++)
            for(int j=;j<=n;j++)
                f[i][j]=f[i][j]|(f[i][k]&f[k][j]);

    for(int i=;i<=n;i++)
        for(int j=;j<=n;j++)
            if(f[i][j] && i!=j) add(i,j+n,);

    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 ;
}