NOIP前 基础动态规划模板

• 背包问题
• • 0/1背包，完全背包，多重背包，分组背包，二维背包，混合背包
• 最长上升子序列（LIS）
• 最长公共子序列（LCS）
• 最长公共上升子序列（？）

背包

0/1背包
（每个物品只能选择一个）

for(i=1;i<=n;++i){
	w[i]=in;c[i]=in;
}
for(i=1;i<=n;++i){
	for(j=m;j>=w[i];--j)//0/1背包要倒序！保证不会将同一个物品多次放入
		f[j]=max(f[j-w[i]]+c[i],f[j]);
}
for(i=1;i<=m;++i) ans=max(ans,f[i]);

完全背包
（每个物品有无穷多个）

for(i=1;i<=n;++i)
	for(j=w[i];j<=m;++j)//正着跑
		f[j]=max(f[j],f[j-w[i]]+c[i]);

多重背包
二进制拆分版
a[i]表示原来一个物品的价值，b[i]表示原来一个物品的体积
c[i]表示这个物品的个数

	for(i=1;i<=n;++i){
		for(j=1;j<=c[i];j<<=1){
			val[++cnt]=a[i]*j;
			w[cnt]=b[i]*j;
			c[i]-=j;
		}
		if(c[i]) val[++cnt]=a[i]*c[i],w[cnt]=b[i]*c[i];
	}

分组背包

for(i=1;i<=t;++i){//枚举每一个组 
	for(j=V;j>=0;--j)//这三个的顺序不能反
		for(k=1;k<=num[i];++k)//枚举这一组里的物品，任选一个
			if(j>=w[i][k])
				f[j]=max(f[j],f[j-w[i][k]]+v[i][k]);
}

混合背包
对于不同的物品，有些只能选一个，有些有选择上线，有些无穷多
分类讨论即可

	for(i=1;i<=n;++i){
		if(!p[i])
			for(j=w[i];j<=m;++j)
				f[j]=max(f[j],f[j-w[i]]+c[i]);
		else
			for(j=1;j<=p[i];++j)
				for(int k=m;k>=w[i];--k)
					f[k]=max(f[k],f[k-w[i]]+c[i]);
	}	

二维背包
（例题：HDU2159）
其实就是把我们之前的一维再增加一维即可
然后其他的就按照0/1背包，多重背包，完全背包的套路用就是了

最长上升子序列

$n\ast log(n)$ 二分
最后最长上升子序列的长度就在k里面
如果要输出路径的话（我只会在n²算法下），就记录一下是从哪个点转移过来的，最后递归输出即可

for(i=1;i<=n;++i){
	if(a[i]>d[k]) d[++k]=a[i];
	else d[lower_bound(d+1,d+k+1,a[i])-d]=a[i];
}

最长公共子序列

放上完整代码：求出$a,b$的最长公共子序列并输出路径

#include<bits/stdc++.h>
using namespace std;
char a[205],b[205];
int f[205][205],g[205][205];
inline void out(int i,int j){
	if(i<0||j<0) return;
	if(g[i][j]==1) out(i-1,j-1),printf("%c",b[j]);
	else if(g[i][j]==0) out(i-1,j);
	else out(i,j-1);
}
int main(){
	scanf("%s%s",a+1,b+1);
	int n=strlen(a+1),m=strlen(b+1);
	for(int i=1;i<=n;++i)
		for(int j=1;j<=m;++j){
			if(a[i]==b[j]) {
				if(f[i-1][j-1]+1>f[i][j]){
					f[i][j]=f[i-1][j-1]+1;
					g[i][j]=1;
				}
			}
			else{
				if(f[i-1][j]>f[i][j]) {
					f[i][j]=f[i-1][j];
					g[i][j]=0;
				}
				if(f[i][j-1]>f[i][j]){
					f[i][j]=f[i][j-1];
					g[i][j]=-1;
				}
			}
		}
	int p,q,ans=0;
	for(int i=1;i<=n;++i)
		for(int j=1;j<=m;++j)
			if(f[i][j]>ans){
				ans=f[i][j];
				p=i;q=j;
			}
	if(!ans)	printf("0");
	else{
		printf("%d\n",ans);
		out(p,q);
	}
	return 0;
}

最长公共上升子序列

	for(i=1;i<=n;++i){
		res=0;
		for(j=1;j<=n;++j){
			if(b[j]<a[i]) res=max(res,f[i-1][j]);
			if(a[i]==b[j])	f[i][j]=res+1;
			else f[i][j]=f[i-1][j];
		}
	}