题目描述 Description
又到暑假了,住在城市A的Car想和朋友一起去城市B旅游。她知道每个城市都有四个飞机场,分别位于一个矩形的四个顶点上,同一个城市中两个机场之间有一条笔直的高速铁路,第I个城市中高速铁路了的单位里程价格为Ti,任意两个不同城市的机场之间均有航线,所有航线单位里程的价格均为t。
那么Car应如何安排到城市B的路线才能尽可能的节省花费呢?她发现这并不是一个简单的问题,于是她来向你请教。
任务
找出一条从城市A到B的旅游路线,出发和到达城市中的机场可以任意选取,要求总的花费最少。
输入描述 Input Description
第一行为一个正整数n(0<=n<=10),表示有n组测试数据。
每组的第一行有四个正整数s,t,A,B。
S(0
输出描述 Output Description
共有n行,每行一个数据对应测试数据。
样例输入 Sample Input
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
样例输出 Sample Output
47.5
思路
预处理飞了飞了。。。
简直就是考代码能力。
数据很弱,floyd都可以过。
代码
#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstring>
#include<algorithm>
using namespace std;
double dis[405][405];
long long caldis(int x,int y)
{
if (x>y) swap(y,x);
return (y-x)*(y-x);
}
struct point
{
int x,y;
}point[405];
int tot=1;
void Floyd(int n)
{
for (int k=1;k<=n;k++)
for (int i=1;i<=n;i++)
for (int j=1;j<=n;j++)
dis[i][j]=min(dis[i][j],dis[i][k]+dis[k][j]);
}
int main()
{
int n;
scanf("%d",&n);
while(n--)
{
int s,t,a,b;
scanf("%d%d%d%d",&s,&t,&a,&b);
int x1,y1,x2,y2,x3,y3,x4,y4,c;
long long pf12,pf13,pf14,pf23,pf24,pf34,p1,p2,p3,p4;
for (int i=1;i<=s;i++)
{
scanf("%d%d%d%d%d%d%d",&x1,&y1,&x2,&y2,&x3,&y3,&c);
pf12=caldis(x2,x1)+caldis(y2,y1);
pf23=caldis(x3,x2)+caldis(y3,y2);
pf13=caldis(x1,x3)+caldis(y1,y3);
if (pf12+pf23==pf13)
{
x4=x1+x3-x2;
y4=y1+y3-y2;
}
if (pf12+pf13==pf23)
{
x4=x2+x3-x1;
y4=y2+y3-y1;
}
if (pf13+pf23==pf12)
{
x4=x1+x2-x3;
y4=y1+y2-y3;
}
point[tot].x=x1;
point[tot].y=y1;
point[tot+1].x=x2;point[tot+1].y=y2;
point[tot+2].x=x3;point[tot+2].y=y3;
point[tot+3].x=x4;point[tot+3].y=y4;
pf14=caldis(x4,x1)+caldis(y4,y1);
pf24=caldis(x4,x2)+caldis(y4,y2);
pf34=caldis(x4,x3)+caldis(y4,y3);
dis[tot][tot+1]=((double)sqrt(pf12))*(double)c;
dis[tot+1][tot]=dis[tot][tot+1];
dis[tot][tot+2]=((double)sqrt(pf13))*(double)c;
dis[tot+2][tot]=dis[tot][tot+2];
dis[tot][tot+3]=((double)sqrt(pf14))*(double)c;
dis[tot+3][tot]=dis[tot][tot+3];
dis[tot+1][tot+2]=((double)sqrt(pf23))*(double)c;
dis[tot+2][tot+1]=dis[tot+1][tot+2];
dis[tot+1][tot+3]=((double)sqrt(pf24))*(double)c;
dis[tot+3][tot+1]=dis[tot+1][tot+3];
dis[tot+2][tot+3]=((double)sqrt(pf34))*(double)c;
dis[tot+3][tot+2]=dis[tot+2][tot+3];
for (int j=1;j<tot;j++)
{
p1=caldis(x1,point[j].x)+caldis(y1,point[j].y);
dis[tot][j]=((double)sqrt(p1))*(double)t;
dis[j][tot]=dis[tot][j];
p2=caldis(x2,point[j].x)+caldis(y2,point[j].y);
dis[tot+1][j]=((double)sqrt(p2))*(double)t;
dis[j][tot+1]=dis[tot+1][j];
p3=caldis(x3,point[j].x)+caldis(y3,point[j].y);
dis[tot+2][j]=((double)sqrt(p3))*(double)t;
dis[j][tot+2]=dis[tot+2][j];
p4=caldis(x4,point[j].x)+caldis(y4,point[j].y);
dis[tot+3][j]=((double)sqrt(p4))*(double)t;
dis[j][tot+3]=dis[tot+3][j];
}
tot+=4;
}
Floyd(s*4);
double ans=214748364;
for (int i=(a-1)*4+1;i<=(a-1)*4+4;i++)
for (int j=(b-1)*4+1;j<=(b-1)*4+4;j++)
ans=min(ans,dis[i][j]);
printf("%0.1lf\n",ans);
}
}