2017.10.17 CF#441 F题 思考记录

。果然没有题解做题感觉是不一样的。

这个F题可能是最像往常B题的题了（往常B题就是dp，还一般都不难）

首先这个是|，所以|的数越多，它一定是不降的

那么合法的一定是一个数到一个数往后的数所构成的所有区间

那么不停地往后|，直到>其中最大的一个值就可以了，但最大值是随|不断变化的，那最大值是多少？

这时可以考虑枚举（化动为静的思想），，当前作为最大值的区间最大是多少（单调栈，前几天刚做过）

然后剩下的就是在左右区间各找一个位置，使 |和 >枚举的最大值

显然其他数位数不可能比最大值大，不然他就成为最大值了

所以考虑把最大值按位取反，然后&这个最大值>0的第一个出现的数

&>0只有最大值这一位是0，而有一个数这一位是1时出现

所以可以用线段树维护这一位是1的最左边的位置和最右边的位置，分别取max和min

因为具体数字满足几个要求无所谓，所以可以直接按位离散，

但线段树会爆空间，所以可以用dp求

 f【这个数是啥】【二进制位】=最靠右的位置  

g【这个数是啥】【二进制位】=最靠左的位置  

半个插数dp。。

然后要注意，数字相同的区间也要算一个，所以可以往后算，

码：

#include<iostream>
#include<cstdio>
using namespace std;
#define zzzz o<<1,l,mid
#define yyyy o<<1|1,mid+1,r
#define N  200005
#define ll long long
int zz[N*32],zy[N*32],c,a,b,er[35],top,n,op,aa[N],sta[N],tong[35],f[N][35],g[N][35],lin[N];
ll zuo[N],you[N],ans,i,j;
bool nozuizuo[N];
int main()
{
	scanf("%d",&n);
	er[0]=1;
	for(i=1;i<=30;i++)
	er[i]=er[i-1]*2;
	for(i=1;i<=n;i++)
	{
		scanf("%d",&aa[i]);
	}
	top=0;
	for(i=1;i<=n;i++)
	{
	zuo[i]=1;
	while(top&&aa[sta[top]]<aa[i])zuo[i]+=zuo[sta[top]],top--;
	sta[++top]=i;	
	}
	top=0;
		for(i=n;i>=1;i--)
	{
	you[i]=1;
	while(top&&aa[sta[top]]<aa[i])you[i]+=you[sta[top]],top--;
	if(aa[sta[top]]==aa[i])nozuizuo[sta[a]]=1;	
	sta[++top]=i;	
	}
	top=0;
		for(i=n;i>=1;i--)
	{
	lin[i]=1;
	while(top&&aa[sta[top]]<=aa[i])
	lin[i]+=you[sta[top]],top--;
	if(nozuizuo[i]!=1)you[i]=lin[i];	
	sta[++top]=i;	
	}
	for(i=1;i<=n;i++)
	{		
		for(j=0;j<=30;j++)
		{
	f[i][j]=tong[j];					
		}		
		for(j=0;j<=30;j++)
		{
	if(er[j]&aa[i])
	tong[j]=i;						
		}
	}
	for(i=0;i<=30;i++)tong[i]=9999999;
		for(i=n;i>=1;i--)
	{		
		for(j=0;j<=30;j++)
		{
	g[i][j]=tong[j];					
		}		
		for(j=0;j<=30;j++)
		{
	if(er[j]&aa[i])
	tong[j]=i;						
		}
	}
	for(i=1;i<=n;i++)
	{    int lin=aa[i],o=0,j=0;	
		ll l1=0,l2=9999999;
		 while(lin)
		 {
		 	if(!(lin&1))
		 	{
		 if(f[i][j]>=i-zuo[i]+1)l1=max(f[i][j],int(l1));	
		 if(g[i][j]<=you[i]+i-1)l2=min(g[i][j],int(l2));	 			 		
		    }
lin/=2;
++j;
		 }
if(l1!=0&&l2!=9999999)
{
	ans+=(l1-(i-zuo[i]+1)+1)*(you[i]);
	ans+=(i-l1)*((you[i]+i-1)-l2+1);
}
if(l1!=0&&l2==9999999)
{
	ans+=(l1-(i-zuo[i]+1)+1)*(you[i]);
}
if(l1==0&&l2!=9999999)
{
	ans+=zuo[i]*((i+you[i]-1)-l2+1);
}	
}
printf("%I64d",ans);	
}