BZOJ4750 密码安全

题目链接：http://www.lydsy.com/JudgeOnline/problem.php?id=4750

题意就是……像题目说的那样，中文题比较好理解。

之前肠胃炎，一个多星期没写题后写的第一道题。

首先我们可以用单调栈把每个数字以它为最大值的区间求出来，然后由于位运算每位之间不互相联系的特点，我们把每个数字分解为二进制，然后再观察一下发现异或有一些可以使用的规律。

比如1 0 1 0 1 0 1 加粗的1是这个区间的最大值，称它为中点吧，我们先计算从中点左边第一个数字往右边一个一个异或，然后算有x个1和y个0，接着从中点开始往右边异或，如果是1的话 ，那么这一位的统计数加上（y+1）(在纸上算一下就可以)，如果是0的话应该加上（x），计算完后这一位要乘上这个区间的最大值再乘上这是第几位（也就是2的多少次方）。最后就是答案了。

#include<cmath>
#include<algorithm>
#include<cstring>
#include<string>
#include<set>
#include<map>
#include<time.h>
#include<cstdio>
#include<vector>
#include<stack>
#include<queue>
#include<iostream>
using namespace std;
#define  LONG long long
const int   INF=0x3f3f3f3f;
const LONG  MOD=1e9+61;
const double PI=acos(-1.0);
#define clrI(x) memset(x,-1,sizeof(x))
#define clr0(x) memset(x,0,sizeof x)
#define clr1(x) memset(x,INF,sizeof x)
#define clr2(x) memset(x,-INF,sizeof x)
#define EPS 1e-10
LONG A[100010] ;
LONG Pre[34][100100] ;
int  p ;
struct Stac
{
    int l , r ;
    LONG val ;
}stac[200010];
int one[200100];
void Push(int x)
{
    while(1)
    {
        if(p <= 0) break ;
        int y = one[p] ;
        if( stac[y].val > stac[x].val )
            break ;
        stac[y].r = x - 1;
        stac[x].l = stac[y].l;
        p -- ;
    }
    one[++p] = x ;
}
int main()
{
    int T;
    cin>>T;
    while(T--)
    {
        p = 0;
        int n ;
        cin>>n;
        for(int i = 1 ; i<= n ;++i)
            scanf("%lld",&A[i]),
            stac[i].val = A[i] ,
            stac[i].l = stac[i].r = i ;
        stac[n+1].val = INF;
        stac[n+1].l = stac[n+1].r = n + 1;
        for(int i = 0 ; i <= 31 ;++ i)
        {
            Pre[i][0] = 0;
            for(int j = 1 ; j <= n ;++ j)
                Pre[i][j] = Pre[i][j-1] ^ ( (A[j]>>i) & 1 ) ;
        }
        for(int i =1 ; i<= n + 1; ++ i)
            Push( i ) ;
        LONG ans = 0;
        for(int i = 1 ; i<= n ; ++ i )
        {
            LONG res = 0;
            for(int k = 0 ; k <= 31 ; ++ k)
            {
                LONG tmp1 = 0, tmp2 = 0;
                LONG cnt = 0;
                for(int j = stac[i].l ; j <= i - 1; ++ j)
                {
                    if(Pre[k][i - 1] ^ Pre[k][j-1])tmp1 ++ ;
                    else tmp2 ++ ;
                }
                for(int j = i ; j<= stac[i].r ; ++ j)
                    if(Pre[k][j]^Pre[k][i-1]) cnt += (tmp2 + 1) ;
                    else cnt += tmp1  ;
                res += ((cnt% MOD ) *(1ll<<k) % MOD ) * stac[i].val;
                res %=MOD;
            }
            ans +=res;
            ans %= MOD;
        }
        cout<<ans<<endl;
    }
}