算法板子--DP总结

01背包问题

详细推导过程：链接

1. 当前背包容量不够（j < v[i]），没得选，因此前 i 个物品最优解即为前 i−1个物品最优解：

   对应代码：f[i][j] = f[i - 1][j]。

2. 当前背包容量够，可以选，因此需要决策选与不选第 i 个物品：

   选：f[i][j] = f[i - 1][j - v[i]] + w[i]。
   不选：f[i][j] = f[i - 1][j] 。

   我们的决策是如何取到最大价值，因此以上两种情况取 max()

朴素版code

#include<iostream>

using namespace std;

const int  maxn=1010;
int n,m;
int v[maxn],w[maxn];
int f[maxn][maxn];

int main() {
    
    cin>>n>>m;
    
    for(int i=1;i<=n;i++) cin>>v[i]>>w[i];
    
    for(int i=1;i<=n;i++)
        for(int j=0;j<=m;j++){
            if(j<v[i])//背包容量小于v[i]，则不能选第i件物品
                f[i][j]=f[i-1][j];
            else
                f[i][j]=max(f[i-1][j],f[i-1][j-v[i]]+w[i]);
        }
            
            
    cout<<f[n][m];
    
    return 0;
}

优化版code

#include<iostream>

using namespace std;

const int  maxn=1010;
int n,m;

int v[maxn],w[maxn];
int f[maxn];
int main() {
    
    cin>>n>>m;
    
    for(int i=1;i<=n;i++) cin>>v[i]>>w[i];
    
    for(int i=1;i<=n;i++)
        for(int j=m;j>=v[i];j--){

            f[j]=max(f[j],f[j-v[i]]+w[i]);
        }
            
            
    cout<<f[m];
    
    return 0;
}

完全背包问题

有 N 种物品和一个容量是 V 的背包，每种物品都有无限件可用。

朴素版code

#include<iostream>

using namespace std;

const int maxn=1010;
int n,m;
int f[maxn][maxn];
int v[maxn],w[maxn];

int main(){
    cin>>n>>m;
    for(int i=1;i<=n;i++) cin>>v[i]>>w[i];
    
    for(int i=1;i<=n;i++)
        for(int j=0;j<=m;j++){
            f[i][j] = f[i-1][j];
            if(j-v[i]>=0)
            	//选了一件第i种物品之后，还可以继续选第i种物品，故仍然停留在第i层
                f[i][j]=max(f[i][j],f[i][j-v[i]]+w[i]);
        }
            
    cout<<f[n][m];
    
    
    return 0;
}

优化版code

#include<iostream>

using namespace std;

const int maxn=1010;
int n,m;
int f[maxn];
int v[maxn],w[maxn];

int main(){
    cin>>n>>m;
    for(int i=1;i<=n;i++) cin>>v[i]>>w[i];
    
    for(int i=1;i<=n;i++)
        for(int j=v[i];j<=m;j++){
            f[j]=max(f[j],f[j-v[i]]+w[i]);
        }
            
    cout<<f[m];
    
    
    return 0;
}

多重背包问题

朴素code

#include<iostream>

using namespace std;

const int maxn=105;
int f[maxn][maxn];
int n,m;
int v[maxn],w[maxn],s[maxn];

int main(){
    cin>>n>>m;
    for(int i=1;i<=n;i++) cin>>v[i]>>w[i]>>s[i];
    
    for(int i=1;i<=n;i++)
        for(int j=0;j<=m;j++)
            for(int k=0;k<=s[i];k++){
                if(j>=k*v[i])
                    f[i][j]=max(f[i][j],f[i-1][j-k*v[i]]+k*w[i]);
            }
    
    cout<<f[n][m];    
    
    return 0;
}

二进制优化方法

#include<iostream>
using namespace std;

const int N = 12010, M = 2010;

int n, m;
int v[N], w[N]; //逐一枚举最大是N*logS
int f[M]; // 体积<M

