当前位置: 首页 > news >正文

动态规划之数列

题目描述

有这样一种数列 a1​,a2​,⋯an​,其中 a1​=0,且对任意一项 ai​ 满足 ∣ai​−ai+1​∣=1(1≤i<n)。设 s=∑i=1n​ai​=a1​+a2​+⋯+an​,表示前 n 项之和。

现在给出数列长度 n 与数列前 n 项之和 s,要求:

输出满足条件的数列的总数对 264 取模后的结果。

输出满足条件的 100 个数列(如果不满 100 个就全部输出)。

输入格式

输入只有一行两个整数,分别表示数列长度 n 与数列之和 s。

输出格式

本题存在 Special Judge

在第一行输出一行一个整数 t,表示满足条件的数列总数对 264 取模后的结果。

接下来若干行,每行 n 个整数,表示一个数列。

输入输出样例

输入 #1复制

4 0

输出 #1复制

2
0 -1 0 1
0 1 0 -1

题解:

假如每一个数字都比前面的一个数字大 1,即数列为 0,1,2,3,4,5,…,n−1,那么这个数列的和为 sum=2n(n−1)​,我们发现 n≤100,那么 sum≤4950,所以如果题目要求的 s>sum,那么一定无解,因为这个数列再大也大不过 sum,更不可能到达 s 了。

反过来,如果一个数字比前面一个数字小 1,即数列为 0,−1,−2,−3,−4,−5,…,−(n−1),那么这个数列的和为 sum=−2n(n−1)​,我们发现 n≤100,那么 sum≥−4950。

所以 −263≤s≤263 是吓你的,真正有用的 −4950≤s≤4950。


接下来,我们来探究一个位子 u 上的数字变化对数组的和 sum 的影响。

如果 au​=au−1​+1,现在改成 au​=au−1​−1,那么从 u 开始的每一个数字都会减去 2。

那么数组的和 sum 就会减去 2(n−u+1)。


那么,我们可以想到让每一个数字都等于前面一个数字 +1,那么和就是 2n(n−1)​。

但是我们想让其变为题目要求的 s,那么要减去的数字就是 k=2n(n−1)​−s。

那么我们只能将有些数字间的 +1,换成 −1 才能达到减去 k 的目的。

那么根据前面对数字间关系的讨论,我们知道,我们要凑出很多个 i,使这些 i 对应的 2(n−i+1) 加起来等于 k 就行了,即让所有的 n−i+1 加起来等于 2k​ 即可。

当然了,如果 k 是奇数,那么一定无解。


我们考虑 DP,设 fi,j​ 表示已经凑到第 i 个数字,和为 j 的方案数。

开始时:f1,0​=1,

目标:fn,2k​​。

我们可以不选择第 i 个数字,fi,j​=fi−1,j​。

我们可以选择第 i 个数字,即在原有基础上加上 (n−i+1),fi,j​=fi−1,j−(n−i+1)​。

代码:

#include <bits/stdc++.h>using namespace std;
using i64 = long long;
using ull = unsigned long long;const int N = 1010, M = 5010;ull f[N][M];int n, k;
i64 s;
int cnt;
int m[N];void dfs(int u, int sum) {if (sum > (k >> 1)) return;if (u > n) {if (sum == (k >> 1)) {cnt++;i64 tmp = 0;for (int i = 1; i <= n; i++) {tmp += m[i];cout << tmp << ' ';}cout << '\n';}if (cnt >= 100) {exit(0);}return;}m[u] = -1;dfs(u + 1, sum + (n - u + 1));m[u] = 1;dfs(u + 1, sum);
}int main() {ios::sync_with_stdio(false);cin.tie(nullptr);cin >> n >> s;if (s > n * (n - 1) / 2) {cout << 0 << '\n';return 0;}k = n * (n - 1) / 2 - s;if (k & 1) {cout << 0 << '\n';return 0;}f[1][0] = 1;for (i64 i = 2; i <= n; i++) {i64 x = (n - i + 1);memcpy(f[i], f[i - 1], sizeof(f[i]));for (int j = x; j < M; j++) {f[i][j] = f[i][j] + f[i - 1][j - x];}}cout << f[n][k >> 1] << '\n';dfs(2, 0);return 0;
}

 

 

相关文章:

  • 免疫浸润分析
  • C语言-8.数组
  • Java大师成长计划之第23天:Spring生态与微服务架构之服务发现与注册中心
  • Manus AI 原理深度解析第三篇:Tools
  • 电总协议调试助手更新-PowerBus-v1.0.5
  • 作业帮Java后台开发面试题及参考答案(下)
  • ACI Fabric 中的各种地址
  • OneNote内容太多插入标记卡死的解决办法
  • 汽配知识(三)|跨境电商平台的汽配类目划分与关键词逻辑
  • Hive PredicatePushDown 谓词下推规则的计算逻辑
  • 嵌入式学习笔记DAY21(双向链表、Makefile)
  • 盲盒:拆开未知的惊喜,收藏生活的仪式感
  • 养生:解锁健康生活的核心密码
  • js在浏览器执行原理
  • golang -- 认识channel底层结构
  • AI软件汇总与功能解析:赋能未来的智能工具库
  • 以项目的方式学QT开发(二)——超详细讲解(120000多字详细讲解,涵盖qt大量知识)逐步更新!
  • mysql 基础复习-安装部署、增删改查 、视图、触发器、存储过程、索引、备份恢复迁移、分库分表
  • 8、SpringBoot集成MinIO
  • 鸿蒙OSUniApp 制作简洁高效的标签云组件#三方框架 #Uniapp
  • 艺术稀缺性和价值坚守如何构筑品牌差异化壁垒?从“心邸”看CINDY CHAO的破局之道
  • 证监会强化上市公司募资监管七要点:超募资金不得补流、还贷
  • MSCI中国指数5月调整:新增5只A股、1只港股
  • 崔登荣任国家游泳队总教练
  • 中国海警舰艇编队5月14日在我钓鱼岛领海巡航
  • 巴方:印度上周导弹袭击造成至少40名平民死亡