动态规划入门_hdu_1421

参考文件：http://pan.baidu.com/s/1ntv7IPF

动态规划是将某个问题分成若干块，每个阶段都在前一阶段的基础上推出最优解，进而得出整体的最优解。

动态规划是典型的空间换时间，通过二维或一维数组记录结果，使得问题不被重复计算。不扯了，直接看题：

链接：http://acm.hdu.edu.cn/showproblem.php?pid=1421

搬寝室

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 16893    Accepted Submission(s): 5727

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
  
  

 

代码：

参考了别人的思路：http://blog.csdn.net/weyuli/article/details/8722940

给定n个物品，每个物品有重量，
   从中选出m对，使得这m对物品重量差的平方和最小。
    疲劳度：m对物品重量差的平方和
   分析与解题思路
   先对n中物品的重量排序
   令dp[i][j]表示前i个物品中选j对的最小疲劳度。
   则dp[i][j]可能含有第i个物品(这种情况下,第i种物品一定是和第i-1个物品配对)，
   则dp[i][j]=dp[i-2][j-1]+w[i-1]
    dp[i][j]的j对也可能不含有第i个物品，此时有
   dp[i][j]=dp[i-1][j]
   状态转移方程
   dp[i][j]=min{dp[i-2][j-1]+w[i-1],dp[i-1][j]}

AC代码：

#include<stdio.h>
#include<stdlib.h>
int cmp(const void*a,const void*b)
{
    return *(int *)a-*(int *)b;
}
int Min(int a,int b)
{
    return a<b?a:b;
}
long long map[2020][2020];
int val[2020];
int main()
{
    int n,k,i,j;
    while(~scanf("%d%d",&n,&k))
    {
        for(i=1;i<=n;i++)scanf("%d",val+i);
        qsort(val+1,n,sizeof(val[0]),cmp);
        for(i=1;i<n;i++)val[i]=(val[i+1]-val[i])*(val[i+1]-val[i]);
        for(i=0;i<=n;i++)
            for(j=0;j<=k;j++)
            map[i][j]=1000000000;
        for(i=0;i<=n;i++)map[i][0]=0;
        for(i=2;i<=n;i++)
            for(j=1;j*2<=i;j++)
            map[i][j]=Min(map[i-2][j-1]+val[i-1],map[i-1][j]);
        printf("%I64d\n",map[n][k]);
    }
    return 0;
}

提交的时候WA了6次，一直找不到哪错了，比着上面链接上的代码改了好久，

最后把map[i][j]=100000;改为map[i][j]=1000000000;通过了

（开始为int后来改为long long ）

分享下队友的模板：

#include<stdio.h>
#include<string.h>
#define max(a,b)(a>b ? a:b)
int v[100000];
int w[100000];
int num[100000];
int bag[100000];
int Weight,n;
int Max(int a,int b)
{
    return a>b?a:b;
}
int zeroOnePack()//01背包
{
    int i,j;
    for(i=1; i<=n; i++)
    {
        for(j=Weight; j>=w[i]; j--)
        {
            bag[j] = max(bag[j], v[i] + bag[j-w[i]]);
        }
    }
    return bag[Weight];
}
int completePack()//完全背包
{
    int i,j;
    for(i=1; i<=n; i++)
    {
        for(j=w[i]; j<=Weight; j++)
        {
            bag[j] = max(bag[j], v[i] + bag[j-w[i]]);
        }
    }
    return bag[Weight];
}
int multiPack()//多重背包
{
    int i,j,k;
    memset(bag,0,sizeof(bag));
    for(i=1; i<=n; i++)
    {
        for(j=w[i]; j<=Weight; j++)
        {
            if(num[i]*w[i]>Weight)
                completePack(n);
            else
            {
                k = 1;
                while(k<num[i])
                {
                    zeroOnePack(n);
                    num[i] -= k;
                    k += k;
                }
                zeroOnePack(n);
            }
        }
    }
    return bag[Weight];
}
int main()
{
    int i;
    while(scanf("%d%d",&n,&Weight),n)
    {
       // for(i=1;i<=n;i++)
       //     scanf("%d",num+i);
      // for(i=1;i<=n;i++)
        //    scanf("%d",w+i);
        for(i=1;i<=n;i++)
            scanf("%d %d",w+i,v+i);
       // int ans = multiPack(n);
        printf("%d\n",zeroOnePack());
    }
    return 0;
}