int main()
{
    cin >> n >> m;
    int cnt = 0; //分组的组别
    for(int i = 1;i <= n;i ++)
    {
        int a,b,s;
        cin >> a >> b >> s;
        int k = 1; // 组别里面的个数
        while(k<=s)
        {
            cnt ++ ; //组别先增加
            v[cnt] = a * k ; //整体体积
            w[cnt] = b * k; // 整体价值
            s -= k; // s要减小
            k *= 2; // 组别里的个数增加
        }
        //剩余的一组
        if(s>0)
        {
            cnt ++ ;
            v[cnt] = a*s; 
            w[cnt] = b*s;
        }
    }

    n = cnt ; //枚举次数正式由个数变成组别数

    //01背包一维优化
    for(int i = 1;i <= n ;i ++)
        for(int j = m ;j >= v[i];j --)
            f[j] = max(f[j],f[j-v[i]] + w[i]);

    cout << f[m] << endl;
    return 0;
}

分组背包问题

#include<iostream>

using namespace std;

int n,m;
const int maxn=110;
int v[maxn][maxn],w[maxn][maxn],s[maxn];
int f[maxn];


int main(){
    cin>>n>>m;
    for(int i=1;i<=n;i++){
        cin>>s[i];
        for(int j=1;j<=s[i];j++){
            cin>>v[i][j]>>w[i][j];
        }
    }
    
    for(int i=1;i<=n;i++)
        for(int j=m;j>=1;j--)
            for(int k=1;k<=s[i];k++)
                if(j>=v[i][k])
                    f[j]=max(f[j],f[j-v[i][k]]+w[i][k]);
    
    
    
    cout<<f[m];
    return 0;
}

最长上升子序列

s[i]的含义是，长度为i的上升子序列末尾值最小是s[i].

#include<iostream>
#include<vector>

using namespace std;

const int maxn=100000;
int n;
int a[maxn];


int main(){
    cin>>n;
    for(int i=1;i<=n;i++){
        cin>>a[i];
    }
    vector<int> s;//模拟堆栈
    s.push_back(a[1]);
    for(int i=2;i<=n;i++){
        
        if(a[i]>s.back())
            s.push_back(a[i]);
        else 
            *lower_bound(s.begin(),s.end(),a[i])=a[i];
        
    }
    cout<<s.size();
    
    
    return 0;
}

/*
例 n: 7
arr : 3 1 2 1 8 5 6

stk : 3

1 比 3 小
stk : 1

2 比 1 大
stk : 1 2

1 比 2 小
stk : 1 2

8 比 2 大
stk : 1 2 8

5 比 8 小
stk : 1 2 5

6 比 5 大
stk : 1 2 5 6

stk 的长度就是最长递增子序列的长度

*/

最长公共子序列

讲解：链接（这个思路比yxc的还要简单）

#include<iostream>

using namespace std;

const int maxn=1010;
int n,m;
char a[maxn],b[maxn];
int f[maxn][maxn];

int main(){
    
    cin>>n>>m>>a+1>>b+1;
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            if(a[i]==b[j])
                f[i][j]=f[i-1][j-1]+1;
            else 
                f[i][j]=max(f[i-1][j],f[i][j-1]);

    cout<<f[n][m]; 
    return 0;
}

最短编辑距离

#include<iostream>

using namespace std;


int n,m;
const int maxn=1010;
char a[maxn],b[maxn];
int f[maxn][maxn];

int main(){
    scanf("%d%s", &n, a + 1);
    scanf("%d%s", &m, b + 1);
    
    for(int i=0;i<=n;i++) f[i][0]=i;
    for(int i=0;i<=m;i++) f[0][i]=i;

    for (int i = 1; i <= n; i ++ )
        for (int j = 1; j <= m; j ++ )
        {
            f[i][j] = min(f[i - 1][j] + 1, f[i][j - 1] + 1);
            if (a[i] == b[j]) f[i][j] = min(f[i][j], f[i - 1][j - 1]);
            else f[i][j] = min(f[i][j], f[i - 1][j - 1] + 1);
        }
            
                
    printf("%d\n", f[n][m]);
    return 0;
}