动态规划之数字三角形模型+最长上升子序列模型

首先,我们从集合角度重新看待DP:

直接看题:https://www.acwing.com/problem/content/1029/

就是取纸条的原题,我们令f[i1,j1,i2,j2]表示从(1,1),(1,1)分别走到(i1,j1),(i2,j2)的路径的max

i1+j1==i2+j2,于是我们可以把状态变为f[k,i1,i2]表示走到(i1,k-i1),(i2,k-i2)的最大值。

对于最后一个状态,我们可以按照上下分成4中,同时转移时重合的化加w[i][j],否则就分别加一个

下面是AC代码:

cpp 复制代码
#include <iostream>
#include <algorithm>

using namespace std;

const int N = 15;

int n;
int w[N][N];
int f[N * 2][N][N];

int main()
{
    scanf("%d", &n);

    int a, b, c;
    while (cin >> a >> b >> c, a || b || c) w[a][b] = c;

    for (int k = 2; k <= n + n; k ++ )
        for (int i1 = 1; i1 <= n; i1 ++ )
            for (int i2 = 1; i2 <= n; i2 ++ )
            {
                int j1 = k - i1, j2 = k - i2;
                if (j1 >= 1 && j1 <= n && j2 >= 1 && j2 <= n)
                {
                    int t = w[i1][j1];
                    if (i1 != i2) t += w[i2][j2];
                    int &x = f[k][i1][i2];
                    x = max(x, f[k - 1][i1 - 1][i2 - 1] + t);
                    x = max(x, f[k - 1][i1 - 1][i2] + t);
                    x = max(x, f[k - 1][i1][i2 - 1] + t);
                    x = max(x, f[k - 1][i1][i2] + t);
                }
            }

    printf("%d\n", f[2*n][n][n]);
    return 0;
}

接题:https://www.acwing.com/problem/content/189/

看到数据范围可以确定爆搜(BFS因为一层层存储量太大不推荐),这样就转换成导弹拦截的思路了:从前往后,对于一个导弹找到比他高的min插入即可,若没有就再开一个(注意防御系统一定是递增的)

下面是AC代码:

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
int n;
int a[10001];
int ans=50;
int up[100010],down[10000];
void dfs(int x,int u,int d)
{
    if(u+d>=ans) return;
    if(x==n+1){
        ans=min(ans,u+d);
        return;
    }
    //先上
    int f=-1;
    for(int i=1;i<=u;i++)
    {
        if(up[i]>a[x]){
            f=i;
            break;
        }
        if(i==u) f=-1;
    }
    if(f==-1){
        up[u+1]=a[x];
        dfs(x+1,u+1,d);
        up[u+1]=0;
    }
    else{
        int t=up[f];
        up[f]=a[x];
        dfs(x+1,u,d);
        up[f]=t;
    }
    f=-1;
    for(int i=1;i<=d;i++)
    {       if(down[i]<a[x]){
            f=i;
            break;
        }
        if(i==d) f=-1;
    }
    if(f==-1){
        down[d+1]=a[x];
        dfs(x+1,u,d+1);
        down[d+1]=0;
    }
    else{
        int t=down[f];
        down[f]=a[x];
        dfs(x+1,u,d);
        down[f]=t;
    }
}
int main()
{
    while(cin>>n)
    {
        if(!n) break;
        for(int i=1;i<=n;i++) cin>>a[i];
        memset(up,0,sizeof(up));
        memset(down,0,sizeof(down));
        ans=50;
        dfs(1,0,0);
        cout<<ans<<endl;
    }
}

接题:https://www.acwing.com/problem/content/1018/

确定f[i][j]的集合:由第一个序列前i个字母,第二个前j个并以b[j]的公共上升子序列。

表示max,然后考虑计算:先按是否包含a[i]分成两类:不包含:f[i-1][j],包含:a[i]==b[j]

我们又可以把b分为以b[1]...b[j-1]以及空,对于b[k]就是f[i-1][k]+1

下面是AC代码:

cpp 复制代码
#include<bits/stdc++.h>
using namespace std;
int n,a[100010],b[100010];
int dp[3010][3005];
int main()
{
    cin>>n;
    for(int i=1;i<=n;i++) cin>>a[i];
    for(int i=1;i<=n;i++) cin>>b[i];
    for(int i=1;i<=n;i++)
    {
        int maxv=0;
        for(int j=1;j<=n;j++)
        {
             
             dp[i][j]=dp[i-1][j];
             if(a[i]==b[j])
             {
                 dp[i][j]=max(dp[i][j],1);
                 dp[i][j]=max(dp[i][j],maxv);
                 
             }
             if(b[j]<a[i]) maxv=max(maxv,dp[i-1][j]+1);
        }
    }
    int res=0;
    for(int i=1;i<=n;i++) res=max(res,dp[n][i]);
    cout<<res;
}
相关推荐
蒙奇D索大18 分钟前
【数据结构】图论最短路圣器:Floyd算法如何用双矩阵征服负权图?
数据结构·算法·矩阵·图论·图搜索算法
白熊1881 小时前
【推荐算法】注意力机制与兴趣演化:推荐系统如何抓住用户的心?
算法·php·推荐算法
Sun_light1 小时前
LeetCode 59.「螺旋矩阵」
javascript·算法·面试
数据与人工智能律师2 小时前
当机床开始“思考”,传统“制造”到“智造”升级路上的法律暗礁
大数据·网络·算法·云计算·区块链
IC 见路不走2 小时前
LeetCode 第73题:矩阵置零
算法·leetcode·矩阵
byte轻骑兵2 小时前
【C++特殊工具与技术】优化内存分配(五):显式析构函数的调用
开发语言·c++
黑听人3 小时前
【力扣 简单 C】141. 环形链表
c语言·开发语言·数据结构·算法·leetcode
谷雨不太卷3 小时前
AVL树的实现
数据结构·c++·算法
别来无恙1495 小时前
岛屿周长问题的三种解法:直接计数法、数学计算法与深度优先搜索
java·c++·算法·深度优先·dfs
周振超的6 小时前
c++编译第三方项目报错# pragma warning( disable: 4273)
开发语言·c++