动态规划

递推：

Tri Tiling

• 定义$f[i][j]$表示填满前$i$列还多$j$个的方案数
• 边界:$f[0][0]=f[0][2]=f[1][1]=1$
• $f[i][0]=f[i-2][0]+f[i-1][1]+f[i-2][2]$
  $f[i][1]=f[i-1][2]$
  $f[i][2]=f[i][0]+f[i-1][1]$
• 注意下标越界问题，从$2$开始推

#include <bits/stdc++.h>
using namespace std;
#define maxn 40

int n;long long f[maxn][5];

void readda_() {
   
	while(scanf("%d",&n)!=EOF) {
   
		if(n==-1) return ;
		f[0][0]=f[0][2]=f[1][1]=1;
		for(int i=2;i<=n;++i) {
   
			f[i][0]=f[i-2][0]+f[i-1][1]+f[i-2][2];
			f[i][1]=f[i-1][2];
			f[i][2]=f[i][0]+f[i-1][1];
		}
		printf("%lld\n",f[n][0]);
	}
}

记忆化搜索：

Function Run Fun

• 为了避免重复的子问题，提高效率，我们将遍历过的子问题记录下来

#include <iostream>
#include <cstdio>
#include <cstring>
#include <algorithm>
using namespace std;
#define maxn 30

int a,b,c;
long long f[maxn][maxn][maxn];

long long work_(int x,int y,int z) {
   
   if(x<=0||y<=0||z<=0) return 1;
   if(x>20||y>20||z>20) return work_(20,20,20);
   if(~f[x][y][z]) return f[x][y][z];
   if(x<y&&y<z) return f[x][y][z]=work_(x,y,z-1)+work_(x,y-1,z-1)-work_(x,y-1,z);
   else {
   
   	return f[x][y][z]=work_(x-1,y,z)+work_(x-1,y-1,z)+work_(x-1,y,z-1)-work_(x-1,y-1,z-1);
   }
}

void readda_() {
   
   while(scanf("%d%d%d",&a,&b,&c)!=EOF) {
   
   	if(a==-1&&b==-1&&c==-1) return ;
   	memset(f,-1,sizeof(f));
   	printf("w(%d, %d, %d) = %lld\n",a,b,c,work_(a,b,c));
   }
}

[SHOI2002]滑雪

#include <bits/stdc++.h>
using namespace std;
#define maxn 105
int ans,n,m,a[maxn][maxn],f[maxn][maxn];
int dx[]={
   0,1,-1,0,0};int dy[]={
   0,0,0,1,-1};

void dfs_(int x,int y) {
   
	if(f[x][y]) {
   
		ans=max(ans,f[x][y]);
		return;
	}
	f[x][y]=1;
	for(int i=1;i<=4;++i) {
   
		int px=x+dx[i];
		int py=y+dy[i];
		if(px<1||py<1||px>n||py>m) continue;
		if(a[px][py]>=a[x][y]) continue;
		if(!f[px][py]) dfs_(px,py);
		f[x][y]=max(f[x][y],f[px][py]+1);
	}
	ans=max(ans,f[x][y]);
}

void readda_() {
   
	n=read_();m=read_();
	for(int i=1;i<=n;++i) {
   
		for(int j=1;j<=m;++j) {
   
			a[i][j]=read_();
		}
	}
	memset(f,0,sizeof(f));
	ans=0;
	for(int i=1;i<=n;++i) {
   
		for(int j=1;j<=m;++j) {
   
			dfs_(i,j);
		}
	}
	printf("%d",ans);
}

最长递增路径

• 先想$dfs$,再记忆化，不难的进一步推出线性方程

线性DP：

[USACO11FEB]属牛的抗议

• 定义$f[i]$表示到$i$的最大分组
• 边界$f[i]=0,if(s[i]>=0)f[i]=1$
• 状态转移：
  $if(f[j]$&&$s[i]-s[j]>=0)$
  $f[i]=max(f[i],f[j]+1)$

#include <bits/stdc++.h>
using namespace std;
#define maxn 1010
int n,s[maxn],f[maxn];
void readda_() {
   
	n=read_();s[0]=0;
	for(int i=1;i<=n;++i) {
   
		s[i]=read_();s[i]+=s[i-1];
		if(s[i]>=0) f[i]=1;
	}
	for(int i=1;i<=n;++i) {
   
		for(int j=1;j<i;++j) {
   
			if(f[j]&&s[i]-s[j]>=0)   
			f[i]=max(f[i],f[j]+1);
		}
	}
	if(!f[n]) printf("Impossible");
	else printf("%d",f[n]);
}

文件排版

• 定义$f[i]$表示将第$i$个单词放在当前行的最小难看程度

• 边界$f[i]=INF,f[0]=0$

• 状态转移：
  $f[i]=min(f[i],f[j]+c(i,j))$

• 注意特判一个单词的情况

三个串的最长公共子序列

• 显然不是先求两个串的最长公共子序列再与第三个求
  反例：$aaaabbbbcefg$和$aabbefcg$和$aabbc$

• 定义$f[i][j][k]$表示三个串分别到$(i,j,k)$位置的最长公共子序列

• 边界$f[i][j][k]=0$

• 状态转移：
  如果$a[i]!=b[j]!=c[k]$
  $f[i][j][k]=max(f[i-1][j][k],f[i][j-1][k],f[i][j][k-1])$
  如果$a[i]==b[j]==c[k]$
  f [ i ] [ j ] [ k ] = m a x ( f [ i