AtCoder Beginner Contest 222 D - Between Two Arrays(前缀和优化dp)

linkkkkk

题意:

给出两个长度为 n n n的序列 a , b a,b a,b,构造一个序列 c c c满足 a i < = c i < = b i a_i<=c_i<=b_i ai<=ci<=bi并且 c c c非单调递减。求方案数。 n < = 3000 , 0 < = a i , b i < = 3000 n<=3000,0<=a_i,b_i<=3000 n<=3000,0<=ai,bi<=3000

思路:

数据范围比较小,考虑 d p dp dp
d p [ i ] [ j ] dp[i][j] dp[i][j]表示选完了前 i i i个数并且第 i i i个数的值为 j j j的方案数。
转移有 d p [ i ] [ j ] = ∑ d p [ i − 1 ] [ k ] , a i − 1 < = k < = j dp[i][j]=\sum{dp[i-1][k]},a_{i-1}<=k<=j dp[i][j]=dp[i1][k],ai1<=k<=j
用前缀和优化一下,记 s u m [ i ] [ j ] sum[i][j] sum[i][j]表示选了 i i i个数并且第 i i i个数的值 < = j <=j <=j的方案数。
转移就变成了 d p [ i ] [ j ] = s u m [ i − 1 ] [ j ] − s u m [ i − 1 ] [ a [ i − 1 ] − 1 ] dp[i][j]=sum[i-1][j]-sum[i-1][a[i-1]-1] dp[i][j]=sum[i1][j]sum[i1][a[i1]1]
有几个小细节:
i − 1 i-1 i1层的 s u m sum sum只维护到了 b i − 1 b_{i-1} bi1,所以如果 j > b i − 1 j>b_{i-1} j>bi1的话,是没有值的,应该取 m i n ( j , b i − 1 ) min(j,b_{i-1}) min(j,bi1)
a i , b i a_i,b_i ai,bi都是可以为 0 0 0的话,如果为 0 0 0的话, a [ i − 1 ] − 1 a[i-1]-1 a[i1]1就会变成 − 1 -1 1,但是不能跟 0 0 0 m a x max max,不然就会漏掉 0 0 0的部分。我的处理方式是在输入的时候让所有的 a , b a,b a,b都增加 1 1 1,这样不会影响相对大小,下标也不是 − 1 -1 1

代码:

// Problem: D - Between Two Arrays
// Contest: AtCoder - Exawizards Programming Contest 2021(AtCoder Beginner Contest 222)
// URL: https://atcoder.jp/contests/abc222/tasks/abc222_d
// Memory Limit: 1024 MB
// Time Limit: 2000 ms
// 
// Powered by CP Editor (https://cpeditor.org)

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;typedef unsigned long long ull;
typedef pair<ll,ll>PLL;typedef pair<int,int>PII;typedef pair<double,double>PDD;
#define I_int ll
inline ll read(){ll x=0,f=1;char ch=getchar();while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}while(ch>='0'&&ch<='9'){x=x*10+ch-'0';ch=getchar();}return x*f;}
#define read read()
#define rep(i, a, b) for(int i=(a);i<=(b);++i)
#define dep(i, a, b) for(int i=(a);i>=(b);--i)
ll ksm(ll a,ll b,ll p){ll res=1;while(b){if(b&1)res=res*a%p;a=a*a%p;b>>=1;}return res;}
const int maxn=4e5+7,maxm=1e6+7;
const ll mod=998244353;

ll dp[3100][3100],n,a[3100],b[3100];
ll sum[3100][3100];

int main(){
	n=read;
	rep(i,1,n) a[i]=read+1;
	rep(i,1,n) b[i]=read+1;
	for(int i=a[1];i<=b[1];i++) dp[1][i]=1,sum[1][i]=sum[1][i-1]+dp[1][i];
	for(int i=2;i<=n;i++)
		for(int j=a[i];j<=b[i];j++){
			if(j>=a[i-1])
				dp[i][j]=(dp[i][j]+sum[i-1][min(j*1ll,b[i-1])]-sum[i-1][max(0ll,a[i-1]-1)]+mod)%mod;
			sum[i][j]=(sum[i][j-1]+dp[i][j])%mod;
		}
	ll ans=0;
	for(int i=a[n];i<=b[n];i++)
		ans=(ans+dp[n][i])%mod;
	cout<<ans;
	return 0;
}
  • 0
    点赞
  • 1
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

豆沙睡不醒

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值