CodeForces 406 D.Hill Climbing （凸包-Graham扫描法+在线倍增LCA）

Description
有n个山坡，第i个山坡顶端坐标为(x[i],y[i])，一个山坡可以到另一个山坡当且仅当从另一个山坡可以看到这个山坡，现在给出q次查询，每次查询给出两个人所在山坡编号，如果左边的人可以到达右边人的山坡那右边人就不动等着左边人来，否则右边人往右到一个可以达到的最右边的山坡继续上面的步骤吗，问两个人最后停留的位置
Input
第一行一整数n表示山坡数量，之后n行每行两个整数x[i],y[i]表示第i个山坡的山顶坐标，之后输入一整数q表示查询数，每次查询输入两个整数a和b表示初始状态两个人所处的山坡编号(1<=n<=1e5,1<=x[i]<=1e7,1<=y[i]<=1e11,1<=q<=1e5,1<=a,b<=n)
Output
对于每次查询，输出两个人最终停留的位置
Sample Input
6
1 4
2 1
3 2
4 3
6 4
7 4
3
3 1
5 6
2 3
Sample Output
5 6 3
Solution
一个人往右走到可达的最右边等价于到达该点所处凸包的下一个顶点，从右往左维护一个凸包，对凸包上相邻顶点建边可以得到一棵树，每次查询其实是求a和b两点在树上的LCA，在线倍增法即可
Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
#define INF 0x3f3f3f3f
#define maxn 111111
int n,q,s[maxn];
struct point
{
    ll x,y;
    point operator -(const point &b)const
    {
        point c;
        c.x=x-b.x,c.y=y-b.y;
        return c;
    }
}a[maxn];
bool check(point a,point b)
{
    return a.x*b.y-b.x*a.y>0;
}
vector<int>g[maxn];
int p[maxn][22],deep[maxn],size[maxn];
void init(int u)
{
    size[u]=1;
    for(int i=0;i<g[u].size();i++)
    {
        int v=g[u][i];
        deep[v]=deep[u]+1;
        p[v][0]=u;
        init(v);
        size[u]+=size[v];
    }
}
void lca_init(int n,int root)
{
    deep[root]=1;
    memset(p,-1,sizeof(p));
    init(root);
    for(int j=1;(1<<j)<=n;j++)
        for(int i=1;i<=n;i++)
            if(~p[i][j-1])
                p[i][j]=p[p[i][j-1]][j-1];
}
int lca(int a,int b)
{
    int i,j;
    if(deep[a]<deep[b])swap(a,b);
    for(i=0;(1<<i)<=deep[a];i++);
    i--;
    for(j=i;j>=0;j--)
        if(deep[a]-(1<<j)>=deep[b])
            a=p[a][j];
    if(a==b) return a;
    for(j=i;j>=0;j--)
    {
        if(p[a][j]!=-1&&p[a][j]!=p[b][j])
        {
            a=p[a][j];
            b=p[b][j];
        }
    }
    return p[a][0];
}
int main()
{
    while(~scanf("%d",&n))
    {
        for(int i=1;i<=n;i++)scanf("%I64d%I64d",&a[i].x,&a[i].y),g[i].clear();
        int t=0;
        for(int i=n;i>=1;i--)
        {
            while(t>=2&&check(a[s[t]]-a[i],a[s[t-1]]-a[i])>0)t--;
            s[++t]=i;
            if(t>=2)g[s[t-1]].push_back(s[t]);
        }
        lca_init(n,n);
        scanf("%d",&q);
        while(q--)
        {
            int x,y;
            scanf("%d%d",&x,&y);
            printf("%d\n",lca(x,y));
        }
    }
    return 0;
}