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. 《Cracking the Coding Interview》——第8章：面向对象设计——题目3

    2014-04-23 18:10 题目:设计一个点唱机. 解法:英文叫Musical Jukebox.这是点唱机么?卡拉OK么?这种题目实在是云里雾里,又没有交流的余地,我索性用一个vector来表示 ...

  2. java中多态的概念

    概念: 简答来说,多态(Polymorphism)是具有表现多种形态的能力的特征.    更专业的说法是:同一个实现接口,使用不同的实例而执行不同的操作. 好处: 通过多态可以减少类中代码量,可以提高 ...

  3. finally在return之后还是之前运行

    finally在运行前打印出来是return的数据,finally是最后修改的数据,如果finally存在对返回值的修改,则以finally修改的值为准. 综上所述,finally最后运行.

  4. 孤荷凌寒自学python第六十二天学习mongoDB的基本操作并进行简单封装1

    孤荷凌寒自学python第六十二天学习mongoDB的基本操作并进行简单封装1 (完整学习过程屏幕记录视频地址在文末) 今天是学习mongoDB数据库的第八天. 今天开始学习mongoDB的简单操作, ...

  5. (原)DirectX11 深度测试(有点另类)

    (问到自己清楚就可) @Author: 白袍小道 @说明:转载随缘,评论随缘,询问建议看书和源码会得到更准确的答案 深度测试的来源.目的.做法 一.问题询问 我们带着一些问题去浏览一番 1.深度测试发 ...

  6. Spring 学习笔记（五）—— Bean之间的关系、作用域、自动装配

    继承 Spring提供了配置信息的继承机制,可以通过为<bean>元素指定parent值重用已有的<bean>元素的配置信息. <?xml version="1 ...

  7. MyBatis:SQL语句中的foreach标签的详细介绍

    foreach 也就是遍历迭代,在SQL中通常用在 in 这个关键词的后面 foreach元素的属性主要有 item,index,collection,open,separator,close. 分别 ...

  8. Linq 聚合函数

    var WReserve = (from m in _db.W_RESERVE select m).ToList().LastOrDefault(); 必须ToList() 之后才能 last fir ...

  9. 使用pdb模块调试Python

    在Python中,我们需要debug时,有三种方式: 加log语句.最简单的方式是添加print()语句来输出我们想要获知的状态或者变量,好处是简单容易操作,坏处是debug完了之后,还需要将prin ...

  10. [洛谷P4389]付公主的背包

    题目大意:有$n(n\leqslant10^5)$种物品,第$i$个物品体积为$v_i$,都有$10^5$件.给定$m(m\leqslant10^5)$,对于$s\in [1,m]$,请你回答用这些商 ...