2017.10.17 CF#441 F题 思考记录

来源:互联网 发布:网络词ip是什么意思 编辑:程序博客网 时间:2024/06/05 20:20

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

这个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 longint 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);}





原创粉丝点击