最近被四区题暴虐。。。
题意:lyk拥有一个区间。
它规定一个区间的价值为这个区间中所有数and起来的值与这个区间所有数or起来的值的乘积。
例如3个数2,3,6。它们and起来的值为2,or起来的值为7,这个区间对答案的贡献为2*7=14。
现在lyk有一个n个数的序列,它想知道所有n*(n+1)/2个区间的贡献的和对1000000007取模后的结果是多少。
区间的and值和区间的or值相乘,实际上等于将and值分解为2的幂次和的形式与or值分解成2的幂次和的形式相乘。
所以对于同一段区间来说,是可以按位来统计的。
首先将数列按位分解成32个位数列。
枚举区间的左端点,,区间的and值要对答案有贡献必须为1,随着右端点右移,and值为不递增的,而区间的or值要对答案有贡献必须为1,随着右端点右移,or值为不递减的。
找到这两个满足条件的最大区间的交集,即可统计出这段区间对答案的贡献。
先用尺取法进行预处理,再依此统计答案。
时间复杂度O(n*loga*loga)。
![](https://images.cnblogs.com/OutliningIndicators/ContractedBlock.gif)
# include <cstdio> # include <cstring> # include <cstdlib> # include <iostream> # include <vector> # include <queue> # include <stack> # include <map> # include <bitset> # include <set> # include <cmath> # include <algorithm> using namespace std; # define lowbit(x) ((x)&(-x)) # define pi acos(-1.0) # define eps 1e-8 # define MOD 1000000007 # define INF 1000000000 # define mem(a,b) memset(a,b,sizeof(a)) # define FOR(i,a,n) for(int i=a; i<=n; ++i) # define FO(i,a,n) for(int i=a; i<n; ++i) # define bug puts("H"); # define lch p<<1,l,mid # define rch p<<1|1,mid+1,r # define mp make_pair # define pb push_back typedef pair<int,int> PII; typedef vector<int> VI; # pragma comment(linker, "/STACK:1024000000,1024000000") typedef long long LL; inline int Scan() { int x=0,f=1; char ch=getchar(); while(ch<'0'||ch>'9'){if(ch=='-') f=-1; ch=getchar();} while(ch>='0'&&ch<='9'){x=x*10+ch-'0'; ch=getchar();} return x*f; } inline void Out(int a) { if(a<0) {putchar('-'); a=-a;} if(a>=10) Out(a/10); putchar(a%10+'0'); } const int N=100005; //Code begin... int R[32][N], L[32][N], n, mod[65]; LL ans; bool wei[32][N]; void init(){ mod[0]=1; FOR(i,1,64) mod[i]=(mod[i-1]<<1)%MOD; FO(i,0,32) { int l=1, r=1; while (l<=n) { if (!wei[i][l]) L[i][l]=0; else { r=max(r,l); while (r<n&&wei[i][r+1]) ++r; L[i][l]=r; } ++l; } } FO(i,0,32) { int l=1, r=1; while (l<=n) { if (wei[i][l]) R[i][l]=l; else { r=max(l,r); while (r<=n&&!wei[i][r]) ++r; R[i][l]=r; } ++l; } } } void sol(){ FO(i,0,32) FOR(j,1,n) { if (!wei[i][j]) continue; FO(k,0,32) { if (R[k][j]>L[i][j]) continue; ans=(ans+(LL)(L[i][j]-R[k][j]+1)*mod[i+k])%MOD; } } } int main () { int x; n=Scan(); FOR(i,1,n) { x=Scan(); FO(j,0,32) wei[j][i]=x%2, x/=2; } init(); sol(); printf("%lld ",ans); return 0; }