JZOJ 1293. 气象牛

Posted zjzjzj

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了JZOJ 1293. 气象牛相关的知识,希望对你有一定的参考价值。

题目

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.
 

Input

  第一行: 两个空格分离的数: N 和 E
  第2..N+1行: 第i+1行包含一次测量记录:M_i

Output

  第一行: 两个空格分开的数: 最少能达到误差小于等于E的测量数目和使用那个测量数目能达到的最小误差.
 

Sample Input

4 20
10
3
20
40

Sample Output

2 17
 

Data Constraint

 
 

Hint

【样例说明】
  选择第二和第四次测量结果能达到最小误差17. 第一次结果的误差是2*|10-3| = 14;第三次结果的误差是|2*20 - (3+40)|=3.

 

分析

 

  • 我们可以用DP
  • 设f[i][j]为前i个用了j个,所产生的最小误差
  • 直接上n^3
  • 预处理出每个点在另一个点前 后 中的值

 

 

代码

 

 

 1 #include <cmath>
 2 #include <cstdio>
 3 #include <cstring>
 4 #include <iostream>
 5 #include <algorithm>
 6 #define ll long long
 7 using namespace std;
 8 int a[201];
 9 int dis[201][201];
10 int f[201][201];
11 int ans[201];
12 int main()
13 {
14     int n,e;
15     scanf("%d%d",&n,&e);
16     for (int i=1;i<=n;i++) scanf("%d",&a[i]);
17     for (int i=0;i<=n+1;i++)
18      for (int j=i+1;j<=n+1;j++)
19       for (int k=i+1;k<=j-1;k++)
20       if (i==0) dis[i+1][j-1]+=abs(2*a[k]-2*a[j]); else 
21       if (j==n+1) dis[i+1][j-1]+=abs(2*a[k]-2*a[i]); else 
22       dis[i+1][j-1]+=abs(2*a[k]-a[i]-a[j]); 
23     memset(f,0x3f,sizeof(f));
24     f[0][0]=0;
25     for (int i=0;i<=n;i++)
26       for (int j=0;j<=n;j++)
27       {
28            if (f[i][j]<=e)
29            {
30                for (int k=i+1;k<=n+1;k++)
31               f[k][j+1]=min(f[k][j+1],f[i][j]+dis[i+1][k-1]);
32          }
33       }
34     for (int i=2;i<=n+1;i++)
35     {
36         if (f[n+1][i]<=e) 
37         {
38             cout<<i-1<<" "<<f[n+1][i]; return 0;
39         }
40         
41     }
42     
43 }

 

以上是关于JZOJ 1293. 气象牛的主要内容,如果未能解决你的问题,请参考以下文章

[usaco jan 09] 气象牛 baric [dp]

[BZOJ1575] [Usaco2009 Jan]气象牛Baric(DP)

动态规划bzoj1575: [Usaco2009 Jan]气象牛Baric

Floyd | | jzoj[1218] | | 我也不知道该怎么写

JZOJ 1252. 天平

JZOJ 1292. 公牛和母牛