layout: post

title: 「kuangbin带你飞」专题二十二 区间DP

author: "luowentaoaa"

catalog: true

tags:

- kuangbin

- 区间DP

- 动态规划

传送门

B.LightOJ - 1422 Halloween Costumes

题意

按顺序参加舞会，参加一个舞会要穿一种衣服，可以在参加完一个舞会后套上另一个衣服再去参加舞会，也可以在参加一个舞会的时候把外面的衣服脱了，脱到合适的衣服，但是脱掉的衣服不能再穿，参加完全部舞会最少穿多少次衣服

题解

区间dp，dp【i】【j】表示区间i->j之间的衣服数 起始情况i->i要穿一件衣服

可以从后往前推

\[dp[i][j]=dp[i+1][j]+1
\]

特别的如果第i 和第j个衣服是一样的

\[dp[i][j]=dp[i][j-1]
\]

然后开始在区间{I,J}之间选择K；

如果K的衣服和I的衣服一样代表 到K的时候可以把衣服脱到剩I的时候 这样就是K前面的衣服加上K->J的衣服

\[dp[i][j]=min(dp[i][j],dp[i][k]+dp[k+1][j])
\]

#include<bits/stdc++.h>

using namespace std;

typedef long long ll;

#define pp pair<int,int>

const ll mod=998244353;

const int maxn=1e6+50;

const ll inf=0x3f3f3f3f3f3f3f3fLL;

int gcd(int a,int b){while(b){int t=a%b;a=b;b=t;}return a;}

int lcm(int a,int b){return a*b/gcd(a,b);}

int dp[110][110];

int a[110];

int main()

{

    std::ios::sync_with_stdio(false);

    std::cin.tie(0);

    std::cout.tie(0);

    int t;

    cin>>t;

    int cnt=1;

    while(t--){

        int n;

        cin>>n;

        memset(dp,0,sizeof(dp));

        for(int i=1;i<=n;i++)cin>>a[i],dp[i][i]=1;

        for(int i=n-1;i>=1;i--){

            for(int j=i+1;j<=n;j++){

                dp[i][j]=dp[i+1][j]+1;

                if(a[i]==a[j])dp[i][j]=dp[i][j-1];

                for(int k=i+1;k<j;k++){

                    if(a[i]==a[k]){

                        dp[i][j]=min(dp[i][j],dp[i][k]+dp[k+1][j]);

                       // cout<<"i=="<<i<<"   j=="<<" dp[i][j]="<<dp[i][j]<<endl;

                    }

                }

            }

        }

        cout<<"Case "<<cnt++<<": ";

        cout<<dp[1][n]<<endl;

    }

    return 0;

}

C.POJ - 2955 Brackets

题意

给出一个的只有'(',')','[',']'四种括号组成的字符串，求最多有多少个括号满足题目里所描述的完全匹配。

题解

用dp[i][j]表示区间[i,j]里最大完全匹配数。

只要得到了dp[i][j]，那么就可以得到dp【i-1][j+1]

dp【i-1]【j+1]=dp 【i]【j]+(s[i-1]与[j+1]匹配 ? 2 : 0)

然后利用状态转移方程更新一下区间最优解即可。

\[dp[i][j]=max(dp[i][j],dp[i][k]+dp[k+1][j])
\]

#include<cstring>

#include<iostream>

#include<algorithm>

using namespace std;

typedef long long ll;

#define pp pair<int,int>

const ll mod=998244353;

const int maxn=1e3+50;

const ll inf=0x3f3f3f3f3f3f3f3fLL;

int gcd(int a,int b){while(b){int t=a%b;a=b;b=t;}return a;}

int lcm(int a,int b){return a*b/gcd(a,b);}

char s[maxn];

int dp[maxn][maxn];

int main()

