SRM 562 Div1 500 CheckerFreeness

/*
通过枚举一对黑点(l,r)。
我们可以计算出其他所有白点和他们的相对位置(用叉积表示)
我们把剩余的白点以和直线(l,r)的左右位置分成两部分
我们用pos表示在线段l -> r左侧的点，其相对于l，相对于r，分别是顺时针第pos[i].fr, pos[i].sc个点。
lim表示在l -> r右侧的点，与它，选中的l,r，能形成凸多边形的左侧点，需要满足fr  >= lim[i].fr, sc <= lim[i].sc
然后我们对pos和lim进行排序
之后我们利用树状数组,把pos中fr>=lim[i].fr的不断sc加入树状数组后,如果从树状数组中找到了对于lim[i]合法的sc时就得到了答案
*/
#include<bits/stdc++.h>
#define cl(x) memset(x, 0, sizeof(x));
#define pii pair<int, int>
#define fr first
#define sc second
#define M 250
typedef long long ll;
using namespace std;
int n,m;
bool side[M]; // side[i] = 1,表示第i个点在线段l -> r的左侧
struct poi{
int x, y, id;
poi(){}
poi(int _x, int _y) { x = _x; y = _y; }
poi operator + (const poi &b) { return poi(x + b.x, y + b.y); }
poi operator - (const poi &b) { return poi(x - b.x, y - b.y); }
ll  operator * (const poi &b) { return 1ll * x * b.y - 1ll * y * b.x; } // 叉积
}W[M],B[M],tmp[M]; // 存白点相对于被选中的两个黑点的坐标;
bool operator < (poi a, poi b) {
return a * b < 0;
}
pii pos[M], lim[M];
//pos表示在线段l -> r左侧的点，其相对于l，相对于r，分别是顺时针第pos[i].fr, pos[i].sc个点。
//lim表示在l -> r右侧的点，与它，选中的l,r，能形成凸多边形的左侧点，需要满足fr  >= lim[i].fr, sc <= lim[i].sc
int c[M]; // 树状数组
void add(int x) { // 单点加
for (int i = x; i <= n; i += (i & (-i)))++c[i];
}
int get(int x) { // 前缀和
int ret = 0;
for (int i = x; i; i -= (i & (-i))) ret += c[i];
return ret;
}
bool cmp(const pii &a, const pii &b) {
return a.fr == b.fr ? a.sc < b.sc : a.fr > b.fr;
}
bool check(poi l, poi r) { // 选择l,r两个黑点为对角上的黑点
int k = 0, p = 0;
for (int i = 1; i <= n; ++i)
side[i]=((W[i] - r) * (l - r) > 0ll); // 判断白点在两个黑点的哪一侧
for (int i = 1; i <= n; ++i)
if (side[i]) tmp[++k] = W[i] - l, tmp[k].id = k; // 计算相对l的极角
sort(tmp + 1, tmp + k + 1); // 因为都在pi的范围内，所以可以直接用叉积来排序
for (int i = 1; i <= k; ++i) pos[tmp[i].id].fr = i; // 确定下相对l极角序的排名
for (int i = 1; i <= n; ++i)
if (!side[i])lim[++p].fr = upper_bound(tmp + 1, tmp + k + 1, l - W[i]) - tmp; // 确定下对应的左侧点，相对l的极角序至少要排在第几位

k = 0; p = 0; // 同理计算相对r的极角序
for (int i = 1; i <= n; ++i)
if (side[i]) tmp[++k] = W[i] - r, tmp[k].id = k;
sort(tmp + 1, tmp + k + 1);
for (int i = 1; i <= k; ++i) pos[tmp[i].id].sc = i;
for (int i = 1; i <= n; ++i)
if (!side[i])lim[++p].sc = lower_bound(tmp + 1, tmp + k + 1, r - W[i]) - tmp - 1;

cl(c);
sort(pos + 1, pos + k + 1, cmp); int tp = 1;
sort(lim + 1, lim + p + 1, cmp); // 都根据两维关键字排序
for (int i = 1; i <= p; ++i) {
while (tp <= k && pos[tp].fr >= lim[i].fr) // 把fr这维满足条件的左侧点加入树状数组的sc位置
add(pos[tp++].sc);
if (get(lim[i].sc)) return 1; // 存在满足条件的点
}return 0;
}
int main(){
scanf("%d",&n);
for (int i = 1; i <= n; ++i) scanf("%d",&W[i].x);
for (int i = 1; i <= n; ++i) scanf("%d",&W[i].y);
scanf("%d",&m);
for (int i = 1; i <= m; ++i) scanf("%d",&B[i].x);
for (int i = 1; i <= m; ++i) scanf("%d",&B[i].y);
for (int i = 1; i <= m; ++i)
for (int j = i + 1; j <= m; ++j)
if (check(B[i], B[j])){puts("NO");return 0;}
puts("YES");
return 0;
}