「网络流 24 题」[3] 最小路径覆盖

题意

求一张有向图的最小不可相交路径覆盖

分析

似乎之前做某一道TC题的时候就顺便记了一下最小不可相交路径覆盖怎么做...

http://note.youdao.com/noteshare?id=624eddad1e2ded71c43939f1a4df0a65&sub=8090A05D94F94F42B2A4794C1153F4DC

然后这里只是把二分图匹配换成网络流了。

但是建图方面还是要重新意思意思说一下的。

code

#include<bits/stdc++.h>
#define M 605
#define E 12005
#define inf 100000000
using namespace std;
void read(int &x){
	x=0; char c=getchar();
	for (;c<48;c=getchar());
	for (;c>47;c=getchar())x=(x<<1)+(x<<3)+(c^48);
}
struct ed{
	int x,cap,nx;
}e[E];
int nx[M],n,ecnt,to[M];
bool mark[M];
void add(int x,int y,int cap){
	e[ecnt]=(ed){y,cap,nx[x]};
	nx[x]=ecnt++;
	e[ecnt]=(ed){x,0,nx[y]};
	nx[y]=ecnt++;
}
struct Dinic{
	int Q[M],nnx[M],level[M],s,t;
	bool bfs(int x){
		int l=0,r=0;
		memset(level,0,sizeof(level));
		Q[r++]=x; level[x]=1;
		for (;l<r;){
			x=Q[l++];
			for (int i=nx[x];~i;i=e[i].nx)if (e[i].cap>0&&!level[e[i].x]){
				level[e[i].x]=level[x]+1;
				Q[r++]=e[i].x;
				if (e[i].x==t)return 1;
			}
		}
		return level[t]>0;
	}
	int dfs(int x,int f){
		if (x==t)return f;
		int sum=0;
		for (int &i=nnx[x];~i;i=e[i].nx)if (e[i].cap>0&&level[e[i].x]==level[x]+1){
			int d=dfs(e[i].x,min(e[i].cap,f-sum));
			if (d){
				to[x]=e[i].x;
				if (e[i].x>n)mark[e[i].x-n]=1;
			}
			e[i].cap-=d;
			e[i^1].cap+=d;
			sum+=d;
			if (sum==f)return f;
		}
		if (!sum)level[x]=0;
		return sum;
	}
	int solve(int S,int T){
		s=S; t=T;
		int res=0;
		for (;bfs(s);){
			memcpy(nnx,nx,sizeof(nx));
			res+=dfs(s,inf);
		}
		return res;
	}
}dinic;

void pt(int x){
	printf("%d",x);
	for (;to[x];){
		x=to[x]-n;
		printf(" %d",x);
	}
	printf("\n");
}
int main(){
	int m,i,x,y;
	memset(nx,-1,sizeof(nx));
	read(n); read(m);
	for (i=1;i<=n;i++)add(0,i,1),add(n+i,2*n+1,1);
	for (i=1;i<=m;i++){
		read(x); read(y);
		add(x,n+y,inf);
	}
	int res=dinic.solve(0,2*n+1);
	for (i=1;i<=n;i++)if (!mark[i])pt(i);
	printf("%d\n",n-res);
	return 0;
}