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

Luogu P7914 [CSP-S 2021] 括号序列 题解 [ 蓝 ] [ 区间 DP ] [ 前缀和优化 ] [ 调试技巧 ]

括号序列:无聊,感觉做过类似的拼接类区间 DP 就直接秒了。

注意到这个超级括号序列定义很复杂,除了两边没有 \(\texttt{*}\) 没有啥很好的性质。于是直接考虑暴力区间 DP:定义 \(dp_{l, r}\) 表示 \(l\sim r\) 的合并方案数。

因为是拼接类区间 DP,所以还需要定义一个辅助转移数组 \(g_{l, r}\),表示将 \(l\sim r\) 合成到一个括号内的方案数。

剩下的就是些无聊的转移了,我们根据题面上的定义分别计数:

  • 规则 \(1\)
    • 直接枚举区间预处理 \(dp, g\) 的初值即可。
  • 规则 \(3\):需要特判区间两边是否能被设为括号。
    • (A):从 \([l + 1, r - 1]\) 处转移 \(dp, g\) 即可。
    • (SA):枚举左侧的 S 转移 \(dp, g\)
    • (AS):同理。
  • 规则 \(2\)
    • AB:拼接类区间 DP 板子,将转移划分为拼接左侧的整段,利用 \(g\) 数组转移即可。
    • ASB:同理,只是多了个后缀和优化而已。

直接 DP,时间复杂度 \(O(n^3)\)

一个计数题调试小技巧 from Lucyna_Kushinada,让我们膜拜他!!!1:

在 DP 转移的时候记录每个 DP 值存的方案,重载运算符并利用 vector 转移。可以比较直观地调试重复和漏掉的计数。

#include <bits/stdc++.h>
#define fi first
#define se second
using namespace std;
typedef long long ll;
using pi = pair<int, int>;
const int N = 505;
const ll mod = 1e9 + 7;
ll n, m, dp[N][N], suf[N], mxlen[N], g[N][N];
char s[N];
void add(ll &x, ll v)
{x += v;if(x >= mod) x -= mod;
}
int main()
{// freopen("bracket.in", "r", stdin);// freopen("bracket.out", "w", stdout);ios::sync_with_stdio(0);cin.tie(0);cout.tie(0);cin >> n >> m >> s + 1;// 预处理每个 DP 的初始值、每个位置的最长延伸值,即第一条规则for(int i = 1; i <= n; i++){for(int j = i + 1; j <= n; j++){bool legal = 1;if(!(s[i] == '(' || s[i] == '?')) legal = 0;if(!(s[j] == ')' || s[j] == '?')) legal = 0;for(int k = i + 1; k < j; k++)if(!(s[k] == '*' || s[k] == '?'))legal = 0;if(j - i + 1 - 2 > m) legal = 0;dp[i][j] = legal;g[i][j] = legal;}mxlen[i] = i - 1;for(int j = i; j <= n; j++){if(!(s[j] == '*' || s[j] == '?')) break;if(j - i + 1 > m) break;mxlen[i] = j;}}// 进行转移for(int len = 3; len <= n; len++){for(int l = 1; l + len - 1 <= n; l++){int r = l + len - 1;if((s[l] == '(' || s[l] == '?') && (s[r] == ')' || s[r] == '?')){// (A)add(dp[l][r], dp[l + 1][r - 1]);add(g[l][r], dp[l + 1][r - 1]);// (SA)bool legal = 1;for(int i = l + 1; i <= r - 2; i++){legal &= (s[i] == '*' || s[i] == '?');legal &= (i - l <= m);if(!legal) break;add(dp[l][r], dp[i + 1][r - 1]);add(g[l][r], dp[i + 1][r - 1]);}// (AS)legal = 1;for(int i = r - 1; i >= l + 2; i--){legal &= (s[i] == '*' || s[i] == '?');legal &= (r - i <= m);if(!legal) break;add(dp[l][r], dp[l + 1][i - 1]);add(g[l][r], dp[l + 1][i - 1]);}}// AB,要求 A 为整段for(int i = l + 1; i <= r - 1; i++)add(dp[l][r], g[l][i] * dp[i + 1][r] % mod);// ASB,要求 A 为整段memset(suf, 0, sizeof(suf));for(int i = r; i >= l; i--){suf[i] = suf[i + 1];add(suf[i], dp[i][r]);}for(int i = l + 1; i <= r - 2; i++){int rx = mxlen[i + 1];if(rx == i) continue;rx = min(rx, r - 1);ll tmp = g[l][i] * (suf[i + 2] - suf[rx + 2]) % mod;tmp = ((tmp % mod) + mod) % mod;add(dp[l][r], tmp);}}}cout << dp[1][n];return 0;
}
http://www.jsqmd.com/news/24923/

相关文章:

  • 扩展BaseMapper类 - 详解
  • 《程序员修炼之道:从小工到专家》前五分之二观后感
  • 矩阵快速幂章节笔记(这里主要介绍的是我的错题)
  • 实验二 现代C++编程初体验
  • P5322 [BJOI2019] 排兵布阵
  • 题解:P9292 [ROI 2018] Robomarathon
  • [题解]P5322 [BJOI2019] 排兵布阵
  • 考前打印
  • 申威服务器安装Nacos 2.0.3 RPM包详细步骤(Kylin V10 sw_64架构)​附安装包
  • ZKY精选冲刺省选国赛仿真训练题
  • MySQL 查询与更新语句执行过程深度解析:从原理到实践​ - 指南
  • ZKY精选冲刺省选国赛技巧训练题
  • 逆向基础--编码(001)
  • 20251027 - 倍增 ST表
  • 周康阳精选冲刺省选国赛思维训练题
  • Luogu P7913 [CSP-S 2021] 廊桥分配 题解 [ 绿 ] [ 贪心 ] [ 前缀和 ] [ STL ]
  • 10-27 CSP 赛前比赛记录
  • P3939 数颜色
  • 完整教程:Docker 搭建 Nginx 并启用 HTTPS 具体部署流程
  • AI开发微信小程序-有感
  • 价值流智能时代:DevOps平台如何成为企业高效交付的核心引擎? - 教程
  • 2025年压力容器品牌综合实力排行榜
  • 2025年压力容器厂家综合评测与选择指南
  • 2025年口碑好的压力容器工厂/厂家前十强
  • 科幻——面包
  • 2025年中国钢结构码垛机制造商Top 5排名解析
  • 2025年钢结构码垛机品牌前十强权威盘点:江苏众利达引领智能制造新浪潮
  • 处理django.db.utils.OperationalError: attempt to write a readonly database错误
  • 10.28代码大全2
  • [GESP202509 二级] 菱形