浅谈区间动态规划

围绕几道题说起。。石子归并、涂色、括号序列

啥是区间动态规划呢,我觉得似乎是指在一段区间上的dp,通过枚举左右子区间来求出解。

那么问题来了,如何去枚举左右子区间呢?

一般来说都是循环一个变量len,表示区间长度,然后循环左区间从开始到结尾,一般来说是1~n。

对于区间dp的话,我大致理解就是先求出小区间(部分)最优解,然后一个又一个小区间合并成稍微大点的大区间,最后合成答案——即总区间。

所以代码就这玩意:

for(int len=1;len<=n;len++)
{
    for(int l=1,r;(r=l+len)<=n;l++)
    {
        //do something
        for(int k=l;k<r;k++)
        {
            //update dp array,such as'min(dp[l][r],dp[l][i]+dp[i+1][r])'
        }
    }
}

所以就是这样咯,循环一个长度,然后枚举区间。

区间的话一定要注意是l<=k<r,不然的话[i+1][r]直接6了。。。

好像。。。区间dp就这些内容了吧,哦对了还有平行四边形优化!

普通的区间dp的时间复杂度大约是O(n^3),从三个嵌套的for就能看出来。

不过有一些区间dp可以优化成O(n^2),要用一个叫做“四边形不等式”的东西进行优化。

大致思想就是保存枚举的k中最优子区间,每次可以将枚举k时的时间复杂度去掉,所以就只剩下了长度与左区间的枚举。

四边形不等式优化

下面就是三道入门题。。。


//石子归并
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cstdlib>
using namespace std;
char in[1000];
int dp[1000][1000],n,w[1000],sum[1000];
int main()
{
    scanf("%d",&n);
    for(int i=1;i<=n;i++)
    {
        scanf("%d",&w[i]);
        sum[i]=sum[i-1]+w[i];//由于合并的是一个区间,所以记录一下前缀和
    }
    for(int len=1;len<=n;len++)//循环长度
    {
        for(int l=1,r;(r=l+len)<=n;l++)
        {
            dp[l][r]=0x3fffffff;//初始化为INF
            for(int i=l;i<r;i++)//枚举中点
                dp[l][r]=min(dp[l][r],dp[l][i]+dp[i+1][r]+sum[r]-sum[l-1]);//判断区间l,r分割成两个小区间后是否更优
        }
    }
    printf("%d\n",dp[1][n]);//答案
}

//bzoj1260涂色
#include <iostream>
#include <algorithm>
#include <cstdio>
#include <cstring>
#include <cstdlib>
using namespace std;
char in[1000];
int dp[1000][1000];
int main()
{
    scanf("%s",in+1);//方便下面dp
    int n=strlen(in+1);
    for(int i=1;i<=n;i++)for(int j=1;j<=n;j++)dp[i][j]=i==j?1:0x3fffffff;//对于一个点的涂色次数,只能是1次
    for(int len=1;len<=n;len++)//循环长度
    {
        for(int l=1,r;(r=l+len)<=n;l++)//生成左右区间
        {
            if(in[l]==in[r])//如果相等那么就可以如下这么转移咯
            {
                if(len==1)dp[l][r]=1;//如果区间长度为1,也就是说l,r是相邻的两个格子,所以只能一笔涂上
                else dp[l][r]=min(dp[l+1][r-1]+1,min(dp[l][r-1],dp[l+1][r]));
                /*
                    否则的话说明可以从区间l,r中进行转移。
                    判断dp[l][r]所包含的三个子区间
                    然后就是dp[l+1][r]跟dp[l][r-1]了。
                    先把最右端/最左端为起点一笔涂到另一头,应该是这意思吧?
                */
            }
            else//否则的话只能将区间l,r分割成两个小区间然后min咯
                for(int i=l;i<r;i++)
                    dp[l][r]=min(dp[l][r],dp[l][i]+dp[i+1][r]);
        }
    }
    printf("%d\n",dp[1][n]);//答案
}

