bzoj 1113: [Poi2008]海报PLA 分治

题意

N个矩形,排成一排. 现在希望用尽量少的矩形海报Cover住它们.
1<=n<=250000

分析

看到题我的第一想法就是分治。。。
每次在区间[l,r]中找到一个最小的，然后左右递归即可。
复杂度由于要求rmq，就是nlogn的。
还有O(n)的做法，就是用单调栈来做。

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
#include<cmath>
#define N 250005
using namespace std;

int n,a[N],a1,rmq[N][25];

void get_rmq()
{
    int lg=log(n)/log(2);
    for (int i=1;i<=n;i++) rmq[i][0]=i;
    for (int j=1;j<=lg;j++)
        for (int i=1;i<=n-(1<<j)+1;i++)
            rmq[i][j]=a[rmq[i][j-1]]<a[rmq[i+(1<<(j-1))][j-1]]?rmq[i][j-1]:rmq[i+(1<<(j-1))][j-1];
}

int get_min(int l,int r)
{
    int lg=log(r-l+1)/log(2);
    return a[rmq[l][lg]]<a[rmq[r-(1<<lg)+1][lg]]?rmq[l][lg]:rmq[r-(1<<lg)+1][lg];
}

int solve(int l,int r,int h)
{
    if (l>r) return 0;
    if (l==r) return 1;
    int x=get_min(l,r);
    if (a[x]==h) return solve(l,x-1,h)+solve(x+1,r,h);
    else return solve(l,x-1,a[x])+solve(x+1,r,a[x])+1;
}

int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++)
    {
        int x;
        scanf("%d%d",&x,&a[++a1]);
        if (a[a1]==a[a1-1]) a1--;
    }
    n=a1;
    get_rmq();
    printf("%d",solve(1,n,0));
    return 0;
}