codevs 1041

又到暑假了,住在城市A的Car想和朋友一起去城市B旅游。她知道每个城市都有四个飞机场,分别位于一个矩形的四个顶点上,同一个城市中两个机场之间有一条笔直的高速铁路,第I个城市中高速铁路了的单位里程价格为Ti,任意两个不同城市的机场之间均有航线,所有航线单位里程的价格均为t。

那么Car应如何安排到城市B的路线才能尽可能的节省花费呢?她发现这并不是一个简单的问题,于是她来向你请教。 任务 找出一条从城市A到B的旅游路线,出发和到达城市中的机场可以任意选取,要求总的花费最少。

输入描述                 Input Description

第一行为一个正整数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个城市高速铁路单位里程的价格。

                输出描述                 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

#include<iostream>

#include<cstring>

#include<cstdio>

#include<cmath>

#define inf 0x1f1f1f1f

using namespace std;

double an[110][4][2];

double tn[110];

double t;

int s,A,B;

struct Edge

{

int c,d,next;

double w;

}edge[100000];

int head[10000][4],e;

void erase()

{

memset(head,-1,sizeof(head));

e = 0;

}

void increase(int a,int b,int c,int d,double w)

{

edge[e].c = c;

edge[e].d = d;

edge[e].w = w;

edge[e].next = head[a][b];

head[a][b] = e++;

}

void fq(int i)

{

double a = an[i][0][0],b = an[i][0][1],c = an[i][1][0],d = an[i][1][1],e = an[i][2][0],f = an[i][2][1];

double l11,l12,l21,l22;

l11 = c-a,l12 = d-b;

l21 = e-a,l22 = f-b;

if(l11*l21+l12*l22 == 0)

{

an[i][3][0] = c+l21;

an[i][3][1] = d+l22;

}

l11 = a-c,l12 = b-d;

l21 = e-c,l22 = f-d;

if(l11*l21+l12*l22 == 0)

{

an[i][3][0] = a+l21;

an[i][3][1] = b+l22;

}

l11 = a-e,l12 = b-f;

l21 = c-e,l22 = d-f;

if(l11*l21+l12*l22 == 0)

{

an[i][3][0] = a+l21;

an[i][3][1] = b+l22;

}

}

double fl(int a,int b,int c,int d)

{

return sqrt((an[a][b][0]-an[c][d][0])*(an[a][b][0]-an[c][d][0])+(an[a][b][1]-an[c][d][1])*(an[a][b][1]-an[c][d][1]));

}

double dis[110][4];

int q[100000][2],l,r;

bool inque[110][4];

void spfa(int a,int b)

{

int i,j;

for(i = 0;i<110;i++)

for(j = 0;j<4;j++)

dis[i][j] = 10000000000;

l = r = 0;

memset(inque,0,sizeof(inque));

dis[a][b] = 0;

q[r][0] = a,q[r++][1] = b;

inque[a][b] = 1;

while(l!=r)

{

int x = q[l][0],y = q[l++][1];

if(l == 100000) l = 0;

for(int t = head[x][y];t!=-1;t = edge[t].next)

{

double w = edge[t].w;

int xx = edge[t].c,yy = edge[t].d;

if(dis[xx][yy]>dis[x][y]+w)

{

dis[xx][yy] = dis[x][y]+w;

if(!inque[xx][yy])

{

inque[xx][yy] = 1;

q[r][0] = xx,q[r++][1] = yy;

if(r==100000)r = 0;

}

}

}

inque[x][y] = 0;

}

}

int main()

{

int z,i,j,k,ii,jj,kk;

cin>>z;

while(z--)

{

erase();

cin>>s>>t>>A>>B;

for(i = 1;i<=s;i++)

{

for(j = 0;j<3;j++)

{

cin>>an[i][j][0]>>an[i][j][1];

}

cin>>tn[i];

fq(i);

for(j = 0;j<4;j++)

{

for(k = j+1;k<4;k++)

{

double l = fl(i,j,i,k);

double w = l*tn[i];

increase(i,j,i,k,w);

increase(i,k,i,j,w);

}

}

}

for(i = 1;i<=s;i++)

{

for(j = i+1;j<=s;j++)

{

for(ii = 0;ii<4;ii++)

{

for(jj = 0;jj<4;jj++)

{

double l = fl(i,ii,j,jj);

double w = l*t;

increase(i,ii,j,jj,w);

increase(j,jj,i,ii,w);

}

}

}

}

double mi = 10000000000;

for(i = 0;i<4;i++)

{

spfa(A,i);

for(j = 0;j<4;j++)

{

if(mi>dis[B][j]) mi = dis[B][j];

}

}

printf("%.1lf\n",mi);

}

return 0;

}

codevs1041的更多相关文章

随机推荐

  1. cookie 和 session 基本使用 以及 封装

    Cookie: 是一小段文本信息,用户请求页面的时候,在浏览器和服务器之间传递.用户每次访问的时候都会记录cookie,cookie里可以包含用户信息,浏览的历史记录等等:Cookie是由服务器端生成 ...

  2. log4net使用简介

    平常我们在开发网站时,有一些比较重要的地方需要添加日志记录.一般日志记录分为两种:1)在数据库中添加一张日志表,用来记录用户操作并给用户提醒(用户可以看到). 2)在系统中添加一个日志文件,用来记录一 ...

  3. 2014年1月24日 Oracle 事务导读

    形象举例:  从 A 账户向 B 账户转账 10000 元 步骤: 1. A - 10000 2. B + 10000 事务的作用就是确保这两步无误执行后提交,若有一个执行不成功则失败. 结束事务的5 ...

  4. Asp.net MVC razor语法参考

    Razor语法的快捷参考http://haacked.com/archive/2011/01/06/razor-syntax-quick-reference.aspx/ 只是copy下来便于查阅! I ...

  5. Aandroid Error之 新导入工程报Unable to resolve target 'android-18'和R cannot be resolved

    有段时间没有写安卓了,今天导入以前的项目,结果看到控制台打印出了这样一句:Unable to resolve target 'android-18', 解决方法: 项目->属性->Andr ...

  6. JSP基础之 C标签中的 varStatues属性

    变量状态参数,该属性有5个常用值count 表示当前遍历集合的元素个数index 表示当前遍历到集合的第几个元素current 表示当前的集合元素first 表示集合的第一个元素last 表示集合的最 ...

  7. 在Yii框架中使用PHPExcel

    PHPExcel是一个比较好用的php读取excel文件的类库,今天遇到了在yii中如何加载PHPExcel类文件的问题,因为Yii的autoload机制是安装类名去找文件,即文件名就是相应的类名,而 ...

  8. MySQL的表分区(转载)

    MySQL的表分区(转载) 一.什么是表分区 通俗地讲表分区是将一大表,根据条件分割成若干个小表.mysql5.1开始支持数据表分区了. 如:某用户表的记录超过了600万条,那么就可以根据入库日期将表 ...

  9. angularjs factory,service,provider 自定义服务的不同

    angularjs框架学了有一段时间了,感觉很好用.可以把angularjs的app理解成php的class,controller是控制器,而内置服务和自定义服务就可以理解成models了.angul ...

  10. 【ecos学习5】redboot 加载运行hello world

    背景: 从主机 192.168.2.14 IP,下载bin文件hello到ecos. redboot>load -v -h 192.168.2.14 hello Using default pr ...