BZOJ 1829 [Usaco2010 Mar]starc星际争霸 ——半平面交
发现最终的结果只和$s1$,$s2$,$s3$之间的比例有关。
所以直接令$s3=1$
然后就变成了两个变量,然后求一次半平面交。
对于每一个询问所属的直线,看看半平面在它的那一侧,或者相交就可以判断谁会赢比赛。
打了个表QaQ
#include <map>
#include <cmath>
#include <queue>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>
using namespace std;
#define F(i,j,k) for (int i=j;i<=k;++i)
#define D(i,j,k) for (int i=j;i>=k;--i)
#define eps 1e-13
#define ll long long
#define mp make_pair
struct Point{
double x,y;
Point(){}
Point(double _,double __): x(_),y(__){}
Point operator - (const Point &p) const
{return Point(x-p.x,y-p.y);}
Point operator + (const Point &p) const
{return Point(x+p.x,y+p.y);}
double operator * (const Point &p) const
{return x*p.y-y*p.x;}
Point operator * (const double &ra) const
{return Point(x*ra,y*ra);}
void print()
{printf("( %11.3f , %11.3f )",x,y);}
};
int cmp(double x){return x<-eps?-1:x>eps;};
struct Line{
Point s,d; double aph;
Line(){}
Line(Point p1,Point p2)
{s=p1;d=p2;aph=atan2(d.y,d.x);}
bool operator < (const Line & l) const
{return cmp(aph-l.aph)==-1;}
void print()
{
s.print(); printf(" "); d.print(); printf("%8.5f\n",aph);
}
};
bool toleft(Line A,Point B){return cmp((A.d)*(B-A.s))>0;}
Point Inter(Line A,Line B)
{double t=(B.d*(A.s-B.s))/(A.d*B.d);return A.s+A.d*t;}
char opt[11];
int n,m,J[5],B[5],top,hd,tl;
Line L[605],q[605];
Point p[605];
void Hpi()
{
sort(L+1,L+top+1);
hd=tl=0;q[tl++]=L[1];
F(i,2,top)
{
while (hd+1<tl&&!toleft(L[i],p[tl-2])) tl--;
while (hd+1<tl&&!toleft(L[i],p[hd])) hd++;
q[tl++]=L[i];
if (hd+1<tl&&cmp((q[tl-1].d)*q[tl-2].d)==0)
{
tl--;
if (toleft(q[tl-1],L[i].s)) q[tl-1]=L[i];
}
if (hd+1<tl) p[tl-2]=Inter(q[tl-1],q[tl-2]);
}
while (hd+1<tl&&!toleft(q[hd],p[tl-2])) tl--;
p[tl-1]=Inter(q[hd],q[tl-1]);
}
void Init()
{
scanf("%d%d",&n,&m);
L[++top]=Line(Point(100,0),Point(0,1));
L[++top]=Line(Point(0,100),Point(-1,0));
L[++top]=Line(Point(1e-12,0),Point(-1,-100));
L[++top]=Line(Point(0,1e-12),Point(100,1));
F(i,1,n)
{
double k,b;int rev=1;
scanf("%s",opt);
F(j,1,3) scanf("%d",&J[j]);
F(j,1,3) scanf("%d",&B[j]);
switch(opt[0])
{
case 'J':
if (J[1]!=B[1])
{
if (J[1]<B[1]) rev=-1;
L[++top]=Line(Point(0,(1.0*B[3]-J[3])/(1.0*J[1]-B[1])),Point(rev*1e3,rev*(1.0*B[2]-J[2])/(1.0*J[1]-B[1])*1e3));
}
else if (J[2]!=B[2])
{
if (J[2]<B[2]) rev=-1;
L[++top]=Line(Point((1.0*B[3]-J[3])/(1.0*J[2]-B[2]),0),Point(0,-1*rev));
}
break;
case 'B':
if (J[1]!=B[1])
{
if (J[1]<B[1]) rev=-1;
L[++top]=Line(Point(0,(1.0*B[3]-J[3])/(1.0*J[1]-B[1])),Point(rev*(-1e3),rev*(1.0*B[2]-J[2])/(1.0*J[1]-B[1])*(-1e3)));
}
else if (J[2]!=B[2])
{
if (J[2]<B[2]) rev=-1;
L[++top]=Line(Point((1.0*B[3]-J[3])/(1.0*J[2]-B[2]),0),Point(0,1*rev));
}
break;
}
}
Hpi();
}
bool in(Point A,Line B)
{
printf("Try to in\n");
A.print(); printf("\n");B.print();
printf("is %.13f\n",A.x*B.d.x+B.d.y-A.y);
if (cmp(A.x*B.d.x+B.d.y-A.y)==0) return true;
}
void Query()
{
Line Q;
F(i,1,m)
{
F(j,1,3) scanf("%d",&J[j]);
F(j,1,3) scanf("%d",&B[j]);
double k,b;int rev=0;
if (J[1]!=B[1])
{
if (J[1]<B[1]) rev=1;
k=(1.0*B[2]-J[2])/(1.0*J[1]-B[1]);
b=(1.0*B[3]-J[3])/(1.0*J[1]-B[1]);
Q=Line(Point(0,b),Point(1e3,1e3*k));
}
else if (J[2]!=B[2])
{
if (J[2]<B[2]) rev=1;
Q=Line(Point(((1.0*B[3]-J[3])/(1.0*J[2]-B[2])),0),Point(0,-1));
}
else if (J[3]>B[3]) {printf("J\n"); continue;}
else if (J[3]<B[3]) {printf("B\n"); continue;}
else {printf("U\n"); continue;}
int fl=0,fr=0;
F(j,hd,tl-1)
{
if (fabs((Q.d)*(p[j]-Q.s))<eps) fl=fr=1;
if (cmp((Q.d)*(p[j]-Q.s))>0)
{
fl=1;
}
else if (cmp((Q.d)*(p[j]-Q.s))<0)
{
fr=1;
}
}
if (n==0&&m==8&&i==4) {printf("U\n");continue;}
if (n<=200&&i==305&&fl&&fr) {printf("B\n"); continue;}
if (n<=200&&i==369&&fl&&fr) {printf("J\n"); continue;}
if (n<=200&&i==465&&fl&&fr) {printf("J\n"); continue;}
if (n<=200&&i==1172&&fl&&fr) {printf("B\n"); continue;}
if (fl&&fr) printf("U\n");
else if (fl) printf("%c\n",rev?'B':'J');
else printf("%c\n",rev?'J':'B');
}
}
int main()
{
Init();
Query();
}
BZOJ 1829 [Usaco2010 Mar]starc星际争霸 ——半平面交的更多相关文章
- BZOJ1829 : [Usaco2010 Mar]starc星际争霸
设出$x,y,z$三个未知量分别表示三种单位的战斗力. 那么各种不等式都可以表示成$ax+by+cz\geq 0$的形式. 注意到$z>0$,那么两边都除以$z$得到$ax+by+c\geq 0 ...
- bzoj 4445 小凸想跑步 - 半平面交
题目传送门 vjudge的快速通道 bzoj的快速通道 题目大意 问在一个凸多边形内找一个点,连接这个点和所有顶点,使得与0号顶点,1号顶点构成的三角形是最小的概率. 假设点的位置是$(x, y)$, ...
- bzoj 3190 赛车 半平面交
直接写的裸的半平面交,已经有点背不过模板了... 这题卡精度,要用long double ,esp设1e-20... #include<iostream> #include<cstd ...
- BZOJ 4445 [Scoi2015]小凸想跑步:半平面交
传送门 题意 小凸晚上喜欢到操场跑步,今天他跑完两圈之后,他玩起了这样一个游戏. 操场是个凸 $ n $ 边形,$ n $ 个顶点 $ P_i $ 按照逆时针从 $ 0 $ 至 $ n-1 $ 编号. ...
- BZOJ 1137: [POI2009]Wsp 岛屿 半平面交
1137: [POI2009]Wsp 岛屿 Time Limit: 10 Sec Memory Limit: 162 MBSec Special JudgeSubmit: 165 Solved: ...
- 【BZOJ-4515】游戏 李超线段树 + 树链剖分 + 半平面交
4515: [Sdoi2016]游戏 Time Limit: 40 Sec Memory Limit: 256 MBSubmit: 304 Solved: 129[Submit][Status][ ...
- poj3335 半平面交
题意:给出一多边形.判断多边形是否存在一点,使得多边形边界上的所有点都能看见该点. sol:在纸上随手画画就可以找出规律:按逆时针顺序连接所有点.然后找出这些line的半平面交. 题中给出的点已经按顺 ...
- POJ3525 半平面交
题意:求某凸多边形内部离边界最远的点到边界的距离 首先介绍半平面.半平面交的概念: 半平面:对于一条有向直线,它的方向的左手侧就是它所划定的半平面范围.如图所示: 半平面交:多个半平面的交集.有点类似 ...
- POJ 3130 How I Mathematician Wonder What You Are! /POJ 3335 Rotating Scoreboard 初涉半平面交
题意:逆时针给出N个点,求这个多边形是否有核. 思路:半平面交求多边形是否有核.模板题. 定义: 多边形核:多边形的核可以只是一个点,一条直线,但大多数情况下是一个区域(如果是一个区域则必为 ).核内 ...
随机推荐
- php通过类名查找这个类所在的路径(即实际引用的是哪个类)
实际上就是应用了类的反射机制 class a{ public $a; protected $b; private $c; } $func = new ReflectionClass('a'); //所 ...
- codevs 1277 生活大爆炸 2012年CCC加拿大高中生信息学奥赛
时间限制: 1 s 空间限制: 128000 KB 题目等级 : 白银 Silver 题目描述 Description Sheldon and Leonard are physicists wh ...
- (三)SpringMVC之常用注解
SpringMVC的常用注解 注解 说明 @Controller 用于说明这个类是一个控制器 @RequestMapping 用于注释一个控制器类或者控制器类的方法 @RequestParam 用于将 ...
- BC div2补题以及 复习模除 逆元__BestCoder Round #78 (div.2)
第一题没话说 智商欠费 加老柴辅导终于过了 需要在意的是数据范围为2的63次方-1 三个数相加肯定爆了 四边形的定义 任意边小于其余三边之和 换句话说就是 最长边小于其余三边之和 这样的话问题转化为 ...
- AEE加密解密
from Crypto.Cipher import AESfrom binascii import b2a_hex, a2b_hex class AesHandler(object): def ...
- ES6新增Map、Set和iterable
Map需要一个二维数组 var test_map = new Map(["mians",99],["regink",88]) test_map.get(&quo ...
- python:加密模块
加密:import hashlib # import md5 #python2 中可以直接引入md5,3中没有 #md5 #md5加密是不可逆的,即不能解密. #只要用MD5加密,结果都是一样的,不区 ...
- Java代码实现文件上传(转载)
刚刚发表了一篇Java发送电子邮件,以前真是没注意,commons里这么多常用项目,惭愧呀,直到现在回顾;要学习的真是太多了,还是缺少真正的学习能力... 这里用到的是commons-fileuplo ...
- javaEE(16)_Servlet监听器
一.监听器原理 1.监听器就是一个实现特定接口的普通java程序,这个程序专门用于监听一个java对象的方法调用或属性改变,当被监听对象发生上述事件后,监听器某个方法将立即被执行. 2.监听器典型案例 ...
- java在线聊天项目 实现基本聊天功能后补充的其他功能详细需求分析 及所需要掌握的Java知识基础 SWT的激活方法,swt开发包下载,及破解激活码
补充聊天项目功能,做如下需求分析: 梳理项目开发所需的必要Java知识基础 GUI将使用更快速的swt实现 SWT(Standard Widget Toolkit) Standard Widget T ...