POJ 1734 图论入门第七题

     题目大意：在桑给巴尔岛的Adelton城镇上有一个旅游机构。它们决定在提供许多的其它旅游点之外，再向客人们提供旅游本镇的服务。 为了从提供的吸引服务中尽可能地获利，这个旅游机构接收了一个精明决定：在相同的起点与终点之间找出一最短路线。你的任务是编写一条程序来找类似的的一条路线。在这个镇上，有N个十字路口(编号1至N)，两个十字路口之间可以有多条道路连接，有M条双向行驶的道路（编号为1至M)。但没有一条道路从一个十字路口出发又回到同一个路口。每一条观光路线都是由一 些路组成的，这些道路序号是：y_1, ..., y_k,且k>2。第y_i(1<=i<=k-1)号路是连接第x_i号十字路口和第x_{i+1}号十字路口的；其中第y_k号路是连接第x_k号十字路口和第x_1号十字路口。而且所有的这些x_1,...,x_k分别代表不同路口的序号。在 某一条观光路线上所有道路的长度的和就是这条观光路线的总长度。换言之L(y_1)+L(y_2)+...+L(y_k) 的和 L(y_i)就是第y_i 号观光路线的长度。你的程序必须找出类似的一条路线：长度必须最小，或者说明在这个城镇上不存在这条观光路线。

     题目分析：判断并求解无向图中的最小环。求图中最小环的算法是扩展floyd。算法思想如下：

 

 

dist[i][j]=i,j之间的边长 dist:=cur; for k:=1 to n do begin for i:=1 to k-1 do for j:=i+1 to k-1 do answer:=min(answer,dist[i][j]+cur[i][k]+cur[k][j]); for i:=1 to n do for j:=1 to n do dist[i][j]:=min(dist[i][j],dist[i][k]+dist[k][j]); end;

 

 

     1734代码：

 

#include<cstdio> #include<cstring> #include<algorithm> using namespace std; #define read freopen("zx.in","r",stdin) #define write freopen("zx.out","w",stdout) const int INF=666666666;//注意INF可能溢出，慎用！ const int MAXN=100+10; int nnum,mnum,cnt,minn; int path[MAXN],fa[MAXN][MAXN],cur[MAXN][MAXN],dist[MAXN][MAXN]; void init() { for(int i=1;i<=MAXN;i++) for(int j=1;j<=MAXN;j++) { if(i==j) dist[i][j]=cur[i][j]=0; else dist[i][j]=cur[i][j]=INF; fa[i][j]=-1; } minn=INF; } void findroute(int i,int j)//求路径 { if(fa[i][j]!=-1) findroute(i,fa[i][j]); path[cnt++]=j; } void extend_floyd() { for(int k=1;k<=nnum;k++) { //扩展folyd，记录环的权值并找出最小环 for(int i=1;i<=k-1;i++) for(int j=i+1;j<=k-1;j++) { int temp=dist[i][j]+cur[i][k]+cur[k][j]; if(minn>temp) { minn=temp; cnt=0; findroute(i,j); path[cnt++]=k; } } //一般的floyd求两点最短路 for(int i=1;i<=nnum;i++) for(int j=1;j<=nnum;j++) if(dist[i][j]>dist[i][k]+dist[k][j]) { dist[i][j]=dist[i][k]+dist[k][j]; fa[i][j]=fa[k][j]; } } } int main() { read, write; while(scanf("%d%d",&nnum,&mnum)!=EOF) { init(); for(int i=1;i<=mnum;i++) { int a,b,cost; scanf("%d%d%d",&a,&b,&cost); dist[a][b]=dist[b][a]=min(dist[a][b],cost); cur[a][b]=cur[b][a]=min(cur[a][b],cost); fa[a][b]=a, fa[b][a]=b; } } extend_floyd(); if(minn==INF) printf("No solution./n"); else { for(int i=0;i<cnt;i++) printf(i==cnt-1 ? "%d/n":"%d ",path[i]); } return 0; }