区间DP

对于某个求区间最优解的问题，将大区间划分为一个个小区间通过求解小区间最优解从而得到大区间最优解的问题。
核心代码：

    for(int len=2;len<=n;len++){//枚举区间长度
        for(int i=1;i<n && i+len-1<=n;i++){//枚举区间左端点
            int j=i+len-1;
            for(int k=i;k<j;k++){//枚举断点位置 在K和k+1之间断开
                dp[i][j]=max(dp[i][j],dp[i][k]+dp[k+1][j]+valu[k][k+1]);//dp[i][j]表示从i到j的最优解
            }
        }
    }

时间复杂度O(n^3)

---

一些例题

Light-Oj-1422
问题描述：有n个舞会要去参加，每个舞会都规定了要穿衣服的编号，你可以把某件衣服套在别的衣服上，也可以脱掉，求最小穿衣服的次数。

/*********************************************************************/
/**********************************cyf*******************************/
/*******************************************************************/
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<string>
#include<vector>
#include<set>
#include<map>
#include<queue>
#include<stack>
#include<utility>
#include<time.h>
using namespace std;
#define ll long long
#define repp(i,x) for(int i=1;i<=x;i++)
#define repx(i,n,j) for(int i=j;i<n;i++)
#define rep(i,begin,end) for(int i=begin;i<=end;i++)
#define mst(a,b) memset(a,b,sizeof(a))
#define drep(i,a,b) for(int i=a;i>=b;i--)

int gcd(int a,int b){
    int t;
    if(a<b){
        t=a;
        a=b;
        b=t;
    }
    while(a%b)
    {
        t=a%b;
        a=b;
        b=t;
    }
    
    return b;
}


const int maxn=1e2+5;
const int MAXN=1e4+5;
const ll mod = 1e9+7;
const int INF=0x3f3f3f3f;
const ll LLINF=0x3f3f3f3f3f3f3f;
int arr[maxn];
int dp[maxn][maxn];
int main(){
    int t,kase=0,n;
    scanf("%d",&t);
    while(t--){
        kase++;
        scanf("%d",&n);
        rep(i,1,n){
            scanf("%d",&arr[i]);
            rep(j,i,n) dp[i][j]=j-i+1;//初始化操作只穿不脱
            dp[i][i]=1;
        }

        rep(len,2,n){//枚举区间长度
            rep(i,1,n-1){//枚举左端点
                int j=i+len-1;
                if(arr[j]==arr[i]) dp[i][j]=dp[i][j-1];//假如某点和左端点相等，则可直接脱至第一件衣服
                for(int k=i;k<j && k<n;k++) dp[i][j]=min(dp[i][j],dp[i][k]+dp[k+1][j]);//枚举断点
            }
        }
        printf("Case %d: %d\n",kase,dp[1][n]);
    }
    return 0;
}

POJ2955

问题描述：给定一串包含 ‘(’ ‘)’ ‘[’ ‘]’ 的字符串求最大匹配数。
套用区间dp，但像([])这样的,用模板的状态转移方程就不行了，因为它没考虑区间包含的情况，所以在枚举断点前要判断 arr[i]和arr[j]是否匹配，匹配则dp[i][j]=dp[i-1][j-1]+1.

/*********************************************************************/
/**********************************cyf*******************************/
/*******************************************************************/
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<string>
#include<vector>
#include<set>
#include<map>
#include<queue>
#include<stack>
#include<utility>
#include<time.h>
using namespace std;
#define ll long long
#define repp(i,x) for(int i=1;i<=x;i++)
#define repx(i,n,j) for(int i=j;i<n;i++)
#define rep(i,begin,end) for(int i=begin;i<=end;i++)
#define mst(a,b) memset(a,b,sizeof(a))
#define drep(i,a,b) for(int i=a;i>=b;i--)
const int maxn=1e2+5;
const int MAXN=1e4+5;
const ll mod = 1e9+7;
const int INF=0x3f3f3f3f;
const ll LLINF=0x3f3f3f3f3f3f3f;
char s[maxn];
int dp[maxn][maxn];
int main(){
    while(~scanf("%s",s+1)){
        getchar();
        if(s[1]=='e') break;

        int n=strlen(s+1);
        rep(i,1,n){
            rep(j,1,n) dp[i][j]=0;//初始化
        }

        rep(len,2,n){
            for (int i = 1; i <= n - 1 && i + len - 1 <= n; i++ ){
                int j=i+len-1;
                if( ( s[i] == '[' && s[j] == ']' ) || (s[i] == '(' && s[j] == ')' ) ) dp[i][j]=dp[i+1][j-1]+1;//判断区间包含的情况
                for(int k = i; k<j; k++){
                    dp[i][j]=max(dp[i][j],dp[i][k]+dp[k+1][j]);//断点枚举
                }
            }
        }
        printf("%d\n",dp[1][n]*2);
    }
    return 0;
}

POJ1651

问题描述：n个数，每次选一个拿出来（不能选两端的），单次得分为这个数和他相邻数的乘积，最后只剩下2个数，总得分为单次得分之和，求最小总得分。

可知对于一个区间，我们可以枚举它最后一个取的点，则有状态转移方程dp[i][j]=min(dp[i][j],dp[i][k]+dp[k][j]+arr[i]*arr[k]*arr[j])
但是要注意，初始化时，dp[i][i+1]要化为0，其它为INF。

/*********************************************************************/
/**********************************cyf*******************************/
/*******************************************************************/
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<cmath>
#include<string>
#include<vector>
#include<set>
#include<map>
#include<queue>
#include<stack>
#include<utility>
#include<time.h>
using namespace std;
#define ll long long
#define repp(i,x) for(int i=1;i<=x;i++)
#define repx(i,n,j) for(int i=j;i<n;i++)
#define rep(i,begin,end) for(int i=begin;i<=end;i++)
#define mst(a,b) memset(a,b,sizeof(a))
#define drep(i,a,b) for(int i=a;i>=b;i--)

int gcd(int a,int b){
    int t;
    if(a<b){
        t=a;
        a=b;
        b=t;
    }
    while(a%b)
    {
        t=a%b;
        a=b;
        b=t;
    }
    
    return b;
}


const int maxn=1e2+5;
const int MAXN=1e4+5;
const ll mod = 1e9+7;
const int INF=0x3f3f3f3f;
const ll LLINF=0x3f3f3f3f3f3f3f;
int arr[maxn];
int dp[maxn][maxn];
int main(){
    int n;
    while(~scanf("%d",&n)){
        rep(i,1,n){
            scanf("%d",&arr[i]);
            rep(j,1,n) dp[i][j]=INF;
        }
        for(int i=1;i+1<=n;i++) dp[i][i+1]=0;
        rep(len,3,n){
            for(int i=1;i<n && i+len-1<=n;i++){
                int j=i+len-1;
                for(int k=i+1;k<j;k++){
                    dp[i][j]=min(dp[i][j],dp[i][k]+dp[k][j]+arr[i]*arr[k]*arr[j]);
                } 
            }
        }
        printf("%d\n",dp[1][n]);
    }
    return 0;
}