关闭

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

标签: OI凸壳单调栈
788人阅读 评论(0) 收藏 举报
分类:

题目大意

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

1n200000,xi106,yi106


题目分析

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

时间复杂度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;
}
0
0

查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:106189次
    • 积分:3481
    • 等级:
    • 排名:第10000名
    • 原创:228篇
    • 转载:0篇
    • 译文:0篇
    • 评论:60条
    博客专栏
    OI

    文章:14篇

    阅读:13981
    最新评论
    文章分类