又到暑假了,住在城市 A 的 Car 想和朋友一起去城市旅游。
她知道每个城市都有 44 个飞机场,分别位于一个矩形的 44 个顶点上,同一个城市中两个机场之间有一条笔直的高速铁路,第 ii 个城市中高速铁路了的单位里程价格为 T_iT
i
,任意两个不同城市的机场之间均有航线,所有航线单位里程的价格均为 tt。
图例(从上而下)
机场
高速铁路
飞机航线
注意:图中并没有标出所有的铁路与航线。
那么 Car 应如何安排到城市B的路线才能尽可能的节省花费呢?她发现这并不是一个简单的问题,于是她来向你请教。
找出一条从城市 A 到 B 的旅游路线,出发和到达城市中的机场可以任意选取,要求总的花费最少。
输入格式
第一行为一个正整数 nn,表示有 nn 组测试数据。
每组的第一行有 44 个正整数 s,t,A,Bs,t,A,B。
SS 表示城市的个数,tt 表示飞机单位里程的价格,AA,BB 分别为城市A,B 的序号。
接下来有 SS 行,其中第 ii 行均有 77 个正整数
)分别是第 ii 个城市中任意 33 个机场的坐标,T_iT
i
为第 ii 个城市高速铁路单位里程的价格。
原题链接:https://www.luogu.com.cn/problem/solution/P1027
#include<bits/stdc++.h>
using namespace std;
const int maxn=100005,maxm=500005,inf=0x3f3f3f3f;
double e[1005][1005];
int wz[1005][11];
int n,tf;
double dist(int a,int b,int c,int d){
return sqrt((a-b)*(a-b)*1.0+1.0*(c-d)*(c-d));
}
double dist1(int a,int b,int c,int d){
return (a-b)*(a-b)*1.0+1.0*(c-d)*(c-d);
}
int main()
{
int t;
cin>>t;
while(t--)
{
int a,b;
scanf("%d%d%d%d",&n,&tf,&a,&b);
for(int i=1;i<=4*n;i++)
for(int j=1;j<=4*n;j++)
e[i][j]=inf;
for(int i=1;i<=n;i++)
{
for(int j=1;j<=6;j++)
scanf("%d",wz[i]+j);
for(int j=1;j<=5;j+=2)
{
wz[i][7]+=wz[i][j];
wz[i][8]+=wz[i][j+1];
}
scanf("%d",wz[i]+9);
double tp[3];
int tp2=inf,tp3;
tp[0]=dist1(wz[i][4],wz[i][6],wz[i][3],wz[i][5]);
tp[1]=dist1(wz[i][2],wz[i][6],wz[i][1],wz[i][5]);
tp[2]=dist1(wz[i][2],wz[i][4],wz[i][1],wz[i][3]);
if(tp[0]+tp[1]==tp[2]){
wz[i][7]-=2*wz[i][5];wz[i][8]-=2*wz[i][6];
}
else if(tp[1]+tp[2]==tp[0]){
wz[i][7]-=2*wz[i][1];wz[i][8]-=2*wz[i][2];
}
else if(tp[0]+tp[2]==tp[1]){
wz[i][7]-=2*wz[i][3];wz[i][8]-=2*wz[i][4];
}
}
for(int i=1;i<=n;i++)
{
for(int j=1;j<=4;j++)
for(int k=j;k<=4;k++)
{
int u=(i-1)*4+j,v=(i-1)*4+k;
double dis=dist(wz[i][j*2-1],wz[i][k*2-1],wz[i][j*2],wz[i][k*2]);
e[u][v]=e[v][u]=dis*wz[i][9];
}
}
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
if(i!=j){
for(int k=1;k<=4;k++)
for(int l=1;l<=4;l++)
{
int u=(i-1)*4+k,v=(j-1)*4+l;
double dis=dist(wz[i][k*2-1],wz[j][l*2-1],wz[i][k*2],wz[j][l*2]);
e[u][v]=dis*tf;
}
}
for(int k=1;k<=n*4;k++)
for(int i=1;i<=n*k;i++)
for(int j=1;j<=n*4;j++)
if(e[i][j]>e[i][k]+e[k][j])
e[i][j]=e[i][k]+e[k][j];
double ans=inf;
for(int i=(a-1)*4+1;i<=a*4;i++)
for(int j=(b-1)*4+1;j<=b*4;j++)
ans=min(ans,e[i][j]);
printf("%.1lf",ans);
}
return 0;
}