区间DP

最新推荐文章于 2022-07-25 23:39:22 发布
最新推荐文章于 2022-07-25 23:39:22 发布
阅读量142 收藏
点赞数
原文链接：http://www.cnblogs.com/HHHG/p/11203890.html
版权

概念

以区间长度作为阶段，使用区间的端点描述每个维度（其实就是dp数组的维度）。在区间DP中，一个状态由若干个比它更小且包含于它的区间所代表的区间转移而来，因此区间DP的决策常常是划分区间的方法。


例题

合并石子

难度：普及+/提高

题目描述

将 n 堆石子绕圆形操场排放，现要将石子有序地合并成一堆。规定每次只能选相邻的两堆合并成新的一堆，并将新的一堆的石子数记做该次合并的得分。
请编写一个程序，读入堆数 nnn 及每堆的石子数，并进行如下计算：

  1. 选择一种合并石子的方案，使得做 n−1n-1n−1 次合并得分总和最大。
  2. 选择一种合并石子的方案，使得做 n−1n-1n−1 次合并得分总和最小。

输入

输入第一行一个整数 n，表示有 n 堆石子。
第二行 n 个整数，表示每堆石子的数量。

 

输出

输出共两行：
第一行为合并得分总和最小值，
第二行为合并得分总和最大值。

提示

对于100%的数据，有1<=n<=200。

思路

其实这就是一道简单的模板题，只用处理一下环的问题，将环剪成链，再分别求出MAX与MIN即可。

#include<bits/stdc++.h>
#define fu(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
int read()
{
    int x=0,k=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')
        k=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        x=(x<<1)+(x<<3)+ch-'0';
        ch=getchar();
    }
    return x*k;
}
int n,se[405],qe[405],fx[405][405],fn[405][405],l,r,maxx,minn;
int main()
{
    n=read();
    memset(qe,0,sizeof(qe));
    memset(fx,0,sizeof(fx));
    memset(fn,0,sizeof(fn));
    fu(i,1,n)
    {
        se[i]=read();
        se[n+i]=se[i];
        qe[i]=qe[i-1]+se[i];
        fx[i][i]=0;
        fn[i][i]=0;
    }
    fu(i,n+1,2*n)
    {
        qe[i]=qe[i-1]+se[i];
        fx[i][i]=fn[i][i]=0;
    }
    fu(jie,2,n)
        fu(l,1,2*n-jie+1)
        {
            r=l+jie-1;
            fx[l][r]=-1;
            fn[l][r]=0x7fffffff;
            fu(k,l,r-1)
            {
                fx[l][r]=max(fx[l][r],fx[l][k]+fx[k+1][r]);
                fn[l][r]=min(fn[l][r],fn[l][k]+fn[k+1][r]);
            }
            fx[l][r]+=qe[r]-qe[l-1];
            fn[l][r]+=qe[r]-qe[l-1];
        }
    maxx=-1;minn=0x7fffffff;
    fu(i,1,n)
    {
        maxx=max(maxx,fx[i][i+n-1]);
        minn=min(minn,fn[i][i+n-1]);
    }
    cout<<minn<<" "<<maxx;
    return 0;
}

 

能量项链 

难度：普及+/提高

题目描述  

在Mars星球上，每个Mars人都随身佩带着一串能量项链。在项链上有N颗能量珠。能量珠是一颗有头标记与尾标记的珠子，这些标记对应着某个正整数。并且，对于相邻的两颗珠子，前一颗珠子的尾标记一定等于后一颗珠子的头标记。因为只有这样，通过吸盘（吸盘是Mars人吸收能量的一种器官）的作用，这两颗珠子才能聚合成一颗珠子，同时释放出可以被吸盘吸收的能量。如果前一颗能量珠的头标记为m，尾标记为r，后一颗能量珠的头标记为r，尾标记为n，则聚合后释放的能量为（Mars单位），新产生的珠子的头标记为m，尾标记为n。
需要时，Mars人就用吸盘夹住相邻的两颗珠子，通过聚合得到能量，直到项链上只剩下一颗珠子为止。显然，不同的聚合顺序得到的总能量是不同的，请你设计一个聚合顺序，使一串项链释放出的总能量最大。
例如：设N=4，4颗珠子的头标记与尾标记依次为(2，3) (3，5) (5，10) (10，2)。我们用记号⊕表示两颗珠子的聚合操作，(j⊕k)表示第j，k两颗珠子聚合后所释放的能量。则第4、1两颗珠子聚合后释放的能量为：
(4⊕1)=10*2*3=60。
这一串项链可以得到最优值的一个聚合顺序所释放的总能量为：
((4⊕1)⊕2)⊕3）=10*2*3+10*3*5+10*5*10=710。

