区间dp及优化

看了下感觉区间dp就是一种套路，直接上的板子代码就好了。

基础题ac代码：石子归并

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
int dir[][]={{,},{,},{,},{,-},{-,},{-,-},{,-},{-,}};
#define pi acos(-1)
#define ls rt<<1
#define rs rt<<1|1
#define me0(s) memset(s,0,sizeof(s))
#define me1(s) memset(s,1,sizeof(s))
#define mef(s) memset(s,-1,sizeof(s))
#define meinf(s) memset(s,inf,sizeof(s))
#define inf 0x3f3f3f
const int N=1e6+;
inline int read() {
    char c=getchar(); int x=, f=;
    while(c<''|c>'') {if(c=='-') f=-;c=getchar();}
    while(c>=''&&c<='') x=x*+c-'',c=getchar();
    return x*f;
}
ll exgcd(ll a,ll b){
    if(b==) return a;
    exgcd(b,a%b);
}
ll q_pow(ll a,ll b,ll mod){
    ll anss=;
    while(b){
        if(b&) anss=anss*a%mod;
        a=a*a%mod;
        b>>=;
    }
    return anss;
}
ll q_mul(ll a,ll b,ll mod){
    ll anss=;
    while(b){
        if(b&) anss=(anss+a)%mod;
        a=(a+a)%mod;
        b>>=;
    }
    return anss;
}
int dp[][];
int sum[];
int stone[];
int main(int argc, char * argv[]){
    ios::sync_with_stdio(false);
    int n;
    cin>>n;
    me0(sum);
    meinf(dp);
    for(int i=;i<=n;i++){
        cin>>stone[i];
        sum[i]=sum[i-]+stone[i];
        dp[i][i]=;
    }
    for(int len=;len<=n;len++){//枚举长度
        for(int j=;j+len<=n+;j++){//枚举起点，ends<=n
            int ends=j+len-;
            for(int i=j;i<ends;i++){//枚举分割点，更新小区间最优解
                dp[j][ends]=min(dp[j][ends],dp[j][i]+dp[i+][ends]+sum[ends]-sum[j-]);
            }
        }
    }
    cout<<dp[][n]<<endl;
    return ;
}

但是这样一眼就看出来了复杂度是n³的复杂度，这个复杂度数据稍稍大点就爆了，所以还是要用到四边形不等式优化。

但是由于个人感觉很复杂，看了不是很懂，直接贴个链接：四边形不等式优化。

优化过的AC的代码：

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
typedef unsigned long long ull;
int dir[][]={{,},{,},{,},{,-},{-,},{-,-},{,-},{-,}};
#define pi acos(-1)
#define ls rt<<1
#define rs rt<<1|1
#define me0(s) memset(s,0,sizeof(s))
#define me1(s) memset(s,1,sizeof(s))
#define mef(s) memset(s,-1,sizeof(s))
#define meinf(s) memset(s,inf,sizeof(s))
#define inf 0x3f3f3f
const int N=1e6+;
inline int read() {
    char c=getchar(); int x=, f=;
    while(c<''|c>'') {if(c=='-') f=-;c=getchar();}
    while(c>=''&&c<='') x=x*+c-'',c=getchar();
    return x*f;
}
ll exgcd(ll a,ll b){
    if(b==) return a;
    exgcd(b,a%b);
}
ll q_pow(ll a,ll b,ll mod){
    ll anss=;
    while(b){
        if(b&) anss=anss*a%mod;
        a=a*a%mod;
        b>>=;
    }
    return anss;
}
ll q_mul(ll a,ll b,ll mod){
    ll anss=;
    while(b){
        if(b&) anss=(anss+a)%mod;
        a=(a+a)%mod;
        b>>=;
    }
    return anss;
}
int dp[][];
int sum[];
int stone[];
int main(int argc, char * argv[]){
    ios::sync_with_stdio(false);
    int n;
    cin>>n;
    me0(sum);
    meinf(dp);
       int s[][];
    for(int i=;i<=n;i++){
        cin>>stone[i];
        sum[i]=sum[i-]+stone[i];
        dp[i][i]=;
        s[i][i]=i;
    }
    for(int len=;len<=n;len++){//枚举长度
        for(int j=;j+len<=n+;j++){//枚举起点，ends<=n
            int ends=j+len-;
            for(int k=s[j][ends-];k<=s[j+][ends];k++){
                if(dp[j][ends]>dp[j][k]+dp[k+][ends]+sum[ends]-sum[j-]){
                    dp[j][ends]=dp[j][k]+dp[k+][ends]+sum[ends]-sum[j-];
                    s[j][ends]=k;
                }
            }
        }
    }
    cout<<dp[][n]<<endl;
    return ;
}