51nod 1674 区间的价值 V2

1674 区间的价值 V2

首先我们要知道对一个数一直and或者or操作减少或增加的次数是log级别的，这个看成二进制之后就很好思考。

这样我们马上就能想到枚举一个端点，然后维护log段相同的线段(and和or显然需要分开维护)。

考虑加入一个新的节点。
只需要把每个线段加入一个节点，合并值相同的线段。

答案的统计也是非常简单。
效率O(n∗2logX)

#include<cstdio>
#include<algorithm>
using namespace std;
const int N=100005,M=25,P=1000000007;
int a[N],ans,rt,tt,n;
struct node{
    int x,R;
}d[M],r[M];
inline int read(int& a)
{
    a=0;char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch>='0'&&ch<='9') a=a*10+ch-48,ch=getchar();
}
void work(int x)
{
    node z[N];
    int num;
    num=tt;for (int i=1;i<=tt;i++) z[i]=d[i];tt=1;d[1].x=a[x];d[1].R=x;
    for (int i=1;i<=num;i++)
    {
        if ((z[i].x&a[x])!=d[tt].x) d[++tt].x=z[i].x&a[x];
        d[tt].R=z[i].R;
    }
    num=rt;
    for (int i=1;i<=rt;i++) z[i]=r[i];rt=1;r[1].x=a[x];r[1].R=x;
    for (int i=1;i<=num;i++)
    {
        if ((z[i].x|a[x])!=r[rt].x) r[++rt].x=z[i].x|a[x];
        r[rt].R=z[i].R;
    }
    int i=1,j=1,last=x;
    while (i<=tt)
    {
        int now=min(d[i].R,r[j].R);
        ans=(ans+(long long)d[i].x*(now-last+1)%P*r[j].x%P)%P;
        i+=(d[i].R==now);j+=(r[j].R==now);last=now+1;
    }
}
int main()
{
    read(n);
    for (int i=1;i<=n;i++) read(a[i]);
    for (int i=n;i>=1;i--) work(i);
    printf("%d\n",ans);
    return 0;
}