近似有序区间 51Nod - 1249

https://www.51nod.com/onlineJudge/questionCode.html#!problemId=1249

单调栈找出所有最大值区间以及最小值区间 所谓近似有序区间一定是一个最大一个最小区间相互包含构成的

假设某近似有序区间为[i,j] 以a[i]为最小值的区间为[ minlef[i],minrgt[i] ],以a[j]为最大值的区间为[ maxlef[j],maxrgt[j] ] 则必然有maxlef[j]<=i<=j<=minrgt[i] 如果说maxlef[j]>i 则[ i,maxlef[j]-1 ]之间必然有一个值是严格大于a[j]的 因为单调栈就在这划得界限了 

然后就是线段树离线处理一下线段 主席树在线也可以

#include <bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=5e4+10;

struct node
{
    int l,r,m;
};

node seg[maxn];
stack <int> stk;
ll sum[4*maxn];
int ary[maxn],maxlef[maxn],maxrgt[maxn],minlef[maxn],minrgt[maxn];
int n;

bool cmp(node n1,node n2)
{
    return n1.r<n2.r;
}

void init()
{
    int i;
    for(i=1;i<=n;i++){
        while(!stk.empty()&&ary[stk.top()]<=ary[i]) stk.pop();
        if(stk.empty()) maxlef[i]=1;
        else maxlef[i]=stk.top()+1;
        stk.push(i);
    }
    while(!stk.empty()) stk.pop();

    for(i=n;i>=1;i--){
        while(!stk.empty()&&ary[stk.top()]<=ary[i]) stk.pop();
        if(stk.empty()) maxrgt[i]=n;
        else maxrgt[i]=stk.top()-1;
        stk.push(i);
    }
    while(!stk.empty()) stk.pop();

    for(i=1;i<=n;i++){
        while(!stk.empty()&&ary[stk.top()]>=ary[i]) stk.pop();
        if(stk.empty()) minlef[i]=1;
        else minlef[i]=stk.top()+1;
        stk.push(i);
    }
    while(!stk.empty()) stk.pop();

    for(i=n;i>=1;i--){
        while(!stk.empty()&&ary[stk.top()]>=ary[i]) stk.pop();
        if(stk.empty()) minrgt[i]=n;
        else minrgt[i]=stk.top()-1;
        stk.push(i);
    }
}

void update(int tar,int l,int r,int cur)
{
    int m;
    sum[cur]++;
    if(l==r) return;
    m=(l+r)/2;
    if(tar<=m) update(tar,l,m,2*cur);
    else update(tar,m+1,r,2*cur+1);
}

ll query(int pl,int pr,int l,int r,int cur)
{
    ll res;
    int m;
    if(pl<=l&&r<=pr) return sum[cur];
    res=0,m=(l+r)/2;
    if(pl<=m) res+=query(pl,pr,l,m,2*cur);
    if(pr>m) res+=query(pl,pr,m+1,r,2*cur+1);
    return res;
}

int main()
{
    ll ans;
    int i,j;
    scanf("%d",&n);
    for(i=1;i<=n;i++) scanf("%d",&ary[i]);
    init();
    for(i=1;i<=n;i++){
        seg[i].l=minlef[i],seg[i].r=minrgt[i],seg[i].m=i;
    }
    sort(seg+1,seg+n+1,cmp);
    ans=0;
    for(i=n,j=n;i>=1;i--){
        while(j>=1&&seg[j].r>=i){
            update(seg[j].m,1,n,1);
            j--;
        }
        ans+=query(maxlef[i],i,1,n,1);
    }
    printf("%lld\n",ans);
    return 0;
}

 

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值