子序列（混合递推+前缀和 || 分块）

问题

• 序列 $a_1，a_2，\cdots ，a_n$，求子序列中长度大于1的传奇序列的个数。$1\le n\le 3\ast 10^5,0\le a_i<2^{18}$，结果对100000007取余。
• 若序列 $b_1，b_2，\cdots ，b_k$ 满足 $b_i\&b_j=b_j且i<j$ ，则称为传奇序列。

分析

• 遍历子序列的一般方式：将子序列按结尾元素分成子空间，子空间数量 $O(n)$；
• $\&$ 运算性质
  • $a\&b=b$，例如：$a=11101010，b=11000010$
  • 若 $a\&b=b，b\&c=c$，那么 $a\&c=c$
  • 小于 $a$ 的最大 $b=(a-1)\&a$，递推式 $b_{i+1}=(b_i-1)\&a(b_i>b_{i+1})$
  • 大于 $b$ 的最小 $a=(b+1)\Vert a$，递推式 $a_{i+1}=(a_i+1)\Vert a(a_i<a_{i+1})$

暴力枚举

分析

• 设 $f[i]$ 是以元素 $a_i$ 结尾的的满足条件的子序列的数量（元素下标：$i$ ）
• 设 $g[a]$ 是以元素 $a$ 结尾的的满足条件的子序列的数量（元素值：$a$ ）
• 则递推过程
  • 到哪里去？
    • $g[a]$ 定义：能与值为 $a$ 的元素构成传奇序列的子序列数量
    • $f[i]=g[a_i]+1$
    • $g[a_j]=g[a_j]+f[i]$，其中 j ∈ { i ， i + k 1 ， i + k 2 ， ⋯   } ， a i    &    a j = a j j \in \{i，i+k_1，i+k_2，\cdots \}，a_i \;\&\;a_j=a_j j∈{i，i+k1​，i+k2​，⋯}，ai​&aj​=aj​，复杂度$O(n)$
    • 可优化：序列的后续元素中可能不出现值为 $a_j$ 的元素，能延迟或不操作更佳；依据查询的“组合”策略，进行组合内修改
  • 从哪里来?
    • $g[a]$ 定义：尾元素值为 $a$ 的传奇数列数量
    • $f[i]={\sum }_{k}g[a_k]$，其中 $1\le k\le i-1，a_k\&a_i=a_i$，复杂度$O(n)$
    • $g[a_i]=g[a_i]+f[i]$
    • 可优化：某些数据可能被反复查询而且是可以成为一个固定的组合（前缀和就是一种组合）
• 无论是哪种递推方式，复杂度都为：$O(n^2)$

混合递推+前缀和

分析

• 对暴力枚举进行优化
• 递推方法：" 到哪里去？“和"从哪里来?” 相结合
• $f[i]$ 定义不变，即以元素 $a_i$ 结尾的的满足条件的子序列的数量（元素下标：$i$ ）
• $g[a]\Rightarrow g[A][B]$ ，$A、B$ 分别等于元素值 $a$ 的高9位和低9位
• 例如：153629的高9位是 $300=153629>>9$，低9位是 $29=153629\&511$
• $g[A][B]$ 的含义：
  • 已有的序列中尾元素的高9位是A且能与值为 $a$ 的元素构成传奇序列的子序列数量（是一个组合）
  • $g[A][B]=\sum f[i]$，其中 $a[i]=(A<<9)+B_k，B_k\&B=B$
  • 类似于是一个前缀和或后缀和， 以153629为例，则 B k ∈ { 29 ， 31 ， ⋯ ， 509 ， 511 } B_k \in \{29，31，\cdots，509，511\} Bk​∈{29，31，⋯，509，511}

   29    31    61    63    93    95   125   127   157   159   189   191   221   223   253   255
  285   287   317   319   349   351   381   383   413   415   445   447   477   479   509   511

• 递推过程：
  • 设 $A=a_i>>9，B=a_i\&511$
  • $f[i]={\sum }_{j=1}^{k}g[A_j][B]+1$，其中 $A=A_1<\cdots <A_k\le 511，A_j\&A=A$，复杂度 $O(\sqrt{n})$，【从哪里来？】
  • $g[A][B_j]=g[A][B_j]+f[i]$，其中 $B=B_1>\cdots >B_k\ge 0，B\&B_j=B_j$，复杂度 $O(\sqrt{n})$，【到哪里去？】
• 复杂度： $O(n\sqrt{n})$

代码

#include<bits/stdc++.h>
using namespace std;
#define MD 1000000007
int n, a, g[520][520], ans;
int main(){
	int t, x, y, cur;
	scanf("%d", &t);
	while(t--){
		scanf("%d", &n);
		memset(g, 0, sizeof g);
        ans = 0;
		for(int i = 1; i <= n; ++i){
			scanf("%d", &a);
			x = a>>9, y = a&511, cur = 1;
			for(int j = x; j <= 511; j = (j+1)|x)
                cur = (cur + g[j][y]) % MD;  // cur: f[i]
			ans = (ans + cur)%MD;
			for(int j = y; ; j = (j-1)&y){
				g[x][j] = (g[x][j] + cur)%MD;
                if(j == 0) break;
            }
		}
		printf("%d\n", ans-n>0?(ans-n):(ans-n+MD));
	}	
    return 0;
}

分块

分析

• 按尾元素下标可以递推

代码