[CSP-S2021] 括号序列
链接:
题意:
有一堆规则,然后判断给定字符串有多少种填法符合规则。
分析:
一眼区间dp,状态数 \(n^2\),我们来分析这些规则。
把这些规则分成三类,第一类可以预处理出区间是否能表达成全部都是 *
的情况并且其长度小于等于 \(k\),后 \(O(1)\) 判断。第二类可以考虑枚举 \(S\) 的长度,\(O(n)\) 处理。第三类如果枚举分割点显然会算重,依套路我们枚举第一个括号的位置且要保证第一个括号无法被继续分割,考虑在每个区间维护一个无法被分割的方案数,也就是不计第三类的方案,记做 \(dp[l][r][0]\),总方案数记做 \(dp[l][r][1]\),然后再枚举点的长度,\(O(n^2)\) 处理,时间复杂度 \(O(n^4)\)。
优化:
对于第三类的处理,假如我们固定第一个括号的位置,那么需要被计算的其实是一个后缀和。这里我固定除去第一个括号的后面括号的位置,那么就是一个前缀和,两者没有什么区别。大概长这样:
于是我们考虑做一个前缀和,就是 \(sum[l][r]=\sum\limits_{k=l}^r dp[l][i][0]\)
那么第三类的转移只需要处理出前缀和左边的边界就可以 \(O(n)\) 转移了,这个边界会受到三个限制,一是至少要大于等于 \(l\),其次与 \(r\) 的距离会受到 \(k\) 的限制,再是如果前面有一个位置不能表示为 *
,显然从这之前的都不能算,所以只需要预处理出一个位置之前最近的不能表示为 *
的位置就好了。然后还有一个细节就是这样做会多减一个没有 *
的情况,加上就好了。
时间复杂度 \(O(n^3)\)
算法:
先预处理出区间是否能表达成全部都是 *
的情况并且其长度小于等于 \(k\),以及一个位置之前最近的不能表示为 *
的位置,然后根据规则 dp 同时维护前缀和就好了。
代码:
#include<bits/stdc++.h>
using namespace std;
#define int long long
#define in read()
inline int read(){
int p=0,f=1;
char c=getchar();
while(!isdigit(c)){if(c=='-')f=-1;c=getchar();}
while(isdigit(c)){p=p*10+c-'0';c=getchar();}
return p*f;
}
const int N=505;
const int mod=1e9+7;
int n,k,sum[N],dp[N][N][2],h[N][N],ex[N][N],last[N];
char s[N];
inline int add(int x,int y){int z=x+y;return (z>=mod)?(z-mod):z;}
inline int mul(int x,int y){int z=x*y;return (z>=mod)?(z%mod):z;}
inline void Add(int &x,int y){x=add(x,y);}
inline void Mul(int &x,int y){x=mul(x,y);}
signed main(){
freopen("bracket.in","r",stdin);
freopen("bracket.out","w",stdout);
cin>>n>>k>>(s+1);
for(int i=1;i<=n;i++)
if(s[i]=='*'||s[i]=='?')
sum[i]=sum[i-1]+1,last[i]=last[i-1];
else
last[i]=i,sum[i]=sum[i-1];
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
if((sum[j]-sum[i-1]==j-i+1&&j-i+1<=k)||j<i)
h[i][j]=1;
else h[i][j]=0;
for(int i=1;i<n;i++)
if((s[i]=='('||s[i]=='?')&&(s[i+1]==')'||s[i+1]=='?'))
dp[i][i+1][0]=dp[i][i+1][1]=ex[i][i+1]=1;
for(int i=3;i<=n;i++)
for(int l=1,r=l+i-1,temp=min(r-l-3,k);r<=n;l++,r=l+i-1){
if((s[l]=='('||s[l]=='?')&&(s[r]==')'||s[r]=='?')){
if(h[l+1][r-1])dp[l][r][0]++;
if(l+1<r-1)Add(dp[l][r][0],dp[l+1][r-1][1]);//1
if(s[l+1]=='('||s[l+1]=='?')
for(int t=1;t<=temp;t++)
Add(dp[l][r][0],dp[l+1][r-1-t][1]*h[r-t][r-1]);
if(s[r-1]==')'||s[r-1]=='?')
for(int t=1;t<=temp;t++)
Add(dp[l][r][0],dp[l+1+t][r-1][1]*h[l+1][l+t]);//2
dp[l][r][1]=dp[l][r][0];
for(int t=r-1;t>=l+2;t--){
int ll=max(t-1-k,l);
ll=max(ll,last[t-1]);
Add(dp[l][r][1],mul(dp[t][r][1],add(add(ex[l][t-1]-ex[l][ll],dp[l][ll][0]),mod)));
}//3
}
ex[l][r]=add(ex[l][r-1],dp[l][r][0]);//前缀和
}
cout<<dp[1][n][1];
return 0;
}
题外话:
区间dp做少了。
[CSP-S2021] 括号序列的更多相关文章
- 上午小测3 T1 括号序列 && luogu P5658 [CSP/S 2019 D1T2] 括号树 题解
前 言: 一直很想写这道括号树..毕竟是在去年折磨了我4个小时的题.... 上午小测3 T1 括号序列 前言: 原来这题是个dp啊...这几天出了好几道dp,我都没看出来,我竟然折磨菜. 考试的时候先 ...
- BZOJ4350: 括号序列再战猪猪侠
Description 括号序列与猪猪侠又大战了起来. 众所周知,括号序列是一个只有(和)组成的序列,我们称一个括号 序列S合法,当且仅当: 1.( )是一个合法的括号序列. 2.若A是合法的括号序列 ...
- DP专题——括号序列
毕竟是个渣,写完一遍之后又按LRJ的写了一遍,再写了一遍递归版,最终加上输出解部分 括号序列 定义如下规则序列(字符串): 空序列是规则序列: 如果S是规则序列,那么(S)和[S]也是规则序列: 如果 ...
- 【BZOJ】2209: [Jsoi2011]括号序列(splay)
http://www.lydsy.com/JudgeOnline/problem.php?id=2209 splay又犯逗........upd1那里的sum忘记赋值反............. 本题 ...
- 51nod1476 括号序列的最小代价
这题应该可以用费用流写吧?不过我想不出贪心来TAT.其实还是单调队列乱搞啊T_T //ÍøÉϵÄÌ°ÐÄËã·¨ºÃÉñ°¡¡£¡£¡£ÎÒÖ»»áÓÃ×îС·ÑÓÃ×î´óÁ÷ÅÜTAT #in ...
- lintcode: 有效的括号序列
题目: 有效的括号序列 给定一个字符串所表示的括号序列,包含以下字符: '(', ')', '{', '}', '[' and']', 判定是否是有效的括号序列. 样例 括号必须依照 "() ...
- uoj #31. 【UR #2】猪猪侠再战括号序列 贪心
#31. [UR #2]猪猪侠再战括号序列 Time Limit: 20 Sec Memory Limit: 256 MB 题目连接 http://uoj.ac/problem/31 Descript ...
- bzoj 1095 [ZJOI2007]Hide 捉迷藏(括号序列+线段树)
[题目链接] http://www.lydsy.com/JudgeOnline/problem.php?id=1095 [题意] 给定一棵树,树上颜色或白或黑而且可以更改,多个询问求最远黑点之间的距离 ...
- CODEVS 3657 括号序列
[问题描述] 我们用以下规则定义一个合法的括号序列: (1)空序列是合法的 (2)假如S是一个合法的序列,则 (S) 和[S]都是合法的 (3)假如A 和 B 都是合法的,那么AB和BA也是合法的 例 ...
- bzoj1095: [ZJOI2007]Hide 捉迷藏 线段树维护括号序列 点分治 链分治
这题真是十分难写啊 不管是点分治还是括号序列都有一堆细节.. 点分治:时空复杂度$O(n\log^2n)$,常数巨大 主要就是3个堆的初始状态 C堆:每个节点一个,为子树中的点到它父亲的距离的堆. B ...
随机推荐
- Django学习day08随堂笔记
今日考题 """ 今日考题 1.聚合查询,分组查询的关键字各是什么,各有什么特点或者注意事项 2.F与Q查询的功能,他们的导入语句是什么,针对Q有没有其他用法 3.列举常 ...
- Django学习day04随堂笔记
每日测验 """ 今日考题 1.列举你知道的orm数据的增删改查方法 2.表关系如何判定,django orm中如何建立表关系,有什么特点和注意事项 3.请画出完整的dj ...
- CentOS7下编译安装redis-5.0.9
CentOS7下编译安装redis-5.0.9 本文地址http://yangjianyong.cn/?p=171转载无需经过作者本人授权 下载redis #code start wget https ...
- 在PHP中灵活使用foreach+list处理多维数组
先抛出问题,有时候我们接收到的参数是多维数组,我们需要将他们转成普通的数组,比如: $arr = [ [1, 2, [3, 4]], [5, 6, [7, 8]], ]; 我们需要的结果是元素1变成1 ...
- Orchard Core Framework Samples
解决方案包含内容 多租户应用 一个ASP.NET Core MVC应用程序,它引用模块项目,并为两个启用了不同模块的租户提供支持. 此Web应用程序的主页提供了更多信息,并链接到两个租户和模块端点.租 ...
- Groovy系列(4)- Groovy集合操作
Groovy集合操作 Lists List 字面值 您可以按如下所示创建列表. 请注意,[]是空列表表达式 def list = [5, 6, 7, 8] assert list.get(2) == ...
- 查看Docker启动jenkins的管理员密码
Docker启动docker后,第一次方法jenkins,需要输入管理员密码. 其实查看启动时候的日志可以看到密码,也可以按照以下方法找到密码. 1.查看docker容器ID:docker ps -a ...
- 低差异序列 (low-discrepancy sequences)之Hammerysley在半球中采样点方法的介绍
半球上的Hammersley 源作者:Holger Dammertz 一组关于如何在2D中使用Hammersley点集以在着色器程序中快速实用地生成半球方向的笔记.如果你发现任何错误或有意见,不要犹豫 ...
- Jmeter压测学习1—第一个项目:登录
我用的是我们公司的测试环境练习的 网址:http://****:9000/login.html 账号是admin 密码:123456 一.打开Jmeter 找到安装目录:bin->Jmeter ...
- React-高阶函数_函数柯里化
高阶函数_函数柯里化 高阶函数(定义) 如果一个函数符合下面两个规范,就是高阶函数: 如果A函数,接收的参数是一个函数,那么A就是一个高阶函数(比如数组方法arr.map()接收的就是一个处理item ...