//tyvj P1193 括号序列
/*
    至于为什么这个dp是对的,小谈一下我的个人理解。
    由于区间长度len是从小到的枚举的,所以每一轮都是从很小的单位区间开始更新,似乎可以看作bfs(雾)?
    因为len是由小到大进行枚举,所以每一次min的时候都是从将已经判断好的子区间进行min。
    或许用滚雪球形容会好点?从一个小区间慢慢滚成了一个大区间?
    以上就是我对于初级区间dp的理解。。。
*/
#include <cstring>
#include <iostream>
#include <algorithm>
#include <cstdio>
using namespace std;
char in[1000];
int dp[1000][1000],n;
bool test(char a,char b)//判断a,b是否匹配
{
    if(a=='(')return b==')';
    if(a=='[')return b==']';
    if(a=='<')return b=='>';
    if(a=='{')return b=='}';
    return 0;
}
int main()
{
    /*
        dp[i][i]=1 若是单个括号的话只能是添加所对应的括号完成匹配,所以为1
        dp[l][r]:区间l,r的最小添加括号数量
        dp[l][r]=min{n,dp[l+1][r-1](match(s[l],s[r]),dp[l][k]+dp[k+1][r]|l<=k<r}
    */
    gets(in+1);//为了能够顺手的写代码,所以下标从1开始
    n=strlen(in+1);
    while(in[n]=='\n'||in[n]=='\r')n--;
    //不知道为啥gets会莫名其妙的读入一个换行/回车,所以第一次就WA了。。。所以加上这个while判一下是否有换行,不过似乎最多就循环一次?
    for(int i=1;i<=n;i++)dp[i][i]=1;
    for(int len=1;len<=n;len++)
    {
        for(int l=1,r;(r=l+len)<=n;l++)
        {
            dp[l][r]=r-l+1;//边界,对于区间l,r一共有r-l+1个括号,所以最坏情况下需要添加r-l+1个括号,不过似乎写成n也能A掉?好吧写成dp[l][r]=n会快点。。。毕竟避免了加减法运算
            if(test(in[l],in[r]))//如果左右两端匹配,那么就可以进行一次转移
                dp[l][r]=min(dp[l][r],dp[l+1][r-1]);//如果当前区间左右都是匹配的括号的话,状态转移为dp[l+1][r-1],也就是说转移到去掉括号后的序列
            for(int i=l;i<r;i++)//枚举区间l,r中的所有区间
                dp[l][r]=min(dp[l][r],dp[l][i]+dp[i+1][r]);
        }
    }
    printf("%d\n",dp[1][n]);//最后的结果就是dp[1][n]
}
最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 215,012评论 6 497
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,628评论 3 389
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 160,653评论 0 350
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,485评论 1 288
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,574评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,590评论 1 293
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,596评论 3 414
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,340评论 0 270
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,794评论 1 307
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 37,102评论 2 330
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,276评论 1 344
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,940评论 5 339
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,583评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,201评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,441评论 1 268
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,173评论 2 366
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,136评论 2 352

推荐阅读更多精彩内容

  • 背景 一年多以前我在知乎上答了有关LeetCode的问题, 分享了一些自己做题目的经验。 张土汪:刷leetcod...
    土汪阅读 12,743评论 0 33
  • 动态规划(Dynamic Programming) 本文包括: 动态规划定义 状态转移方程 动态规划算法步骤 最长...
    廖少少阅读 3,275评论 0 18
  • 归去来兮。 1.1 说明 本篇为《挑战程序设计竞赛(第2版)》[http://www.ituring.com.cn...
    尤汐Yogy阅读 14,321评论 0 160
  • 以前在linux下看源码一直习惯用source insight,然而最近发现sublime text的插件好多好强...
    少阁主_enfj阅读 9,081评论 1 2
  • 在听完陈奕迅的歌 又换到了twins 时间一点一点 歌声辗转离合 一个一个字符音符 像极了奔跑的雪花 在雪地翻滚打...
    言由衷阅读 237评论 0 0