bzoj 1598: [Usaco2008 Mar]牛跑步

1598: [Usaco2008 Mar]牛跑步

Time Limit: 10 Sec   Memory Limit: 162 MB
Submit: 595   Solved: 351
[ 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

A*求最短路裸题

关于A*可以去看我的另一篇博客

里面有推荐两篇很好的博客，看完就懂啦

NOIP2012普及组第四题 文化之旅

这道题选取的估价函数是从终点到该点的最短路

是小于或等于实际值的

所以可以找到正确的解

#include<cstdio>
#include<cstring>
#include<queue>
const int N=1e3+7,M=1e4+7,K=1e2+7;
struct node
{
	int to,next,c;
}e[M],f[M];
int cnt,ans[K],first[N],dis[N],last[N];
int qu[N],visit[N];
int n,m,k;
inline void insert(int u,int v,int c)
{
	e[++cnt]=(node){v,first[u],c};first[u]=cnt;
	f[cnt]=(node){u,last[v],c};last[v]=cnt;
}
void spfa()
{
	memset(dis,127,sizeof(dis));
	int head=1,tail=2;
	qu[1]=1;visit[1]=1;dis[1]=0;
	while(head!=tail)
	{
		int rr=qu[head++];if(head==N)	head=0;
		for(int k=last[rr];k;k=f[k].next)
		{
			if(dis[f[k].to]>dis[rr]+f[k].c)
			{
				dis[f[k].to]=dis[rr]+f[k].c;
				if(!visit[f[k].to])
				{
					if(dis[f[k].to]<dis[qu[head]])
					{
						head--;if(head<0)	head=N-1;
						qu[head]=f[k].to;
					}
					else {
						qu[tail++]=f[k].to;if(tail==N)	tail=0;
					}
				}
			}
		}
		visit[rr]=0;
	}
}
struct edgt
{
	int dis,x;
	bool operator < (const edgt &y)
	const
	{
		return dis>y.dis;
	}
};
std::priority_queue<edgt>q;
void astar()
{
	memset(ans,-1,sizeof(ans));
	q.push((edgt){dis[n],n});
	cnt=0;
	while(!q.empty())
	{
		edgt t=q.top();q.pop();
		if(t.x==1)
		{
			ans[++cnt]=t.dis;
			if(cnt==k)	return;
		}
		for(int k=first[t.x];k;k=e[k].next)
		q.push((edgt){e[k].c+t.dis-dis[t.x]+dis[e[k].to],e[k].to});
	}
}
int main()
{
	scanf("%d %d %d",&n,&m,&k);
	for(int i=1;i<=m;i++)
	{
		int u,v,q;
		scanf("%d %d %d",&u,&v,&q);
		insert(u,v,q);
	}
	spfa();
	astar();
	for(int i=1;i<=k;i++)	printf("%d\n",ans[i]);
	return 0;
}