消失的运算符（和脑子）

最新推荐文章于 2024-07-25 21:08:08 发布

OneInDark

最新推荐文章于 2024-07-25 21:08:08 发布

阅读量173

点赞数 1

分类专栏： C++ 动态规划文章标签：算法 c语言

本文链接：https://blog.csdn.net/qq_42101694/article/details/120598262

版权

C++ 同时被 2 个专栏收录

547 篇文章 3 订阅

订阅专栏

动态规划

220 篇文章 2 订阅

订阅专栏

题目

题目描述
定义合法的表达式：

数字 $1,2,3,\dots,9$ 是合法表达式。
如果 $A, B$ 是合法表达式，那么 $(A)$ 和 $A + B$ 和 $A\times B$ 是合法的表达式。
以上规则可以构造出所有合法表达式。

现在有一个合法表达式，它的运算符都被拿掉了。你需要填入恰好 $k$ 个加号（剩下的就都是乘号），并计算出表达式的值；对于所有可能的还原后的表达式，输出其求和。

数据范围与提示
字符串长度 $n\leqslant 10^5$ ，很可能有括号。运算符数量 $m\leqslant 2.5\times 10^3$ ，都是你需要填入的。

思路

我知道了， $\sf OneInDark$ 是傻子！而且是傻的无可救药的那种！

我很快地写出了一个完全错误、八竿子打不着的 $\tt dp$ 式，然后把它优化。优化之后，我开始 犹豫代码怎么写，是用栈处理括号匹配吗？还是用函数做递归？用 $\rm vector$ 直接封装 $\rm struct$ 多项式吗？还是用数组存储下来？

想了很久才开始打。~~一边打一边骂娘~~。打到最后，发现样例过不了，输出中间变量调试，仔细想想——我求出了个寂寞啊？什么睿智 $\rm nt$ 可以写出这个 $j 8$ $\tt dp$ 式啊？

以上种种证据充分说明了这一点： $\sf OneInDark$ 就是个垃圾，还是有害垃圾。在此我衷心祝愿大家，写代码前想想思路的正确性，不要重蹈覆辙……

从头讲一下吧。括号可以递归处理，因为求 $\cdots+a\times (X)\times b+\cdots$ 的所有情况的和，根据分配率可知，就是 $\cdots+a\times (\sum X)\times b+\cdots$ 的结果。

如果递归，就要额外记录一维，表示用了多少个加号。粗略一算， $m$ 次合并，每次是 $\mathcal O(m^2)$ 的卷积，竟然是 $\mathcal O(m^3)$ 的，该死！尝试优化，失败；试着敲这个 $\mathcal O(m^3)$ 的暴力，没敲出来。

仔细想想，真的吗？一定要注意这个问题：事实上这是一个树上背包（总量是 $m$ ，任意合并），复杂度是 $\mathcal O(m^2)$ 的。~~也就是说，暴力敲出来，我就能过？~~

代码

#include <cstdio>
#include <iostream>
#include <cstring>
#include <vector>
#include <algorithm>
using namespace std;
typedef long long int_;
# define rep(i,a,b) for(int i=(a); i<=(b); ++i)
# define drep(i,a,b) for(int i=(a); i>=(b); --i)
inline int readint(){
	int a = 0; char c = getchar(), f = 1;
	for(; c<'0'||c>'9'; c=getchar())
		if(c == '-') f = -f;
	for(; '0'<=c&&c<='9'; c=getchar())
		a = (a<<3)+(a<<1)+(c^48);
	return a*f;
}
inline void writeint(int x){
	if(x > 9) writeint(x/10);
	putchar((x-x/10*10)^48);
}

const int MaxM = 2505;
int dp[MaxM<<1][MaxM];
int len[MaxM<<1];
vector< pair<int,int> > num; // <id,position>
vector< int > bracket;

const int Mod = 1e9+7;
int tmp[MaxM], tmp_len;
inline void copyTo(int p){
	len[p] = tmp_len;
	memcpy(dp[p],tmp,tmp_len<<2);
}
void mul(int a,int b){
	tmp_len = len[a]+len[b]-1;
	memset(tmp,0,tmp_len<<2);
	rep(i,0,len[a]-1) rep(j,0,len[b]-1)
		tmp[i+j] = (tmp[i+j]+int_(dp[a][i])*dp[b][j])%Mod;
}
void add(int a,int b){
	if(len[a] < len[b]) swap(a,b);
	tmp_len = len[a]; // max
	memcpy(tmp,dp[a],tmp_len<<2);
	rep(i,0,len[b]-1) // add up
		tmp[i] = (tmp[i]+dp[b][i])%Mod;
}
void shift(const int &a){
	drep(i,len[a],1)
		dp[a][i] = dp[a][i-1];
	++ len[a]; *dp[a] = 0;
}

const int MaxN = 100005;
char str[MaxN];
int sta[MaxM], top;
int main(){
	int n = readint(), k = readint();
	dp[MaxM][0] = len[MaxM] = 1;
	rep(i,MaxM+1,MaxM+(MaxM-5)){
		len[i] = len[i-1]+1;
		rep(j,dp[i][0]=1,len[i]-1)
			dp[i][j] = (dp[i-1][j-1]+dp[i-1][j])%Mod;
	}
	scanf("%s",str+1);
	str[0] = '(', str[n+1] = ')';
	for(int i=0,xyx=0; i<=n+1; ++i){
		if(str[i] == '-') continue;
		else if(str[i] == '(')
			bracket.push_back(i);
		else if(str[i] == ')'){
			const int bound = bracket.back();
			bracket.pop_back(); // match
			const int p = num.back().first;
			const int loca = num.back().second;
			num.pop_back(); // special processed
			const int sum = (MaxM<<1)-1; // sigma
			const int ans = sum-1; // result
			len[sum] = len[ans] = 1;
			*dp[sum] = *dp[ans] = 0; // init = 0
			add(sum,p), copyTo(sum); // sum = p
			int len_ans = len[p]; // for accuracy
			while(!num.empty() && num.back().second > bound){
				const int cur = num.back().first;
				num.pop_back(); // computed
				const int len_cur = len[cur];
				mul(ans,MaxM+len_cur), copyTo(ans);
				mul(sum,MaxM+len_cur-1), copyTo(p);
				shift(p); // the plus '+' on the left
				add(ans,p), copyTo(ans); // (l_0,r]
				mul(sum,cur), copyTo(sum);
				mul(cur,MaxM+len_ans-1), copyTo(cur);
				shift(cur); // the operator on the right
				add(sum,cur), copyTo(sum);
				len_ans += len_cur; // total operators
			}
			add(ans,sum), copyTo(p);
			num.push_back(make_pair(p,loca));
		}
		else{
			const int &p = ++ xyx;
			dp[p][0] = str[i]-'0'; len[p] = 1;
			num.push_back(make_pair(p,i));
		}
	}
	printf("%d\n",dp[num.back().first][k]);
	return 0;
}