#LOJ2264. 「CTSC2017」吉夫特 DP计数

最新推荐文章于 2019-11-02 17:47:55 发布

lunch__

最新推荐文章于 2019-11-02 17:47:55 发布

阅读量202

点赞数

分类专栏：数学相关

本文链接：https://blog.csdn.net/lunch__/article/details/82777248

版权

数学相关专栏收录该内容

4 篇文章 0 订阅

订阅专栏

这篇博客探讨了LOJ2264竞赛题的解题思路，主要涉及一个长度为n的序列中不上升子序列的计数问题。通过卢卡斯定理转化题目的条件，发现满足条件的子序列对应于前一个数的二进制子集。作者提出使用动态规划（DP）的方法，以序列结尾的数字i来定义状态f[i]，并枚举i的子集进行转移，复杂度为O(3^log2maxa)。尽管看起来复杂度较高，但实际运行时由于常数项小，能在规定时间内完成。此外，还讨论了使用meet in middle优化技术，将复杂度降低到O(6^9)，使得解法更加高效，适用于解决类似题目。

摘要由CSDN通过智能技术生成

题意

输入一个长度为 $n$ 的序列 $a_i$ ，求有多少个长度大于等于 $2$ 的不上升子序列满足 $\displaystyle\prod_{i = 2}^{k} \binom{a_{b_{i - 1}}}{a_{b_i}} \bmod 2 > 0$ .

用卢卡斯定理转化一下性质

发现这个东西模 $2$ 为 $0$ ，当且仅当前一个数存在一个二进制位为 $0$ 而后一个数那一位是 $1$

再转化一下那后一个数就是前一个数的子集，那么不妨设 $f [i]$ 为以 $i$ 结尾的序列数

直接枚举 $i$ 的子集转移即可复杂度 $O(3^{log_2^{maxa}})$

大概是 $O(3^{18})$ 看起来这个复杂度比较虚但实际上

这个 $d p$ 常数很小而且又有两秒而且复杂度又没有卡满就过去了

Codes

#include<bits/stdc++.h>

using namespace std;

const int N = 3e5 + 10;
const int mod = 1e9 + 7;

int a[N], n, f[N];

int add(int x, int y) {
	return x + y >= mod ? x + y - mod : x + y;
}

int main() {
#ifdef ylsakioi
	freopen("3773.in", "r", stdin);
	freopen("3773.out", "w", stdout);
#endif
	int ans = 0;
	scanf("%d", &n);
	for(int i = 1; i <= n; ++ i) {
		scanf("%d", &a[i]); ++ f[a[i]];
		for(int sub = a[i]; sub; (sub -= 1) &= a[i])
			if(sub != a[i])
				f[sub] = add(f[sub], f[a[i]]);
		ans = add(ans, f[a[i]] - 1);
	}
	cout << ans << endl;
	return 0;
}

还是学一学复杂度比较正的正解做法吧

毕竟 $C C F$ 的老爷机是不会有 $L O J$ 这样的速度

一般二进制集合相关的题目都可以采用 $\ in \ middle$ 来优化

这样子复杂度直接变成了根号

像这道题我们记 $f [i] [j]$ 为前 $9$ 位是 $i$ ，后 $9$ 位为 $j$ 子集的方案数，那么每次枚举 $i$ 的超集记下来之后再转移到 $j$ 的子集上，那么这个复杂度就是 $O(6^9)$ ，就可以过了

Codes

#include<bits/stdc++.h>

using namespace std;

const int N = sqrt(3e5);
const int mod = 1e9 + 7;
const int all = (1 << 9) - 1;

int f[N][N], n;

int main() {
#ifdef ylsakioi
	freopen("4903.in", "r", stdin);
	freopen("4903.out", "w", stdout);
#endif
	int x, y, ans = 0, z;
	scanf("%d", &n);
	for(int i = 1; i <= n; ++ i) {
		int now = 0;
		scanf("%d", &z);
		x = (z >> 9), y = z & all;
		for(int j = x; j <= all; j = (j + 1) | x)
			(now += f[j][y]) %= mod;
		(ans += now) %= mod; ++ now;
		for(int j = y; ; j = (j - 1) & y) {
			(f[x][j] += now) %= mod;
			if(!j) break;
		}
	}
	cout << ans << endl;
	return 0;
}