[JZOJ3401]Pty爬山/[2013年集训队论文答辩]登顶计划

题目大意

一座山，有$n$个山峰，用整点坐标$(x_i,y_i)$表示。
相邻两个山峰连起来形成一段山路。
一个人要爬山，他从某一个点出发，然后找到他能看到的最高的山峰（三点共线算作能看到，高度相同$x_i$大的算作更高），然后向那个山峰走去。每到一座山峰，他都会重新观察，如果当前看到更高的山峰，就转向那里走去。显然最后他会到达最高峰。
题目要求对于每一座山峰，求出以此为起始点，走多少段山路能到顶点。

$1\le n\le 200000,x_i\le10^6,y_i\le 10^6$

---

题目分析

首先要明确一点，我们从某座山峰出发，到达某一山峰，然后更改了目标，从这以后的路径，和开始时就从当前山峰出发是一样的。
用正反两次上凸壳可以求出任意一座山峰能看到的最高山峰，然后我们再用正反两次单调栈得到我们从每个点出发向可见最高山峰前进过程中会被哪个点打断（即找到更高峰，如果没有就设为目标山峰），将这个转折点设为开始点的父亲，就可以发现它形成了一个树结构。自上而下扫一遍即可出解（注意$\mathrm{DFS}$会爆栈）。

时间复杂度$\mathrm O(n)$。

---

代码实现

#include <iostream>
#include <cstdio>
#include <cctype>
#include <queue>
#include <cmath>

using namespace std;

typedef long long LL;

int read()
{
    int x=0,f=1;
    char ch=getchar();
    while (!isdigit(ch)) ch=='-'?f=-1:f,ch=getchar();
    while (isdigit(ch)) x=x*10+ch-'0',ch=getchar();
    return x*f;
}

const int N=200050;

struct P
{
    int x,y;

    P (int x0=0,int y0=0){x=x0,y=y0;}
}p[N];

P operator+(P a,P b){return P(a.x+b.x,a.y+b.y);}
P operator-(P a,P b){return P(a.x-b.x,a.y-b.y);}

LL product(P v1,P v2){return v1.x*1ll*v2.y-v1.y*1ll*v2.x;}

int st[N],aim[N],f[N],fa[N],last[N],next[N],tov[N];
queue<int> q;
int n,tot,rt;
LL dis[N];

void insert(int x,int y){tov[++tot]=y,next[tot]=last[x],last[x]=tot;}

void scan()
{
    st[0]=0;
    for (int i=1;i<=n;i++)
    {
        while ((st[0]>1&&product(p[st[st[0]]]-p[i],p[st[st[0]-1]]-p[i])<=0)||(st[0]&&p[st[st[0]]].y<=p[i].y)) st[0]--;
        aim[i]=st[st[0]],st[++st[0]]=i;
    }
    st[0]=0;
    for (int i=n;i>=1;i--)
    {
        while ((st[0]>1&&product(p[st[st[0]]]-p[i],p[st[st[0]-1]]-p[i])>=0)||(st[0]&&p[st[st[0]]].y<p[i].y)) st[0]--;
        if (p[aim[i]].y<=p[st[st[0]]].y) aim[i]=st[st[0]];
        st[++st[0]]=i;
    }
    for (int i=1;i<=n;i++) f[i]=p[aim[i]].y;
}

bool cmp(int i,int j)
{
    return f[i]<f[j]||f[i]==f[j]&&aim[i]<aim[j];
}

void pre()
{
    st[0]=0;
    for (int i=1;i<=n;i++)
    {
        if (!aim[i]) continue;
        while (st[0]&&cmp(st[st[0]],i)) st[0]--;
        if (aim[i]<i) fa[i]=(st[0]?(aim[i]<st[st[0]]?st[st[0]]:aim[i]):aim[i]);
        st[++st[0]]=i;
    }
    st[0]=0;
    for (int i=n;i>=1;i--)
    {
        if (!aim[i]) continue;
        while (st[0]&&cmp(st[st[0]],i)) st[0]--;
        if (aim[i]>i) fa[i]=(st[0]?(st[st[0]]<aim[i]?st[st[0]]:aim[i]):aim[i]);
        st[++st[0]]=i;
    }
    for (int i=1;i<=n;i++)
        if (fa[i]) insert(fa[i],i);
        else rt=i;
}

void solve()
{
    q.push(rt);
    int x;
    while (!q.empty())
    {
        x=q.front(),q.pop();
        for (int i=last[x],y;i;i=next[i])
            q.push(y=tov[i]),dis[y]=dis[x]+abs(x-y);
    }
}

int main()
{
    freopen("mountain.in","r",stdin),freopen("mountain.out","w",stdout);
    n=read();
    for (int i=1;i<=n;i++) p[i].x=read(),p[i].y=read();
    scan(),pre(),solve();
    for (int i=1;i<=n;i++) printf("%lld\n",dis[i]);
    fclose(stdin),fclose(stdout);
    return 0;
}