【算法】动态规划DP

时间复杂度：状态数（f）×转移复杂度O(1)

背包

体积v，价值w
背包装得下，总价值最大

01背包

$f(v,j)$表示最大的价值（v个物体，j剩余）
集合划分为两部分：
包含v，不包含v
$f(v,j)=max(f(v-1,j),f(v-1,j-v_v))+w_v$
Dp优化：
空间上可以优化：
二维——一维：f(v)只用到了f(v-1)：滚动数组
特点：每件物品只用一次

#include<bits/stdc++.h>

using namespace std;

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

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])
                f[i][j]=max(f[i][j],f[i-1][j-v[i]]+w[i]);
        }
    }
    
    cout<<f[n][m];
}

完全背包问题

特点：每件物品无限次

#include<bits/stdc++.h>

using namespace std;

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

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++)
        {
            for(int k=0; k*v[i]<=j;k++)
                f[i][j]=max(f[i-1][j],f[i-1][j-k*v[i]]+k*w[i]);
        }
    }
    
    cout<<f[n][m];
}

可以优化：

#include<bits/stdc++.h>

using namespace std;

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

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])
                f[i][j]=max(f[i][j],f[i][j-v[i]]+w[i]);//
        }
    }
    
    cout<<f[n][m];
}

多重背包问题

特点：每件物品个数不同

#include<bits/stdc++.h>

using namespace std;

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

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*v[i]<=j&&k<=s[i] ;k++)
                f[i][j]=max(f[i][j],f[i-1][j-k*v[i]]+k*w[i]);//变成了ij
        }
    }
    
    cout<<f[n][m];
}

分组背包问题

特点：分组，每组只能选一个

#include<bits/stdc++.h>

using namespace std;
const int N = 110;
int n,m;
int v[N][N], w[N][N], s[N];
int f[N];

int main()
{
    cin>>n>>m;
    for(int i = 1;i<=n;i++)
    {
        cin>>s[i];
        for(int j =0;j<s[i];j++)
            cin>>v[i][j]>>w[i][j];
    }
    
    for(int i=1;i<=n;i++)
        for(int j = m;j>=0;j--)//体积从大到小枚举
            for(int k=0;k<s[i];k++)
                if(v[i][k]<=j) //第i组第k个
                    f[j] = max(f[j],f[j-v[i][k]]+w[i][k]);
                    
    cout<<f[m]<<endl;
    return 0;          
}

线性dp

数字三角形

我觉得我现在比较理解了，dp永远是一个矩阵（二维的话）这里的维数指的是每一种方式的性质，也代表每一种选择方式。
状态表示：如何定义[i][j]的值
状态计算：找不同方式之间的方程，用到划分的思想
通过这两步，就可以确定方程：

#include<bits/stdc++.h>

using namespace std;

const int N=510, INF=1e9;
int a[N][N];
int n;
int f[N][N];

int main()
{
    cin>>n;
    for(int i = 1; i <= n; i++)
        for(int j =1;j<=i;j++)
            cin>>a[i][j];

    for(int i=0;i<=n;i++)
        for(int j=0;j<=i+1;j++)
            f[i][j]=-INF;   //初始化第0行

    f[1][1]  = a[1][1];     //初始化第1行

    for(int i = 2; i <= n; i++) 
        for(int j =1;j<=i;j++)
            f[i][j]=max(f[i-1][j-1]+a[i][j],f[i-1][j]+a[i][j]);

    int res=-INF;
    for(int i=1;i<=n;i++) res=max(res,f[n][i]);

    cout<<res;

}

最长上升子序列

#include<bits/stdc++.h>

using namespace std;
const int N=1010;
int f[N];
int a[N];
int n;

int main()
{
    cin>>n;
    for(int i=1;i<=n;i++) cin>>a[i];
    
    for(int i=1;i<=n;i++)
    {
        f[i]=1;
        for(int j =1;j<n;j++)
        {
            if(a[i]>a[j])
            f[i]=max(f[i],f[j]+1);
        }
    }
    int res=0;
    for(int i=1;i<=n;i++) res=max(res,f[i]);
    cout<<res;
}

最长公共子序列

$f[i-1,j]$包含了00和01两种情况

整数划分

#include <iostream>  //完全背包
#include <algorithm>

using namespace std;

const int N = 1010, mod = 1e9 + 7;

int n;
int f[N];

int main()
{
    cin >> n;

    f[0] = 1;
    for (int i = 1; i <= n; i ++ )
        for (int j = i; j <= n; j ++ )
            f[j] = (f[j] + f[j - i]) % mod;

    cout << f[n] << endl;

    return 0;
}

区间dp

石子合并

区间划分感觉这样可以，找到中间线

S可以前缀和来做

枚举：区间长度从小到大

#include<bits/stdc++.h>

using namespace std;
const int N=310;
int n;
int s[N];
int f[N][N];

int main()
{
    cin>>n;
    for(int i=1;i<=n;i++) cin>>s[i];
    for(int i=1;i<=n;i++) s[i]+=s[i-1];
    
    for(int len=2;len<=n;len++)
        for(int i=1;i+len-1<=n;i++)
        {
        int l=i,r=i+len-1;
        f[l][r]=1e8;
        for(int k=l;k<r;k++)
            {f[l][r]=min(f[l][r],f[l][k]+ f[k+1][r] + s[r] - s[l-1]);}
            
        
        }
    cout<<f[1][n];
}

数位统计dp

#include <iostream>
#include <algorithm>
#include <vector>

using namespace std;

const int N = 10;

/*

001~abc-1, 999

abc
    1. num[i] < x, 0
    2. num[i] == x, 0~efg
    3. num[i] > x, 0~999

*/

int get(vector<int> num, int l, int r)
{
    int res = 0;
    for (int i = l; i >= r; i -- ) res = res * 10 + num[i];
    return res;
}

int power10(int x)
{
    int res = 1;
    while (x -- ) res *= 10;
    return res;
}

int count(int n, int x)
{
    if (!n) return 0;

    vector<int> num;
    while (n)
    {
        num.push_back(n % 10);
        n /= 10;
    }
    n = num.size();

    int res = 0;
    for (int i = n - 1 - !x; i >= 0; i -- )
    {
        if (i < n - 1)
        {
            res += get(num, n - 1, i + 1) * power10(i);
            if (!x) res -= power10(i);
        }

        if (num[i] == x) res += get(num, i - 1, 0) + 1;
        else if (num[i] > x) res += power10(i);
    }

    return res;
}

int main()
{
    int a, b;
    while (cin >> a >> b , a)
    {
        if (a > b) swap(a, b);

        for (int i = 0; i <= 9; i ++ )
            cout << count(b, i) - count(a - 1, i) << ' ';
        cout << endl;
    }

    return 0;
}