[CTSC2017]吉夫特

最新推荐文章于 2022-12-21 08:52:00 发布

OneInDark

最新推荐文章于 2022-12-21 08:52:00 发布

阅读量148

点赞数

分类专栏： C++ 分块动态规划

本文链接：https://blog.csdn.net/qq_42101694/article/details/113817231

版权

C++ 同时被 3 个专栏收录

547 篇文章 3 订阅

订阅专栏

动态规划

220 篇文章 2 订阅

订阅专栏

分块

13 篇文章 0 订阅

订阅专栏

题目

传送门 to luogu

思路

根据卢卡斯定理，问题转化为 $a_{b_i}\subseteqq a_{b_{i-1}}$ ，即二进制下，后一个数是前一个数的子集。

从后往前 $\tt dp$ ，用 $f (x)$ 表示 $x$ 作为开头的数量。显然
$f(a_i)=\sum_{x\subseteqq a_i}f(x)$

于是我们有了一个 $\mathcal O(na)$ 的算法。这不足以通过此题（理论上）。

利用一个很简单的技巧，可以做到 $\mathcal O(n\sqrt{a})$ 。每个数存储后 $\log a\over 2$ 位可以选择子集的 $f$ 求和。转移时，枚举前 $\log a\over 2$ 位，再查这张表即可。更新这张表时，枚举后 $\log a\over 2$ 位即可。

代码

#include <cstdio>
#include <iostream>
#include <cstring>
using namespace std;
typedef long long int_;
inline int readint(){
	int a = 0; char c = getchar(), f = 1;
	for(; c<'0'||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}
inline void writeint(int_ x){
	if(x > 9) writeint(x/10);
	putchar((x-x/10*10)^48);
}

const int MaxN = 300000;
const int Mod = 1000000007;

int f[1<<9][1<<9], dp[MaxN];
const int all = ((1<<9)-1);
void modify(int x,int v){
	const int S = x&all;
	f[x>>9][S] = (f[x>>9][S]+v)%Mod;
	for(int i=all^S; i; i=(i-1)&(all^S))
		f[x>>9][S^i] = (v+f[x>>9][S^i])%Mod;
}

int a[MaxN];
int main(){
	int n = readint();
	for(int i=1; i<=n; ++i)
		a[i] = readint();
	modify(0,1); // initial
	int ans = 0; dp[0] = 1;
	for(int i=n; i>=1; --i){
		int x = a[i]>>9;
		dp[a[i]] = f[0][a[i]&all];
		for(int j=x; j; j=(j-1)&x)
			dp[a[i]] = (dp[a[i]]
				+f[j][a[i]&all])%Mod;
		modify(a[i],dp[a[i]]);
		ans = (ans+dp[a[i]]-1)%Mod;
	}
	printf("%d\n",(ans+Mod)%Mod);
	return 0;
}