URAL 1750 Pakhom and the Gully 计算几何+floyd
题目链接:点击打开链接
gg。。
。
#pragma comment(linker, "/STACK:1024000000,1024000000")
#include <cstdio>
#include <iostream>
#include <algorithm>
#include <string.h>
#include <vector>
#include <queue>
#include <math.h>
using namespace std;
#define ll int
#define point Point
const double eps = 1e-8;
const double PI = acos(-1.0);
double ABS(double x){return x>0?x:-x;}
int sgn(double x){
if(fabs(x) < eps)return 0;
if(x<0)return -1;
return 1;
}
struct Point
{
double x, y;
void put(){printf("%.0f,%.0f\n",x,y);}
Point(){}
Point(double _x, double _y){
x = _x; y = _y;
}
Point operator - (const Point & b)const{
return Point(x-b.x, y-b.y);
}
double operator ^ (const Point &b) const{
return x*b.y - y*b.x;
}
double operator *(const Point &b) const{
return x*b.x + y*b.y;
}
void transXY(double B){
double tx = x, ty = y;
x = tx*cos(B) - ty*sin(B);
y = tx*sin(B) + ty*cos(B);
}
};
double dist(point a,point b){ return sqrt((a.x-b.x)*(a.x-b.x)+(a.y-b.y)*(a.y-b.y));}
struct Line
{
Point s,e;
double len(){return dist(s,e);}
Line(){}
Line(Point _s, Point _e)
{
s = _s; e = _e;
}
pair<int,Point> operator & (const Line &b) const{
Point res = s;
if(sgn((s-e) ^ (b.s-b.e)) == 0)
{
if(sgn((s-b.e) ^ (b.s-b.e)) == 0)
return make_pair(0,res);
else return make_pair(1,res);
}
double t = ((s-b.s) ^ (b.s-b.e)) / ((s-e)^(b.s-b.e));
res.x += (e.x - s.x) *t;
res.y += (e.y - s.y) *t;
return make_pair(2, res);
}
};
double cross(point a, point b, point c) {
return (a.x - c.x) * (b.y - c.y) - (a.y - c.y) * (b.x - c.x);
}
bool intersect(Line l1, Line l2) {
return
max(l1.s.x, l1.e.x) > (min(l2.s.x, l2.e.x)-eps) &&
max(l2.s.x, l2.e.x) > (min(l1.s.x, l1.e.x)-eps) &&
max(l1.s.y, l1.e.y) > (min(l2.s.y, l2.e.y)-eps) &&
max(l2.s.y, l2.e.y) > (min(l1.s.y, l1.e.y)-eps) &&
(cross(l2.s, l1.e, l1.s) * cross(l1.e, l2.e, l1.s) > -eps) &&
(cross(l1.s, l2.e, l2.s) * cross(l2.e, l1.e, l2.s) > -eps);
}
bool online(Line l, point P) {
return (fabs(cross(l.e, P, l.s)) < eps) &&
((((P.x > l.s.x - eps) && (P.x < l.e.x + eps)) || ((P.x > l.e.x - eps) && (P.x < l.s.x+eps)))&&
(((P.y > l.s.y - eps) && (P.y < l.e.y + eps)) || ((P.y > l.e.y - eps) && (P.y < l.s.y +eps))));
}
bool inter(Line u, Line v){
return ((intersect(u, v)) && (!online(u, v.s)) && (!online(u, v.e)) && (!online(v, u.s)) && (!online(v, u.e)));
}
point S, T, A, B, C;
Line AB, BC, ST, AC, SA, SB, SC, TA, TB, TC;
double ans;
bool ok(Line x){
if(inter(x, AB))return false;
if(inter(x, BC))return false;
return true;
}
double d[5][5];
void addedge()
{
for (int i=0;i<5;i++)
for (int j=0;j<5;j++)
{
if (i==j) d[i][j]=0;
else d[i][j]=1e8;
} if (!inter(SA, BC))
{
d[0][1]=d[1][0]=min(d[0][1],dist(S,A));
}
//d[0][2]=d[2][0]=min(d[0][2],dis(S,B)); d[1][2]=d[2][1]=min(d[1][2],dist(A,B));
d[1][3]=d[3][1]=min(d[1][3],dist(A,C));
d[2][3]=d[3][2]=min(d[2][3],dist(B,C)); if (!inter(SC,AB))
{
d[0][3]=d[3][0]=min(d[0][3],dist(S,C));
} if (!intersect(ST,AB) && !intersect(ST,BC))
{
d[0][4]=d[4][0]=min(d[0][4],dist(S,T));
}
//printf("1->4 %d\n",is_inter(A,T,B,C));
if (!inter(TA,BC))
{
d[1][4]=d[4][1]=min(d[1][4],dist(A,T));
} if (cross(B,A,T)*cross(B,C,T)>-eps)
{
d[2][4]=d[4][2]=min(d[2][4],dist(B,T));
} if (cross(S,A,B)*cross(S,C,B)>-eps)
{
d[0][2]=d[2][0]=min(d[0][2],dist(S,B));
} if (!inter(TC,AB))
{
d[3][4]=d[4][3]=min(d[3][4],dist(C,T));
} }
int main(){
int i, j, k, cas;scanf("%d",&cas);
while(cas--) {
scanf("%lf %lf", &S.x, &S.y);
scanf("%lf %lf", &T.x, &T.y);
scanf("%lf %lf", &A.x, &A.y);
scanf("%lf %lf", &B.x, &B.y);
scanf("%lf %lf", &C.x, &C.y);
AB = Line(A,B); BC = Line(B,C); ST = Line(S,T); AC = Line(A,C);
SA = Line(S,A); SB = Line(S,B); SC = Line(S,C);
TA = Line(T,A); TB = Line(T,B); TC = Line(T,C);
addedge();
for(i = 0; i < 5; i++)
for(j = 0; j < 5; j++)
for(k = 0; k < 5; k++)
d[i][j] = min(d[i][j], d[i][k]+d[k][j]);
printf("%.8f\n", d[0][4]);
}
return 0;
}
/* */
URAL 1750 Pakhom and the Gully 计算几何+floyd的更多相关文章
- BZOJ 1027 JSOI2007 合金 计算几何+Floyd
题目大意:给定一些合金,选择最少的合金,使这些合金能够按比例合成要求的合金 首先这题的想法特别奇异 看这题干怎么会想到计算几何 并且计算几何又怎么会跟Floyd挂边 好强大 首先因为a+b+c=1 所 ...
- hdu 4606 简单计算几何+floyd+最小路径覆盖
思路:将所有的直线的两个端点和城市混在一起,将能直接到达的两个点连线,求一次floyd最短路径.二分枚举bag容量,然后按给的要先后占领的城市由前向后,把能到一步到达的建一条边.然后求一次最小路径覆盖 ...
- BZOJ_1027_[JSOI2007]_合金_(计算几何+Floyd求最小环)
描述 http://www.lydsy.com/JudgeOnline/problem.php?id=1027 共三种金属,\(m\)种材料,给出每种材料中三种金属的占比. 给出\(n\)种合金的三种 ...
- bzoj 1027 [JSOI2007]合金(计算几何+floyd最小环)
1027: [JSOI2007]合金 Time Limit: 4 Sec Memory Limit: 162 MBSubmit: 2970 Solved: 787[Submit][Status][ ...
- BZOJ1027 [JSOI2007]合金 【计算几何 + floyd】
题目 某公司加工一种由铁.铝.锡组成的合金.他们的工作很简单.首先进口一些铁铝锡合金原材料,不同种类的 原材料中铁铝锡的比重不同.然后,将每种原材料取出一定量,经过融解.混合,得到新的合金.新的合金的 ...
- 【JZOJ5094】【GDSOI2017第四轮模拟day3】鸽子 计算几何+floyd
题面 养鸽人要监视他的鸽子,有n只鸽子站在平面上,他可以在m个给定的点上设置监视器,如果一只鸽子在某个监视器上或者在两个监视器所连直线上或者在三个监视器所连直线的三角形内则其就咕咕咕了,现在养鸽人要让 ...
- URAL 2099 Space Invader题解 (计算几何)
啥也不说了,直接看图吧…… 代码如下: #include<stdio.h> #include<iostream> #include<math.h> using na ...
- BZOJ 1027: [JSOI2007]合金 (计算几何+Floyd求最小环)
题解就看这位仁兄的吧-不过代码还是别看他的了- 同样的方法-我200ms,他2000ms. 常数的幽怨- CODE #include <bits/stdc++.h> using names ...
- 杭电ACM分类
杭电ACM分类: 1001 整数求和 水题1002 C语言实验题——两个数比较 水题1003 1.2.3.4.5... 简单题1004 渊子赛马 排序+贪心的方法归并1005 Hero In Maze ...
随机推荐
- golang笔记:net/smtp
跟go语言的net/smtp斗争了一天,记录下历程. 先用最标准的例子 host := net.JoinHostPort(hostname, port) auth := smtp.PlainAut ...
- SPOJ BGSHOOT - Shoot and kill (线段树 区间修改 区间查询)
BGSHOOT - Shoot and kill no tags The problem is about Mr.BG who is a great hunter. Today he has gon ...
- 循环节(BFS)
循环节 时间限制: 1 Sec 内存限制: 64 MB提交: 56 解决: 16[提交][状态][讨论版] 题目描述 第一节是英语课.今天,老师又教了桐桐很多单词.桐桐发现所有单词都有循环节(大写 ...
- 素数路(prime) (BFS)
问题 C: 素数路(prime) 时间限制: 1 Sec 内存限制: 64 MB提交: 8 解决: 5[提交][状态][讨论版] 题目描述 已知一个四位的素数,要求每次修改其中的一位,并且要保证修 ...
- Scrum生命周期
Recently while cleaning up my photo albums I found some interesting old pictures which were captured ...
- 【可持久化Trie】bzoj3261 最大异或和
对原序列取前缀异或值,变成pre[1...N],然后询问等价于求max{a[N]^x^pre[i]}(l-1<=i<=r-1). #include<cstdio> #defin ...
- Mac sublime 编译Python UnicodeEncodeError: 'ascii' codec can't encode characters in position 6-8: ordinal not in range(128)
刚学Python,想打印个“hello 张林峰”,代码如下: #!/usr/bin/env python3 # -*- coding: utf-8 -*- print('hello 张林峰') 用su ...
- spark checkpoint机制
首先rdd.checkpoint()本身并没有执行任何的写操作,只是做checkpointDir是否为空,然后生成一个ReliableRDDCheckpointData对象checkpointData ...
- Linux查找某个时间点后生成的文件(转)
需要找到某天(例如2017-04-13)以及这之后生成的空文件.那么这个要怎么处理呢?这个当然是用find命令来解决.如下所示, -mtime -5表示查找距现在5*24H内修改过的文件 -type ...
- 常见的 HTTP 状态代码及原因
代码 说明 备注 200 确定 IIS 7.0.IIS 7.5 和 IIS 8.0 成功处理了请求. 304 未修改 客户端浏览器请求已处于缓存中的文档,并且自从该文档被缓存后,未修改此文档.客户端浏 ...