POJ--1847[Tram] 最小费用流

题意：

给你N个点，告诉你每个点u可以到达的点v，但是可以到达的点又分两类：一类可以直接到达；另一类要经过一次“转换”后可以到达。现在问你点A到点B至少要“转换”几次能到达。如果A不能到B则输出“-1”，否则输出要“转换”的最少次数。

 

思路：

最小费用最大流

建图：
(1):每个点向他可以到的第一个点连<INF,0>,和他可以到达的其余点连<INF,1>
(2):scr向A连<1,0>
(3):B向sink连<1,0>
(4):然后从scr向sink找增广路
(5):找到则输出dis[sink];否则输出-1

 

CODE：

/*最小费用最大流*/
/*
建图：
(1):每个点向他可以到的第一个点连<INF,0>,和他可以到达的其余点连<INF,1>
(2):scr向A连<1,0>
(3):B向sink连<1,0>
(4):然后从scr向sink找增广路
(5):找到则输出dis[sink];否则输出-1
*/
/*AC代码：16ms*/
#include <iostream>
#include <cstdio>
#include <memory.h>
#include <algorithm>
#include <queue>
#define MAXN 105
#define INF 1e8
#define max(a,b) (a>b?a:b)
#define min(a,b) (a<b?a:b)
using namespace std;
struct edge
{
	int u,v,w,c,next; 
}E[200000];
int head[MAXN],ecnt;
int N,A,B,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,j,v,m;
	memset(head,-1,sizeof(head));ecnt=0;
	scr=0;sink=N+1;vn=sink+1;
	for(i=1;i<=N;i++)
	{
		scanf("%d",&m);
		for(j=1;j<=m;j++)
		{
			scanf("%d",&v);
			if(j==1) Insert(i,v,INF,0);
			else Insert(i,v,INF,1); 
		} 
	}
	Insert(scr,A,1,0);
	Insert(B,sink,1,0);
}
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()
{
	if(SPFA(scr,sink,vn))
		printf("%d\n",dis[sink]);
	else 
		printf("-1\n");
}
int main()
{
	while(scanf("%d%d%d",&N,&A,&B)!=EOF)
	{
		Init();
		Solve();
	}
return 0;
}