DP

一、 HDU 1081
  求最大和的子矩阵  100*100矩阵
解题思路：这道题是个DP问题，如果用穷举的方法的话，时间复杂度为O（N^4），会超时，所以我们就考虑用DP的方法。
 #include <stdio.h>
#include <string.h>
int main()
{
const int size = 101;
int Max,sum;
int N,i,j,k,val,a[size][size];
while (scanf("%d",&N)!=EOF)
{
memset(a,0,sizeof(a));
Max = -128;
for(i=1;i<=N;i++)
for(j=1;j<=N;j++)
{
scanf("%d",&val);
a[i][j]+=a[i][j-1]+val; //a[i][j]表示第i行前j个数之和
}
for(i=1;i<=N;i++)
for (j=i;j<=N;j++)
for(sum = 0,k=1;k<=N;k++)
{
//a[k][j]-a[k][i-1]表示第k行第j列与第i列之间的数
sum=(sum>0?sum:0)+a[k][j]-a[k][i-1];
if(sum>Max)
Max = sum;
}
printf("%d/n",Max);
}
}




无语了，自己大意阿扁，漏洞百出！！！
#include<cstdio>
#include<algorithm>
#include<cstring>
int main()
{
    int n;
    int a[101][101];
    while(scanf("%d",&n)!=EOF)
    {
        memset(a,0,sizeof(a));//还有这里一定要初始化为0，因为后面要直接求和
        int max=-128;   //果然出现了问题！！这里要注意和可能出现负值，所以，一定不要初始化为0
        for(int i=1; i<=n; i++)
            for(int j=1; j<=n; j++)  //要从1到N，也是因为后面求和 有j-1，不然越栈
            {
                int val;
                scanf("%d",&val);
                a[i][j]+=a[i][j-1]+val; //这里之前也写错了，，这样可以节省开销，自己好好理解一下
            }
        int ans=0,k;
        for(int i=1; i<=n; i++)
            for(int j=i; j<=n; j++)  //这里也要注意，从i行开始，刚开始写成了，从0开始是不对的！！！！因为他是从左往右扫描的！！！！
                for(ans=0,k=1; k<=n; k++)   //这里尤其要注意！！这里是以确定了矩阵的左右两边为基础，当再重新进入这个K的时候，ans还要初始化为0
                {
                    ans=(ans>0?ans:0)+a[k][j]-a[k][i-1];   //第K行i和j之间的和
                    if(ans>max) max=ans;
                }
        printf("%d\n",max); //这里又错了，之前ans默认为是答案，其实答案是最大值！！！
    }
}


二、 HDU1421

解题思路：一道DP题
   给定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]+(val[i]-val[i-1])*(val[i]-val[i-1])
    dp[i][j]的j对也可能不含有第i个物品，此时有
   dp[i][j]=dp[i-1][j]
   状态转移方程
   dp[i][j]=min{dp[i-2][j-1]+(val[i]-val[i-1])*(val[i]-val[i-1]),dp[i-1][j]


#include <stdio.h>
#include <stdlib.h>
#define size 2005
#define INIT 2147483646

int cmp(const void *a,const void *b)
{
return *(int *)a-*(int *)b;
}

int Min(int a,int b)
{
return a<b?a:b;
}

int dp[size][1005];

int main()
{
int n,k,i,j;
int val[size] = {0};
dp[0][0] = 0;
while (scanf("%d%d",&n,&k)!=EOF)
{
val[0] = 0;
for(i=1;i<=n;i++)
scanf("%d",&val[i]);
qsort(val+1,n,sizeof(val[0]),cmp);
for (i=0;i<=n;i++)
{
for(j=1;j<=k;j++)
dp[i][j] = INIT;
}
for (i=2;i<=n;i++)
{
for (j=1;j*2<=i;j++)
dp[i][j] = Min(dp[i-2][j-1]+(val[i]-val[i-1])*(val[i]-val[i-1]),dp[i-1][j]);
}
printf("%d/n",dp[n][k]);
}
return 0;
}