{

    std::ios::sync_with_stdio(false);

    std::cin.tie(0);

    std::cout.tie(0);

    while(cin>>s+1){

        if(s[1]=='e')break;

        memset(dp,0,sizeof(dp));

        int n=strlen(s+1);

        for(int len=2;len<=n;len++)

        for(int i=1;i<=n;i++){

            int j=i+len-1;

            if(j>n)break;

            if(s[i]=='('&&s[j]==')'||s[i]=='['&&s[j]==']'){

                dp[i][j]=dp[i+1][j-1]+2;

            }

            for(int k=i;k<j;k++){

                dp[i][j]=max(dp[i][j],dp[i][k]+dp[k][j]);

            }

        }

        cout<<dp[1][n]<<endl;

    }

    return 0;

}

D.CodeForces - 149D Coloring Brackets

题意

给你一个只有括号的字符串，

1,每个括号只有三种选择：涂红色，涂蓝色，不涂色。

2,每对括号有且仅有其中一个被涂色。

3,相邻的括号不能涂相同的颜色，但是相邻的括号可以同时不涂色。

题解

因为是一个合法的括号序列。

所以每个括号与之匹配的位置是一定的。

那么就可以将这个序列分成两个区间。 (L - match[L] ) (match[L]+1, R)

用递归先处理小区间，再转移大区间。

因为条件的限制，所以记录区间的同时，还要记录区间端点的颜色。

然后就是一个递归的过程。

#include<bits/stdc++.h>

using namespace std;

typedef long long ll;

#define pp pair<int,int>

const ll mod=1e9+7;

const int maxn=1e6+50;

const ll inf=0x3f3f3f3f3f3f3f3fLL;

int gcd(int a,int b){while(b){int t=a%b;a=b;b=t;}return a;}

int lcm(int a,int b){return a*b/gcd(a,b);}

ll dp[800][800][3][3];

int a[110];

ll ans;

char s[800];

int match[800];

int stk[800];

void dfs(int l,int r){

    if(l+1==r){

        dp[l][r][0][1]=dp[l][r][0][2]=dp[l][r][1][0]=dp[l][r][2][0]=1;

        return;

    }

    if(match[l]==r){

        dfs(l+1,r-1);

        for(int i=0;i<3;i++){

            for(int j=0;j<3;j++){

                if(j!=1)dp[l][r][0][1]=(dp[l][r][0][1]+dp[l+1][r-1][i][j])%mod;

                if(j!=2)dp[l][r][0][2]=(dp[l][r][0][2]+dp[l+1][r-1][i][j])%mod;

                if(i!=1)dp[l][r][1][0]=(dp[l][r][1][0]+dp[l+1][r-1][i][j])%mod;

                if(i!=2)dp[l][r][2][0]=(dp[l][r][2][0]+dp[l+1][r-1][i][j])%mod;

            }

        }

    }

    else{

        dfs(l,match[l]);

        dfs(match[l]+1,r);

        for(int i=0;i<3;i++){

            for(int j=0;j<3;j++){

                for(int k=0;k<3;k++){

                    for(int p=0;p<3;p++){

                        if(j&&j==k)continue;

                        dp[l][r][i][p]=(dp[l][r][i][p]+(dp[l][match[l]][i][j]*dp[match[l]+1][r][k][p]%mod)%mod)%mod;

                    }

                }

            }

        }

    }

}

int main()

{

    std::ios::sync_with_stdio(false);

    std::cin.tie(0);

    std::cout.tie(0);

    int t;

    int cnt=1;

    cin>>s+1;

    int n=strlen(s+1);

    memset(match,0,sizeof(match));

    int top=-1;

    for(int i=1;i<=n;i++){

        if(s[i]=='(')stk[++top]=i;

        else match[stk[top--]]=i;

    }

    memset(dp,0,sizeof(dp));

    dfs(1,n);

    ans=0;

    for(int i=0;i<3;i++)

    for(int j=0;j<3;j++){

        ans=(ans+dp[1][n][i][j])%mod;

    }

    cout<<ans<<endl;

    return 0;

}

E.POJ - 1651 Multiplication Puzzle

题意

