pku3638 Moogle (动态规划)

 题目地址：http://acm.pku.edu.cn/JudgeOnline/problem?id=3638

 题目大意：一个地图软件需要需要存储h个房子的坐标信息，但为节约存储空间可以只存储其中的c个房子的坐标信息，其他房子的坐标信息可以由插值得到。问怎样选择这个c个房子，使得所有房子的插值误差之和最小， 输出最小的插值误差均值。

 题目分析：动态规划，

状态定义为： dp[i][j]表示第i号房子被定为第j个房子需要存储的房子，此时总的最小误差。

 

#include <cstdio> #include <cmath> double dp[201][201]; double sum[201][201]; int house[201]; #define min(a,b) (((a) < (b)) ? (a) : (b)) int main() { int t, h, c, i, j, k, cindex, mid, p; double err, tmp; //freopen("data.txt", "r", stdin); scanf("%d", &t); while (t--) { scanf("%d%d", &h, &c); for (i=1; i<=h; i++) scanf("%d", &house[i]); for (i=1; i<=h; i++) for (j=i+2; j<=h; j++) { err = 0; for (k=i+1; k<j; k++) err += abs((double)house[i] + (house[j]-house[i])*(k-i)/(double)(j-i) - house[k]); sum[i][j] = err; } tmp = 1000000000; for (i=1; i<=h; i++) for (j=1; j<=c; j++) dp[i][j] = tmp; dp[1][1] = 0; for (i=2; i<=h; i++) for (j=2; j<=c; j++) { if (j>i) continue; dp[i][j] = min(dp[i][j], dp[i-1][j-1]); for (p=i-2; p>=j-1; p--) //p...i的之间的误差计算 if (dp[p][j-1]<tmp) //p...i dp[i][j] = min(dp[i][j], dp[p][j-1]+sum[p][i]); } printf("%.4f/n", dp[h][c]/h); } return 0; }