最佳路线【SGOI-7】

最新推荐文章于 2024-05-21 23:16:53 发布

cgj0128_1

最新推荐文章于 2024-05-21 23:16:53 发布

阅读量123

点赞数 1

分类专栏：图论文章标签：算法

本文链接：https://blog.csdn.net/cgj0128_1/article/details/129103698

版权

图论专栏收录该内容

9 篇文章 0 订阅

订阅专栏

题目：最佳路线【SGOI-7】（标号法）

每题提交次数限制：无限制

题目描述

　　塞内加尔是非洲的一个小国家，你也许很难在世界地图上找到它，甚至你有可能从未听说过它--它实在是个太小、太贫穷的国家了。可是，就是这个人口不足900万、全国仅有2个标准足球场地的小国，在2002韩日世界杯的非洲区预选赛中脱颖而出，取得了世界杯决赛圈的入场券（幸好，去年中国队也进入了世界杯决赛圈，不然可就丢脸了）。

　　在塞内加尔全国球迷欣喜若狂，世界足球行家大跌眼镜的同时，塞内加尔足协却发现自己面临着一个颇为尴尬的问题--说起来令人不可思议，由于打非洲区预选赛时四处征战，加上足协经营不力，现在足协的预算以几近赤字--也就是说，塞内加尔足协支付不起从本国乘飞机到达韩国参加世界杯的费用！经过三思，塞内加尔足协向非洲足联递交了一份《关于减免球队旅行费用》的申请；可是--众所周知的，非洲足联也是惨淡经营，幸好非洲足联秘书长神通广大，弄来了M张优惠乘机券：每张优惠券可以作用于一条航线，使全队通过此航线的费用减半；多张优惠券用于同一条线路，其效果叠加--即在一条航线上用两张优惠券，其费用降为原费用的1/4，依此类推。

【问题】

　　塞内加尔足球队要从塞内加尔国家机场出发，途经一些中转机场，最后要到达韩国釜山机场。为了合理地分配各张优惠券，使得所需费用最少，塞内加尔足协找到了你，请你编程解决这个问题。

输入输出格式

输入格式：

　　第1行有两个数N、M（0 < N <= 70，0 <= M <= 20 ）并用空格隔开，分别表示包括起点（塞内加尔国家机场）、终点（韩国釜山机场）的机场数，以及塞内加尔足协现有的优惠券数量。
　　从第2行到第N+1行起，每行有N个数，其中第I行的第J个数代表从机场I到机场J所需费用；为零的数代表两机场无航线。
　　假设起点标号为1，终点标号为N。

输出格式：

　　1行，仅有一个数（保留两位小数），代表所求得的从机场1到机场N的最小费用；
　　输入数据保证从塞内加尔机场可达釜山机场。

输入输出样例

输入样例#1：

5 2
0 0 80 96 0
70 0 72 54 0
18 0 0 99 82
72 18 71 0 0
69 0 0 70 0

输出样例#1：

81.00

提示信息

样例的一个具体方案：
　　1->3->5

1.算法分析

代表所求得的从机场1到机场N的最小费用

由1到n，即单源最短路

采用dijkstra算法（数据量很小，N、M（0 < N <= 70，0 <= M <= 20 ））

但是，题目中加入了优惠劵这个东西，因此要对狄克斯特拉算法进行改变

考虑在数组中加入多一个维度：dis[i][j]表示由源点到i点途中使用j张优惠券的最小花费

则对于任意dis[i][j],有

$dis[j][k]=min(dis[j][k],dis[kl[k]][l]+a[kl[k]][j][k-l]);$

其中，kl[k]表示在当前循环中，使用k张优惠券的使用最小花费的点，a[i][j][k]表示由i点到j点使用k张优惠券的初始最小花费：

$a[i][j][k]=a[i][j][k-1]/2;$

2.代码实现：

    for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
		    int x;
			cin>>x;
			if(x!=0)
			{
				a[i][j][0]=x;
				if(i==1) dis[j][0]=x;
				for(int k=1;k<=you;k++)
				{
					if(i==1) dis[j][k]=dis[j][k-1]/2;
					a[i][j][k]=a[i][j][k-1]/2;
				}
			}
		}
	}

1.初始化

    for(int i=1;i<=n-1;i++)
	{
		memset(kl,0,sizeof(kl));
		double minn[110];
		for(int j=0;j<=100;j++) minn[j]=1e9;
		for(int j=1;j<=n;j++)
		{
			for(int k=0;k<=you;k++)
			{
				if(vis[j][k]==0&&minn[k]>dis[j][k]) 
				{
					minn[k]=dis[j][k];
					kl[k]=j;
				}
			}
		}//找最小
		int p=0;
		for(int j=0;j<=you;j++)
		{
			if(kl[j]==0) //当前连通分量找不到没有遍历过的点
			{
                p=1;
				break;
				
			}
			vis[kl[j]][j]=1;//标记
		}
		if(p) break;
		for(int j=1;j<=n;j++)
		{
			for(int k=0;k<=you;k++)
			{
				for(int l=0;l<=k;l++)//l表示给dis[kl[k]]使用l张优惠券，k-l则是给a[kl[k]][j]用剩下的优惠券。
				{
					dis[j][k]=min(dis[j][k],dis[kl[k]][l]+a[kl[k]][j][k-l]);
				}
			}
		}
	}

2.代码实现

3.完整代码

#include<bits/stdc++.h>
using namespace std;
double a[110][110][30],you,vis[110][30];int kl[110];
double dis[110][30];
int main()
{
	int n,you;
	cin>>n>>you;
	for(int i=1;i<=100;i++)
	{
		for(int j=1;j<=100;j++)
		{
			for(int k=0;k<=27;k++)
			{
				a[i][j][k]=1e9;
				dis[i][k]=1e9; 
			}
		}
	}
	
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<=n;j++)
		{
		int x;
			cin>>x;
			if(x!=0)
			{
				a[i][j][0]=x;
				if(i==1)dis[j][0]=x;
				for(int k=1;k<=you;k++)
				{
					if(i==1)dis[j][k]=dis[j][k-1]/2;
					a[i][j][k]=a[i][j][k-1]/2;
				}
			}
		}
	}
	for(int i=0;i<=you;i++)
	{
		vis[1][i]=1;
	}
	for(int i=1;i<=n-1;i++)
	{
		memset(kl,0,sizeof(kl));
		double minn[110];
		for(int j=0;j<=100;j++) minn[j]=1e9;
		for(int j=1;j<=n;j++)
		{
			for(int k=0;k<=you;k++)
			{
				if(vis[j][k]==0&&minn[k]>dis[j][k]) 
				{
					minn[k]=dis[j][k];
					kl[k]=j;
				}
			}
		}
		int p=0;
		for(int j=0;j<=you;j++)
		{
			if(kl[j]==0) 
			{p=1;
				break;
				
			}
			vis[kl[j]][j]=1;
		}
		if(p) break;
		for(int j=1;j<=n;j++)
		{
			for(int k=0;k<=you;k++)
			{
				for(int l=0;l<=k;l++)
				{
					dis[j][k]=min(dis[j][k],dis[kl[k]][l]+a[kl[k]][j][k-l]);
				}
			}
		}
	}
	long long x=dis[n][you]*10000;
	if(x%100==50) printf("%.2lf",dis[n][you]-0.005);
	else printf("%.2lf",dis[n][you]);
	return 0;
}