Greater and Greater

题目

传送门

题目大意

给定大小为n的序列A和大小为m的序列B,计算A中所有大小为m的子区间S,满足
在这里插入图片描述

分析

本题使用了一个special的STL:bitset
考虑bitset,对每个A求一个长为m的bitset Si,其中Si[j]=1当且仅当Ai≥Bj。注意到本质只有O(m)种不同的bitset,具体就是把m个数排完序之后,第i个bitset就在第i-1个bitset的基础上在第i大的数对应的位置多一个1,所以预处理这些bitset复杂度是O(m2/w)的。
再设n个长为m的bitset,记为curi,其中 curi[j]=1当且仅当∀k∈[ j,m],Ai+k-j ≥Bk
可知有:
curi= ((curi+1 >>1)|Im)& Si
其中Im表示一个只有第m位为1的bitset。
所以只要统计有多少个curi的curi[1]=1即可,这就是最后的答案。写代码的时候可以
不用把每个curi都存起来,只弄一个cur然后不断滚动就可以了。
时间复杂度:O(nm/W)

代码

#include

using namespace std;

const int MAXN=2e5+10,MAXM=4e4+10;
int n,m,ans,a[MAXN],b[MAXM],pos[MAXN];
bitset<MAXM> s[MAXM],cur,w;

bool cmp(int x,int y){return b[x]<b[y];}

int main()
{
    scanf("%d%d",&n,&m);w[m]=1;
    for(int i=1;i<=n;i++) scanf("%d",a+i);
    for(int i=1;i<=m;i++) scanf("%d",b+i),pos[i]=i;
    sort(pos+1,pos+1+m,cmp);
    sort(b+1,b+1+m);
    for(int i=1;i<=m;i++) s[i]=s[i-1],s[i][pos[i]]=1;
    for(int i=n;i>=1;i--)
	{
        int d=upper_bound(b+1,b+1+m,a[i])-b-1;
        cur=(((cur>>1)|w)&s[d]);
        if(cur[1]) ans++;
    }
    printf("%d\n",ans);
}

你可能感兴趣的:(数据结构,算法)