DP问题汇总，快乐背包背一背，小学算数数一数

参考自アリ本(プログラミングコンテスト チャレンジブック)

---

01背包

$n$个重$w_i$值$v_i$的物品，选总重$\le W$的组合的最大价值。
（略）

---

最长公共子序列(LCS)

给两个字符串（s，t），求最长公共子序列。

$O(nm)$

int n,m;
char s[MAX_N],t[MAX_M]; // inputs

int dp[MAX_N + 1][MAX_M + 1]; // DP
void solve(){
	for(int i=0;i<n;i++){
		for(int j=0;j<0;j++){
			if(s[i]==y[j])
				dp[i+1][j+1]=dp[i][j]+1;
			else
				dp[i+1][j+1]=max(dp[i][j+1],dp[i+1][j]);
		}
	}
	printf("%d\n",dp[n][m]);
}

最长上升子序列(LIS)

长$n$数列$a_n$。求最长上升子序列。

$O(n\log n)$

$dp[i]:=$所有长度为$i+1$的上升子列之中末尾元素最小值。
注意：最终算出来的$dp$数组并没有存这个LIS是什么，只能用来算LIS的长度。

int dp[MAXN];

void solve(){
	fill(dp,dp+n,INF);
	for(int i=0;i<n;i++){
		*lower_bound(dp,dp+n,a[i])=a[i];
	}
	printf("%d\n",lower_bound(dp,dp+n,INF)-dp);
}

---

多重部分和

$n$个不同数字$a_i$，每种$m_i$个，能否选出正好和为$K$的组合。

$O(nK)$

$dp[i+1][j]:=$区间$[0,i]$加到$K$的$a_i$的剩余，无法得到$K$的话设为$-1$

// inputs
int dp[MAXK + 1];

void solve(){
	memset(dp,-1,sizeof(dp));
	dp[0]=0;
	for(int i=0;i<n;i++){
		for(int j=0;j<=k;j++){
			if(dp[j]>=0)
				dp[j]=m[i];
			else if(j<a[i] || dp[j-a[i]]<=0)
				dp[j]=-1;
			else
				dp[j]=dp[j-a[i]]-1;
		}
	}
	if(dp[k] >= 0) printf("Yes\n");
	else printf("No\n");
}

---

多重集组合数

$n$个不同数字，每种$a_i$个，从中取$m$个的组合数。

$O(nm)$

$dp[i+1][j]:=$从前$i$钟物品中取$j$个的组合数
$dp[i+1][j]=dp[i+1][j+1]+dp[i][j]-dp[i][j-1-a_i]$

int n,m;
int a[MAXN];

int dp[MAXN+1][MAXN+1];

void solve(){
	for(int i=0;i<=n;i++) dp[i][0]=1;
	for(int i=0;i<n;i++) 
	for(int j=1;j<=m;j++)
	if(j-1-a[i]>=0) dp[i+1][j]=dp[i+1][j-1]+dp[i][j]-dp[i][j-1-a[i]]; //Add MOD! otherwise it will <0
	else dp[i+1][j]=dp[i+1][j-1]+dp[i][j];
}

---

划分数&狭义划分数

把数m划分成不超过n个数的和，求划分组合。

定义
划分数：$H_n^m$
狭义划分数：$h_n^m$
定义狭义划分数为数m划分成正好n个数的和的组合数

这里我们有
$$H_n^m=\sum _{k=1}^n{h}_k^m$$
$$h_1^m=H_1^m$$
$$h_n^m=H_n^m-H_{n-1}^m,(n\ge 2)$$

且
$$H_n^m=H_n^{m-n}+H_{n-1}^m(m-n\ge 0,n\ge 1)$$
$$H_n^m=H_{n-1}^m(m-n<0,n\ge 1)$$
$$H_n^0=1,H_0^m=0(n\ge 0,m\ge 1)$$

同时
$$h_n^m=H_n^{m-n}$$

$$H_n^m=\sum _{k=1}^m{h}_k^m=H_n^m+\sum _{k=1}^{m-1}{h}_k^m=H_n^{m-n}+H_{n-1}^m$$

$O(mn)$

// inputs

int n,m;

int dp[MAXM+1][MAXN+1];

void solve(){
	dp[0][0]=1;
	for(int i=1;i<=m;i++){
		for(int j=0;j<=n;j++){
			if(j-i>=0)
				dp[i][j]=dp[i-1][j]+dp[i][j-i];
			else
				dp[i][j]=dp[i-1][j];
		}
	}
	printf("%d\n",dp[m][n]);
	
}