BZOJ 1575 [Usaco2009 Jan]气象牛Baric

Description
为了研究农场的气候,Betsy帮助农夫John做了N(1 <= N <= 100)次气压测量并按顺序记录了结果M_1…M_N(1 <= M_i <= 1,000,000). Betsy想找出一部分测量结果来总结整天的气压分布. 她想用K(1 <= K <= N)个数s_j (1 <= s_1 < s_2 < … < s_K <= N)来概括所有测量结果. 她想限制如下的误差: 对于任何测量结果子集,每一个非此子集中的结果都会产生误差.总误差是所有测量结果的误差之和.更明确第说, 对于每一个和所有s_j都不同的i: * 如果 i 小于 s_1, 误差是: 2 * | M_i - M_(s_1) | * 如果i在s_j和s_(j+1)之间,误差是: | 2 * M_i - Sum(s_j, s_(j+1)) | 注:Sum(x, y) = M_x + M_y; (M_x 和 M_y 之和) * 如果i大于s_K,误差为: 2 * | M_i - M_(s_K) | Besty给了最大允许的误差E (1 <= E <= 1,000,000),找出最小的一部分结果史得误差最多为E.

---

【题目分析】
动态规划。

---

【代码】

#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
int n,e,f[101][101];
int ans[101][101],lst[101],pre[101],a[101],minn=0x3f3f3f3f;
int main()
{
    scanf("%d%d",&n,&e);
    for (int i=1;i<=n;++i) scanf("%d",&a[i]);
    for (int i=1;i<=n;++i)
        for (int j=1;j<i;++j)
            pre[i]+=2*abs(a[j]-a[i]);
    for (int i=n;i;--i)
        for (int j=n;j>i;--j)
            lst[i]+=2*abs(a[i]-a[j]);
    for (int i=1;i<n;++i)
        for (int j=i+1;j<=n;++j)
            for (int k=i+1;k<j;++k)
                ans[i][j]+=abs(2*a[k]-(a[i]+a[j]));
    memset(f,0x3f,sizeof f);
    for (int i=1;i<=n;++i) f[i][1]=pre[i];
    for (int i=1;i<=n;++i)
        for (int j=1;j<=n;++j)
            for (int k=1;k<i;++k)
            {
                f[i][j]=min(f[i][j],f[k][j-1]+ans[k][i]);
                if (f[i][j]+lst[i]<=e) minn=min(minn,j);
            }
    int tmp=0x3f3f3f3f;
    for (int i=1;i<=n;++i) tmp=min(tmp,f[i][minn]+lst[i]);
    printf("%d %d\n",minn,tmp);
}