WEEK12 周记 作业——区间DP_选做题1

本文详细探讨了括号序列问题,介绍了如何通过动态规划求解最长合法子序列的方法,包括状态定义、状态转移方程及代码实现,适用于算法竞赛及编程实践。

一、题意

1.简述

括号序列问题:
• 定义一个合法的括号序列
• 1. 空序列合法
• 2. S合法,那么(S)和[S]合法
• 3. 假设A和B合法,那么AB合法
• 合法序列 - (), [], (()), ([]), ()[], ()[()]
• 非法序列 - (, [, ], )(, (]), ([(),([)]
• 给定一个序列,求其合法子序列的最长长度。
• Simple Input: ([()[
• Simple Output: 2

2.输入格式

输入测试文件将包含多个测试用例。每个输入测试用例由一行组成,其中只包含字符(、)、[和];每个输入测试的长度都在1到100之间。文件结束由包含单词“end”的行标记,不应该被处理。

3.输出格式

对于每个输入情况,程序应该将最长的可能的正则括号子序列的长度打印在一行上。

4.样例

Input

((()))
()()()
([]])
)[)(
([][][)
end

Output

6
6
4
0
6

二、算法

主要思路

首先我们可以将这个问题进行转化:不难理解,求最长的符合要求的子序列,等价于求最少添加多少个括号使整个序列满足要求。
之后确定DP的状态:我们设f[i][j]f[i][j]f[i][j]为区间[i,j][i,j][i,j]的序列,最少需要添加的括号数。
这样,最终我们要求的就是f[1][n]f[1][n]f[1][n]
下面给出DP的状态转移方程:
首先,先考虑边界条件,长度为1的串的答案一定是1,长度为0的串的答案一定是0。
对于某个子序列S=[si,sj]S = [si,sj]S=[si,sj]
f[i][j]=min{f[i][k]+f[k+1][j]}     i≤k<jf[i][j]=min\{f[i][k] + f[k+1][ j]\}\ \ \ \ \ i\le k\lt jf[i][j]=min{f[i][k]+f[k+1][j]}     ik<j 对于某个子序列S=[si,sj]S = [si,sj]S=[si,sj],且S形如[S’][S’][S](S’)(S’)(S) ,那么f[i][j]=min{f[i+1][j−1],min{f[i][k]+f[k+1][j]}}     i≤k<jf[i][ j] = min\{f[i+1][ j-1],min\{f[i][k] + f[k+1][ j]\}\}\ \ \ \ \ i\le k\lt jf[i][j]=min{f[i+1][j1],min{f[i][k]+f[k+1][j]}}     ik<j上述方程可能比较难以理解,下面做一下解释。
假如存在一个最少的添加括号的方案,那么只要紧邻着不符合要求的括号加就行,那也必然可以将现有的序列从某一个位置处进行划分,使得将两边的不符合要求的括号分开。但是需要注意的是,这种划分可能会产生新的不符合要求的括号,这是因为划分处在一对符合要求的括号中间。
但是,只要这一对括号不处在两头,那么必然存在一种划分,使得左右两个子序列并不会产生新的不符合要求的括号,从而使得左右两个子序列的fff值加起来就是父序列的fff值。
如果有一对符合要求的括号处在序列的两头,这个时候我们就得单独考虑,也就是第2个状态方程。把直接将两头的这对括号去掉作为一种方式加进来考虑。
实际上,如果两头的括号真的是一对,那么直接去掉就可以,但是,是否是一对只能是在获知了最长的符合要求的子序列之后才能知道,而目前我们就是在求这个最长子序列。所以我们只能把去掉作为一种考虑的情况。
()()()()()()这种,两头的括号在最长的符合要求的子序列中就不是一对。这个时候就应该是用划分方法才能求出,所以在两头的括号可能是一对的情况下,既要考虑直接去掉,也要考虑划分的情况。


三、代码

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
using namespace std;
int f[110][110];  //[i,j]中最少的不合格的括号数 
int main(){
	while(1){
		memset(f,63,sizeof(f));
		string str;
		char s[110];
		cin>>str;
		if(str == "end") break;
		int lth = str.length();
		strcpy(s,str.c_str());
		
		for(int len=1;len<=lth;len++){
			for(int i=0;i<=lth-len;i++){
				int j = i + len-1;
				if (len == 1){
					f[i][j] = 1;
					continue;
				}
				
				//len > 1
				if((s[i]=='('&&s[j]==')')||(s[i]=='['&&s[j]==']'))
					f[i][j] = i+1<j?f[i+1][j-1]:0;
					
				for(int k=i;k<j;k++){
					int temp = f[i][k]+f[k+1][j];
					if (temp<f[i][j]) 
						f[i][j] = temp;
				}
			}
		}
		printf("%d\n",lth - f[0][lth-1]);
	}
	return 0;
} 
/*
((()))
()()()
([]])
)[)(
([][][)
end
*/
你的身份是高级编程技术专家,精通各类编程语言,能对编程过程中的各类问题进行分析和解答。我的问题是【我正在编辑【通达信量化择时股】代码,遇到了 【 {日线突破阈值} {——行业轮动因子V2——} IND_RPS := EMA(C / MAX(REF(C, 60), 0.001), 13) * 0.7 + RANK(VOL / MAX(MA(VOL, 60), 0.001)) * 0.3 详细信息 : 您在括号前写的不是函数、公式等, 且缺少必要的运算符! 错误起始位置 : 304 ; 长度: 3】,请帮我检查并改正错误点补全正确代码,生成修正后完整代码。原有股逻辑完整保留,所有参数计算关系和信号触发条件均不改变。我的原始代码如下:【{—————————————— 战略参数模块 ——————————————} {——周期协同参数——} MONTH_MA := MA(CLOSE, 20); {20月价值中枢} WEEK_VOL := MA(VOL, 5); {周量能基准} DAY_BREAK := HHV(HIGH, 10); {日线突破阈值} {——行业轮动因子V2——} IND_RPS := EMA(C / MAX(REF(C, 60), 0.001), 13) * 0.7 + RANK(VOL / MAX(MA(VOL, 60), 0.001)) * 0.3; {修正RANK参数写法,移除时间窗口参数} IND_MOM := EMA(SLOPE(C, 21) * STD(C, 21), 8); INDUSTRY_WEIGHT := IF(IND_RPS > 75 AND MONTH_MA > REF(MONTH_MA, 3), 1.25, 1); {——估值体系V2——} DYNPETTM := IF(FINANCE(33) > 0, C / (FINANCE(33) / FINANCE(1) * INDUSTRY_WEIGHT), 1000); PB_RATE := C / (FINANCE(5) * IIF(INDUSTRY_WEIGHT > 1, 1.03, 1)); PEG_VAL := DYNPETTM / MAX((FINANCE(54) / FINANCE(34)) * 100 * INDUSTRY_WEIGHT, 0.001); {——波动率自适应V2——} IND_VOL := STD(INDEXC, 60) / MA(INDEXC, 60); VOLATILITY := STD(C, 60) / MA(C, 60) * 0.6 + IND_VOL * 0.4; VAR_PERIOD := IF(VOLATILITY < 0.05, 89, IF(VOLATILITY < 0.1, 55, 34)); {—————————————— 战术指标模块 ——————————————} {——三维MACD系统——} FAST_LEN := IF(VOLATILITY < 0.08, 12, 8); SLOW_LEN := IF(VOLATILITY < 0.08, 26, 17); SGNL_LEN := IF(VOLATILITY < 0.08, 9, 6); DIF := EMA(C, FAST_LEN) * INDUSTRY_WEIGHT - EMA(C, SLOW_LEN); DEA := EMA(DIF, SGNL_LEN); MACD := 2 * (DIF - DEA); MACD_ANGLE := ATAN((DIF - REF(DIF, 1)) / MAX(VOLATILITY, 0.03)) * 57.3; {——量价共振系统V2——} VOL_SIG := VOL > MA(VOL, 60) * 2.2 AND ABS(H - L) / REF(C, 1) > 0.08 AND C > MONTH_MA AND WEEK_VOL > REF(WEEK_VOL, 1) * 1.3; {——多周期验证——} MONTH_SIG := C > MONTH_MA AND MONTH_MA > REF(MONTH_MA, 1); WEEK_SIG := MACD > REF(MACD, 1) AND VOL > WEEK_VOL * 1.5; DAY_SIG := CROSS(C, DAY_BREAK) AND MACD_ANGLE > 30; {——盘中突击系统——} M60_BREAK := BARSLAST(CROSS(C, HHV(H, 40))) < 3; M60_RSI := SMA(MAX(C - REF(C, 1), 0), 14, 1) / SMA(ABS(C - REF(C, 1)), 14, 1) * 100; INTRA_SIG := TIME >= 93000 AND TIME <= 145500 AND VOL > DYNAINFO(8) * 0.45 AND (C - L) / (H - L + 0.0001) > 0.75 AND M60_BREAK AND M60_RSI > 62; {—————————————— 信号合成模块 ——————————————} STRATEGIC_SIGNAL := MONTH_SIG AND WEEK_SIG; TACTICAL_SIGNAL := DAY_SIG AND VOL_SIG; FINAL_SIGNAL := (STRATEGIC_SIGNAL AND TACTICAL_SIGNAL) OR INTRA_SIG;】】
最新发布
03-24
经过分析您的代码及报错提示,发现主要的问题出在 `{日线突破阈值}` 的部分。以下是具体的修复步骤: --- ### **问题分析** 1. 报错指出您在 `DAY_BREAK` 定义处存在问题:**"括号前写的不是函数、公式等,且缺少必要的运算符!"** 具体来看: ```plaintext DAY_BREAK := HHV(HIGH, 10); ``` 这里的语法实际上是正确的,但可能是某些环境对其解析有特殊限制。 2. 确保其他变量如 `HHV`, `HIGH` 和后续依赖项都已正确定义,并确认无遗漏的运算符或括号匹配错误。 --- ### **修正后的完整代码** ```plaintext {—————————————— 战略参数模块 ——————————————} {——周期协同参数——} MONTH_MA := MA(CLOSE, 20); {20月价值中枢} WEEK_VOL := MA(VOL, 5); {周量能基准} DAY_BREAK := HHV(HIGH, 10); {日线突破阈值} {——行业轮动因子V2——} IND_RPS := EMA(C / MAX(REF(C, 60), 0.001), 13) * 0.7 + RANK((VOL / MAX(MA(VOL, 60), 0.001)), COUNT()) * 0.3; {修正RANK参数写法,添加时间窗口参数COUNT()} IND_MOM := EMA(SLOPE(C, 21) * STD(C, 21), 8); INDUSTRY_WEIGHT := IF(IND_RPS > 75 AND MONTH_MA > REF(MONTH_MA, 3), 1.25, 1); {——估值体系V2——} DYNPETTM := IF(FINANCE(33) > 0, C / (FINANCE(33) / FINANCE(1) * INDUSTRY_WEIGHT), 1000); PB_RATE := C / (FINANCE(5) * IIF(INDUSTRY_WEIGHT > 1, 1.03, 1)); PEG_VAL := DYNPETTM / MAX((FINANCE(54) / FINANCE(34)) * 100 * INDUSTRY_WEIGHT, 0.001); {——波动率自适应V2——} IND_VOL := STD(INDEXC, 60) / MA(INDEXC, 60); VOLATILITY := STD(C, 60) / MA(C, 60) * 0.6 + IND_VOL * 0.4; VAR_PERIOD := IF(VOLATILITY < 0.05, 89, IF(VOLATILITY < 0.1, 55, 34)); {—————————————— 战术指标模块 ——————————————} {——三维MACD系统——} FAST_LEN := IF(VOLATILITY < 0.08, 12, 8); SLOW_LEN := IF(VOLATILITY < 0.08, 26, 17); SGNL_LEN := IF(VOLATILITY < 0.08, 9, 6); DIF := EMA(C, FAST_LEN) * INDUSTRY_WEIGHT - EMA(C, SLOW_LEN); DEA := EMA(DIF, SGNL_LEN); MACD := 2 * (DIF - DEA); MACD_ANGLE := ATAN((DIF - REF(DIF, 1)) / MAX(VOLATILITY, 0.03)) * 57.3; {——量价共振系统V2——} VOL_SIG := VOL > MA(VOL, 60) * 2.2 AND ABS(H - L) / REF(C, 1) > 0.08 AND C > MONTH_MA AND WEEK_VOL > REF(WEEK_VOL, 1) * 1.3; {——多周期验证——} MONTH_SIG := C > MONTH_MA AND MONTH_MA > REF(MONTH_MA, 1); WEEK_SIG := MACD > REF(MACD, 1) AND VOL > WEEK_VOL * 1.5; DAY_SIG := CROSS(C, DAY_BREAK) AND MACD_ANGLE > 30; {——盘中突击系统——} M60_BREAK := BARSLAST(CROSS(C, HHV(H, 40))) < 3; M60_RSI := SMA(MAX(C - REF(C, 1), 0), 14, 1) / SMA(ABS(C - REF(C, 1)), 14, 1) * 100; INTRA_SIG := TIME >= 93000 AND TIME <= 145500 AND VOL > DYNAINFO(8) * 0.45 AND (C - L) / (H - L + 0.0001) > 0.75 AND M60_BREAK AND M60_RSI > 62; {—————————————— 信号合成模块 ——————————————} STRATEGIC_SIGNAL := MONTH_SIG AND WEEK_SIG; TACTICAL_SIGNAL := DAY_SIG AND VOL_SIG; FINAL_SIGNAL := (STRATEGIC_SIGNAL AND TACTICAL_SIGNAL) OR INTRA_SIG; ``` --- ### **关键改动说明** 1. **针对 `DAY_BREAK` 的调整** 原来的定义已经基本正确,但由于环境中可能存在隐式约束,因此建议明确其上下文逻辑是否满足特定需求。 2. **`RANK` 函数的时间窗口补充** 在原代码中,`RANK` 缺少显式的排名范围(即需要指定统计区间),通过增加 `COUNT()` 参数解决此问题。 --- ### **测试与验证** 建议将上述修改应用到实际交易软件中运行调试,观察是否存在进一步异常情况。如果仍有问题,则需逐一排查各组件间的关联性以及数据输入的有效性。
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值