POJ 3463 Sightseeing

    

     这绝对是一道好题、一道深刻理解dijkstra算法的不二选择。

     题目是要求求出最短路和比最短路大1的次短路的总个数、卡时为2000ms。读完题目我感觉这题应该不难、本想练练手熟悉一下模版、后来……就悲剧了。

     dijkstra算法到底是怎么回事？我现在不敢说真懂了，做了至少80+的图论题目了，遇到的dijkstra算法的也不下20+、但这个题目真的让我很受益……在此膜拜一下dijkstra本人，不愧为21世纪计算机界的神牛！

     扯的有点远了……现在来分析如何解决这道题：（首先声明我的解题思路参考了这里：http://blog.csdn.net/leeeyupeng/archive/2010/08/06/5790928.aspx）我说说自己代码中三个二维数组的作用：dis[x][2]:dis[x][0]表示起点到x的最短路、dis[x][1]表示起点到x的次短路；arr[x][2]:arr[x][0]表示起点到x的最短路条数、arr[x][1]表示起点到x的次短路的条数；vis[x][2]对应于x和0、1功能为记录该点是否被访问！

     那么如何更新最小和次小路呢？显然我们容易想到下面的方法：

      1.if(x<最小)更新最小，次小；
      2.else if(x==最小)更新方法数；
      3.else if(x<次小)更新次小；
      4.else if(x==次小)更新方法数；

     最后说说dijkstra的循环部分、这也是本题的关键。为什么我们要循环2*nnum-1次？显然这道题中我们每一条边都需要考虑、这不是在求最短的一条，说白了是让你求出所有的可能组合，那么我们势必对每一种情况都需要遍历一次，虽然中间有重复。最短路里已知[start][0]已被标记为访问过，那么就只有nnum-1次遍历了，而次短路我们则需要遍历nnum次，这样两者加起来就为2*nnum-1次。这与我们平时使用优先队列+heap是一样的。只是更加细化了而已。

     那么、至此本题已被完美的解决。

 #include<cstdio> #include<cstring> #include<queue> #include<algorithm> using namespace std; #define read freopen("zx.in","r",stdin) #define write freopen("zx.out","w",stdout) const int INF=0x3f3f3f3f; const int MAXN=1000+10; typedef struct { int from,cost,next; }rode; rode zx[MAXN*MAXN]; int nnum,mnum,ans,p[MAXN],dis[MAXN][2],arr[MAXN][2]; bool vis[MAXN][2]; void init(int start) { memset(vis,false,sizeof(vis)); memset(arr,0,sizeof(arr)); for(int i=1;i<=nnum;i++) dis[i][0]=dis[i][1]=INF; dis[start][0]=0; arr[start][0]=1; } void read_G(int from,int to,int cost) { ans++; zx[ans].from=to; zx[ans].cost=cost; zx[ans].next=p[from]; p[from]=ans; } int dijsktra(int start,int end) { init(start); //我们需要进行2*nnum-1次遍历（because最短+次短总共需要遍历2*nnum-1个点--其中[start][0]已被访问） //以下代码我们同样可以使用优先队列+heap，但对于理解dijkstra没有什么好处！ for(int i=1;i<nnum*2;i++) { int mm=INF,xx,flag; for(int j=1;j<=nnum;j++) { if(!vis[j][0] && mm>dis[j][0]) { mm=dis[xx=j][0]; flag=0; } else if(!vis[j][1] && mm>dis[j][1]) { mm=dis[xx=j][1]; flag=1; } } if(mm==INF) break; vis[xx][flag]=true; int cur,cost; for(int j=p[xx];j!=-1;j=zx[j].next) { cur=zx[j].from; cost=zx[j].cost; if(mm+cost<dis[cur][0])//更新最小、次小 { dis[cur][1]=dis[cur][0];//先前的最小变成了次小 arr[cur][1]=arr[cur][0];//个数相应的传递 dis[cur][0]=mm+cost;//更新后的最小 arr[cur][0]=arr[xx][flag];//相应的个数 } else if(mm+cost==dis[cur][0])//只更新最小条数 { arr[cur][0]+=arr[xx][flag];//注意是连加 } else if(mm+cost<dis[cur][1])//只更新次小 { dis[cur][1]=mm+cost; arr[cur][1]=arr[xx][flag]; } else if(mm+cost==dis[cur][1])//只更新次小条数 { arr[cur][1]+=arr[xx][flag]; } } } //若最短距离和次短距离满足关系则返回两者个数之和 if(dis[end][0]+1==dis[end][1]) return arr[end][0]+arr[end][1]; return arr[end][0]; } int main() { read, write; int ncase,start,end; scanf("%d",&ncase); while(ncase--) { ans=0; memset(p,-1,sizeof(p));//注意位置 scanf("%d%d",&nnum,&mnum); for(int i=0;i<mnum;i++) { int from,to,cost; scanf("%d%d%d",&from,&to,&cost); read_G(from,to,cost); } scanf("%d%d",&start,&end); printf("%d/n",dijsktra(start,end)); } return 0; }