Codevs1041 Car的旅行路線
阿新 • • 發佈:2019-01-31
- 題目大意:給定一張無向圖,其中的某些點構成某個集合,求從一個集合至另一個集合的最短路。
- 思路:根據題目中的條件求出兩兩點間的花費,用floyd就能過。但計算時需分類討論,在同一個城市裡的點用鐵路費用計算,不在的用航空費用計算。
- 在程式中:利用vector記錄表示城市的節點,getmap( )用於建圖。getanother( )用於得到第四個飛機場的座標,其中利用其它三個點的座標計算時,公式為 xi=ct[1].x+ct[3].x-ct[2].x與yi=ct[1].y+ct[3].y-ct[2].y。需要注意的是,ct[2]是處於直角頂上的點,所以在程式中不斷交換ct[ ]的三個元素以保證這一條件的成立。判斷條件使用(ct[1].x-ct[2].x)(ct[3].x-ct[2].x)+(ct[1].y-ct[2].y)(ct[3].y-ct[2].y)=0,不能直接使用k1*k2=-1,因為可能存在斜率不存在的情況。
- 程式碼如下:
#include<iostream>
#include<cstdio>
#include<cmath>
#include<vector>
struct node
{
int x,y,id;
node(int x=0,int y=0,int id=0):x(x),y(y),id(id){}
};
using namespace std;
const int oo=1000000000;
int n,s,t,a,b,tot=0;
double map[303][303];
int ti[105];
vector<node> city[105 ];
double dist(int x1,int y1,int x2,int y2)
{
double ans=(x1-x2)*(x1-x2)+(y1-y2)*(y1-y2);
return sqrt(ans);
}
void getmap()
{
int i,j,k;
vector<node>::iterator a,b;
for (i=1;i<=s;++i)
for (a=city[i].begin();a!=city[i].end();++a)
for (b=city[i].begin();b!=city[i].end();++b)
if (a!=b)
{
map[(*a).id][(*b).id]=map[(*b).id][(*a).id]=dist((*a).x,(*a).y,(*b).x,(*b).y)*ti[i];
}
for (i=1;i<=s;++i)
for (j=1;j<=s;++j)
if (i!=j)
{
for (a=city[i].begin();a!=city[i].end();++a)
for (b=city[j].begin();b!=city[j].end();++b)
map[(*a).id][(*b).id]=map[(*b).id][(*a).id]=dist((*a).x,(*a).y,(*b).x,(*b).y)*t;
}
}
void floyd()
{
int i,j,k,ui=city[a].front().id,vi=city[b].front().id;
for (k=1;k<=tot;++k)
for (i=1;i<=tot;++i)
for (j=1;j<=tot;++j)
if (i!=j && i!=k && j!=k)
if (map[i][j]>map[i][k]+map[k][j])
map[i][j]=map[i][k]+map[k][j];
double ans=oo;
for (i=ui;i<=ui+3;++i)
for (j=vi;j<=vi+3;++j)
{
if (map[i][j]<ans)
ans=map[i][j];
}
printf("%.1f\n",ans);
}
void getanother(int i)
{
int num=0;
node ct[4];
vector<node>::iterator it;
for (it=city[i].begin();it!=city[i].end();++it)
{
ct[++num].x=(*it).x;
ct[num].y=(*it).y;
}
while ( (ct[1].x-ct[2].x)*(ct[3].x-ct[2].x)+(ct[1].y-ct[2].y)*(ct[3].y-ct[2].y) )
{
node tmp=ct[1];
ct[1]=ct[2];
ct[2]=ct[3];
ct[3]=tmp;
}
int xi=ct[1].x+ct[3].x-ct[2].x;
int yi=ct[1].y+ct[3].y-ct[2].y;
city[i].push_back(node(xi,yi,++tot));
}
void init()
{
int xi,yi,i,j;
cin>>n;
while (n--)
{
scanf("%d%d%d%d",&s,&t,&a,&b);
for (i=1;i<=s;++i)
for (j=1;j<=s;++j)
map[i][j]=oo;
for (i=1;i<=s;++i)
{
for (j=1;j<=3;++j)
{
scanf("%d%d",&xi,&yi);
city[i].push_back(node(xi,yi,++tot));
}
getanother(i);
scanf("%d",&ti[i]);
}
getmap();
floyd();
}
}
int main()
{
init();
return 0;
}