首页 > 其他分享 >P7961 [NOIP2021] 数列

P7961 [NOIP2021] 数列

时间:2024-04-24 21:11:24浏览次数:26  
标签:P7961 数列 read ll times 进位 权值 NOIP2021 dp

P7961 [NOIP2021] 数列

这题想了一半,后面有点不敢想结果直接看题解了。

思考后发现,对于 \(a_i\le x\),也就是二进制中第 \(x\) 位前的部分,它们都可能会影响到二进制中第 \(x\) 位后的进位,而 \(a_i>x\) 的部分是不会影响到 \(x\) 位前的进位的。

所以为了满足无后效性,我们从低位向高位考虑,这样之后的答案对之前的答案没有影响。

设 \(dp(i,j,k,p)\) 为已经考虑完了前 \(i-1\) 位,已经用了 \(j\) 个 \(a_i\),前 \(i-1\) 位中有 \(k\) 个 \(1\),并且进位到第 \(i\) 为 \(p\) 个 \(1\) 的合法序列权值和。

考虑向后转移更简单,由于乘法分配律,

\((g_1+g_2+\cdots+g_l)\times v_i=g_1\times v_i+g_2\times v_i+\cdots+g_l\times v_i\)

\(g_i\) 表示当前的合法方案,所以当前权值和是可以通过先前权值和成当前位的权值得来。

考虑再枚举 \(t\) 表示第 \(i\) 位用了 \(t\) 个 \(a_i\),所以

\(dp(i+1,j+t,k+(t+p)\bmod2,\lfloor\frac{t+p}{2}\rfloor)\leftarrow dp(i,j,k,p)\times v_i^t\times \tbinom{n-j}{r}\)

组合数是因为 \(t\) 个 \(a_i\) 可以放在剩下 \(n-j\) 个位置的任意位置。

答案即为 \(\sum dp(m+1,n,k,p)\),满足 \(k+popcnt(p)\le K\),\(popcnt(p)\) 计算的是第 \(m\) 位之后的进位。

复杂度为 \(O(n^4m)\)。

总结:无后效性的转移,状态表示,不同方案之间的转移满足分配律,贡献的计算。

#include <bits/stdc++.h>
typedef long long ll;
ll read() {
	ll x = 0, f = 1;
	char c = getchar();
	while(!isdigit(c)) {
		if(c == '-') f = -1;
		c = getchar();
	}
	while(isdigit(c)) {
		x = (x << 3) + (x << 1) + (c - '0');
		c = getchar();
	} 
	return x * f;
}
ll n, m, K, mod = 998244353;
ll v[100010], c[40][40], pw[110][40], dp[110][32][32][55];
void init() {
	c[0][0] = 1;
	for(int i = 1; i <= n; i++) {
		c[i][0] = 1;
		for(int j = 1; j <= i; j++) {
			c[i][j] = (c[i - 1][j] + c[i - 1][j - 1]) % mod;
		}
	}
}
ll popcnt(ll x) {
	ll cnt = 0;
	while(x) {
		cnt += (x & 1);
		x >>= 1;
	}
	return cnt;
}
void Solve() {
	n = read(), m = read(), K = read();
	init();
	for(int i = 0; i <= m; i++) {
		v[i] = read();
		pw[i][0] = 1;
		for(int j = 1; j <= n; j++) pw[i][j] = pw[i][j - 1] * v[i] % mod;
	}
	dp[0][0][0][0] = 1;
	for(int i = 0; i <= m; i++) {
		for(int j = 0; j <= n; j++) {
			for(int k = 0; k <= K; k++) {
				for(int p = 0; p <= (n >> 1); p++) {
					for(int t = 0; t <= n - j; t++) {
						dp[i + 1][j + t][k + (t + p) % 2][(t + p) / 2] = (dp[i + 1][j + t][k + (t + p) % 2][(t + p) / 2] + dp[i][j][k][p] * pw[i][t] % mod * c[n - j][t] % mod) % mod;
						// std::cout << dp[i + 1][j + t][k + (t + p) % 2][(t + p) / 2] << "\n";
					}
				}
			}
		}
	}
	ll ans = 0;
	for(int k = 0; k <= K; k++) {
		for(int p = 0; p <= (n >> 1); p++) {
			if(k + popcnt(p) <= K) ans = (ans + dp[m + 1][n][k][p]) % mod;
		}
	}
	std::cout << ans << "\n";
}

int main() {

	Solve();

	return 0;
}