输入

每组输入数据的第一行是一个正整数N（4≤N≤100），表示项链上珠子的个数。第二行是N个用空格隔开的正整数，所有的数均不超过1000。第i个数为第i颗珠子的头标记（1≤i≤N），当i<N时，第i颗珠子的尾标记应该等于第i+1颗珠子的头标记。第N颗珠子的尾标记应该等于第1颗珠子的头标记。
至于珠子的顺序，你可以这样确定：将项链放到桌面上，不要出现交叉，随意指定第一颗珠子，然后按顺时针方向确定其他珠子的顺序。

输出

每组输出只有一行，是一个正整数E（E≤2.1*109），为一个最优聚合顺序所释放的总能量。

思路

这道题表述成链就一定有它的道理，但它实际上还是环，我们用tou存头标记，用wei存尾标记，截环为链，即可。

#include<bits/stdc++.h>
#define fu(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
int read()
{
    int x=0,k=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')
        k=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        x=(x<<1)+(x<<3)+ch-'0';
        ch=getchar();
    }
    return x*k;
}
int tou[205],wei[205],f[205][205];
int E=0,n;
int main()
{
    n=read();
    fu(i,1,n)
    {
        tou[i]=read();
        tou[n+i]=tou[i];
    }
    fu(i,1,2*n-1)
    {
        wei[i]=tou[i+1];
        f[i][i]=0;
    }
    wei[2*n]=tou[1];
    fu(jie,1,n-1)
        fu(l,1,2*n-jie)
        {
            int r=l+jie;
            fu(k,l,r-1)
                f[l][r]=max(f[l][r],f[l][k]+f[k+1][r]+tou[l]*wei[k]*wei[r]);
        }
    fu(i,1,n)
        E=max(E,f[i][i+n-1]);
    cout<<E;
    return 0;
 }

 

凸多边形的划分

难度：提高+/省选-

题目描述

给定一个具有 N 个顶点的凸多边形，将顶点从 1 至 N 标号，每个顶点的权值都是一个正整数。将这个凸多边形划分成 N-2 个互不相交的三角形，试求这些三角形顶点的权值乘积和至少为多少。 

输入

输入第一行为顶点数 N
第二行依次为顶点 1 至顶点 N 的权值。

输出

输出仅一行，为这些三角形顶点的权值乘积和的最小值。

提示

对于 100% 的数据，有 N<= 50，每个点权值小于 10^9。

思路

#include<bits/stdc++.h>
#define fu(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
#define ll long long 
using namespace std;
int read()
{
    int x=0,k=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')
        k=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        x=(x<<1)+(x<<3)+ch-'0';
        ch=getchar();
    }
    return x*k;
}
ll f[105][105][105],a[105];
ll s1[105],s2[105],s3[105];
int n;
void mark(ll c[])
{
    fu(i,1,c[0])
    {
        c[i+1]+=c[i]/10000;
        c[i]%=10000;
    }
    while(c[c[0]+1])
    {
        c[0]++;
        c[c[0]+1]+=c[c[0]]/10000;
        c[c[0]]%=10000;
    }
}
void mul(ll a1,ll a2,ll a3,ll c[])
{
    c[0]=c[1]=1;
    fu(i,1,c[0])
    {
        c[i]*=a1;
    }
    mark(c);
    fu(i,1,c[0])
    {
        c[i]*=a2;
    }
    mark(c);
    fu(i,1,c[0])
    {
        c[i]*=a3;
    }
    mark(c); 
}
void add(ll a[],ll b[],ll c[])
{
    if(a[0]>b[0])
        c[0]=a[0];
    else
        c[0]=b[0];
    fu(i,1,c[0])
        c[i]=a[i]+b[i];
    mark(c);
}
int ce(ll a[],ll b[])
{
    if(a[0]<b[0]) return 0;
    if(a[0]>b[0]) return 1;
    fd(i,a[0],1)
        if(a[i]<b[i])
            return 0;
        else if(a[i]>b[i])
            return 1;
    return 0;
}
void p()
{
    int i;
    cout<<f[1][n][f[1][n][0]];
    fd(i,f[1][n][0]-1,1)
    {
        cout<<f[1][n][i]/1000;
        cout<<f[1][n][i]/100%10;
        cout<<f[1][n][i]/10%10;
        cout<<f[1][n][i]%10;
    }
    cout<<endl;
}
int main()
{
    n=read();
    fu(i,1,n)
        a[i]=read();
    fu(i,1,n)
        fu(j,1,n)
            f[i][j][0]=0;
    fu(jie,2,n-1)
        fu(i,1,n-jie)
        {
            int j=i+jie;
            f[i][j][0]=60;
            fu(k,i+1,j-1)
            {
                memset(s1,0,sizeof(s1));
                memset(s2,0,sizeof(s2));
                memset(s3,0,sizeof(s3));
                mul(a[i],a[k],a[j],s1);
                add(f[i][k],f[k][j],s2);
                add(s1,s2,s3);
                if(ce(f[i][j],s3))
                    memcpy(f[i][j],s3,sizeof(s3));
            }
        }
    p();
    return 0;
}
OPEN

 

