一、区间dp:
区间dp概念
1.区间dp就是在区间上进行动态规划,求解一段区间上的最优解。主要是通过合并小区间的 最优解进而得出整个大区间上最优解的dp算法。
2.他和线性dp在列状态方程的区别是前者考虑的是从第一个到最后一个
3.解题思路:
求解在一个区间上的最优解,那么我把这个区间分割成一个个小区间,求解每个小区间的最优解,再合并小区间得到大区间即可。所以在代码实现上,我可以枚举区间长度len为每次分割成的小区间长度(由短到长不断合并),内层枚举该长度下可以的起点,自然终点也就明了了。然后在这个起点终点之间枚举分割点,求解这段小区间在某个分割点下的最优解。
4.常用模板
for(int len = 1;len<=n;len++){//枚举长度
for(int j = 1;j+len<=n+1;j++){//枚举起点,ends<=n
int ends = j+len - 1;
for(int i = j;i<ends;i++){//枚举分割点,更新小区间最优解
dp[j][ends] = min(dp[j][ends],dp[j][i]+dp[i+1][ends]+something);
}
}
}
5.第一种类型
石子合并1:有N堆石子,现要将石子有序的合并成一堆,规定如下:每次只能移动任意的2堆石子合并,合并花费为新合成的一堆石子的数量。求将这N堆石子合并成
#include <iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
#define INF 0x3f3f3f
int stone[105];
int dp[105][105];
int sum[105];
int main()
{
int n;
scanf("%d",&n);
memset(sum,0,sizeof(sum));
memset(dp,INF,sizeof(dp));
for(int i =1;i<=n;i++){
scanf("%d",&stone[i]);
sum[i] = sum[i - 1] + stone[i];//重量
dp[i][i] = 0;
}
for(int len = 1;len<=n;len++){//枚举长度
for(int j = 1;j+len<=n+1;j++){//枚举起点,ends<=n
int ends = j+len - 1;//枚举终点
for(int i = j;i<ends;i++){//枚举分割点
dp[j][ends] = min(dp[j][ends],dp[j][i]+dp[i+1][ends]+sum[ends]-sum[j-1]);//更新状态
}
}
}
cout<<dp[1][n]<<endl;
return 0;
}
前缀和:sum[i]=sum[i-1]+w[i];
区间和:sum[i]-sum[k-1] (k<i)(注意这里是k-1)
石子合并2:原题与上面相同,但是石子排列由线性排列变成环状排列
思路:环状以后合并区间的情况就可以从后往前合并,最后合并完成可能是1n,2n1,3n~2…这种n个石子合并的情况。所以我们可以破环成链,将前n-1各元素也放到n后面构成一个线性的环状序列,在对这个序列dp即可
#include <iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
using namespace std;
#define INF 0x3f3f3f
int stone[105];
int dpmin[205][205];//最小
int dpmax[205][205];//最大
int sum[205];
int main()
{
int n;
scanf("%d",&n);
memset(sum,0,sizeof(sum));
memset(dpmin,INF,sizeof(dpmin));
memset(dpmax,-1,sizeof(dpmax));
for(int i =1;i<=n;i++){
scanf("%d",&stone[i]);
sum[i] = sum[i - 1] + stone[i];
dpmin[i][i] = 0;
dpmax[i][i] = 0;
}
for(int i = 1;i<=n;i++){
sum[i+n] = sum[i+n-1]+stone[i];//展开的n后面的n-1~1重量
dpmin[i+n][i+n] = 0;
dpmax[i+n][i+n] = 0;
}
for(int len = 1;len<=n;len++){//长度还是最大n
for(int j = 1;j+len<=2*n;j++){//起点枚举最大到2*n-1,ends<=2*n-1
int ends = j+len - 1;
for(int i = j;i<ends;i++){//注意!i<ends!!!因为i=ends时,dp[ends+1][ends]是不成立的!
dpmin[j][ends] = min(dpmin[j][ends],dpmin[j][i]+dpmin[i+1][ends]+sum[ends]-sum[j-1]);
dpmax[j][ends] = max(dpmax[j][ends],dpmax[j][i]+dpmax[i+1][ends]+sum[ends]-sum[j-1]);
}
}
}
int ansmin = 0xfffffff;
int ansmax = -1;
for(int i = 1;i<=n;i++){
ansmin = min(ansmin,dpmin[i][i+n-1]);//找1~n,2~n~1,3~n~2....的合并n个堆的中最大和最小的值
ansmax = max(ansmax,dpmax[i][i+n-1]);
}
cout<<ansmin<<endl;
cout<<ansmax<<endl;
return 0;
}
转换成环方法:
for(int i = 1;i<=n;i++)
{
a[i+n]=a[i];
dp[i+n] [i+n]=dp[i] [i]=???;
}
6.第二种题型
题意:给出一个的只有’(’,’)’,’[’,’]'四种括号组成的字符串,求 最多 有多少个括号满足题目里所描述的完全匹配
#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int dp[105][105];
int main()
{
char s[105];
while(scanf("%s",s+1)!=EOF)
{
memset(dp,0,sizeof(dp));//dp初始化为0,因为一方面是找最大之,一方面初始匹配数为0
int len = strlen(s+1);//dp[i][i]不用处理,因为自己和自己不匹配就是0
if(s[1]=='e')break;
for(int l = 1;l<=len;l++){
for(int i = 1;i+l<=len+1;i++){
int j= i+l-1;
if((s[i]=='('&&s[j]==')')||(s[i]=='['&&s[j]==']')){//如果匹配,先更新
dp[i][j] = dp[i+1][j-1]+2;
}
for(int k = i;k<j;k++){//k<j
dp[i][j] = max(dp[i][j],dp[i][k]+dp[k+1][j]);
}
}
}
cout<<dp[1][len]<<endl;
}
return 0;
}
题意:给你n个数字,要求不能删除两端点的数字,然后删除其他数字的代价是该数字和左右相邻数字的乘积,问把数字(除端点)删完后的最小总代价。
#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
int dp[105][105];
int num[105];
int main()
{
int N;
cin>>N;
memset(dp,0x3f,sizeof(dp));//dp初始化为inf,因为找最小值
for(int i = 1;i<=N;i++){
cin>>num[i];
dp[i][i] = 0;//dp[i][i]要初始化,作为由小到大累计的初始条件,自己取出自己就是为0
}
for(int len = 1;len<=N;len++){
for(int i = 2;i+len<=N+1;i++){//起点从2开始,因为不包括两端点
int j = i+len-1;
for(int k = i;k<j;k++){
dp[i][j] = min(dp[i][j],dp[i][k]+dp[k+1][j]+num[i-1]*num[k]*num[j]);
}
}
}
cout<<dp[2][N]<<endl;//由dp[i][j]的定义,输出合并第2~N-1个卡片的结果
return 0;
}
题意:给出两个整数 n , m ,要求在 n 中加入m - 1 个乘号,将n分成m段,求出这m段的最大乘积
#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
ll dp[50][50];
ll num[50][50];
int main()
{
int T;
scanf("%d",&T);
char s[100];
while(T--)
{
int m;
scanf("%s%d",s+1,&m);
int len = strlen(s+1);
memset(dp,0,sizeof(dp));//初始化dp
memset(num,0,sizeof(num));
for(int i = 1; i<=len; i++)
{
for(int j = i; j<=len; j++)
{
for(int k = i;k<=j;k++){
num[i][j]*=10;
num[i][j]+=(s[k]-'0');
}
}
dp[i][0] = num[1][i];//初始插入0个乘号时是自己,由小到大累计的基础
}
for(int j = 1;j<m;j++){//乘号个数由1~m-1个
for(int i = 1;i<=len;i++){//结尾坐标(表示长度)
for(int k = 1;k<i;k++){//分割点插入第j个乘号
dp[i][j] = max(dp[i][j],dp[k][j-1]*num[k+1][i]);
}
}
}
cout<<dp[len][m-1]<<endl;//输出在1~len插入m-1个乘号的结果
}
return 0;
}
题意:给你一个字符串,求出其最多的可构成的回文字串(不要求连续),注:这里不同的回文字串只要求位置不同即可视为不同,如:aaaaa 的最多回文子串数目是 31.
思路:用dp[i][j]表示状态,表示i~j里最多的回文字串数目
#include <iostream>
#include<cstdio>
#include<cstring>
using namespace std;
typedef long long ll;
char s[1005];
int dp[1005][1005];
int main()
{
int T;
scanf("%d",&T);
int t = 0;
while(T--){
t++;
scanf("%s",s+1);
memset(dp,0,sizeof(dp));
int len = strlen(s+1);
for(int i = 1;i<=len;i++)dp[i][i] = 1;//自己是长度为1的
for(int l = 1;l<=len;l++){
for(int i = 1;i+l<=len+1;i++){
int j = i+l-1;//容斥原理 dp[i][j] = dp[i+1][j]并dp[i][j-1] (因为是自序列可以不连续)
dp[i][j] = (dp[i+1][j] + dp[i][j-1] - dp[i+1][j-1]+10007)%10007;
if(s[i] == s[j])dp[i][j]=(dp[i][j]+dp[i+1][j-1] + 1)%10007;//如果两端相等,dp[i][j] = 原来的 + 两端与中间每一个回文也可以构成回文(dp[i+1][j-1]) + 两端单独构成一个回文(1)!!!
}
}
printf("Case %d: %d\n",t,dp[1][len]);
}
return 0;
}
7.第三种类型
能量链球
#include<stdio.h>
#include<string.h>
#define Max(a,b) a>b?a:b
int main()
{
long long int dp[202][202];
long long int s[202];
memset(dp,0,sizeof(dp));
long long int sum=0;
int n,i,j,k,len;
scanf("%d",&n);
for(i=1;i<=n;i++)
scanf("%lld",&s[i]),s[i+n]=s[i];
for(len=2;len<=n;len++)
for(i=1;i+len-1<2*n;i++)
{ j=i+len-1;
for(k=i;k<j;k++)
dp[i][j]=Max(dp[i][j],dp[i][k]+dp[k+1][j]+s[i]*s[k+1]*s[j+1]);
}
for(i=1;i<=n;i++)
if(sum<dp[i][i+n-1])sum=dp[i][i+n-1];
printf("%lld\n",sum);
return 0;
}
小屌丝
#include <iostream>
#include <cstdio>
#include <cstring>
using namespace std;
const int INF=0x3f3f3f3f;
int n,D[120],sum[120];
int dp[120][120];
int main()
{
int T_T,cas=1;
scanf("%d",&T_T);
while(T_T--)
{
scanf("%d",&n);
sum[0]=0;
memset(dp,0,sizeof(dp));
for(int i=1;i<=n;i++)
{
scanf("%d",D+i);
sum[i]=sum[i-1]+D[i];
}
for(int i=0;i<120;i++) for(int j=i+1;j<120;j++) dp[i][j]=INF;
for(int len=1;len<=n;len++)
{
for(int i=1;i+len-1<=n;i++)
{
int j=i+len-1;
for(int k=1;k<=len;k++)
{
dp[i][j]=min(dp[i][j],dp[i+1][i+k-1]+(k-1)*D[i]+dp[i+k][j]+k*(sum[j]-sum[i+k-1]));
}
}
}
printf("Case #%d: %d\n",cas++,dp[1][n]);
}
return 0;
}