POJ--3921[Destroying the bus stations] 最小费用流+拆点

题目大意：
一个有向图有N(N<=50)个点,M(M<=4000)条边。现在要求删掉最少的点,使得不存在从1号点到N号点的长度<=K(K<1000)的路径。当然,不能直接删掉1号点或N号点。输出最少需要删掉的点数(删掉一个点的时候会将与之相连的边也删掉)。

PS.不存在一条边直接连1-N。

 

思路：
最小费用流+拆点

 

构图：
(1):把[2..N-1]个点拆成i->i+N连一条容量为1费用为0的边<1,0>;//保证每个点被破坏一次
(2):把1和N拆点,连一条<INF,0>;//保证不会破坏
(3):如果存在一条边(u->v)怎连一条u+N->v容量为INF费用为1的边<INF,1>;//确定每条边的费用
(4):最后从1点到2*N点不停地找最短路增广，知道最短路>K为止

PS.也可以用深搜+枚举删点的方案：

CODE:

 

/*最小费用流+拆点*/
/*AC代码：32ms*/
#include <iostream>
#include <cstdio>
#include <memory.h>
#include <cstdio>
#include <cstdlib>
#include <queue>
#define MAXN 200
#define INF 1e8
#define min(a,b) (a<b?a:b)
using namespace std;
struct edge
{
	int u,v,w,c,next; 
}E[20000];
int head[MAXN],ecnt;
int N,M,K,scr,sink,vn;
bool vis[MAXN];
int dis[MAXN],pre[MAXN];

void Insert(int u,int v,int w,int c)
{
	E[ecnt].u=u;
	E[ecnt].v=v;
	E[ecnt].w=w;
	E[ecnt].c=c;
	E[ecnt].next=head[u];
	head[u]=ecnt++;
	E[ecnt].u=v;
	E[ecnt].v=u;
	E[ecnt].w=0;
	E[ecnt].c=-c;
	E[ecnt].next=head[v];
	head[v]=ecnt++;
}
void Init()
{
	int i,u,v;
	memset(head,-1,sizeof(head));ecnt=0;
	scr=1;sink=2*N;vn=sink;
	for(i=2;i<=N-1;i++)
		Insert(i,i+N,1,0);
	Insert(1,1+N,INF,0);
	Insert(N,N+N,INF,0);
	for(i=1;i<=M;i++)
	{
		scanf("%d%d",&u,&v);
		Insert(u+N,v,INF,1);
	}
}
queue<int>Q;
bool SPFA(int s,int t,int n)
{
	int i,u,v;
	int c;
	while(!Q.empty()) Q.pop();
	memset(vis,false,sizeof(vis));
	for(i=0;i<=n;i++)//找最长路
		dis[i]=INF;
	Q.push(s);
	pre[s]=-1;
	dis[s]=0;
	vis[s]=true;
	while(!Q.empty())
	{
		u=Q.front();Q.pop();
		vis[u]=false;
		for(i=head[u];i!=-1;i=E[i].next)
		{
			v=E[i].v;c=E[i].c;
			if(E[i].w>0&&dis[v]>dis[u]+c)
			{
				dis[v]=dis[u]+c;
				pre[v]=i;
				if(!vis[v])
				{
					vis[v]=true;
					Q.push(v);
				}
			}
		}
	}
	if(dis[t]!=INF) return true;
	return false;
}
void Solve()
{
	int i,ans=0;
	while(SPFA(scr,sink,vn))
	{
		
		if(dis[sink]>K) break;
		ans++;
		for(i=pre[sink];i!=-1;i=pre[E[i].u])//更新
		{
			E[i].w-=1;
			E[i^1].w+=1; 
		} 
	}
	printf("%d\n",ans);
}
int main()
{
	while(scanf("%d%d%d",&N,&M,&K)!=EOF)
	{
		if(N==0&&M==0&&K==0) break;
		Init();
		Solve();
	}
return 0;
}