2020牛客多校#2 G - Greater and Greater 动态规划+bitset优化

题目链接：
https://ac.nowcoder.com/acm/contest/5667/G

  思路与官方题解相同，只是描述方式不同。

  设$d{p}_{i,j}$表示子序列$[A_i,A_{i+m-j}]$与子序列$[B_j,B_m]$是否匹配，即$\forall k\in [0,m-j]，A_{i+k}\ge B_{j+k}$匹配则为$1$，不匹配则为$0$。可以想到答案就是$j=1$时子序列$[A_i,A_{i+m-1}]$与子序列$[B_1,B_m]$匹配的总数，即$$ans=\sum _{i=1}^{n-m+1}d{p}_{i,1}$$

  考虑$dp$的转移方程$$d{p}_{i,j}=d{p}_{i+1,j+1}\&(A_i\ge B_j)$$即若想子序列$[A_i,A_{i+m-j}]$与子序列$[B_j,B_m]$匹配，则需要子序列$[A_{i+1},A_{i+m-j}]$与子序列$[B_{j+1},B_m]$匹配，且$A_i\ge B_j$

  可以想到，这个转移的计算复杂度是$O(nm)$，约$6e9$，超时，也超空间。

  由于所有$dp$值都是0或1，可以采用bitset优化，这样复杂度就变为$O(\frac{nm}{w})$，$w=32或64$，处于可以接受的范围。

  设$bitd{p}_i[j]=d{p}_{i,j}$，这样有关$j$的操作，都可以用bitset优化。可以用右移将$j$向前移动一位，使得$bitd{p}_{i+1}[j+1]$转移到$bitd{p}_i[j]$。对于$A_i\ge B_j$的判断，可以预处理一个长度为$m$的bitset表$S_i$，$S_i[j]=(A_i\ge B_j)$，这样，转移式就优化为$$bitd{p}_i=(bitd{p}_{i+1}>>1|(1<<m))\&S_i$$
  计算$S_i$时，把$m$个$B_j$排序后，那么每个$A_i$只会插在这些$B_j$中间或两边，因此最多只有$m+1$个$S_i$（从全是$0$的情况，到全是$1$的情况）。第$j$种$S_i$，就是在前一种$S_i$的基础上在第$j$大的$B_j$位置上从$0$变为$1$。这样就可以在$O(\frac{m^2}{w})$的复杂度内预处理完$S_i$。

代码如下（不是我写的）：

#include<bits/stdc++.h>

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);
}