codevs1041
2015-07-19 11:02
344 查看
题目地址:http://codevs.cn/problem/1041/
分析:
首先求出第四个点坐标
先根据的读入的三个点确定哪个是直角顶点,设为A B C 三点其中A为直角顶点
然后设另一点为D,连BC取中点E,则AD的中点也为E再由线段中点坐标公式((x1+x2)/2,(y1+y2)/2)可推知x1+x4=x2+x3 y1+y4=y2+y3 这就求出了第四点
将每个飞机场看成一个节点,则有4*s个节点
构建两个虚节点,一个作为0节点,另一个作为4*s+1,然后求出任两点的距离,然后就赤裸裸的最短路了
代码:
program travel;
var i,j,k,l,t1,s,c,d,n,c12,c23,c13,j1,j2,k1,k2:longint;
min1:real;
t:array[1..100]of longint;
f:array[0..401,0..401]of real;
cost:array[0..401]of real;
b:array[0..401]of boolean;
a:array[1..100,1..4,1..2]of longint;
procedure init;
begin
readln(s,t1,c,d);
for j:=1 to s do
begin
for k:=1 to 3 do
for l:=1 to 2 do
read(a[j,k,l]);
readln(t[j]);
end;
end;
procedure swap(var a,b:longint);
var x:longint;
begin
x:=a;a:=b;b:=x;
end;
function dis(var x1,y1,x2,y2:longint):real;
begin
exit(sqrt(sqr(x1-x2)+sqr(y1-y2)));
end;
procedure fourth;
begin
for j:=1 to s do
begin
c12:=sqr(a[j,1,1]-a[j,2,1])+sqr(a[j,1,2]-a[j,2,2]);
c13:=sqr(a[j,1,1]-a[j,3,1])+sqr(a[j,1,2]-a[j,3,2]);
c23:=sqr(a[j,2,1]-a[j,3,1])+sqr(a[j,2,2]-a[j,3,2]);
if c12+c23=c13 then begin swap(a[j,1,1],a[j,2,1]);swap(a[j,1,2],a[j,2,2]);end
else if c13+c23=c12 then begin swap(a[j,1,1],a[j,3,1]);swap(a[j,1,2],a[j,3,2]);end;
a[j,4,1]:=a[j,2,1]+a[j,3,1]-a[j,1,1];a[j,4,2]:=a[j,2,2]+a[j,3,2]-a[j,1,2];
end;
end;
function min(a,b:real):real;
begin
if a<b then exit(a) else exit(b);
end;
procedure count;
begin
for j:=1 to 4*s do
for k:=1 to 4*s do
if j<>k then
begin
j1:=(j-1)div 4+1;
j2:=(j-1)mod 4+1;
k1:=(k-1)div 4+1;
k2:=(k-1)mod 4+1;
if j1=k1 then
begin
f[j,k]:=dis(a[j1,j2,1],a[j1,j2,2],a[k1,k2,1],a[k1,k2,2])*t[j1];
if f[j,k]>dis(a[j1,j2,1],a[j1,j2,2],a[k1,k2,1],a[k1,k2,2])*t1 then
f[j,k]:=dis(a[j1,j2,1],a[j1,j2,2],a[k1,k2,1],a[k1,k2,2])*t1;
f[k,j]:=f[j,k];
end
else
begin
f[j,k]:=dis(a[j1,j2,1],a[j1,j2,2],a[k1,k2,1],a[k1,k2,2])*t1;
f[k,j]:=f[j,k];
end;
end;
end;
procedure dijkstra;
begin
for j:=1 to 4*s+1 do f[0,j]:=maxlongint;
for j:=1 to 4 do
f[0,c*4-4+j]:=0;
for j:=1 to 4*s do f[j,4*s+1]:=maxlongint;
for j:=1 to 4 do
f[d*4-4+j,4*s+1]:=0;
for j:=1 to 4*s+1 do
cost[j]:=f[0,j];
for j:=1 to 4*s+1 do
begin
min1:=maxlongint;l:=-1;
for k:=1 to 4*s+1 do
if(not b[k])and(cost[k]<min1) then
begin
min1:=cost[k];
l:=k;
end;
if l=-1 then break;
b[l]:=true;
for k:=0 to 4*s+1 do
if cost[l]+f[l,k]<cost[k] then
cost[k]:=cost[l]+f[l,k];
end;
end;
procedure main;
begin
readln(n);
for i:=1 to n do
begin
init;
fourth;
count;
dijkstra;
writeln(cost[4*s+1]:0:1);
end;
end;
begin
main;
end.
分析:
首先求出第四个点坐标
先根据的读入的三个点确定哪个是直角顶点,设为A B C 三点其中A为直角顶点
然后设另一点为D,连BC取中点E,则AD的中点也为E再由线段中点坐标公式((x1+x2)/2,(y1+y2)/2)可推知x1+x4=x2+x3 y1+y4=y2+y3 这就求出了第四点
将每个飞机场看成一个节点,则有4*s个节点
构建两个虚节点,一个作为0节点,另一个作为4*s+1,然后求出任两点的距离,然后就赤裸裸的最短路了
代码:
program travel;
var i,j,k,l,t1,s,c,d,n,c12,c23,c13,j1,j2,k1,k2:longint;
min1:real;
t:array[1..100]of longint;
f:array[0..401,0..401]of real;
cost:array[0..401]of real;
b:array[0..401]of boolean;
a:array[1..100,1..4,1..2]of longint;
procedure init;
begin
readln(s,t1,c,d);
for j:=1 to s do
begin
for k:=1 to 3 do
for l:=1 to 2 do
read(a[j,k,l]);
readln(t[j]);
end;
end;
procedure swap(var a,b:longint);
var x:longint;
begin
x:=a;a:=b;b:=x;
end;
function dis(var x1,y1,x2,y2:longint):real;
begin
exit(sqrt(sqr(x1-x2)+sqr(y1-y2)));
end;
procedure fourth;
begin
for j:=1 to s do
begin
c12:=sqr(a[j,1,1]-a[j,2,1])+sqr(a[j,1,2]-a[j,2,2]);
c13:=sqr(a[j,1,1]-a[j,3,1])+sqr(a[j,1,2]-a[j,3,2]);
c23:=sqr(a[j,2,1]-a[j,3,1])+sqr(a[j,2,2]-a[j,3,2]);
if c12+c23=c13 then begin swap(a[j,1,1],a[j,2,1]);swap(a[j,1,2],a[j,2,2]);end
else if c13+c23=c12 then begin swap(a[j,1,1],a[j,3,1]);swap(a[j,1,2],a[j,3,2]);end;
a[j,4,1]:=a[j,2,1]+a[j,3,1]-a[j,1,1];a[j,4,2]:=a[j,2,2]+a[j,3,2]-a[j,1,2];
end;
end;
function min(a,b:real):real;
begin
if a<b then exit(a) else exit(b);
end;
procedure count;
begin
for j:=1 to 4*s do
for k:=1 to 4*s do
if j<>k then
begin
j1:=(j-1)div 4+1;
j2:=(j-1)mod 4+1;
k1:=(k-1)div 4+1;
k2:=(k-1)mod 4+1;
if j1=k1 then
begin
f[j,k]:=dis(a[j1,j2,1],a[j1,j2,2],a[k1,k2,1],a[k1,k2,2])*t[j1];
if f[j,k]>dis(a[j1,j2,1],a[j1,j2,2],a[k1,k2,1],a[k1,k2,2])*t1 then
f[j,k]:=dis(a[j1,j2,1],a[j1,j2,2],a[k1,k2,1],a[k1,k2,2])*t1;
f[k,j]:=f[j,k];
end
else
begin
f[j,k]:=dis(a[j1,j2,1],a[j1,j2,2],a[k1,k2,1],a[k1,k2,2])*t1;
f[k,j]:=f[j,k];
end;
end;
end;
procedure dijkstra;
begin
for j:=1 to 4*s+1 do f[0,j]:=maxlongint;
for j:=1 to 4 do
f[0,c*4-4+j]:=0;
for j:=1 to 4*s do f[j,4*s+1]:=maxlongint;
for j:=1 to 4 do
f[d*4-4+j,4*s+1]:=0;
for j:=1 to 4*s+1 do
cost[j]:=f[0,j];
for j:=1 to 4*s+1 do
begin
min1:=maxlongint;l:=-1;
for k:=1 to 4*s+1 do
if(not b[k])and(cost[k]<min1) then
begin
min1:=cost[k];
l:=k;
end;
if l=-1 then break;
b[l]:=true;
for k:=0 to 4*s+1 do
if cost[l]+f[l,k]<cost[k] then
cost[k]:=cost[l]+f[l,k];
end;
end;
procedure main;
begin
readln(n);
for i:=1 to n do
begin
init;
fourth;
count;
dijkstra;
writeln(cost[4*s+1]:0:1);
end;
end;
begin
main;
end.
相关文章推荐
- CoreAnimation的使用小结
- Hadoop—MapReduce进行数据查询和实现推简单荐系统---练习7
- 记录一下cocos2d-x——开发环境安装
- 团购(cell高度相同时自定义cell)
- [LeetCode][Java] Remove Duplicates from Sorted List
- 关于eclipse安装Genymotion插件的方法
- 用矩阵键盘使数码管显示不同的值
- 经典的排序算法
- 专题:01分数规划
- 3.hibernate基础配置
- 关于cocoapods一些错误的解决问
- OpenCV Access pixel value
- hasse数列
- 发送邮件
- POJ 3352 Road Construction (边双连通分量 Tarjan缩点)
- QT 环境的搭建
- busybox的作用和功能
- java学习之旅64--常用类_自动装箱和拆箱_缓存处理
- boost::implicit_cast
- 《JAVA 核心技术 基础知识》第三章 笔记