Educational Codeforces Round #46 D. Yet Another Problem On a Subsequence(动态规划,组合数学)

原题链接:D. Yet Another Problem On a Subsequence
题目大意:给出n个数字组成的数组,规定一个数组[a1,a2,a3, … ,an],如果有a1=n-1且a1>0,这个数组就被叫为好数组,再规定“好序列”是可以分为若干个“好数组”的序列,现在问这n个数字中有多少个子序列是“好序列”
答案对998244353取模
数据范围:1≤n≤103,−109≤ai≤109
解题思路:
对于一个数ai,我们可以考虑从ai后的数组中选择ai个数,与ai一同构成好数组,共有Cnai种选择,考虑到一个“好序列”是可以分为若干个“好数组”的序列,我们可以从ai前面的数组中选择好数组进行拼接,再由乘法原理,即可得到答案
首先是组合数的计算,由于题目数据范围较小,可以由杨辉三角算出:

	c[0][0]=1;
	for(int i=1;i<=n;i++)c[i][0]=c[i][i]=1;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<i;j++)
		{
			c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
		}
	}

可以考虑用动态规划来解决问题,设置dp数组,其中dp[i]表示前i项中好数组的数量,设当前下标为i,对i+a[i]之后的数组进行更新,当前的数组以a[i]开头,可以在i~j的范围内选择a[i]个数组成子序列,然后与前面的好数组进行拼接,由dp[i-1]种方案,若不拼接,方案数为1种,共dp[i-1]+1种拼接方案,根据乘法原理
共有dp[j]+c[j-i][a[i]]*(dp[i-1]+1)种好数组,由此的到状态转移方程为

dp[j]=(dp[j]+1LL*c[j-i][a[i]]*(dp[i-1]+1))%mod;

AC代码:

#include<bits/stdc++.h>
using namespace std;
#define int long long
int mod=998244353;
int a[1005],dp[1005],c[1005][1005];
void init(int n)
{
	c[0][0]=1;
	for(int i=1;i<=n;i++)c[i][0]=c[i][i]=1;
	for(int i=1;i<=n;i++)
	{
		for(int j=1;j<i;j++)
		{
			c[i][j]=(c[i-1][j]+c[i-1][j-1])%mod;
		}
	}
}
signed main(){
	int n;
	cin>>n;
	init(n);//用杨辉三角计算组合数
	for(int i=1;i<=n;i++)cin>>a[i];
	for(int i=1;i<=n;i++)
	{
		if(a[i]>0)
		{
			for(int j=i+a[i];j<=n;j++)
			{
				dp[j]=(dp[j]+1LL*c[j-i][a[i]]*(dp[i-1]+1))%mod;
            	dp[j]%=mod;
			}
		}
	}
    cout<<dp[n]%mod;
	return 0;
}
### Codeforces Educational Round 168 Div. 2 Problem D 解析 #### 题目概述 题目描述涉及给定一个长度为 n 的数组 a 和整数 k,目标是在不超过 k 次操作的情况下使数组中的所有元素相等。每次操作可以选择任意数量的元素并将其增加或减少 1。 为了实现这一目标,需要找到一种策略来最小化所需的操作次数,使得最终所有的数组元素都相同[^1]。 #### 关键思路 核心在于理解到最优解通常是将所有元素变为中位数值。这是因为对于任何一组数字来说,调整至中位数所需的总移动量是最少的。具体而言: - 如果 n 是奇数,则存在唯一的中间值作为最佳目标; - 若 n 为偶数,则两个可能的目标分别是位于中心位置相邻的一对平均值之一; 考虑到最多可以执行 k 次修改动作,在某些情况下即使不是严格意义上的“中位数”,也可以通过适当选取接近它的其他值达到目的,只要满足条件即可。 #### 实现方法 以下是解决该问题的一种算法框架: ```cpp #include <bits/stdc++..h> using namespace std; int main() { int t; cin >> t; while (t--) { long long n, k; cin >> n >> k; vector<long long> a(n); for(auto &i : a) cin >> i; sort(a.begin(), a.end()); map<int,int> m; // 记录每个差异数出现频率 for(int i=0; i<n-1; ++i){ if(abs(a[i]-a[i+1])!=0)m[abs(a[i]-a[i+1])]++; } auto it=m.rbegin(); bool flag=true; while(k && !m.empty()){ if(it->second<=k){ k-=it->second; m.erase((++it).base()); }else{ cout << "NO\n"; flag=false;break; } it=m.rbegin(); } if(flag || m.size()==0)cout<<"YES\n"; } } ``` 此代码片段实现了上述逻辑的一部分,即尝试尽可能多地消除差异最大的一对元素之间的差距直到无法继续为止。需要注意的是这段代码并不完全适用于所有情况下的直接提交测试,而是提供了一个基本的方向指引。
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值