SGU 185 Two shortest

Two shortest

Time Limit: 500ms

Memory Limit: 4096KB

This problem will be judged on SGU. Original ID: 185
64-bit integer IO format: %I64d      Java class name: Solution

 

Yesterday Vasya and Petya quarreled badly, and now they don't want to see each other on their way to school. The problem is that they live in one and the same house, leave the house at the same time and go at the same speed by the shortest road. Neither of them wants to change their principles, that is why they want to find two separate shortest routes, which won't make them go along one road, but still they can meet at any junction. They ask you to help them. They number all the junctions with numbers from 1 to N (home and school are also considered as junctions). So their house has the number 1 and the school has the number N, each road connects two junctions exactly, and there cannot be several roads between any two junctions.

 

Input

The first line contains two integer numbers N and M (2<=N<=400), where M is the number of roads Petya and Vasya noticed. Each of the following M lines contains 3 integers: X, Y and L (1<=X, Y<=N, 1<=L<=10000), where X and Y - numbers of junctions, connected by the road and L is the length of the road.

 

Output

Write to the first line numbers of the junctions in the way they passed them on the first route. Write to the second line numbers of the junctions in the way they passed them on the second route. If it is impossible to help guys, then output "No solution".

 

Sample Input

Sample test(s)

Input

 

 

6 8
1 2 1
3 2 1
3 4 1
1 3 2
4 2 2
4 5 1
5 6 1
4 6 2

 

 

Output

 

 

1 3 4 5 6
1 2 4 6

 

 

Source

ACM International Collegiate Programming Contest 2003-2004 North-Eastern European Region, Southern Subregion

 

解题：妈拉个巴子，写了三遍才过。。

 

 #include <bits/stdc++.h>
 #define pii pair<int,int>
 using namespace std;
 const int INF = 0x3f3f3f3f;
 const int maxn = ;
 struct arc {
     int to,w,next;
     arc(int x = ,int y = ,int z = -) {
         to = x;
         w = y;
         next = z;
     }
 } e[maxn*maxn*];
 int hd[maxn],hd2[maxn],d[maxn],cur[maxn],tot,n,m;
 void add(int *head,int u,int v,int w) {
     e[tot] = arc(v,w,head[u]);
     head[u] = tot++;
 }
 void dijkstra() {
     priority_queue<pii,vector<pii>,greater<pii > >q;
     memset(d,0x3f,sizeof d);
     bool done[maxn] = {false};
     q.push(pii(d[] = ,));
     while(!q.empty()) {
         int u = q.top().second;
         q.pop();
         if(done[u]) continue;
         done[u] = true;
         for(int i = hd[u]; ~i; i = e[i].next) {
             if(d[e[i].to] > d[u] + e[i].w) {
                 d[e[i].to] = d[u] + e[i].w;
                 q.push(pii(d[e[i].to],e[i].to));
             }
         }
     }
 }
 bool bfs() {
     queue<int>q;
     memset(d,-,sizeof d);
     d[] = ;
     q.push();
     while(!q.empty()) {
         int u = q.front();
         q.pop();
         for(int i = hd2[u]; ~i; i = e[i].next) {
             if(e[i].w && d[e[i].to] == -) {
                 d[e[i].to] = d[u] + ;
                 q.push(e[i].to);
             }
         }
     }
     return d[n] > -;
 }
 int dfs(int u,int low) {
     if(u == n) return low;
     int tmp = ,a;
     for(int &i = cur[u]; ~i; i = e[i].next) {
         if(e[i].w && d[e[i].to] == d[u]+&&(a=dfs(e[i].to,min(low,e[i].w)))) {
             e[i].w -= a;
             e[i^].w += a;
             tmp += a;
             low -= a;
             if(!low) break;
         }
     }
     if(!tmp) d[u] = -;
     return tmp;
 }
 int dinic() {
     int ret = ;
     while(bfs()) {
         memcpy(cur,hd2,sizeof cur);
         ret += dfs(,INF);
     }
     return ret;
 }
 void solve(int u) {
     if(u == ) printf("%d",u);
     else printf(" %d",u);
     if(u == n) {
         putchar('\n');
         return;
     }
     for(int i = hd2[u]; ~i; i = e[i].next) {
         if((~i&) && !e[i].w) {
             e[i].w = ;
             solve(e[i].to);
             break;
         }
     }
 }
 int main() {
     int u,v,w;
     while(~scanf("%d%d",&n,&m)) {
         memset(hd,-,sizeof hd);
         memset(hd2,-,sizeof hd2);
         for(int i = tot = ; i < m; ++i) {
             scanf("%d%d%d",&u,&v,&w);
             add(hd,u,v,w);
             add(hd,v,u,w);
         }
         dijkstra();
         for(int i = ; i <= n; ++i) {
             for(int j = hd[i]; ~j; j = e[j].next) {
                 if(d[e[j].to] == d[i] + e[j].w) {
                     add(hd2,i,e[j].to,);
                     add(hd2,e[j].to,i,);
                 }
             }
         }
         if(dinic() >= ) {
             solve();
             solve();
         } else puts("No solution");
     }
     return ;
 }