我的第一个dp程序 哎! 可惜状态方程不是自己想的 有点小小的遗憾!!
搬寝室
Time Limit : 2000/1000ms (Java/Other) Memory Limit : 65536/32768K (Java/Other)
Total Submission(s) : 107 Accepted Submission(s) : 40
Total Submission(s) : 107 Accepted Submission(s) : 40
Problem Description
搬寝室是很累的,xhd深有体会.时间追述2006年7月9号,那天xhd迫于无奈要从27号楼搬到3号楼,因为10号要封楼了.看着寝室里的n件物品,xhd开始发呆,因为n是一个小于2000的整数,实在是太多了,于是xhd决定随便搬2*k件过去就行了.但还是会很累,因为2*k也不小是一个不大于n的整数.幸运的是xhd根据多年的搬东西的经验发现每搬一次的疲劳度是和左右手的物品的重量差的平方成正比(这里补充一句,xhd每次搬两件东西,左手一件右手一件).例如xhd左手拿重量为3的物品,右手拿重量为6的物品,则他搬完这次的疲劳度为(6-3)^2 = 9.现在可怜的xhd希望知道搬完这2*k件物品后的最佳状态是怎样的(也就是最低的疲劳度),请告诉他吧.
Input
每组输入数据有两行,第一行有两个数n,k(2<=2*k<=n<2000).第二行有n个整数分别表示n件物品的重量(重量是一个小于2^15的正整数).
Output
对应每组输入数据,输出数据只有一个表示他的最少的疲劳度,每个一行.
Sample Input
2 1 1 3
Sample Output
4
状态方程
dp[k][i]=min(dp[k-1][i-2]+(a[i]-a[i-1])^2,dp[k][i-1]);
dp[k][i] 表示 k 对物品在前 i 个物品的最小值
代码:
View Code
#include<iostream> #include<stdio.h> #include<algorithm> #define N 2005 using namespace std; int dp[N/2][N]; int a[N]; void Init(int k) { int i; dp[k][2*k]=0; for(i=2;i<=2*k;i+=2) dp[k][2*k]+=(a[i]-a[i-1])*(a[i]-a[i-1]); } int work(int n,int m) { int i,k; sort(a+1,a+n+1); for(k=1;k<=m;k++) Init(k); for(k=1;k<=m;k++) { for(i=2*k+1;i<=n;i++) { if((dp[k-1][i-2]+(a[i]-a[i-1])*(a[i]-a[i-1]))<dp[k][i-1]) dp[k][i]=dp[k-1][i-2]+(a[i]-a[i-1])*(a[i]-a[i-1]); else dp[k][i]=dp[k][i-1]; } } return dp[m][n]; } int main() { int n,k; while(scanf("%d%d",&n,&k)!=EOF) { int i; for(i=1;i<=n;i++) cin>>a[i]; int ans=work(n,k); printf("%d\n",ans); } return 0; }