[JSOI2018]战争
题目描述
九条可怜是一个热爱读书的女孩子。
在她最近正在读的一本小说中,描述了两个敌对部落之间的故事。第一个部落有 nnn 个人,第二个部落有 mmm 个人,每一个人的位置可以抽象成二维平面上坐标为 (xi,yi)(x_i,y_i)(xi,yi) 的点。
在这本书中,人们有很强的领地意识,对于平面上的任何一个点,如果它被三个来自同一部落的人形成的三角形(可能退化成一条线段)包含(包括边界),那么这一个点就属于这一个部落的领地。如果存在一个点同时在两个阵营的领地中,那么这两个部落就会为了争夺这一个点而发生战争。
常年的征战让两个部落不堪重负,因此第二个部落的族长作出了一个英明的决定,他打算选择一个向量 (dx,dy)(dx,dy)(dx,dy) ,让所有的族人都迁徙这个向量的距离,即所有第二阵营的人的坐标都变成 (xi+dx,yi+dy)(x_i+dx,y_i+dy)(xi+dx,yi+dy) 。
现在他计划了 qqq 个迁徙的备选方案,他想要你来帮忙对每一个迁徙方案,计算一下在完成了迁徙之后,两个部落之间还会不会因为争夺领地而发生战争。
输入格式
第一行输入三个整数 n,m,qn,m,qn,m,q,表示两个部落里的人数以及迁徙的备选方案数。
接下来 nnn 行每行两个整数 xi,yix_i,y_ixi,yi 表示第一个部落里的人的坐标。
接下来 mmm 行每行两个整数 xi,yix_i,y_ixi,yi 表示第二个部落里的人的坐标。
接下来 qqq 行每行两个整数 dxi,dyidx_i,dy_idxi,dyi 表示一个迁徙方案。
输入数据保证所有人的坐标两两不同。
输出格式
对于每个迁徙方案,输出一行一个整数,000 表示不会发生冲突,111 表示会发生冲突。
输入输出样例
4 4 3
0 0
1 0
0 1
1 1
-1 0
0 3
0 2
0 -1
0 0
2 3
0 -1
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long lol;
struct Node
{
lol x,y;
Node operator + (const Node &b) const
{
return (Node){x+b.x,y+b.y};
}
Node operator - (const Node &b) const
{
return (Node){x-b.x,y-b.y};
}
}a[],b[],sa[],sb[],sta[],bs,s1[],s2[],s[];
int n,m,q,top;
bool cmp(Node a,Node b)
{
if (a.y==b.y) return a.x<b.x;
return a.y<b.y;
}
lol cross(Node a,Node b)
{
return (a.x*b.y-a.y*b.x);
}
lol dist(Node a)
{
return a.x*a.x+a.y*a.y;
}
bool cmp1(Node A,Node B)
{
lol t=cross((a[]-A),(a[]-B));
if (t==) return dist(a[]-A)<dist(a[]-B);
return t>;
}
bool cmp2(Node A,Node B)
{
lol t=cross((b[]-A),(b[]-B));
if (t==) return dist(b[]-A)<dist(b[]-B);
return t>;
}
int grahama(int N)
{int i;
sort(a+,a+N+,cmp);
sort(a+,a+N+,cmp1);
top=;
sa[++top]=a[];
if (N==) return ;
sa[++top]=a[];
for (i=;i<=N;i++)
{
while (top>&&cross(a[i]-sa[top-],sa[top]-sa[top-])>=) top--;
top++;
sa[top]=a[i];
}
sa[top+]=a[];
return top;
}
int grahamb(int N)
{int i;
sort(b+,b+N+,cmp);
sort(b+,b+N+,cmp2);
top=;
sb[++top]=b[];
if (N==) return ;
sb[++top]=b[];
for (i=;i<=N;i++)
{
while (top>&&cross(b[i]-sb[top-],sb[top]-sb[top-])>=) top--;
top++;
sb[top]=b[i];
}
sb[top+]=b[];
return top;
}
bool cmpp(Node A,Node B)
{
lol t=cross((s[]-A),(s[]-B));
if (t==) return dist(s[]-A)<dist(s[]-B);
return t>;
}
int grahams(int N)
{int i;
sort(s+,s+N+,cmp);
sort(s+,s+N+,cmpp);
top=;
sta[++top]=s[];
if (N==) return ;
sta[++top]=s[];
for (i=;i<=N;i++)
{
while (top>&&cross(s[i]-sta[top-],sta[top]-sta[top-])>=) top--;
top++;
sta[top]=s[i];
}
sta[top+]=s[];
return top;
}
bool cmp3(Node A,Node B)
{
return cross(A,B)>||(cross(A,B)==&&dist(A)<dist(B));
}
lol find(Node A)
{
if(cross(A,sta[])>||cross(sta[top],A)>) return ;
lol ps=lower_bound(sta+,sta+top+,A,cmp3)-sta-;
return cross((A-sta[ps]),(sta[ps%top+]-sta[ps]))<=;
}
void Minkowski()
{
for(lol i=;i<n;i++) s1[i]=sa[i+]-sa[i];s1[n]=sa[]-sa[n];
for(lol i=;i<m;i++) s2[i]=sb[i+]-sb[i];s2[m]=sb[]-sb[m];
top=;
s[top]=sa[]+sb[];
lol i=,j=;
while(i<=n&&j<=m) ++top,s[top]=s[top-]+(cross(s1[i],s2[j])>=?s1[i++]:s2[j++]);
while(i<=n) ++top,s[top]=s[top-]+s1[i++];
while(j<=m) ++top,s[top]=s[top-]+s2[j++];
}
int main()
{int i,j;
lol dx,dy;
cin>>n>>m>>q;
for (i=;i<=n;i++)
scanf("%lld%lld",&a[i].x,&a[i].y);
for (i=;i<=m;i++)
scanf("%lld%lld",&b[i].x,&b[i].y),b[i].x=-b[i].x,b[i].y=-b[i].y;
n=grahama(n);
m=grahamb(m);
Minkowski();
top=grahams(top);
bs=sta[];
for (i=;i<=top;i++)
sta[i]=sta[i]-bs; for (i=;i<=q;i++)
{
scanf("%lld%lld",&dx,&dy);
if (find((Node){dx,dy}-bs))
printf("1\n");
else printf("0\n");
}
}
[JSOI2018]战争的更多相关文章
- P4557 [JSOI2018]战争
首先可以题目描述的两个点集是两个凸包,分别设为A和B. 考虑一个向量w不合法的条件. 即存在b+w=a,其中a属于A,b属于B. 也就是a-b=w. 即对b取反后和a的闵可夫斯基和. 求出闵可夫斯基和 ...
- BZOJ5317:[JSOI2018]战争(闵可夫斯基和)
令 \(a\in A,b\in B\) 则移动向量 \(\omega\) 使得存在 \(b+\omega=a\) 那么 \(\omega\) 需要满足 \(\omega=a−b\) 黑科技:闵可夫斯基 ...
- 洛谷P4557 [JSOI2018]战争(闵可夫斯基和+凸包)
题面 传送门 题解 看出这是个闵可夫斯基和了然而我当初因为见到这词汇是在\(shadowice\)巨巨的\(Ynoi\)题解里所以压根没敢学-- 首先您需要知道这个 首先如果有一个向量\(w\)使得\ ...
- [JSOI2018]战争(闵可夫斯基和)
害怕,可怜几何题 果然不会 题目就是说给你两个凸包,每次询问给你一个向量 \(c\) 问你能不能从两个凸包 \(A\) , \(B\) 里分别找到一个点 \(a\) , \(b\) 满足 \(a+c= ...
- 【LuoguP4557】[JSOI2018]战争
题目链接 题意 给你两个点集. q次询问 , 每次把其中一个点集往一个方向移动 , 问两个点集的凸包还有没有交. Sol 闵可夫斯基和板子题. 把问题做如下转换: 我们本来两个凸包相交是相当于是对于移 ...
- 计算几何细节梳理&模板
点击%XZY巨佬 向量的板子 #include<bits/stdc++.h> #define I inline using namespace std; typedef double DB ...
- HHHOJ #151. 「NOI模拟 #2」Nagisa
计算几何板子题(我才没有拷板子的说--) 众所周知,三角形的重心坐标是\((\frac{x_1+x_2+x_3}{3},\frac{y_1+y_2+y_3}{3})\) 然后我们发现如果我们有一个点集 ...
- 【学习笔记】Minkowski和
这还是个被我咕了N久的玩意 Minkowski和是一个奇怪的玩意 他长这样 $S={a+b \| a \in A , b \in B}$ AB可以是点集也可是向量集(显然) 他可以处理一些奇怪的东西 ...
- JSOI部分题解
JSOI部分题解 JSOI2018 战争 问题转化为给定你两个凸包\(\mathbb S,\mathbb T\),每次独立的询问将\(\mathbb T\)中的每个点移动一个向量,问\(\mathbb ...
随机推荐
- 【C++札记】标准输入与输出
概述 C语言中使用函数scanf和printf作为标准输入和输出,在C++中引入了类cin和cout进行标准输入和输出,所需头文件为<iostream>,命名空间是std.这里所说的标准输 ...
- Python基础系列讲解-自动控制windows桌面
原链接:https://zhuanlan.zhihu.com/p/73001806 在使用PC时与PC交互的主要途径是看屏幕显示.听声音,点击鼠标和敲键盘等等.在自动化办公的趋势下,繁琐的工作可以让程 ...
- tensorboard 拒绝访问解决方法
打开Anaconda Prompt,切换到TensorFlow环境(activate tensorflow) 切换成功之后,输入tensorboard --logdir='路径' 注意:--logdi ...
- quartz2.3.0(六)job任务异常处理方式
Job1类 package org.quartz.examples.example6; import org.quartz.DisallowConcurrentExecution; import or ...
- spring源码解析前瞻
很多人有疑问:为什么要读源码?读源码有什么用?我也一直问自己这些问题,读源码非常枯燥,工作中又用不到,慢慢的自己读源码越发现自己知识的不足,无法把知识串起来,形成知识体系.从单系统中常用的Spring ...
- diverta 2019 Programming Contest 2
A:签到. #include<bits/stdc++.h> using namespace std; #define ll long long #define inf 1000000010 ...
- ADO.NET 一般操作(c#)
ADO.NET五大对象:SqlConnection.SqlCommand.SqlDataReader.SqlDataAdapter .DataSet ,其中SqlDataAdapter 不建议使用 一 ...
- 第一个.NET小程序
一.用户需求 做一个简单的网页版销售合同签核系统 1.业务员需要在手机或者电脑上操作,Key入销售合同 2.业务员填入相应的合同信息,对应主管签核 3.最终签核完,生成PDF版的销售合同,且上面自动加 ...
- MSMQ消息加密
证书实现非对称加密/解密的代码如下 //非对称加密密钥 static byte[] RSAEncrypt(byte[] enkey, X509Certificate2 Certificate) { R ...
- 子进程的LD_PRELOAD
一个指定LD_PRELOAD的进程创建的子进程是否受LD_PRELOAD的影响? 1. fork()后在子进程中执行函数. main.c #include <unistd.h> #incl ...