2989:糖果
2989:糖果
这道题暴力加上一点优化可以过很多数据,正统做法是背包DP,不过若是直接套也不行,这里的糖果总数太多,时间和空间都不够用。所以在这里设f(i,j)=前i个糖果选一些构成j%k的最大值,函数里的j是mod k过的,这样优化后情况就会好很多。方程为f(i,j)=max(f(i,j),f(i-1,(j-a[i]%k+k)%k)+a[i]; 边界为d[i][0]=0。答案即为d[n][0]。
————————————————
两个关卡:一是初始化dp数组(一是初始化所有,二是初始化i=1,i=2这些最前面的有意义的状态,用于之后的递推),二是写出状态转移方程
dp数组初始化为负无穷
恰好装满背包 和 取得的糖果数恰好为k的倍数
初始化负无穷
恰好装满背包,只能由装了若干物品对应容量的背包 递推之后的状态,最后得到的dp【m】才是装满情况下的最大值
恰好模k余j,只能由取到了糖果的余数 递推过来 ,
最后得到的dp【n】【0】才是一步一步由有意义的余数状态推来的
如果初始化为0,假设一个不可能存在的比如dp【2】【10】(k取12,前两包糖果分别数量 3,4,dp【2】【10】不可能有意义),但
dp[i][j]=max(dp[i-1][j],dp[i-1][((j-a[i])%k+k)%k]+a[i]);
假设在状态dp【2】【10】的基础前提下取第i包糖果,dp【2】【10】初始值加上a【i】完全可能大于dp[i-1][j],而导致这种没有意义的取法,
所以把有意义的取法和没有意义的取法区分开来,没有意义的取法初始化为负无穷,有意义的取法初始化为非负数(当然了,只需要初始化一部分有意义的值,比如i=0,i=1,之后的状态由这些基本的有意义的状态递推),这样通过状态转移方程的max可以自动筛出不合意义的取法。
至于不初始化为负无穷,在选取状态时就要加上若干限制判条件,在背包问题例题中有讲,不再赘述。
初始化的两个步骤,第一步全部初始化为负无穷,第二步初始化i=1,i=2这些最前面的有意义的状态,用于之后的递推,要知道状态转移方程用于递推是非常强大的,很多时候,第二步仅仅初始化dp【0】【0】=0(一定有意义)就足以递推出所有
所以我想说的是,不求初始太多有意义的状态,只求不要把没有意义的也从无穷大初始化到了0
这题,取前0包糖,可能的、有意义的状态只能是dp【0】【0】,dp[0][1–k-1]这必定无意义哎
如下,初始化第二步角度众多,只要揪出一个有意义的就行,当然了,由于状态转移方程,从i=1,开始递推,必须要初始化一个i=0的某个状态的值,如果仅仅初始化for(int i=1;i<n;i++)dp[i][a[i]%k]=a[i];
,那必然是不对的,因为这只是有意义的状态,并不是i=1状态的最佳结果,最佳结果还得从i=0推过来
以上说的是最正经的,要从i=1开始推也不是不可以,dp[i][(k+j+a[i])%k]=max(dp[i-1][j]+a[i],max(dp[i-1][(k+j+a[i])%k],dp[i][(k+j+a[i])%k]));
i=1时,dp[i-1]对应的都是负无穷,总要由一个不是负无穷的、有意义的状态转化而来,那就是dp[1]的初始值了,说多了会混乱
#include <iostream>
#include <string.h>
using namespace std;
const int N=105;
int a[N];
int dp[N][N];
//dp[i][j]表示对于前i包糖果 满足取到的总数模k意义下为j 的最大糖果数
int main(){
int n,k;
cin>>n>>k;
for(int i=1;i<=n;i++){
cin>>a[i];//每包的糖果数
}
memset(dp,0xc0,sizeof(dp));
//1、初始化,
//a.对前0包糖(取0包糖)
//for(int i=0;i<k;i++)dp[0][i]=0;//❌
dp[0][0]=0;
//b.前i堆取得模k为a[i]%k的数目,可以包含上面那种啦,a[0]=0
for(int i=1;i<n;i++)dp[i][a[i]%k]=a[i];//这是一定取的到的
//c.前i堆取得模k为0的数目(直接不取嘛)
for(int i=0;i<=n;i++)dp[i][0]=0;
//2、状态转移方程
for(int i=1;i<=n;i++){
for(int j=0;j<k;j++){
// a.前i包的情况可以由 前i-1包哪些情况妆化而来
// 即余数为j的状态可以由哪些余数状态转化而来
//由状态转移方程可见,需要初始化
dp[i][j]=max(dp[i-1][j],dp[i-1][((j-a[i])%k+k)%k]+a[i]);
//b.前i-1包可以转化为前i包的哪些情况
// dp[i][(j+a[i])%k]=max(dp[i-1][j]+a[i],dp[i-1][(j+a[i])%k]);
}
}
cout<<max(0,dp[n][0]);
return 0;
}