给出N个数，每次从中抽出一个数(第一和最后一个不能抽)，每次的得分为抽出的数与相邻两个数的乘积，直到只剩下首尾两个数为止，问最小得分是多少

题解1

枚举起点和终点 ,但是这次长度要从3开始了因为中间必要取一个

#include<cstring>

#include<iostream>

#include<algorithm>

using namespace std;

typedef long long ll;

#define pp pair<int,int>

const ll mod=998244353;

const int maxn=1e2+50;

const int inf=0x3f3f3f3f;

int gcd(int a,int b){while(b){int t=a%b;a=b;b=t;}return a;}

int lcm(int a,int b){return a*b/gcd(a,b);}

int a[maxn];

int dp[maxn][maxn];

int main()

{

    std::ios::sync_with_stdio(false);

    std::cin.tie(0);

    std::cout.tie(0);

    int n;

    while(cin>>n){

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

        memset(dp,0,sizeof(dp));

        for(int len=3;len<=n;len++){

            for(int i=1;i<=n;i++){

                int j=i+len-1;

                if(j>n)continue;

                dp[i][j]=inf;

                for(int k=i+1;k<j;k++){

                    dp[i][j]=min(dp[i][j],dp[i][k]+dp[k][j]+a[i]*a[j]*a[k]);

                    //cout<<"i=="<<i<<"   j="<<j<<"   k=="<<k<<"   dp[i][j]="<<dp[i][j]<<endl;

                }

            }

        }

        cout<<dp[1][n]<<endl;

    }

    return 0;

}

题解2

记忆化搜索

#include<cstring>

#include<iostream>

#include<algorithm>

using namespace std;

typedef long long ll;

#define pp pair<int,int>

const ll mod=998244353;

const int maxn=1e2+50;

const int inf=0x3f3f3f3f;

int gcd(int a,int b){while(b){int t=a%b;a=b;b=t;}return a;}

int lcm(int a,int b){return a*b/gcd(a,b);}

int a[maxn];

int dp[maxn][maxn];

int dfs(int s,int t){

    if(dp[s][t]!=-1)return dp[s][t];

    if(t-s<=1)return 0;

    int ans=inf;

    for(int i=s+1;i<t;i++)

        ans=min(ans,dfs(s,i)+dfs(i,t)+a[s]*a[t]*a[i]);

    dp[s][t]=ans;

    return ans;

}

int main()

{

    std::ios::sync_with_stdio(false);

    std::cin.tie(0);

    std::cout.tie(0);

    int n;

    while(cin>>n){

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

        memset(dp,-1,sizeof(dp));

        cout<<dfs(1,n)<<endl;

    }

    return 0;

}

G.HDU - 4283 You Are the One

题意

有一个队列，每个人有一个愤怒值D，如果他是第K个上场，不开心指数就为(K-1)*D。但是边上有一个小黑屋（其实就是个堆栈），可以一定程度上调整上场程序

题意