标签:P7961,数列,read,ll,times,进位,权值,NOIP2021,dp
From: https://www.cnblogs.com/FireRaku/p/18092159

相关文章

  • 6-1 使用函数求特殊a串数列和
    给定两个均不超过9的正整数a和n,要求编写函数fn(a,n)求a+aa+aaa++⋯+aa⋯aa(n个a)之和,fn须返回的是数列和函数接口定义: fn(a,n)其中a和n都是用户传入的参数。a的值在[1,9]范围;n是[1,9]区间内的个位数。函数须返回级数和裁判测试程序样例: /*请在这里填写答......
  • 「洛谷」题解:P1720 月落乌啼算钱(斐波那契数列)
    题目传送门比较经典的一道斐波那契数列的模版题,原题中给了一个很复杂的公式(也就是下面这个),但是实际上题目跟它毛关系没有……(所以放这个公式干什么)\[F_n=\dfrac{\left(\frac{1+\sqrt{5}}{2}\right)^n-\left(\frac{1-\sqrt{5}}{2}\right)^n}{\sqrt{5}}\]看见题解去了有很多人都......
  • 肚子数列即斐波那契数比列
    问题第一个月,小兔子1没有繁殖能力,所以还是1对。第二个月,小兔子1进入成熟期,仍然是1对。第三个月,兔子1生了1对小兔子2,于是这个月共有2(1+1=2)对兔子第四个月,兔子1又生了1对小兔子3.因此共有3(1+2=3)对兔子第五个月,兔子1又生了1对小兔子4,而在第三个月出生的兔子2也......
  • 力扣-665. 非递减数列
    1.题目题目地址(665.非递减数列-力扣(LeetCode))https://leetcode.cn/problems/non-decreasing-array/题目描述给你一个长度为 n 的整数数组 nums ,请你判断在最多改变 1个元素的情况下,该数组能否变成一个非递减数列。我们是这样定义一个非递减数列的: 对于数组中任......
  • 考点 1:数列极限概念
    考点1:数列极限概念考点点拨:考查数列极限的概念,即数列极限的\(\varepsilon-N\)语客谜述.【试题1-1-1】(江苏大学2006年)设\(p\)为正整数,证明:若\(p\)不是完全平方数,则\(\sqrt{p}\)是无理数.分析:考查实数的性质.证明:反证法.假设\(\sqrt{p}\)是有理......
  • 第二节 数列的极限
    第二节数列的极限数列的概念:如果按照某一法则,对每个\(n\inN\),对应着一个确定的实数\(x_n\),这些实数\(x_n\),按照下标n从小到大排列得到的一个序列\(x₁,x₂,x₃,\cdots,x_n,\cdots,\)就叫做数列,简记为数列\({x_n}\).数列中的每一个数叫做数列的项,第n项\(x_......
  • P10342 [THUSC 2019] 数列 题解
    形式化题面:求\[\sum_{l=1}^{n}\sum_{r=l}^{n}\max_{i=l}^{r}(i-l+1)\timesf(i,r)\]其中\(f(l,r)\)为\(a_l,...,a_r\)中有多少个不同的数字。注意到,除了Sub2,其余数据点都有\(\maxf\le800\),这启发我们考虑\(O(nm)\)的算法。套路地,扫描线枚举右端点,则现在只需要考虑......
  • P3986 斐波那契数列
    题目链接:P3986斐波那契数列推式子观察题目所给的序列:根据题意我们可以知道k=f(i-1)+f(i-2)那么浅浅的推一下就可以发现:k=f(3)时k=a+bk=f(4)时k=a+2bk=f(5)时k=2a+3b......故可以得出k=fib(i)a+fib(i+1)b因为a,b属于正整数故fib(n-2)+fib(n-1)>k时停止那么我们......
  • 1030 完美数列
    直观:(25分)#include<bits/stdc++.h>usingnamespacestd;#definelllonglongvector<ll>vec;map<ll,int>mp;//元素-所在下标intmain(){ intn; llp; cin>>n>>p; for(inti=0;i<n;i++){ intno; cin>>no; vec.push_bac......
  • 猴子吃桃 递归 循环 等比数列
    do-while#include<stdio.h>intmain(){intn=1;//第十天只剩下1个桃子,所以初始值为1intday=9;//第十天是已知条件,所以循环从第九天开始do{n=(n+1)*2;//每天都是前一天的一半加1,所以这里计算后一天的桃子数day--;//天数减1}while(day>=0);//......