优化

倍增思想（注意是思想）

例题：LG-P3147[USACO16OPEN]262144与LG-P3146[USACO16OPEN]248

解析

一维的2048换了个规则：两数合并结果+1。两题的不同在于：数据范围不同；这难道不是调整数组大小就能解决的事么？但是USACO的出题人不会傻到同一场出两道相同的题，

 

于是我们开始做数据范围小的题，直接套板子，若两个区间数值相同就合并，再用个ans取最大值。

#include<bits/stdc++.h>
#define fu(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
int read()
{
    int x=0,k=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')
        k=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        x=(x<<1)+(x<<3)+ch-'0';
        ch=getchar();
    }
    return x*k;
}
int n,f[250][250],ans=-1;
int main()
{
    n=read();
    fu(i,1,n)
    {
        int op;
        op=read();
        if(ans<op) 
           ans=op;
        f[i][i]=op;
    }
    fu(jie,2,n)
    {
        fu(i,1,n-jie+1)
        {
            int j=i+jie-1;
            fu(k,i,j-1)
            {
                if(f[i][k]==f[k+1][j])
                    f[i][j]=max(f[i][j],f[i][k]+1);
                ans=max(ans,f[i][j]);
            }
        }
    }
    cout<<ans;
    return 0;
}
(mod)

 

接下来是n<=262144，咦，怎么跟题目一样，一定有不可告人的秘密，2的18次方等于262144，倍增思想石锤。

 

设 f[ i ][ j ] 表示位置在 j 的数向右合并至 f[ i ][ j ] 得到 i ，那么在 j 处的数合并至 f[ i-1 ][ j ] 得到 i-1 ，便有在 f[ i-1 ][ j ] 处的数合并至 f [ i-1 ][ f[ i-1 ][ j ] ] 得到 i ;

出现了，转移方程！

#include<bits/stdc++.h>
#define fu(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>=b;i--)
using namespace std;
int read()
{
    int x=0,k=1;
    char ch=getchar();
    while(ch<'0'||ch>'9')
    {
        if(ch=='-')
        k=-1;
        ch=getchar();
    }
    while(ch>='0'&&ch<='9')
    {
        x=(x<<1)+(x<<3)+ch-'0';
        ch=getchar();
    }
    return x*k;
}
int n,f[61][270000],ans;
int main()
{
    n=read();
    for(int i=1;i<=n;i++)
    {
        int op;
        op=read();
        f[op][i]=i+1;
    }
    for(int two=2;two<=58;two++)
        for(int t=1;t<=n;t++)
        {
            if(f[two][t]==0)
                f[two][t]=f[two-1][f[two-1][t]];
            if(f[two][t]!=0)
                ans=two;
            //cout<<ans<<endl;
        }
    cout<<ans;
    return 0;
}
INTERESTING

 

转载于:https://www.cnblogs.com/HHHG/p/11203890.html

