区间DP入门+优化

区间DP原理

先对小区间进行操作得到局部最优解， 然后通过小区间的最优解来得到大区间的最优解。利用了最优子结构和无后继影响的特性。

DP模板

#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
ll dp[2000][2000];
int main()
{
	init();  // 初始化DP
	for(int len=1;len<=n;len++)   // 遍历子区间的所有可能的长度，从只有一个数到全部（n个数）
		for(int begin=1;begin+len<limit;begin++) //遍历所有起点，其中limit是题目中理解出的范围，环形可能是2*n,线性可能是n
		{
			int end=begin+len;
			for(int k=begin;k<end;k++)  //k代表着分割点，[begin,k]是第一个子区间，(k,end]是第二个子区间  k包含在前一个区间内，所以是<end   
				dp[begin][end]=max(dp[begin][end],dp[begin][k]+dp[k+1][end]+sth)  //这里的sth代表着对两个子区间合并和产生的代价 
		}
	int ans=0;
	for(int i=0;i<n;i++)
		ans=max(ans,dp[i][i+n]);
	cout<<ans;
}

例题

P1880 [NOI1995]石子合并

传送门

题目描述
在一个圆形操场的四周摆放 N 堆石子,现要将石子有次序地合并成一堆.规定每次只能选相邻的2堆合并成新的一堆，并将新的一堆的石子数，记为该次合并的得分。

试设计出一个算法,计算出将 N 堆石子合并成 1 堆的最小得分和最大得分。

输入格式
数据的第 1 行是正整数 N，表示有 N 堆石子。

第 2 行有 N 个整数，第 i 个整数 a_i表示第 i 堆石子的个数。
输出格式
输出共 2 行，第 1 行为最小得分，第 2 行为最大得分。

输入输出样例
输入
4
4 5 9 4
输出
43
54
说明/提示
1\leq N\leq 1001≤N≤100，0\leq a_i\leq 200≤a
i≤20。

#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio> 
using namespace std;
typedef long long ll;
ll s[500];
ll dpmin[500][500];
ll dpmax[500][500];
ll sum[500];
int main()
{
	int n;
	cin>>n;
	memset(dpmin,0x3f3f3f,sizeof(dpmin));
	for(int i=0;i<n;i++)
	{
		scanf("%d",&s[i]);
		sum[i+1]=sum[i]+s[i];
		dpmin[i][i]=0;
        dpmax[i][i]=0;
	}
	for(int i=0;i<n;i++)
	{
		sum[i+n+1]=sum[i+n]+s[i];
		dpmin[i+n][i+n]=0;
        dpmax[i+n][i+n]=0;
	}//想要获取 区间 i->j 的区间和，sum下标应该为sum[j+1]-sum[i] 
	for(int len=1;len<=n;len++)
	{
		for(int begin=0;begin+len<2*n;begin++)
		{
			int end=begin+len;
			for(int k=begin;k<end;k++)//k包含于前一个区间
			{
				dpmin[begin][end]=min(dpmin[begin][end],dpmin[begin][k]+dpmin[k+1][end]+sum[end+1]-sum[begin]);
				dpmax[begin][end]=max(dpmax[begin][end],dpmax[begin][k]+dpmax[k+1][end]+sum[end+1]-sum[begin]);
			}
		}
	}
	ll ansmin=0xfffffff;
    ll ansmax=-1;
    for(int i=0;i<n;i++)
	{
        ansmin=min(ansmin,dpmin[i][i+n-1]);
        ansmax=max(ansmax,dpmax[i][i+n-1]);
    }
    cout<<ansmin<<endl;
    cout<<ansmax<<endl;
    return 0;
}

P1063 能量项链

传送门

题目描述
在Mars星球上，每个Mars人都随身佩带着一串能量项链。在项链上有N颗能量珠。能量珠是一颗有头标记与尾标记的珠子，这些标记对应着某个正整数。并且，对于相邻的两颗珠子，前一颗珠子的尾标记一定等于后一颗珠子的头标记。因为只有这样，通过吸盘（吸盘是Mars人吸收能量的一种器官）的作用，这两颗珠子才能聚合成一颗珠子，同时释放出可以被吸盘吸收的能量。如果前一颗能量珠的头标记为m，尾标记为r，后一颗能量珠的头标记为r，尾标记为n，则聚合后释放的能量为m * r * n（Mars单位），新产生的珠子的头标记为m，尾标记为n。

需要时，Mars人就用吸盘夹住相邻的两颗珠子，通过聚合得到能量，直到项链上只剩下一颗珠子为止。显然，不同的聚合顺序得到的总能量是不同的，请你设计一个聚合顺序，使一串项链释放出的总能量最大。

例如：设N=4，颗珠子的头标记与尾标记依次为(2,3) (3,5) (5,10) (10,2)我们用记号⊕表示两颗珠子的聚合操作，(jj⊕kk)表示第j,kj,k两颗珠子聚合后所释放的能量。则第4、1两颗珠子聚合后释放的能量为：

(4⊕1)=10×2×3=60。

这一串项链可以得到最优值的一个聚合顺序所释放的总能量为：

((4⊕1)⊕2)⊕3）=10×2×3+10×3×5+10×5×10=710。

输入格式
第一行是一个正整数N(4≤N≤100)，表示项链上珠子的个数。第二行是N个用空格隔开的正整数，所有的数均不超过1000。第i个数为第i颗珠子的头标记(1≤i≤N)，当i<N时，第i颗珠子的尾标记应该等于第i+1颗珠子的头标记。第N颗珠子的尾标记应该等于第1颗珠子的头标记。

至于珠子的顺序，你可以这样确定：将项链放到桌面上，不要出现交叉，随意指定第一颗珠子，然后按顺时针方向确定其他珠子的顺序。

输出格式
一个正整数E(E≤2.1×(10) ^9)，为一个最优聚合顺序所释放的总能量。

输入输出样例
输入 #1复制
4
2 3 5 10
输出 #1复制
710
说明/提示
NOIP 2006 提高组 第一题

#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
ll dp[400][400];
ll stone[400];
int main()
{
	int n;
	cin>>n;
	for(int i=0;i<n;i++)
	{
		cin>>stone[i];
		stone[i+n]=stone[i];
	}
	for(int len=1;len<=n;len++)
	{
		for(int be=0;be+len<2*n;be++)
		{
			int end=be+len;
			for(int k=be;k<end;k++)
				dp[be][end]=max(dp[be][end],dp[be][k]+dp[k+1][end]+stone[be]*stone[k+1]*stone[end+1]);
		}
	}
	ll ans=0;
	for(int i=0;i<n;i++)
		ans=max(ans,dp[i][i+n-1]);
	cout<<ans;
}

平行四边形优化

	init()
	for(len=1;len<=n;len++)
    	for(i=1;i+len<limit;i++)
		{
        	j=i+len;
        	for(k=s[i][j-1];k<s[i+1][j];k++)//在最优分割点范围内枚举分割点
            	if(dp[i][j]>dp[i][k]+dp[k+1][j]+sth)
				{
                	dp[i][j]=dp[i][k]+dp[k+1][j]+sth;
                	s[i][j]=k;//更新最佳分割点
            	}
    	}