P1027 Car的旅行路线

题目描述

又到暑假了，住在城市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

喜闻乐见的Code

#include<cmath>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
#define maxn 405
#define inf 99999999.0
struct City {
    int x,y;
    City(int x=0,int y=0):
        x(x),y(y) {}
}city[maxn];
int timee[maxn];
double dis[maxn][maxn];
inline int calc_dis(int x1,int y1,int x2,int y2) {
    return (x1-x2)*(x1-x2)+(y1-y2)*(y1-y2);
}
int main() {
    int T;
    scanf("%d",&T);
    while(T--) {
        int n,t,A,B;
        scanf("%d%d%d%d",&n,&t,&A,&B);
        for(int i=1;i<=n;i++) {
            int x1,y1,x2,y2,x3,y3,x4,y4;
            scanf("%d%d%d%d%d%d%d",&x1,&y1,&x2,&y2,&x3,&y3,&timee[i]);
            int d1=calc_dis(x2,y2,x3,y3);
            int d2=calc_dis(x1,y1,x3,y3);
            int d3=calc_dis(x1,y1,x2,y2);
            if(d1>d2&&d1>d3) x4=x3-(x1-x2),y4=y3-(y1-y2);
            else if(d2>d1&&d2>d3) x4=x3-(x2-x1),y4=y3-(y2-y1);
            else if(d3>d1&&d3>d2) x4=x1-(x3-x2),y4=y1-(y3-y2);
            city[4*i-3]=City(x1,y1);city[4*i-2]=City(x2,y2);
            city[4*i-1]=City(x3,y3);city[4*i]=City(x4,y4);
        }
        n*=4,A*=4,B*=4;
        for(int i=0;i<=n;i++)
            for(int j=0;j<=n;j++) dis[i][j]=inf*(i!=j);
        for(int i=1;i<=n;i++)
            for(int j=1;j<=n;j++)
                if((i-1)/4==(j-1)/4) {
                    double w=(double)calc_dis(city[i].x,city[i].y,city[j].x,city[j].y);
                    dis[i][j]=(double)sqrt(w)*timee[(i+3)/4];
                } else {
                    double w=(double)calc_dis(city[i].x,city[i].y,city[j].x,city[j].y);
                    dis[i][j]=(double)sqrt(w)*t;
                }
        for(int i=0;i<=3;i++) dis[0][A-i]=0.0;
        for(int k=0;k<=n;k++)
            for(int i=0;i<=n;i++)
                for(int j=0;j<=n;j++)
                    dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
        double ans=inf;
        for(int i=0;i<=3;i++) ans=min(ans,dis[0][B-i]);
        printf("%.1lf",ans);
    }
    return 0;
}