【无聊放个模板系列】HDU 3506 （四边形不等式优化DP-经典石子合并问题[环形]）

#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
#include<queue>
#include<cmath>
using namespace std;
#define Maxn 1010
#define INF 0xfffffff

int a[2*Maxn],sum[2*Maxn];
int f[2*Maxn][2*Maxn],s[2*Maxn][2*Maxn];

int mymin(int x,int y) {return x<y?x:y;}

int main()
{
    int n;
    while(scanf("%d",&n)!=EOF)
    {
        for(int i=1;i<=n;i++) scanf("%d",&a[i]);
        for(int i=1;i<=n;i++) a[i+n]=a[i];
        sum[0]=0;
        for(int i=1;i<=2*n;i++) sum[i]=sum[i-1]+a[i];
        memset(f,63,sizeof(f));
        for(int i=1;i<=2*n;i++) f[i][i]=0,s[i][i]=i;
        for(int i=2*n;i>=1;i--)
         for(int j=i+1;j<=2*n;j++)
         {
             if(j-i+1>n) break;
             for(int k=s[i][j-1];k<=s[i+1][j];k++)
             {
                 if(f[i][j]>f[i][k]+f[k+1][j]+sum[j]-sum[i-1])
                 {
                    f[i][j]=f[i][k]+f[k+1][j]+sum[j]-sum[i-1];
                    s[i][j]=k;
                 }
             }
         }
        int ans=INF;
        for(int i=1;i<=n;i++) ans=mymin(ans,f[i][i+n-1]);
        printf("%d
",ans);
    }
    return 0;
}

其实就是个区间DP的优化，环形石子，把它拆成链就好了【在右面copy一倍】

四边形不等式优化DP

2016-11-18 08:51:01