对于区间 [x, y] 中的数字 a[x]，有可能第1个出栈，也有可能最后一个出栈，如果第k个出栈的话前i+1->k个数字要按照顺序出栈，然后a[x]出栈，[k+1,y]出栈。对于a[x]要加上a[x]×(k-i）的花费，对于[x+1, y]要加上sum(x+1, y) * (i-x+1)的额外花费。

#include<bits/stdc++.h>

using namespace std;

typedef long long ll;

#define pp pair<int,int>

const ll mod=998244353;

const int maxn=1e6+50;

const ll inf=0x3f3f3f3f3f3f3f3fLL;

int gcd(int a,int b){while(b){int t=a%b;a=b;b=t;}return a;}

int lcm(int a,int b){return a*b/gcd(a,b);}

int dp[110][110];

int a[110];

int ans[110];

int sum[110];

int main()

{

    std::ios::sync_with_stdio(false);

    std::cin.tie(0);

    std::cout.tie(0);

    int t;

    int cnt=1;

    cin>>t;

    while(t--){

        int n;

        cin>>n;

        memset(dp,0,sizeof(dp));

        memset(sum,0,sizeof(sum));

        for(int i=1;i<=n;i++)cin>>a[i],dp[i][i]=0,sum[i]=sum[i-1]+a[i];

        for(int len=2;len<=n;len++)

        for(int i=1;i<=n;i++){

            int j=i+len-1;

            if(j>n)break;

            dp[i][j]=inf;

            for(int k=i;k<=j;k++){

                dp[i][j]=min(dp[i][j],a[i]*(k-i)+(k-i+1)*(sum[j]-sum[k])+dp[i+1][k]+dp[k+1][j]);

            }

        }

        cout<<"Case #"<<cnt++<<": ";

        cout<<dp[1][n]<<endl;

    }

    return 0;

}

H.HDU - 2476 String painter

题意

给出两个长度相同的字符串a, b.对a可以进行区间更新(选择一个区间，把这个区间全部更新成相同的字符)，问要最少操作a多少次，a才能等于b？

题解

先处理一下b串 用DP存b串i到j的改变次数 然后再对a串进行状态转移。对于串b，最要是看两个相同字符之间的区间了。我们可以预处理出来改变区间[x, y]里面的字符，操作的最小次数用dp[x, y]来存储。

　　预处理完了以后，我们就可以用dp来匹配a了。ans[i] 表示 a字符区间[0, i] 与 b字符区间[0, i]相同的最小操作数目。状态转移就是：ans[i] = min (dp[0, i], ans[j] + dp【j+1】[i], ans【i-1](a[i] == b[i]))

#include<bits/stdc++.h>

using namespace std;

typedef long long ll;

#define pp pair<int,int>

const ll mod=998244353;

const int maxn=1e6+50;

const ll inf=0x3f3f3f3f3f3f3f3fLL;

int gcd(int a,int b){while(b){int t=a%b;a=b;b=t;}return a;}

int lcm(int a,int b){return a*b/gcd(a,b);}

int dp[110][110];

char a[110],b[110];

int ans[110];

int main()

{

    std::ios::sync_with_stdio(false);

    std::cin.tie(0);

    std::cout.tie(0);

    while(cin>>a>>b){

        memset(dp,0,sizeof(dp));

        int n=strlen(b);

        for(int i=0;i<n;i++)dp[i][i]=1;

        for(int len=2;len<=n;len++)

        for(int i=0;i<n;i++){

            int j=i+len-1;

            if(j>=n)break;

            if(b[i]==b[j])dp[i][j]=dp[i+1][j];

            else

                dp[i][j]=dp[i+1][j]+1;

            for(int k=i+1;k<j;k++){

                if(b[k]==b[i])

                    dp[i][j]=min(dp[i][j],dp[i+1][k]+dp[k+1][j]);

            }

        }

        ans[0]=a[0]==b[0]?0:1;

        for(int i=1;i<n;i++){

            ans[i]=dp[0][i];

            if(a[i]==b[i])ans[i]=min(ans[i],ans[i-1]);

            for(int j=0;j<i;j++)

                ans[i]=min(ans[i],ans[j]+dp[j+1][i]);

        }

        cout<<ans[n-1]<<endl;

    }

    return 0;

}

