洛谷P2764：最小路径覆盖问题

问题描述

      给定有向图G=(V,E)。设P 是G 的一个简单路（顶点不相交）的集合。如果V 中每个顶点恰好在P 的一条路上，则称P是G 的一个路径覆盖。P 中路径可以从V 的任何一个顶点开始，长度也是任意的，特别地，可以为0。G 的最小路径覆盖是G 的所含路径条数最少的路径覆盖。设计一个有效算法求一个有向无环图G 的最小路径覆盖。

 解法

      首先我们要明白一个定理：一张图中，路径数（点不重复）=点数-点之间匹配数（连边且不重复）。这个很好证明。因为当点之间没有边时，路径数=点数-0=点数；当有一条边（x，y）时，相当于把这两个点合在一起，所以路径数=点数-1；以此类推。

      所以我们要求最小路径覆盖，点数不变，那么就要使匹配数最大，就变成了求最大匹配的问题。

      同飞行员匹配问题一样，我们把一个点拆开两个点，然后有边的两个点之间连流量为1的边，begin到每个出点连一条流量为1的边，每个入点到end连一条流量为1的边，跑一遍最大流就可以知道最大匹配了，用总点数减一减就可以得出答案。

      关于方案输出的问题，用last记录一下前一个访问这个点的是哪个，用match记录一下当前这个点访问谁。具体可以到上面的链接中学习，每找到一个路径的结尾（match[i]==i）那么就递归访问last，输出即可。

代码如下：

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<queue>
using namespace std;

int n,m;
struct edge{int y,next,c;};
edge s[120010];
int first[310];
int len=1;
int begin,end;
queue<int> f;
int h[310];
int last[310],match[310];
bool tf[310];

int mmin(int x,int y)
{
	return x<y?x:y;
}

void ins(int x,int y,int c)
{
	len++;
	s[len].y=y;s[len].c=c;s[len].next=first[x];first[x]=len;
	len++;
	s[len].y=x;s[len].c=0;s[len].next=first[y];first[y]=len;
}

bool bfs()
{
	memset(h,0,sizeof(h));
	h[begin]=1;
	f.push(begin);
	while(!f.empty())
	{
		int x=f.front();
		for(int i=first[x];i!=0;i=s[i].next)
		{
			int y=s[i].y;
			if(h[y]==0 && s[i].c>0)
			{
				h[y]=h[x]+1;
				f.push(y);
			}
		}
		f.pop();
	}
	return h[end];
}

int dfs(int x,int t)
{	
	if(x==end) return t;
	int tot=0;
	for(int i=first[x];i!=0;i=s[i].next)
	{
		int y=s[i].y;
		if(tot==t) return t;
		if(h[y]==h[x]+1 && s[i].c>0)
		{
			int now=dfs(y,mmin(s[i].c,t-tot));
			tot+=now;s[i].c-=now;s[i^1].c+=now;
			if(now!=0 && x!=begin && y!=end)
			{
				match[x]=y-n;
				last[y-n]=x;
			}
		}
	}
	if(tot==0) h[x]=0;
	return tot;
}

int maxflow()
{
	int flow=0;
	while(bfs())
	{
		int dx=dfs(begin,1e9);
		while(dx!=0)
		{
			flow+=dx;
			dx=dfs(begin,1e9);
		}
	}
	return flow;
}

void out_put(int x)
{
	if(x==0) return ;
	if(last[x]!=x) out_put(last[x]);
	tf[x]=true;
	printf("%d ",x);
}

int main()
{
	memset(tf,false,sizeof(tf));
	scanf("%d %d",&n,&m);
	begin=0;end=2*n+1;
	for(int i=1;i<=n;i++)
		match[i]=last[i]=i;
	for(int i=1;i<=m;i++)
	{
		int x,y;
		scanf("%d %d",&x,&y);
		ins(x,y+n,1);
	}
	for(int i=1;i<=n;i++)
	{
		ins(begin,i,1);
		ins(i+n,end,1);
	}
	int p=maxflow();
	bool f[320];	
	for(int i=1;i<=n;i++)
	{
		if(match[i]==i && tf[i]==false)
		{
			out_put(i);
			printf("\n");
		}
	}
	printf("%d",n-p);
}