<USACO09JAN>气象测量/气象牛The Baric Bovine

我..莫名其妙搞出来的

随便搞搞..幻想中的dp结果对了

我也很迷茫

 1 #include<cstdio>
 2 #include<cstring>
 3 #include<iostream>
 4 #include<algorithm>
 5 using namespace std;
 6 int n,e;
 7 long long m[110],ero[110][110],dp[110][110];//ero选i和j产生的误差//dp以j结尾的时候,有i个最少测量数目 能达到的最小误差 
 8 long long abss(long long a)
 9 {
10     if(a<0)return -a;
11     return a;
12 }
13 int main()
14 {
15     int i,j,k;
16     scanf("%d%d",&n,&e);
17     for(i=1;i<=n;i++)scanf("%lld",&m[i]);    
18     long long ans=n,anss=999999999;
19     
20     for(k=1;k<=n;k++)
21     {
22         for(i=k+1;i<=n;i++)
23             for(j=k+1;j<=i-1;j++) /*枚举题中的i*/ ero[k][i]+=abss(2*m[j]-m[i]-m[k]);
24             
25         for(i=1;i<k;i++) ero[k][0]+=2*abss(m[i]-m[k]); //取k时,小于k的i值产生的误差 
26         for(i=k+1;i<=n;i++) ero[k][n+1]+=2*abss(m[i]-m[k]);// 取k时 大于k的i值产生的误差 
27     }//关于误差 
28     
29     for(i=1;i<=n;i++) 
30     {
31         dp[1][i]=ero[i][0]+ero[i][n+1]; 
32         if((dp[1][i]<=e) && (dp[1][i]<anss))
33         {
34             anss=dp[1][i];ans=1;
35         }
36     }//处理只取一个数的情况,也就是只取数i的情况 预处理 
37     
38     if(ans==1){printf("%lld %lld\n",ans,anss);return 0;}//如果可以只取一个数emm且满足条件啦 
39     
40     for(k=2;k<=n;k++)//最少能达到误差小于等于E的测量数目
41         for(i=k;i<=n;i++)//枚举结尾的数 
42         {
43             dp[k][i]=9999999999;
44             for(j=k-1;j<i;j++)//枚举上一轮的结尾数 
45             {
46                 long long err=ero[j][i]+ero[i][n+1]-ero[j][n+1]; //相比上一轮多的误差 
47                 dp[k][i]=min(dp[k][i],dp[k-1][j]+err);
48             }
49             if((dp[k][i]<=e) && (dp[k][i]<anss) && (k<=ans)) anss=dp[k][i],ans=k;//找最优解emm 
50         }
51     printf("%lld %lld\n",ans,anss);
52 return 0;
53 }

猜你喜欢

转载自www.cnblogs.com/pile8852/p/9294319.html