洛谷——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~

 #include <algorithm>
 #include <cstring>
 #include <cstdio>
 #include <cmath>

 #define y1 mmp

 using namespace std;

 const int M();
 const int N(<<);
 const double INF(1e15);

 int n,city_num,A,B;
 double dis[N][N],ans;
 int air_cost,road_cost[N];
 int x1,x2,x3,x4,y1,y2,y3,y4,d1,d2,d3;

 struct Node
 {
     int x,y;
     Node(int x=,int y=): x(x),y(y) {}
 }airport[M];

 int Get_dis(int a,int aa,int b,int bb)
 {
     return (a-b)*(a-b)+(aa-bb)*(aa-bb);
 }

 void init()
 {
     memset(dis,,sizeof(dis));
     memset(airport,,sizeof(airport));
     memset(road_cost,,sizeof(road_cost));
 }

 int main()
 {
     scanf("%d",&n);
     for(;n;n--)
     {
         init();
         scanf("%d%d%d%d",&city_num,&air_cost,&A,&B);
         for(int i=;i<=city_num;i++)
         {
             scanf("%d%d%d%d%d%d",&x1,&y1,&x2,&y2,&x3,&y3);
             scanf("%d",road_cost+i);
             d1=Get_dis(x1,y1,x2,y2);
             d2=Get_dis(x2,y2,x3,y3);
             d3=Get_dis(x3,y3,x1,y1);
             if(d1>d2&&d1>d3)
             {
                 x4=x1+x2-x3;
                 y4=y1+y2-y3;
             }
             else if(d2>d1&&d2>d3)
             {
                 x4=x2+x3-x1;
                 y4=y2+y3-y1;
             }
             else if(d3>d1&&d3>d2)
             {
                 x4=x3+x1-x2;
                 y4=y3+y1-y2;
             }
             airport[(i<<)-]=Node(x1,y1);
             airport[(i<<)-]=Node(x2,y2);
             airport[(i<<)-]=Node(x3,y3);
             airport[(i<<)-]=Node(x4,y4);
         }
         for(int i=;i<=(city_num<<);i++)
             for(int j=;j<=(city_num<<);j++)
                 if(i!=j) dis[i][j]=INF;
         for(int i=;i<=(city_num<<);i++)
             for(int j=;j<=(city_num<<);j++)
             {
                 double dist=sqrt(Get_dis(airport[i].x,airport[i].y,airport[j].x,airport[j].y));
                 if(i+>>==j+>>)
                     dis[i][j]=dist*(road_cost[i+>>]);
                 else
                     dis[i][j]=dist*air_cost;
             }
         dis[][(A<<)-]=;
         dis[][(A<<)-]=;
         dis[][(A<<)-]=;
         dis[][(A<<)-]=;
         for(int k=;k<=(city_num<<);k++)
           for(int i=;i<=(city_num<<);i++)
             for(int j=;j<=(city_num<<);j++)
                 dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
         ans=min(min(dis[][(B<<)-],dis[][(B<<)-]),min(dis[][(B<<)-],dis[][B<<]));
         printf("%.1lf\n",ans);
     }
     return ;
 }