[bzoj1598][Usaco2008 Mar]牛跑步

1598: [Usaco2008 Mar]牛跑步

Time Limit: 10 Sec Memory Limit: 162 MB

Submit: 565 Solved: 333

[Submit][Status][Discuss]

Description

BESSIE准备用从牛棚跑到池塘的方法来锻炼. 但是因为她懒,她只准备沿着下坡的路跑到池塘, 然后走回牛棚. BESSIE也不想跑得太远,所以她想走最短的路经. 农场上一共有M (1 <= M <= 10,000)条路, 每条路连接两个用1..N(1 <= N <= 1000)标号的地点. 更方便的是,如果X>Y,则地点X的高度大于地点Y的高度. 地点N是BESSIE的牛棚;地点1是池塘. 很快, BESSIE厌倦了一直走同一条路.所以她想走不同的路,更明确地讲,她想找出K (1 <= K <= 100)条不同的路经.为了避免过度劳累,她想使这K条路经为最短的K条路经. 请帮助BESSIE找出这K条最短路经的长度.你的程序需要读入农场的地图, 一些从X_i到Y_i 的路经和它们的长度(X_i, Y_i, D_i). 所有(X_i, Y_i, D_i)满足(1 <= Y_i < X_i; Y_i < X_i <= N, 1 <= D_i <= 1,000,000).

Input

第1行: 3个数: N, M, 和K

第 2..M+1行: 第 i+1 行包含3个数 X_i, Y_i, 和 D_i, 表示一条下坡的路.

Output

第1..K行: 第i行包含第i最短路经的长度,或-1如果这样的路经不存在.如果多条路经有同样的长度,请注意将这些长度逐一列出.

Sample Input

5 8 7

5 4 1

5 3 1

5 2 1

5 1 1

4 3 4

3 1 1

3 2 1

2 1 1

Sample Output

1

2

2

3

6

7

-1

HINT

输出解释:

路经分别为(5-1), (5-3-1), (5-2-1), (5-3-2-1), (5-4-3-1),

(5-4-3-2-1).

Source

Gold

sol:

求k短路，A*

#include<iostream>
#include<algorithm>
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<cstring>
#include<string>
#include<queue>
using namespace std;
const int N=110000;
int n,m,dis
,fir
,nex
,go
,val
,h
,tot;
bool vis
;
inline void add(int x,int y,int z)
{
nex[++tot]=fir[x];fir[x]=tot;go[tot]=y;val[tot]=z;
}
inline int read()
{
char c;
int res,flag=0;
while((c=getchar())>'9'||c<'0') if(c=='-')flag=1;
res=c-'0';
while((c=getchar())>='0'&&c<='9') res=(res<<3)+(res<<1)+c-'0';
return flag?-res:res;
}
struct cc
{
int point,val;
friend bool operator <(cc a,cc b)
{
return a.val>=b.val;
}
};
priority_queue<cc> q;
struct dd
{
int x,y,z;
}a
;
inline void spfa()
{
int t=0,w=1,u,v,e;
h[1]=1;
memset(dis,127/3,sizeof(dis));
dis[1]=0;
while(t<w)
{
u=h[++t];
for(e=fir[u];v=go[e],e;e=nex[e])
if(dis[v]>dis[u]+val[e])
{
dis[v]=dis[u]+val[e];
if(!vis[v]) vis[v]=1,h[++w]=v;
}
vis[u]=0;
}
}
int ans
,tim
,k;
inline void astar()
{
int v,e;cc u;
if(dis
==707406378) return;
memset(ans,-1,sizeof(ans));
q.push((cc){n,dis
});
while(!q.empty())
{
u=q.top();
++tim[u.point];
q.pop();
if(u.point==1) ans[tim[u.point]]=u.val;
if(tim[u.point]>k) continue;
for(e=fir[u.point];v=go[e],e;e=nex[e])
q.push((cc){v,u.val-dis[u.point]+dis[v]+val[e]});
}
}
int main()
{
//  freopen("1598.txt","r",stdin);
//  freopen("1598.out","w",stdout);
n=read();m=read();k=read();
for(int i=1;i<=m;++i)
{
a[i].x=read();a[i].y=read();a[i].z=read();
add(a[i].y,a[i].x,a[i].z);
}
spfa();
memset(fir,0,sizeof(fir));tot=0;
for(int i=1;i<=m;++i) add(a[i].x,a[i].y,a[i].z);
astar();
for(int i=1;i<=k;++i)
printf("%d\n",ans[i]);
}