[ ST表 ] Codeforces875D High Cry

最新推荐文章于 2024-02-05 23:05:03 发布

gjghfd

最新推荐文章于 2024-02-05 23:05:03 发布

阅读量253

点赞数

分类专栏： ST表文章标签： codeforces

本文链接：https://blog.csdn.net/gjghfd/article/details/78265864

版权

ST表专栏收录该内容

3 篇文章 0 订阅

订阅专栏

枚举左端点，注意到区间或值只会变化 $O(\log n)$ 次，于是可以二分出每个或值相等的区间，再二分出最大值小于或值的区间。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
#define N 200010
inline char nc(){
    static char buf[100000],*p1=buf,*p2=buf;
    return p1==p2&&(p2=(p1=buf)+fread(buf,1,100000,stdin),p1==p2)?EOF:*p1++;
}
inline void Read(int& x){
    char c=nc();
    for(;c<'0'||c>'9';c=nc());
    for(x=0;c>='0'&&c<='9';x=(x<<3)+(x<<1)+c-48,c=nc());
}
int i,j,k,n,m;
int c[N],lg2[N];
int a1[N][35],a2[N][35],b1[N][35],b2[N][35],p[40];
long long Ans;
inline int Get_max(int x,int y){
    int k=lg2[y-x+1];
    return max(a1[x][k],a2[y][k]);
}
inline int Get_or(int x,int y){
    int k=lg2[y-x+1];
    return (b1[x][k]|b2[y][k]);
}
inline int Find(int x,int l){
    int r=n,y=Get_or(x,l);
    while(l<=r){
        int Mid=l+r>>1;
        if(Get_or(x,Mid)==y)l=Mid+1;else r=Mid-1;
    }
    return r;
}
inline int Calc(int x,int l,int r){
    int y=Get_or(x,l),L=l;
    while(l<=r){
        int Mid=l+r>>1;
        if(Get_max(x,Mid)<y)l=Mid+1;else r=Mid-1;
    }
    return r-L+1;
}
int main(){
    Read(n);
    for(i=1;i<=n;i++)Read(c[i]),a1[i][0]=a2[i][0]=b1[i][0]=b2[i][0]=c[i];
    for(i=2;i<=n;i++)lg2[i]=lg2[i>>1]+1;
    for(p[0]=i=1;i<=lg2[n];i++)p[i]=p[i-1]<<1;
    for(j=1;j<=lg2[n];j++)
    for(i=1;i<=n;i++){
        if(i+p[j-1]<=n)a1[i][j]=max(a1[i][j-1],a1[i+p[j-1]][j-1]),b1[i][j]=b1[i][j-1]|b1[i+p[j-1]][j-1];
        if(i-p[j-1]>=1)a2[i][j]=max(a2[i][j-1],a2[i-p[j-1]][j-1]),b2[i][j]=b2[i][j-1]|b2[i-p[j-1]][j-1];
    }
    for(i=1;i<n;i++){
        for(j=i+1;j<=n;j=k+1){
            k=Find(i,j);
            Ans+=Calc(i,j,k);
        }
    }
    cout<<Ans<<endl;
    return 0;
}