codeforces round518 div1A Array Without Local Maximums(dp)

题目链接
感谢大佬myx12345

题意

给出一段长度为n的序列,序列中的值都在1到200之间序列满足

  1. a 1 ≤ a 2 a_1 \leq a2 a1a2
  2. a n ≤ a n − 1 a_n \leq a_{n-1} anan1
  3. a i ≤ m a x { a i − 1 , a i + 1 } a_i \leq max\{a_{i-1},a_{i+1}\} aimax{ai1,ai+1}

此时,小明他打翻了一瓶墨水,使得序列中的某些数被墨水遮挡住了(有没有似曾相识的感觉),我们需要猜测出符合条件的序列有多少种,答案模998244353。

题解

a i − 1 a_{i-1} ai1 a i a_i ai有三种关系, &lt; , = , &gt; &lt;, = ,&gt; <,=,>。 这里设 d p [ i ] [ j ] [ k ] dp[i][j][k] dp[i][j][k]代表第 i i i个位置为 j j j且与前面 a i − 1 a_{i-1} ai1的关系为 k k k的方案数。
那么状态转移方程为

  1. d p [ i ] [ j ] [ 0 ] = ∑ k = 1 j ( d p [ i − 1 ] [ k ] [ 0 ] + d p [ i − 1 ] [ k ] [ 1 ] + d p [ i − 1 ] [ k ] [ 2 ] ) dp[i][j][0] = \sum_{k=1}^{j} (dp[i-1][k][0]+dp[i-1][k][1]+dp[i-1][k][2]) dp[i][j][0]=k=1j(dp[i1][k][0]+dp[i1][k][1]+dp[i1][k][2])
    因为 a i − 1 &lt; a i a_{i-1}&lt;a_i ai1<ai时,那 a i − 1 a_{i-1} ai1可以和 a i − 2 a_{i-2} ai2满足任何条件。
  2. d p [ i ] [ j ] [ 1 ] = d p [ i − 1 ] [ j ] [ 0 ] + d p [ i − 1 ] [ j ] [ 1 ] + d p [ i − 1 ] [ j ] [ 2 ] dp[i][j][1] = dp[i-1][j][0]+dp[i-1][j][1]+dp[i-1][j][2] dp[i][j][1]=dp[i1][j][0]+dp[i1][j][1]+dp[i1][j][2]
    这里与1的不同之处就是 a i − 1 = a i a_{i-1} = a_i ai1=ai
  3. d p [ i ] [ j ] [ 2 ] = ∑ k = j 200 ( d p [ i − 1 ] [ k ] [ 2 ] + d p [ i − 1 ] [ k ] [ 1 ] ) dp[i][j][2] = \sum_{k=j}^{200} (dp[i-1][k][2]+dp[i-1][k][1]) dp[i][j][2]=k=j200(dp[i1][k][2]+dp[i1][k][1])
    因为 a i − 1 &gt; a i a_{i-1}&gt;a_{i} ai1>ai,所以必须要 a i − 1 ≤ a i − 2 a_{i-1} \leq a_{i-2} ai1ai2

如果暴力计算前缀和,后缀和的话时间复杂度为 O ( n ∗ 200 ∗ 200 ) O(n*200*200) O(n200200),可以边计算边记录前缀和和后缀和,那么时间复杂度为 O ( n ∗ 200 ) O(n*200) O(n200)

代码

#include <bits/stdc++.h>
using namespace std;
const int maxn = 1e5+5;
const int mod = 998244353;
int a[maxn];
long long dp[maxn][205][3];
// 0: a_i-1 < a_i   1: =   2: >
int main() {
	int n;
	scanf("%d", &n);
	for(int i = 0; i < n; ++i)
		scanf("%d", &a[i]);
	if(a[0] != -1) 
		dp[0][a[0]][0] = 1;
	else {
		for(int i = 1; i <= 200; ++i)
			dp[0][i][0] = 1;
	}
	for(int i = 1; i < n; ++i) {
		int sum = 0;
		for(int j = 1; j <= 200; ++j) {
			if(a[i] == -1 || j == a[i]) {
				dp[i][j][0] = sum;
			}
			
			sum = (sum+dp[i-1][j][1]+dp[i-1][j][0]+dp[i-1][j][2])%mod;
		}
		sum = 0;
		for(int j = 1; j <= 200; ++j) {
			if(a[i] == -1 || j == a[i]) {
				dp[i][j][1] = (dp[i-1][j][0]+dp[i-1][j][1]+dp[i-1][j][2])%mod;
			}
			
		}
		for(int j = 200; j >= 1; --j) {
			if(a[i] == -1 || j == a[i]) {
				dp[i][j][2] = sum;
			}
			sum = (sum+dp[i-1][j][2]+dp[i-1][j][1])%mod;
		}
	}
	long long ans = 0;
	for(int i = 1; i <= 200; ++i) {
		ans = (ans+dp[n-1][i][1]+dp[n-1][i][2])%mod;
	}
	printf("%lld\n", ans);
	return 0;
}
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值