1011 - 网络流 - 祭祀river（BZOJ 1143）

传送门

 

题意

输出最长反链

 

分析

根据某D定理：最长反链 = 最小链覆盖，这道题就分分钟了

突然发现这类题，翻来覆去，换了件外套还是同一个板子

难点就在于你能不能把外套给别人扒下来，（这么说好像有点粗俗……不过没关系能懂就行）

话说原题其实还要求输出一种方案的，但我太菜了……还是先把BZOJ上的交了再说吧

 

代码

对了，不要觉得这个代码眼熟，你没看错，就是上一道题读入略作改动而已

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<cstdlib>
#include<queue>
#define in read()
#define M 2000000
#define inf 0x3f3f3f3f
#define N 21009
using namespace std;
inline int read(){
	char ch;int f=1,res=0;
	while((ch=getchar())<'0'||ch>'9') if(ch=='-') f=-1;
	while(ch>='0'&&ch<='9'){
		res=(res<<3)+(res<<1)+ch-'0';
		ch=getchar();
	}
	return f==1?res:-res;
}

int S,T;
int nxt[M],to[M],cap[M],head[N],cur[N],cnt=1,lev[N];
int f[505][505];
void add(int x,int y,int z){
	nxt[++cnt]=head[x];head[x]=cnt;to[cnt]=y;cap[cnt]=z;
	nxt[++cnt]=head[y];head[y]=cnt;to[cnt]=x;cap[cnt]=0;
}
bool bfs(){
	for(int i=S;i<=T;++i){	cur[i]=head[i];lev[i]=-1;} 
	queue<int > q;
	q.push(S);lev[S]=0;
	while(!q.empty()){
		int u=q.front();q.pop();
		for(int e=head[u];e;e=nxt[e]){
			int v=to[e];
			if(lev[v]!=-1||cap[e]<=0) continue;
			lev[v]=lev[u]+1;
			if(v==T) return true;
			q.push(v);
 		}
	}
	return false;
}
int dinic(int u,int flow){
	if(u==T)	return flow;
	int delta,res=0;
	for(int &e=cur[u];e;e=nxt[e]){
		int v=to[e];
		if(lev[v]>lev[u]&&cap[e]>0){
			delta=dinic(v,min(cap[e],flow-res));
			if(delta){
				res+=delta;cap[e]-=delta;
				cap[e^1]+=delta;if(res==flow) return flow;
			}
		}
	}
	return res;
}
int tt,m,n;
int main(){
	
		n=in;m=in;
		cnt=1;memset(head,0,sizeof(head));
		memset(f,0,sizeof(f));
		S=0;T=2*n+1;
		int i,j,k;
		for(i=1;i<=m;++i){
			j=in;k=in;
			f[j][k]=1;
		}
		for(k=1;k<=n;++k)
		for(i=1;i<=n;++i)
			for(j=1;j<=n;++j)
				if(f[i][j]==1||(f[i][k]==1&&f[k][j]==1))
					add(i,j+n,1);
		for(i=1;i<=n;++i) add(S,i,1),add(i+n,T,1);
		int maxflow=0;
		while(bfs()) maxflow+=dinic(S,inf);
		printf("%d\n",n-maxflow);
	
	return 0;
}