最短路和次短路问题,dijkstra算法

  1 /* 
  2  *题目大意： 
  3  *在一个有向图中,求从s到t两个点之间的最短路和比最短路长1的次短路的条数之和; 
  4  * 
  5  *算法思想： 
  6  *用A*求第K短路,目测会超时,直接在dijkstra算法上求次短路; 
  7  *将dist数组开成二维的,即dist[v][2],第二维分别用于记录最短路和次短路; 
  8  *再用一个cnt二维数组分别记录最短路和次短路的条数; 
  9  *每次更新路径的条数时,不能直接加1,,应该加上cnt[u][k],k为次短路径或者最短路径的标记; 
 10  *图有重边,不能用邻接矩阵存储; 
 11  *不知道为什么,题目上说的是N and M, separated by a single space, with 2≤N≤1000 and 1 ≤ M ≤ 10000; 
 12  *而我的代码硬是把N开成1W了才过,求解释,RE了无数次,擦; 
 13 **/ 
 14 #include<iostream>  
 15 #include<cstdio>  
 16 #include<cstring>  
 17 #include<string>  
 18 #include<algorithm>  
 19 using namespace std;  
 20    
 21 const int N=11111;  
 22 const int M=111111;  
 23 const int INF=0xffffff;  
 24    
 25 struct node  
 26 {  
 27     int to;  
 28     int w;  
 29     int next;  
 30 };  
 31    
 32 node edge[N];  
 33 int head[N];  
 34    
 35 int dist[N][2],cnt[N][2];  
 36 bool vis[N][2];  
 37 int n,m,s,t,edges;  
 38    
 39 void addedge(int u,int v,int w)  
 40 {  
 41     edge[edges].w=w;  
 42     edge[edges].to=v;  
 43     edge[edges].next=head[u];  
 44     head[u]=edges++;  
 45 }  
 46    
 47 int dijkstra()  
 48 {  
 49     int k;  
 50     for(int i=0; i<=n; i++)  
 51     {  
 52         dist[i][0]=dist[i][1]=INF;  
 53         vis[i][0]=vis[i][1]=0;  
 54         cnt[i][0]=cnt[i][1]=0;  
 55     }  
 56     cnt[s][0]=1,dist[s][0]=0;  
 57    
 58     for(int i=1; i<=n*2; i++)  
 59     {  
 60         int u=-1;  
 61         int min_dist=INF;  
 62         for(int j=1; j<=n; j++)  
 63             for(int flag=0; flag<2; flag++)  
 64                 if(!vis[j][flag]&&dist[j][flag]<min_dist)  
 65                 {  
 66                     min_dist=dist[j][flag];  
 67                     u=j;  
 68                     k=flag;  
 69                 }  
 70         if(u==-1)  
 71             break;  
 72         vis[u][k]=true;  
 73         for(int e=head[u]; e!=-1; e=edge[e].next)  
 74         {  
 75             int j=edge[e].to;  
 76             int tmp=dist[u][k]+edge[e].w;  
 77    
 78             if(tmp<dist[j][0])//tmp小于最短路径长:  
 79             {  
 80                 dist[j][1]=dist[j][0];//次短路径长  
 81                 cnt[j][1]=cnt[j][0];//次短路径计数  
 82                 dist[j][0]=tmp;//最短路径长  
 83                 cnt[j][0]=cnt[u][k];//最短路径计数  
 84             }  
 85    
 86             else if(tmp==dist[j][0])//tmp等于最短路径长：  
 87             {  
 88                 cnt[j][0]+=cnt[u][k];//最短路径计数  
 89             }  
 90    
 91             else if(tmp<dist[j][1])//tmp大于最短路径长且小于次短路径长：  
 92             {  
 93                 dist[j][1]=tmp;//次短路径长  
 94                 cnt[j][1]=cnt[u][k];//次短路径计数  
 95             }  
 96    
 97             else if(tmp==dist[j][1])//tmp等于次短路径长：  
 98             {  
 99                 cnt[j][1]+=cnt[u][k];//次短路径计数  
100             }  
101         }  
102     }  
103    
104     int res=cnt[t][0];  
105     if(dist[t][0]+1==dist[t][1])//判断最短路和次短路是否相差1  
106         res+=cnt[t][1];  
107     return res;  
108 }  
109    
110 int main()  
111 {  
112     //freopen("C:\\Users\\Administrator\\Desktop\\kd.txt","r",stdin);  
113     int tcase;  
114     scanf("%d",&tcase);  
115     while(tcase--)  
116     {  
117         edges=0;  
118         scanf("%d%d",&n,&m);  
119         memset(head,-1,sizeof(head));  
120         int u,v,w;  
121         for(int i=0; i<m; i++)  
122         {  
123             scanf("%d%d%d",&u,&v,&w);  
124             addedge(u,v,w);  
125         }  
126         scanf("%d%d",&s,&t);  
127         printf("%d\n",dijkstra());  
128     }  
129     return 0;  
130 }