看了下感受區間dp就是一種套路，直接上的板子代碼就行了。html
基礎題ac代碼：石子歸併ios
#include<bits/stdc++.h> using namespace std; typedef long long ll; typedef unsigned long long ull; int dir[8][2]={{1,0},{0,1},{1,1},{1,-1},{-1,1},{-1,-1},{0,-1},{-1,0}}; #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+6; inline int read() { char c=getchar(); int x=0, f=1; while(c<'0'|c>'9') {if(c=='-') f=-1;c=getchar();} while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar(); return x*f; } ll exgcd(ll a,ll b){ if(b==0) return a; exgcd(b,a%b); } ll q_pow(ll a,ll b,ll mod){ ll anss=1; while(b){ if(b&1) anss=anss*a%mod; a=a*a%mod; b>>=1; } return anss; } ll q_mul(ll a,ll b,ll mod){ ll anss=0; while(b){ if(b&1) anss=(anss+a)%mod; a=(a+a)%mod; b>>=1; } return anss; } int dp[105][105]; int sum[105]; int stone[105]; int main(int argc, char * argv[]){ ios::sync_with_stdio(false); int n; cin>>n; me0(sum); meinf(dp); for(int i=1;i<=n;i++){ cin>>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; }
可是這樣一眼就看出來了複雜度是n3的複雜度，這個複雜度數據稍稍大點就爆了，因此仍是要用到四邊形不等式優化。c++
可是因爲我的感受很複雜，看了不是很懂，直接貼個連接：四邊形不等式優化。優化
優化過的AC的代碼：url
#include<bits/stdc++.h> using namespace std; typedef long long ll; typedef unsigned long long ull; int dir[8][2]={{1,0},{0,1},{1,1},{1,-1},{-1,1},{-1,-1},{0,-1},{-1,0}}; #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+6; inline int read() { char c=getchar(); int x=0, f=1; while(c<'0'|c>'9') {if(c=='-') f=-1;c=getchar();} while(c>='0'&&c<='9') x=x*10+c-'0',c=getchar(); return x*f; } ll exgcd(ll a,ll b){ if(b==0) return a; exgcd(b,a%b); } ll q_pow(ll a,ll b,ll mod){ ll anss=1; while(b){ if(b&1) anss=anss*a%mod; a=a*a%mod; b>>=1; } return anss; } ll q_mul(ll a,ll b,ll mod){ ll anss=0; while(b){ if(b&1) anss=(anss+a)%mod; a=(a+a)%mod; b>>=1; } return anss; } int dp[105][105]; int sum[105]; int stone[105]; int main(int argc, char * argv[]){ ios::sync_with_stdio(false); int n; cin>>n; me0(sum); meinf(dp); int s[111][111]; for(int i=1;i<=n;i++){ cin>>stone[i]; sum[i]=sum[i-1]+stone[i]; dp[i][i]=0; s[i][i]=i; } 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 k=s[j][ends-1];k<=s[j+1][ends];k++){ if(dp[j][ends]>dp[j][k]+dp[k+1][ends]+sum[ends]-sum[j-1]){ dp[j][ends]=dp[j][k]+dp[k+1][ends]+sum[ends]-sum[j-1]; s[j][ends]=k; } } } } cout<<dp[1][n]<<endl; return 0; }