标签：gravity 程序设计 ati net return com nbsp 递归 term
1.背包问题
（1）01背包
从n个重量和价值分别为wi,vi的物品，从中选出不超过W的物品，每种物品仅有一件，求所有方案中V的最大值。
最朴素最简单也最费时的方法：O(2^n) int rec(int i,int j)//从第i个开始挑选总重小于j的部分
递归 递归终止条件：i==n return 0;
递归分支：① j<wi res=rec(i+1,j); //无法挑选，看下一个
② res=max(rec(i+1,j),rec(i+1,j-w[i])+v[i]))//挑选，不挑选，选其中大的
分析：递归搜索深度→n，每层两次分支（选与不选），有重复计算
优化：记录每次递归的结果（记忆化搜索）
Int dp[MAX_N][MAX_N];
递归 int rec(int i,int j)
初始条件：memset(dp,-1,sizeof(dp));
终止条件：①dp[i][j]>=0 returndp[i][j]//已计算过
②i==n return 0;
递归分支：① j<wi res=rec(i+1,j); //无法挑选，看下一个
② res=max(rec(i+1,j),rec(i+1,j-w[i])+v[i]))//挑选，不挑选，选其中大的
分析及remark：复杂度O(nW)
数组初始化：①memset(type *arrary,figure,sizeof(arrary))
只能填充（0，+-1，0x3f3f3f3f），其他值不可以
memset按照1字节为单位对内存填充，-1的二进制每一位均为1
②fill(type* arrary,type *arrary+n,figure) 可赋值任意值
递归
↓↓↓↓↓↓ for(int i=n-1;i>=0;i--) //
递推(双重循环) for(int j=0;j<=w;j++)
{ if(j<w[i]) dp[i][j]=dp[i+1][j]; //不选
else dp[i][j]=max(dp[i+1][j],dp[i+1][j-w[i]]+v[i]);}
其他3种递推写法：（来源：《挑战程序设计竞赛》）
（2）完全背包
递推关系1：（3重循环，有重复计算）复杂度O(nW^2)
For(int i=0;i<n;i++)
For(int j=0;j<=w;j++)
For(int k=0;k*w[i]<=j;k++)
Dp[i+1][j]=max(dp[i+1][j],dp[i+1][j-k*w[i]]+k*v[i])
优化：(左上→右下 变为 左→右)
Dp数组初始化为0
For(int i=0;i<n;i++)
For(int j=0;j<=w;j++)
{
If(j<w[i]) dp[i+1][j]=dp[i][j];
Else dp[i+1][j]=max(dp[i][j],dp[i+1][j-w[i]]+v[i])
} 只有这里与01背包不同，前j个已更新过，可直接用
进一步优化：用一个数组实现，只需要记录当前最优状态
比较01背包与完全背包（循环方向不同）
2.LCS(Longest common subsequence)
dp[n][m]即为所求
for(int i=0;i<n;i++)
for(int j=0;j<m;j++)
{
if(s[i]==t[i])
dp[i+1][j+1]=dp[i][j]+1;
else
dp[i+1][j+1]=max(dp[i+1][j],dp[i][j+1])
}
3.LIS(Longest Increasing subsequence)
dp[i]:以ai为结尾的最长上升子序列长度
dp[i]=max{1,dp[i]+1|j<i且aj<ai}
O(n^2)
int res;
for(int i=0;i<n;i++)
{
dp[i]=1;
for(int j=0;j<i;j++)
if(a[j]<a[i]) //每存在aj<ai&&j<i,dp[i]更新一次
dp[i]=max{dp[i],dp[j]+1};
}
res=max{dp[i]|0<=i<n}
Remark:
其他方法：
可以用lower_bound();
dp[i]:长度为i+1的上升子列中末尾元素的最小值（不存在的话为inf）
dp[max_n]初始化为inf，按顺序逐个考虑数列的元素，对于每个ai，如果i=0||dp[i-1]<ai,就用dp[i]=min(dp[i],ai)更新，最终找出使得dp[i]<inf的最大的i+1即为结果。DP直接实现，可以在O（n^2）的时间内给出结果，但可以进一步优化，dp数组中除inf之外是单调递增的，对于每个ai最多有一次更新，更新的位置可用二分的方法优化，时间复杂度可以降低到O（nlogn）
int dp[max_n]
void solve（）
{
fill(dp,dp+n,inf);
for(int i=0;i<n;i++)
*lower_bound(dp,dp+n,a[i])=a[i];
res=lower_bound(dp,dp+n,inf)-dp;
}
// lower_bound()可以从已排好序的a中利用二分搜索找出满足ai>=k的ai的最小的指针，类似的还有upper_bound,找出的为ai>k的最小指针
//求n个有序数组a中k的个数，可以用：upper_bound(a,a+n,k)-lower_bound(a,a+n,k);
标签：gravity 程序设计 ati net return com nbsp 递归 term
原文地址：http://www.cnblogs.com/Egoist-/p/7391224.html