「kuangbin带你飞」专题二十二 区间DP的更多相关文章

  1. 「kuangbin带你飞」专题二十 斜率DP

    layout: post title: 「kuangbin带你飞」专题二十 斜率DP author: "luowentaoaa" catalog: true tags: mathj ...

  2. 「kuangbin带你飞」专题十二 基础DP

    layout: post title: 「kuangbin带你飞」专题十二 基础DP author: "luowentaoaa" catalog: true tags: mathj ...

  3. 「kuangbin带你飞」专题十九 矩阵

    layout: post title: 「kuangbin带你飞」专题十九 矩阵 author: "luowentaoaa" catalog: true tags: mathjax ...

  4. 「kuangbin带你飞」专题十八 后缀数组

    layout: post title: 「kuangbin带你飞」专题十八 后缀数组 author: "luowentaoaa" catalog: true tags: - kua ...

  5. 「kuangbin带你飞」专题十四 数论基础

    layout: post title: 「kuangbin带你飞」专题十四 数论基础 author: "luowentaoaa" catalog: true tags: mathj ...

  6. 「kuangbin带你飞」专题十七 AC自动机

    layout: post title: 「kuangbin带你飞」专题十七 AC自动机 author: "luowentaoaa" catalog: true tags: - ku ...

  7. 「kuangbin带你飞」专题十五 数位DP

    传送门 A.CodeForces - 55D Beautiful numbers 题意 一个正整数是 漂亮数 ,当且仅当它能够被自身的各非零数字整除.我们不必与之争辩,只需计算给定范围中有多少个漂亮数 ...

  8. kuangbin带你飞 生成树专题 ： 次小生成树; 最小树形图;生成树计数

    第一个部分 前4题 次小生成树 算法:首先如果生成了最小生成树,那么这些树上的所有的边都进行标记.标记为树边. 接下来进行枚举,枚举任意一条不在MST上的边,如果加入这条边,那么肯定会在这棵树上形成一 ...

  9. KUANGBIN带你飞

    KUANGBIN带你飞 全专题整理 https://www.cnblogs.com/slzk/articles/7402292.html 专题一 简单搜索 POJ 1321 棋盘问题    //201 ...

随机推荐

  1. 【Luogu P2257】YY 的 GCD

    题目 求: \[ \sum_{i = 1}^n \sum_{j = 1}^m [\gcd(i, j) \in \mathbb P] \] 有 \(T\) 组数据, \(T\le 10^4, n, m\ ...

  2. css 之 border-radius属性

    css中给盒子设置圆角可以通过 border-radius 属性来实现(具体原理就不深入探讨了); 在开发过程中都会遇到浏览器兼容问题,这问题其实也不难解决,无非就是加上私有前缀,在这里先忽略掉.   ...

  3. 牛客网暑期ACM多校训练营（第一场）：J-Different Integers（分开区间不同数+树状数组）

    链接:J-Different Integers 题意:给出序列a1, a2, ..., an和区间(l1, r1), (l2, r2), ..., (lq, rq),对每个区间求集合{a1, a2, ...

  4. ExtJs学习之MessAgeBox的使用

    1.Ext.MessageBox.alert() 调用格式: alert( String title, String msg, [Function fn], [Object scope] ) 参数说明 ...

  5. django通用视图之TemplateView和ListView简单介绍

    django支持类视图,与此同时django为我们提供了许多非常好用的通用视图供我们使用,这其中TemplateView.ListView和DetailView是我们经常使用到的,这里就对Templa ...

  6. URAL 1934 spfa算法

    D - Black Spot Time Limit:1000MS     Memory Limit:65536KB     64bit IO Format:%I64d & %I64u Subm ...

  7. StringBuilder_学习笔记

    参考:https://www.jianshu.com/p/160c9be0b132 连接符号 "+" 本质 字符串变量(非final修饰)通过 "+" 进行拼接 ...

  8. [经验分享]WeTouch中使用VueInputCode

    如果你想在WeTouch中使用以下输入框的效果,需要做以下几个步骤: 第一步:安装Node.js 第二步:使用Npm在全局安装vue-input-code ,使用命令提示符,在Node安装根目录输入以 ...

  9. 解决PKIX path building failed

    起因 上周在生产环境部署时,把安全证书加到k8s-ingress中时发现报该错误 解决 找网上解决方案,因为这种问题相对比较少见,也没百度,直接谷歌,找到解决方案如下:https://stackove ...

  10. Action中动态方法的调用 Action中通配符的使用 Result的配置

       Action中动态方法的调用 动态方法调用(Dynamic Method Invocation,DMI) 标识符:! 一.通过以下选中的文件来查看是否禁止调用动态方法