51nod 1022 石子归并 V2

题解:

这是一道石子合并的加强版DP。对于环，我们可以延长到2*n就好了。但是传统的石子合并DP是O(n^3)的，对于养的数据强度显然会超时。
这里就要引入一个四边形不等式的概念。我们定义s[i][j] 表示将i到j合并中间的最小分割点是哪个（即为k）。
我们先写一下DP转移方程：f[i][j]=f[i][k]+f[k+1][j]+sum[i][j]。
我们发现这个DP中的sum[i][j] 是满足:
sum(i,j)+sum(i1,j1) <= sum(i1,j)+sum(i,j1)( i< i1< j< j1)
那么根据四边形不等式的定理，就可以推出转移方程f也是满足这个性质的。
所以s就满足s[i][j-1]< s[i][j]< s[i+1][j]
那么我们枚举k就能变得很快。可以证明，这样做之后的DP复杂度是O(n^2)。

代码如下：

#include<cstdio>
#include<string>
using namespace std;
const int maxn=1005;
int n,ans,a[2*maxn],sum[2*maxn],f[2*maxn][2*maxn],s[2*maxn][2*maxn];
inline int read(){
    int x=0; char ch=getchar();
    while (ch<'0'||ch>'9') ch=getchar();
    while (ch>='0'&&ch<='9') x=x*10+ch-48,ch=getchar();
    return x;
}
int main(){
    n=read();
    for (int i=n+1;i<=2*n;i++) a[i]=read(),a[i-n]=a[i];
    for (int i=1;i<=2*n;i++) sum[i]=sum[i-1]+a[i];
    for (int i=0;i<=2*n;i++) f[i][i]=0,s[i][i]=i;
    for (int l=2;l<=n;l++)
    for (int i=1;i<=2*n-l+1;i++){
        int j=i+l-1; f[i][j]=1<<30;
        for (int k=s[i][j-1];k<=s[i+1][j];k++)
        if (f[i][k]+f[k+1][j]+sum[j]-sum[i-1]<f[i][j])
        f[i][j]=f[i][k]+f[k+1][j]+sum[j]-sum[i-1],s[i][j]=k;
    }
    ans=1<<30;
    for (int i=1;i<=n;i++) ans=min(ans,f[i][i+n-1]);
    printf("%d\n",ans);
    return 0;
}