POJ 3463 最（次）短路条数

Posted 2020-10-07 liylho

Sightseeing

Time Limit: 2000MS		Memory Limit: 65536K
Total Submissions: 9497		Accepted: 3340

Description

Tour operator Your Personal Holiday organises guided bus trips across the Benelux. Every day the bus moves from one city S to another city F. On this way, the tourists in the bus can see the sights alongside the route travelled. Moreover, the bus makes a number of stops (zero or more) at some beautiful cities, where the tourists get out to see the local sights.

Different groups of tourists may have different preferences for the sights they want to see, and thus for the route to be taken from S to F. Therefore, Your Personal Holiday wants to offer its clients a choice from many different routes. As hotels have been booked in advance, the starting city S and the final city F, though, are fixed. Two routes from S to F are considered different if there is at least one road from a city A to a city B which is part of one route, but not of the other route.

There is a restriction on the routes that the tourists may choose from. To leave enough time for the sightseeing at the stops (and to avoid using too much fuel), the bus has to take a short route from S to F. It has to be either a route with minimal distance, or a route which is one distance unit longer than the minimal distance. Indeed, by allowing routes that are one distance unit longer, the tourists may have more choice than by restricting them to exactly the minimal routes. This enhances the impression of a personal holiday.

For example, for the above road map, there are two minimal routes from S = 1 to F = 5: 1 → 2 → 5 and 1 → 3 → 5, both of length 6. There is one route that is one distance unit longer: 1 → 3 → 4 → 5, of length 7.

Now, given a (partial) road map of the Benelux and two cities S and F, tour operator Your Personal Holiday likes to know how many different routes it can offer to its clients, under the above restriction on the route length.

Input

The first line of the input file contains a single number: the number of test cases to follow. Each test case has the following format:

• One line with two integers N and M, separated by a single space, with 2 ≤ N ≤ 1,000 and 1 ≤ M ≤ 10, 000: the number of cities and the number of roads in the road map.

• M lines, each with three integers A, B and L, separated by single spaces, with 1 ≤ A, B ≤ N, A ≠ B and 1 ≤ L ≤ 1,000, describing a road from city A to city B with length L.

  The roads are unidirectional. Hence, if there is a road from A to B, then there is not necessarily also a road from B to A. There may be different roads from a city A to a city B.

• One line with two integers S and F, separated by a single space, with 1 ≤ S, F ≤ N and S ≠ F: the starting city and the final city of the route.

  There will be at least one route from S to F.

Output

For every test case in the input file, the output should contain a single number, on a single line: the number of routes of minimal length or one distance unit longer. Test cases are such, that this number is at most 10⁹ = 1,000,000,000.

Sample Input

2
5 8
1 2 3
1 3 2
1 4 5
2 3 1
2 5 3
3 4 2
3 5 4
4 5 3
1 5
5 6
2 3 1
3 2 1
3 1 10
4 5 2
5 2 7
5 2 7
4 1

Sample Output

3
2

Hint

The first test case above corresponds to the picture in the problem description.

 

题意：旅行团每天固定的从S地出发到达T地，为了省油要求尽量走最短路径或比最短路径长1单位距离的路径，求满足条件的路径条数

算法：最短路和次短路。Dijkstra算法。采用邻接表建图。 

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是一样的。只是更加细化了而已。

代码：

  1 //#include<bits/stdc++.h>
  2 //#include<regex>
  3 #include <vector>
  4 #include <cstdio>
  5 #include "cstring"
  6 #include <algorithm>
  7 #define db double
  8 #define ll long long
  9 #define vec vector<ll>
 10 #define Mt  vector<vec>
 11 #define ci(x) scanf("%d",&x)
 12 #define cd(x) scanf("%lf",&x)
 13 #define cl(x) scanf("%lld",&x)
 14 #define pi(x) printf("%d\n",x)
 15 #define pd(x) printf("%f\n",x)
 16 #define pl(x) printf("%lld\n",x)
 17 #define MP make_pair
 18 #define PB push_back
 19 #define fr(i,a,b) for(int i=a;i<=b;i++)
 20 using namespace std;
 21 const int N=1e6+5;
 22 const int mod=1e9+7;
 23 const int MOD=mod-1;
 24 const db  eps=1e-18;
 25 //const db  pi = acos(-1.0);
 26 const int inf = 0x3f3f3f3f;
 27 const ll  INF = 0x3f3f3f3f3f3f3f3f;
 28 int h[N],d[N][2],c[N][2];
 29 bool vis[N][2];
 30 struct P
 31 {
 32     int v,w,nx;
 33     P(int x,int y,int z):v(x),w(y),nx(z){};
 34     P(){};
 35 };
 36 vector<P> e;
 37 int cnt,n,m,t;
 38 void add(int x,int y,int z)
 39 {
 40 
 41     e.push_back(P(y,z,h[x]));
 42     h[x]=cnt++;
 43 }
 44 void dij(int s,int g)
 45 {
 46     memset(vis,0, sizeof(vis));
 47     memset(c,0, sizeof(c));
 48     for(int i=0;i<=n;i++){
 49         d[i][0]=inf;
 50         d[i][1]=inf;
 51     }
 52     int k,ok;
 53     d[s][0]=0,c[s][0]=1;
 54 
 55     for(int ii=0;ii<2*n;ii++)
 56     {
 57         int ans=inf;
 58         for(int i=1;i<=n;i++){
 59             if(!vis[i][0]&&d[i][0]<ans){//取最近的点
 60                 k=i;
 61                 ok=0;
 62                 ans=d[i][0];
 63             }
 64             else if(!vis[i][1]&&d[i][1]<ans){
 65                 k=i;
 66                 ok=1;
 67                 ans=d[i][1];
 68             }
 69         }
 70         if(ans==inf) break;
 71         vis[k][ok]=1;
 72         for(int j=h[k];j!=-1;j=e[j].nx){
 73             int v=e[j].v,w=e[j].w;
 74             if(d[v][0]>ans+w){//更新最短路
 75                 d[v][1]=d[v][0];
 76                 c[v][1]=c[v][0];
 77                 d[v][0]=ans+w;
 78                 c[v][0]=c[k][ok];
 79             }
 80             else if(d[v][0]==ans+w){//更新最短路条数
 81                 c[v][0]+=c[k][ok];
 82             }
 83             else if(d[v][1]>ans+w)//更新次短路
 84                 d[v][1]=ans+w,c[v][1]=c[k][ok];
 85             else if(d[v][1]==ans+w)//次短路条数
 86                 c[v][1]+=c[k][ok];
 87         }
 88     }
 89     if(d[g][1]==d[g][0]+1)
 90         c[g][0]+=c[g][1];
 91 
 92     pi(c[g][0]);
 93 
 94 }
 95 int main()
 96 {
 97     ci(t);
 98     while(t--)
 99     {
100         ci(n),ci(m);
101         e.clear();
102         cnt=0;
103         memset(h,-1, sizeof(h));
104         for(int i=0;i<m;i++){
105             int x,y,z;
106             ci(x),ci(y),ci(z);
107             add(x,y,z);
108         }
109         int s,g;
110         ci(s),ci(g);
111         dij(s,g);
112     }
113     return 0;
114 
115 }