洛谷——P1027 Car的旅行路线
https://www.luogu.org/problem/show?pid=1027#sub
题目描述
又到暑假了,住在城市A的Car想和朋友一起去城市B旅游。她知道每个城市都有四个飞机场,分别位于一个矩形的四个顶点上,同一个城市中两个机场之间有一条笔直的高速铁路,第I个城市中高速铁路了的单位里程价格为Ti,任意两个不同城市的机场之间均有航线,所有航线单位里程的价格均为t。
图例(从上而下)
机场 高速铁路
飞机航线
注意:图中并没有
标出所有的铁路与航线。
那么Car应如何安排到城市B的路线才能尽可能的节省花费呢?她发现这并不是一个简单的问题,于是她来向你请教。
找出一条从城市A到B的旅游路线,出发和到达城市中的机场可以任意选取,要求总的花费最少。
输入输出格式
输入格式:
第一行为一个正整数n(0<=n<=10),表示有n组测试数据。
每组的第一行有四个正整数s,t,A,B。
S(0<S<=100)表示城市的个数,t表示飞机单位里程的价格,A,B分别为城市A,B的序号,(1<=A,B<=S)。
接下来有S行,其中第I行均有7个正整数xi1,yi1,xi2,yi2,xi3,yi3,Ti,这当中的(xi1,yi1),(xi2,yi2),(xi3,yi3)分别是第I个城市中任意三个机场的坐标,T I为第I个城市高速铁路单位里程的价格。
输出格式:
共有n行,每行一个数据对应测试数据。 保留一位小数
输入输出样例
输入样例#1:
1 3 10 1 3 1 1 1 3 3 1 30 2 5 7 4 5 2 1 8 6 8 8 11 6 3
输出样例#1:
47.5
恶心——————————太恶心了
首先用几何算出第四个点
用dis[0][A] A的四个点作为起点,dis[0][B]的四个点做终点
剩下的就是Floyd O(∩_∩)O~
1 #include <algorithm> 2 #include <cstring> 3 #include <cstdio> 4 #include <cmath> 5 6 #define y1 mmp 7 8 using namespace std; 9 10 const int M(2333); 11 const int N(110<<2); 12 const double INF(1e15); 13 14 int n,city_num,A,B; 15 double dis[N][N],ans; 16 int air_cost,road_cost[N]; 17 int x1,x2,x3,x4,y1,y2,y3,y4,d1,d2,d3; 18 19 struct Node 20 { 21 int x,y; 22 Node(int x=0,int y=0): x(x),y(y) {} 23 }airport[M]; 24 25 int Get_dis(int a,int aa,int b,int bb) 26 { 27 return (a-b)*(a-b)+(aa-bb)*(aa-bb); 28 } 29 30 void init() 31 { 32 memset(dis,0,sizeof(dis)); 33 memset(airport,0,sizeof(airport)); 34 memset(road_cost,0,sizeof(road_cost)); 35 } 36 37 int main() 38 { 39 scanf("%d",&n); 40 for(;n;n--) 41 { 42 init(); 43 scanf("%d%d%d%d",&city_num,&air_cost,&A,&B); 44 for(int i=1;i<=city_num;i++) 45 { 46 scanf("%d%d%d%d%d%d",&x1,&y1,&x2,&y2,&x3,&y3); 47 scanf("%d",road_cost+i); 48 d1=Get_dis(x1,y1,x2,y2); 49 d2=Get_dis(x2,y2,x3,y3); 50 d3=Get_dis(x3,y3,x1,y1); 51 if(d1>d2&&d1>d3) 52 { 53 x4=x1+x2-x3; 54 y4=y1+y2-y3; 55 } 56 else if(d2>d1&&d2>d3) 57 { 58 x4=x2+x3-x1; 59 y4=y2+y3-y1; 60 } 61 else if(d3>d1&&d3>d2) 62 { 63 x4=x3+x1-x2; 64 y4=y3+y1-y2; 65 } 66 airport[(i<<2)-3]=Node(x1,y1); 67 airport[(i<<2)-2]=Node(x2,y2); 68 airport[(i<<2)-1]=Node(x3,y3); 69 airport[(i<<2)-0]=Node(x4,y4); 70 } 71 for(int i=0;i<=(city_num<<2);i++) 72 for(int j=0;j<=(city_num<<2);j++) 73 if(i!=j) dis[i][j]=INF; 74 for(int i=0;i<=(city_num<<2);i++) 75 for(int j=0;j<=(city_num<<2);j++) 76 { 77 double dist=sqrt(Get_dis(airport[i].x,airport[i].y,airport[j].x,airport[j].y)); 78 if(i+3>>2==j+3>>2) 79 dis[i][j]=dist*(road_cost[i+3>>2]); 80 else 81 dis[i][j]=dist*air_cost; 82 } 83 dis[0][(A<<2)-3]=0; 84 dis[0][(A<<2)-2]=0; 85 dis[0][(A<<2)-1]=0; 86 dis[0][(A<<2)-0]=0; 87 for(int k=0;k<=(city_num<<2);k++) 88 for(int i=0;i<=(city_num<<2);i++) 89 for(int j=0;j<=(city_num<<2);j++) 90 dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]); 91 ans=min(min(dis[0][(B<<2)-3],dis[0][(B<<2)-2]),min(dis[0][(B<<2)-1],dis[0][B<<2])); 92 printf("%.1lf\n",ans); 93 } 94 return 0; 95 }
——每当你想要放弃的时候,就想想是为了什么才一路坚持到现在。