做题打卡第七天

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