区间dp

题目列表：P1040 hdu3506 hdu4632 hdu4745 hdu2476 hdu4283 poj2955 poj1651

P1040加分二叉树

#include<bits/stdc++.h>
using namespace std;
#define ll long long
int n;
ll dp[30][30];
int root[30][30];

int flag=0;
void print(int l,int r){
	if(l>r)return ;
	if(flag)putchar(' ');
	flag=1;
	printf("%d",root[l][r]);
		
	if(l==r)return ;
	
	print(l,root[l][r]-1);
	print(root[l][r]+1,r);
}

int main(){
	scanf("%d",&n);
	for(int i=1;i<=n;i++){
		scanf("%lld",&dp[i][i]);
		root[i][i]=i;
	}
	
	for(int i=1;i<=n;i++){//枚举长度 
		for(int l=1;l<=n;l++){//左端点 
			int r=l+i;//右端点 
			if(r>n)break;

			for(int j=l;j<=r;j++){//枚举根节点 
				if(j==l&&dp[l][r]<dp[j][j]+dp[j+1][r]){//从无左孩子 
					root[l][r]=j;
					dp[l][r]=dp[j][j]+dp[j+1][r];
				}
				if(j==r&&dp[l][r]<dp[l][j-1]+dp[j][j]){//枚举至无右孩子 
					root[l][r]=j;
					dp[l][r]=dp[l][j-1]+dp[j][j];
				}
				if(j!=l&&j!=r&&dp[l][r]<dp[l][j-1]*dp[j+1][r]+dp[j][j]){
					root[l][r]=j;
					dp[l][r]=dp[l][j-1]*dp[j+1][r]+dp[j][j];
				}
			}
		}
	}
	
	cout<<dp[1][n]<<endl;
	print(1,n);

	return 0;
} 

hdu-3506

状态转移方程：dp[i][j]=min{dp[i][k],dp[k+1][j]}+w[i][j](i≤k<j)（min也可以改为max）  O(n^3)复杂度

上述的dp[i][j]表示区间[i,j]上的某个最优值。w[i][j]表示在转移时需要额外付出的代价。该方程的时间复杂度为O(N3)。

区间包含的单调性：如果对于 i≤i'<j≤j'，有 dp[i'][j]≤w[i][j']，那么说明w具有区间包含的单调性。（如果小区间包含于大区间中，那么小区间的w值不超过大区间的w值）

四边形不等式：如果对于 i≤i'<j≤j'，有 dp[i][j]+w[i'][j']≤w[i'][j]+w[i][j']，我们称函数w满足四边形不等式。（两个交错区间的w的和不超过小区间与大区间的w的和）

定理一：如果上述的 w 函数同时满足区间包含单调性和四边形不等式性质，那么函数 m 也满足四边形不等式性质

我们再定义 s[i][j]表示 m[i][j]取得最优值时对应的下标（即 i≤k≤j 时，k 处的 w 值最大，则 s[i][j]=k）。此时有如下定理

定理二：假如 m[i][j] 满足四边形不等式，那么 s[i][j] 单调，即 s[i][j]≤s[i][j+1]≤s[i+1][j+1]。

优化后可将转移方程变为：dp[i][j]=min{dp[i][k],dp[k+1][j]}+w[i][j](s[i][j-1]≤k<s[i+1][j])（min也可以改为max）  接近O(n^2)复杂度

#include<bits/stdc++.h>
using namespace std;
#define MAXN 2005
#define inf 0x3f3f3f3f
int n;
int sum[MAXN],num[MAXN];
int dp[MAXN][MAXN];
int s[MAXN][MAXN];
int main(){
    while(cin>>n){
        memset(dp,0x3f,sizeof(dp));
        
        for(int i=1;i<=n;i++){
            scanf("%d",&num[i]);
            num[i+n]=num[i];
        }
        
        for(int i=1;i<=2*n-1;i++){
            dp[i][i]=0;
            s[i][i]=i;
            sum[i]=sum[i-1]+num[i];
        }
        
        for(int len=2;len<=n;len++){
            for(int i=1;i<=n*2-1;i++){
                int l=i,r=i+len-1;
                if(r>n*2-1)break;
                for(int k=s[l][r-1];k<=s[l+1][r];k++){
                    if(dp[l][r]>dp[l][k]+dp[k+1][r]+sum[r]-sum[l-1]){
                        dp[l][r]=dp[l][k]+dp[k+1][r]+sum[r]-sum[l-1];
                        s[l][r]=k;
                    }
                }
            }
        }
        
        int ans=inf;
        for(int i=1;i<=n;i++){
            ans=min(ans,dp[i][i+n-1]);
        }
        cout<<ans<<endl;
    }

    return 0;
}

hdu4632

 

hdu4745

 

hdu2476

 

hdu4283

 

poj2955

 

poj1651