先考虑如何确定对于某个 i i i 的答案。
题目的要求是:对于每一个 i = 1 , 2 , ⋯ , n − 1 i=1,2,\cdots,n-1 i=1,2,⋯,n−1,找一个最小的 j j j 使得 j > i j>i j>i 且线段 P j P j + 1 P_jP_{j+1} PjPj+1(含端点)与将射线 P i P i + 1 P_iP_{i+1} PiPi+1 略向上平移后所得的射线相交。
这个要求可以转化为:拟一条射线 P i P i + 1 P_iP_{i+1} PiPi+1,找到最小的 j j j 使得 P j P_j Pj 在这条射线上方(不含射线),那么 P j − 1 P j P_{j-1}P_j Pj−1Pj 就是我们要找的答案。
原因也很简单,由于 j j j 是最小的,所以对于任意的 k ∈ [ i + 1 , j − 1 ] k\in [i+1,j-1] k∈[i+1,j−1],都有 P k P_{k} Pk 在这条射线下方(含射线),于是就能保证射线穿过线段 P j − 1 P j P_{j-1}P_j Pj−1Pj,在 P j − 1 P j P_{j-1}P_j Pj−1Pj 之前不存在射线能穿过的线段。
那么现在问题就变成了:对于每一个 i = 1 , 2 , ⋯ , n − 1 i=1,2,\cdots,n-1 i=1,2,⋯,n−1,找一个最小的 j j j 使得 j > i + 1 j>i+1 j>i+1 且点 P j P_j Pj 在 P i P i + 1 P_iP_{i+1} PiPi+1 上方(不含射线,下同)。
考虑二分,判断是否存在一个 j ≤ m i d j\leq mid j≤mid 满足条件,那么就需要判断点集 { P i + 2 , P i + 3 , ⋯ , P j } \{P_{i+2},P_{i+3},\cdots,P_j\} {Pi+2,Pi+3,⋯,Pj} 中是否存在点在射线上方。
看似我们并没有把问题简单化,甚至更加复杂化,因为貌似我们还是要枚举所有的点,甚至还加上了个二分。
但实际上我们能利用一个很好的性质:把二分的判断转化为判断点集 { P i + 2 , P i + 3 , ⋯ , P j } \{P_{i+2},P_{i+3},\cdots,P_j\} {Pi+2,Pi+3,⋯,Pj} 的上凸壳中是否存在点在射线上方,我们就能通过三分来 O ( log n ) O(\log n) O(logn) 判断了。
而且一段区间的点的上凸壳是不会随着询问改变而改变的,而且我们可以通过线段树来维护一段区间的点的上凸壳。
于是我们可以直接在线段树上二分:假设我们当前要在区间 [ l , r ] [l,r] [l,r] 寻找最小的 j j j,那我们用三分判断左儿子内有没有符合条件的 j j j,有的话就去左儿子找,否则去右儿子找。
总时间复杂度 O ( T n log 2 n ) O(Tn\log ^2n) O(Tnlog2n)。
#include<bits/stdc++.h>
#define N 100010
#define ll long long
using namespace std;
inline int read()
{
int x=0,f=1;
char ch=getchar();
while(ch<'0'||ch>'9')
{
if(ch=='-') f=-1;
ch=getchar();
}
while(ch>='0'&&ch<='9')
{
x=(x<<1)+(x<<3)+(ch^'0');
ch=getchar();
}
return x*f;
}
struct Point
{
int x,y;
Point(){};
Point(int xx,int yy){x=xx,y=yy;}
}p[N],sta[N],s,t;
Point operator - (Point a,Point b){return Point(a.x-b.x,a.y-b.y);}
ll cross(Point a,Point b){return 1ll*a.x*b.y-1ll*b.x*a.y;}
int n;
int top;
vector<Point>v[N<<2];
void insert(Point u)
{
while(top>1&&cross(sta[top]-sta[top-1],u-sta[top-1])>=0) top--;
sta[++top]=u;
}
void merge(int k)
{
top=0;
int lc=k<<1,rc=k<<1|1;
int l=0,r=0;
while(l<v[lc].size()&&r<v[rc].size())
{
if(v[lc][l].x<v[rc][r].x) insert(v[lc][l]),l++;
else insert(v[rc][r]),r++;
}
while(l<v[lc].size()) insert(v[lc][l]),l++;
while(r<v[rc].size()) insert(v[rc][r]),r++;
for(int i=1;i<=top;i++) v[k].push_back(sta[i]);
}
void build(int k,int l,int r)
{
v[k].clear();
if(l==r)
{
v[k].push_back(p[l]);
return;
}
int mid=(l+r)>>1;
build(k<<1,l,mid);
build(k<<1|1,mid+1,r);
merge(k);
}
bool check(int k)
{
int l=0,r=v[k].size()-1,ans=-1;
while(l<=r)
{
int mid=(l+r)>>1;
if(!mid||cross(t-s,v[k][mid-1]-s)<cross(t-s,v[k][mid]-s)) ans=mid,l=mid+1;
else r=mid-1;
}
return cross(t-s,v[k][ans]-s)>0;
}
int query(int k,int l,int r,int x)
{
if(x<=l)
{
if(!check(k)) return 0;
if(l==r) return l-1;
}
int mid=(l+r)>>1;
if(x<=mid)
{
int tmp=query(k<<1,l,mid,x);
if(tmp) return tmp;
}
return query(k<<1|1,mid+1,r,x);
}
int main()
{
int T=read();
while(T--)
{
n=read();
for(int i=1;i<=n;i++)
p[i].x=read(),p[i].y=read();
build(1,1,n);
for(int i=1;i<n-1;i++)
{
s=p[i],t=p[i+1];
printf("%d ",query(1,1,n,i+2));
}
puts("0");
}
return 0;
}