weixin_30768661
关注
凸多边形的划分（区间dp，高精度）
Ray.C.L的博客
09-29 953
思路：dp[i][j]表示从顶点i到j划分成N-2个互不相交的三角形权值最小。 为了使三角形不相交，我们枚举中间点K，和端点那么dp[i][j]=min(dp[i][j],dp[i][k]+dp[k][r]+w[i]*w[k]*w[r])，因为数据范围过大要用高精度， #pragma GCC optimize(2) #include <cstdio> #include <cstring> #include <algorithm> #include <set&gt.
区间dp
iamcht的博客
11-15 374
区间dp对区间dp的简单介绍：P1880 [NOI1995]石子合并：分析：完整代码如下：明天更新。。。 对区间dp的简单介绍：   区间dp也属于线性dp中的一种，它以“区间长度”作为dp的“阶段”，使用两个坐标（区间的左、右端点）描述每个维度。在区间dp中，一个状态由若干个比它更小且包含于它的区间所代表的的状态转移而来，因此区间dp的决策往往就是划分区间的方法。区间dp的初态一般就又长度为1的“元区间”构成。 P1880 [NOI1995]石子合并：   洛谷：P1880 石子合并 题目描述   在一个
区间DP-----------------凸多边形的划分
qq_45772483的博客
05-17 441
给定一个具有 NN 个顶点的凸多边形，将顶点从 11 至 NN 标号，每个顶点的权值都是一个正整数。 将这个凸多边形划分成 N−2N−2 个互不相交的三角形，对于每个三角形，其三个顶点的权值相乘都可得到一个权值乘积，试求所有三角形的顶点权值乘积之和至少为多少。 输入格式 第一行包含整数 NN，表示顶点数量。 第二行包含 NN 个整数，依次为顶点 11 至顶点 NN 的权值。 输出格式 输出仅一行，为所有三角形的顶点权值乘积之和的最小值。 数据范围 N≤50N≤50, 数据保证所有顶点的权值都小于109109
多边形 （区间dp）
qq_44866969的博客
12-14 159
#include<bits/stdc++.h> using namespace std; const int N=55; int a[2*N]; char b[2*N]; int fmx[2*N][2*N]; int fmn[2*N][2*N]; int main(){ int n;scanf("%d",&n); for(int i=1;i<=2*n;i++){ if(i&1){ char c;cin>>.
Loj10149多边形划分（区间DP+高精度计算）
qq_44647282的博客
03-25 291
题目描述 给定一个具有 N个顶点的凸多边形，将顶点从 1至 N标号，每个顶点的权值都是一个正整数。将这个凸多边形划分成 N-2个互不相交的三角形，试求这些三角形顶点的权值乘积和至少为多少。 输入格式 输入第一行为顶点数 第二行依次为顶点 1至顶点 N的权值。 输出格式 输出仅一行，为这些三角形顶点的权值乘积和的最小值。 题解 将一个n边形划分为n-2个三角形，我们可以看成由两个分别具有k个顶点和n-k个顶点的已经划分完成的多边形再划分一次得到。此时我们只要枚举k，使得k取到最优位置。显然，我们也可以用这种方
凸多边形的划分 区间dp
行码棋
07-29 222
原题链接 状态表示：f[i][j][k]f[i][j][k]f[i][j][k]表示选择以第i个点开头，以第j个点结尾的围成凸多边形的最小值，k是为了写高精度对第三维做的扩充，表示数的每一位。 转移方程：f[i][j]=min(f[i][j],f[i][k]+f[k][j]+w[i]∗w[j]∗w[k])f[i][j] = min(f[i][j], f[i][k] + f[k][j] + w[i] * w[j] * w[k])f[i][j]=min(f[i][j],f[i][k]+f[k][j]+w[i]∗
CodeForces – 1312E Array Shrinking(区间dp)
01-03
题目分析：最优性问题，且是对区间操作的，而且数据范围满足 n^3 的时间复杂度，综上可以考虑区间dp，因为题目已经明确了需要求什么，所以我们不妨设 dp[ i ][ j ] 为区间 [ i , j ] 合并后的最短数列的长度，因为...
【AcWing】线性DP、区间DP、计数类DP
m0_70948728的博客
07-25 473
线性DP、区间DP、计数类DP
3区间DP1
08-08
在这个题目中，我们通过C++语言实现了一个基础的3区间DP求解器。 首先，我们需要理解动态规划的核心思想：将复杂问题分解为更小的子问题，并存储子问题的解，避免重复计算。3区间DP在此基础上增加了对3个连续子数组...
区间DP概率DP树形DP插头DP
09-05
区间DP概率DP树形DP插头DP，每种DP一道典型例题，有助于初学者
AcWing 1069. 凸多边形的划分（区间 dp）
做一只熊猫不好吗？
11-25 465
代码 #include <bits/stdc++.h> using namespace std; #define db double #define ll long long #define Pir pair<int, int> #define fi first #define se second #define pb push_back #define m_p make_pair #define inf 0x3f3f3f3f #define INF 0x3f3f3f3f.
经典问题四. 【区间dp】 凸多边形最优三角形划分
学无纸巾
11-02 1634
（区间dp） 凸多边形最优三角形划分 问题描述： 思路： 将凸多边形的点数组化。发现三角形划分是满足区间叠加的。 dp[i][j]，1<=i<=j<=N,代表凸子多边形{vi-1,vi,…,vj}的最优三角剖分所对应的权函数值，即其最优值。#include <cstdio> #include <cstring> #include <algorithm> using namespace std
划分凸多边形 区间DP + 高精度运算
正月看雪花的博客
04-28 206
思路 对于每个三角形， 我们的len从3开始枚举， f[l, r]代表区间[l, r]的最小值， 然后我们就从一般的区间DP角度去做， 就可以得到样例的结果 #include <iostream> #include <cstring> #include <algorithm> using namespace std; const int N = 100...
凸多边形的划分(区间dp,高精,记忆化搜索实现)
weixin_45590210的博客
03-03 285
题目链接 区间dp练习题,不过这题除了区间dp外…还有高精度的问题,先放一放…不用高精度的话在oj上面只能拿50分. 因为写区间dp我都习惯了写记忆化搜索.可能会不太好…毕竟可能会爆栈.速度也相对慢一点. 思路分析:由于是自顶向下的方法,所以我一般都先从最终状态出发,来把主框架先写好.最终状态是由一个一个的小三角形所组成的,区间范围是[1,n].那么在大区间之中选一个点和大区间的头尾组成一个三角形...
AcWing1069.凸多边形的划分(区间DP)题解
逐梦er的博客
06-07 3630
凸多边形的划分 题目传送门 题目描述 给定一个具有 N 个顶点的凸多边形，将顶点从 1 至 N 标号，每个顶点的权值都是一个正整数。 将这个凸多边形划分成 N−2 个互不相交的三角形，对于每个三角形，其三个顶点的权值相乘都可得到一个权值乘积，试求所有三角形的顶点权值乘积之和至少为多少。 输入格式 第一行包含整数 N，表示顶点数量。 第二行包含 N 个整数，依次为顶点 1 至顶点 N 的权值。 输出格式 输出仅一行，为所有三角形的顶点权值乘积之和的最小值。 数据范围 N≤50 数据保证所有顶点的权值都小于1
1069. 凸多边形的划分（思维，区间dp）
qq12323qweeqwe的博客
02-15 571
限制条件： 分为N-2个三角形 互不相交 分析： 由于一个三角形是需要三个顶点，我们可以先枚举三个顶点，然后观察一下有什么性质 如图，当我们枚举出i，j，k三个顶点的时候，当前可以划分为三个集合 左边的多边形集合+右边的多边形集合+当前的三角形 这样子看就和能量项链很像了，通过枚举三个点，然后计算集合，由于需要先计算出所有的单个三角形的集合 然后 再递推到一整个多边形中的权值（枚举两个端点，至少中间间隔一个点，所以len为3，才能保证中间间隔至少一个点） 所以len至多从3开始枚举...
AcWing 1069. 凸多边形的划分(区间dp+高精度)
qq_43543086的博客
07-08 189
const int N = 55; ll w[N], f[N][N][N],n; void add(ll a[], ll b[]) { static ll c[N]; memset(c, 0, sizeof c); for (int i = 0, t = 0;i < N;i++) { t += b[i] + a[i]; c[i] = t % 10; t = t / 10; } memcpy(a, c, sizeof c); } void pri(ll a[]) { int i
凸多边形的划分
qq_61185410的博客
05-31 949
题目： 给定一个具有 NN 个顶点的凸多边形，将顶点从 11 至 NN 标号，每个顶点的权值都是一个正整数。 将这个凸多边形划分成 N−2N−2 个互不相交的三角形，对于每个三角形，其三个顶点的权值相乘都可得到一个权值乘积，试求所有三角形的顶点权值乘积之和至少为多少。 输入格式 第一行包含整数 NN，表示顶点数量。 第二行包含 NN 个整数，依次为顶点 11 至顶点 NN 的权值。 输出格式 输出仅一行，为所有三角形的顶点权值乘积之和的最小值。 数据范围 N≤50N≤50, 数据保证所有顶
codeforces 区间dp
最新发布
10-04
区间DP是一种动态规划的方法，用于解决区间范围内的问题。在Codeforces竞赛中，区间DP经常被用于解决一些复杂的字符串或序列相关的问题。 在区间DP中，dp[i][j]表示第一个序列前i个元素和第二个序列前j个元素的最优...
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值