破译密码 动态规划 重庆一中高2018级竞赛班第十次测试 2016.9.16 Problem 3

【问题描述】
Lucky得到总部发来的消息，我军特种部队已经截获敌人的一个密码本，但是这个密码本本身是由密码写成的。为了给敌人造成沉重的打击，Lucky必须尽快破译密码。经过一天一夜的探索，Lucky发现敌人密码本实际上记载着一个数列，而最终密码由这个数列经过下列运算得到：
1、把数列从小到大排序。
2、在排好序的数列中，任选一个数，这个数将把原数列分成左右两个数列（选出的数不在新数列中，并且新数列有可能为空）。
3、每个新数列进行第2步操作，直到最后得到的数列长度都为1，即全部变成单个数。
4、将第3步得到的每个数*（得到它所需的分拆次数），并累加得到一个和。
5、重复2，3，4操作，直到遍历所有的分拆可能，这时，所得的和当中最小的一个就是日军的最终密码。
现在Lucky请求你帮助，尽快破译这段密码！

【输入格式】
第一行为N（N<=1000），表示密码本记录的数列的长度。
下一行共有N个数，即敌人密码本记载的数列。

【输出格式】
一个整数，即敌人最终密码。

【输入样例】

3
1 3 2

【输出样例】

10

【数据范围】
N<=1000，数列元素的值不超过1000000。

思路：此题与石子归并相似，运用动态规划完成。
f[i][j]:拆分i…j的结果的最小值。
sum[i]:1…i的总和（前缀和优化）。
……………………………………………………………………………………
边界：f[i][i]=a[i]; d[i][i]=i;
状态转移方程：f[i][j]=f[i][k-1]+f[k+1][j]+sum[j]-sum[i-1];
优化：使用四边形不等式减小范围。

#include<cstdio>
#include<algorithm>
#include<iostream>
#define inf 100000000000000ll
using namespace std;

int n,a[1005];
long long sum[1005],f[1005][1005],d[1005][1005];

int main()
{
    scanf("%d",&n);
    sum[0]=0;
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&a[i]);
    }
    sort(a+1,a+1+n);
    for(int i=1;i<=n;i++)
    {
        sum[i]=sum[i-1]+a[i];
    }
    for(int i=1;i<=n;i++)
    {
        f[i][i]=a[i];
        d[i][i]=i;
    }

    for(int i=n;i>=1;i--)
    {
        for(int j=i+1;j<=n;j++)
        {
            f[i][j]=inf;
            for(int k=d[i][j-1];k<=d[i+1][j];k++)
            {
                if(f[i][k-1]+f[k+1][j]+sum[j]-sum[i-1]<f[i][j])
                {
                    d[i][j]=k;
                    f[i][j]=f[i][k-1]+f[k+1][j]+sum[j]-sum[i-1];
                }
            }
        }
    }

    cout<<f[1][n]<<endl;
    return 0;
}