码迷,mamicode.com
首页 > 其他好文 > 详细

[题解+总结]动态规划大合集II

时间:2015-08-07 00:03:04      阅读:138      评论:0      收藏:0      [点我收藏+]

标签:

1、前言

        大合集总共14道题,出自江哥之手(这就没什么好戏了),做得让人花枝乱颤。虽说大部分是NOIP难度,也有简单的几道题目,但是还是做的很辛苦,有几道题几乎没思路,下面一道道边看边分析一下。

 

2、lis 最长上升子序列

        唯一一道裸题,但是O(n^2)过不了,临时看了看O(n log n)的二分做法和线段树做法。先来讲讲简单的二分做法,其本质就是在O(n^2)上进行优化,需要证明一个结论。设当前处理数列第k位,存在:

  (1)a[i]<a[j]<a[k];

  (2)i<j<k;

  (3)f[i]=f[j]。

  此时我们选择a[i]好还是a[j]好?显然是a[i]好。为什么?倘若存在a[i]<a[x]<a[j],i<x<j,则选择a[i]可以得到更长的子序列——a[i],a[x],a[k]……这样,对于f[t]=k,我们只需要保存a[t]值最小的那个t了,可以证明得f数组将会单调递增,故我们每次不需要向前每一位都进行查找,而可以进行二分查找,故时间复杂度为O(n log n)。

 

代码:

------------------------------------------------------------------------------------------------------

 

#include<cstdio>

int min(int a,int b) { return (a<b)?a:b; }

int v[200010],n,x,ans;

 

 

int main()
{
  freopen("lis.in","r",stdin);
  freopen("lis.out","w",stdout);
  scanf("%d", &n);
  for (int i=1;i<=n;i++)
  {
    scanf("%d",&x);
    int l=1,r=ans,mid=(l+r)>>1,loc=0;
    while (l<=r)
    {
      if (v[mid]<=x) loc=mid,l=mid+1;
      else r=mid-1;
      mid=(l+r)>>1;
    }
    if (loc==ans) v[++ans]=x;
    else v[loc+1]=min(v[loc+1],x);
  }
  printf("%d\n",ans);

}

------------------------------------------------------------------------------------------------------

 

 

3、chess 黑白棋盘

        根据题目分析,我们需要维护该列是否存在棋子;该行是否存在棋子;由于n<=15,我们可以考虑二进制状态压缩进行转移,最多有2^(n+1)种状态。设f[i][j]为当前第i行,状态为j(压缩的二进制状态 )时的方案数。今天才刚刚学会用&、^、|等位运算符,上次写状压的时候还是手写的。。。方程:f[i][j]=sum{f[i-1][j^(1<<(k-1))]} (k∈[1,n],j&(1<<(k-1)),!a[i][k])

 

代码:

------------------------------------------------------------------------------------------------------

#include <cstdio>
#define MOD 1e9+7

int n,a[20][20],f[20][100000],ans;
char ch;

int main()
{
  freopen("chess.in","r",stdin);
  freopen("chess.out","w",stdout);
  scanf("%d",&n);
  for (int i=1;i<=n;i++)
    for (int j=1;j<=n;j++)
    {
      do ch=getchar(); while ((ch!=‘0‘)&&(ch!=‘1‘));
      a[i][j]=ch-‘0‘;
    }
  f[0][0]=1;
  for (int i=1;i<=n;i++)
    for (int j=0;j<1<<n;j++)
    {
      f[i][j]=f[i-1][j];
      for (int k=1;k<=n;k++)
        if ((j&(1<<(k-1))) && (!a[i][k])) // 如果状态j和
        {
          f[i][j]+=f[i-1][j^(1<<(k-1))];
          if (f[i][j]>=MOD) f[i][j]-=MOD;
        }
    }
  for (int i=0;i<1<<n;i++) { ans+=f[n][i]; if (ans>=MOD) ans-=MOD; }
  printf("%d",ans);
  return 0;
}

------------------------------------------------------------------------------------------------------

 

4、bound 多背包问题

        相信这道题应该考场上没几个人能够做出来啊。。。做出来了也很难证明。。。实在是令人感到意外。题目本身难度似乎并不高,但是我们注意到数据范围<=10^9,这是连O(n)都难以接受的长度。此时可以考虑到将路径进行一些修改以满足我们的范围要求。首先求出1-10里面任意两个数的最小公倍数,然后取最大值max,可以证明当两石块之间的距离大于max的时候,那么大于它的部分的每一个点都可以通过这两个数的某一种组合跳到,所以当两个石块间的距离大于这个最小公倍数,那么就把它们的距离缩小到这个最小公倍数。我仍然不清楚这种说法是否正确,因为网上还有的说可以>=72均压缩,而取最大的最小公倍数只能是>=100压缩。正确性有待考证。用f[i]表示跳到i时最少踩到的石子数。方程:f[i]=min{f[i-k]+x} (k∈[s,t],x={0,1},用以表示i是否有石子)

 

5、NOIP2006 能量项链

        石子合并加强版,由于是一个环,所以我们要考虑将其事先切割一次。有O(n^4)的做法——枚举切割点,在枚举左端点和右端点的基础上枚举中间点。由于n<=100,10^8的范围勉强能过。对于切割点的枚举,能不能考虑优化呢?我们不需要在实际上切割,设置数组时可以将数据复制一次,然后直接进入合并过程。用f[i][j]表示从第i个数到第j个数全部合并的最大得分,则ans=max{f[k][n+k+1]} (k∈[1,n-1]),f[k][n+k+1]表示切割点在k-1到k之间。方程:f[i][j]=max{f[i][k]+f[k+1][j]+a[i]*a[k+1]*a[j+1]} (k∈[i,j-1])

 

6、NOIP2006 金明的预算方案

        背包问题加强版,我没有想到这种做法,太弱了。。。同样地,f[i][j]表示前i个物品用了不超过j元所得到的重要度乘价格的最大值,由于主件存在附件,购买附件的前提是主件要购买,在动态规划中我们不可能进行类似于搜索中的条件判断,这一考虑将这一要求进行压缩——枚举物品时,若该物品为附件,直接跳过;若该物品为主件,则进行三次转移——单独买主件;买主件和一个附件;买主件和两个附件。这样,存在方程:f[i][j]=f[i-1][j-w[i]]+p[i]*v[i]

 

7、NOIP2007 矩阵取数

        n行m列,由于行与行之间没有任何联系,所以我们对于每一行进行动态规划,然后计算总和即可。f[i][j]表示当前行剩余第i个数到第j个数没有被取走时的最大得分,存在取走左边和取走右边两种情况,故f[i][j]=max{f[i][j-1]+a[j]*2^(m-j+i),f[i+1][j]+a[i]*2^(m-i+j)}。由于数据很大,需要用到高精度加法。

 

8、NOIP2008 传纸条

        待补充。

 

9、NOIP2010 乌龟棋

        待补充。

 

10、NOIP2013 花匠

        这道题的正解貌似是贪心?找到拐点就行了。。。但是动态规划同样可行。首先容易想到的是一种O(n^2)的算法。因为题目要求数列保持波浪形,存在两种情况,一高一矮或一矮一高(奇偶性讨论),f[i][0]表示以第i朵花为终点,且其为高花时的最多剩余花数;f[i][1]表示以第i朵花为终点,且其为矮花时的最多剩余花数。故根据波浪形态进行交替转移,可得方程:f[i][0]=max{f[j][1]}+1(j∈[1,i-1],h[i]>h[j]); f[i][1]=max{f[j][0]}+1(j∈[1,i-1],h[i]<h[j])。

        但是,n<=100000,O(n^2)显然过不了。我们能否省去一个维度?看了题解之后豁然开朗,我们在转移的时候根本不需要用到之前的所有状态,只需要考虑前面一朵花与当前花的关系。为什么?根据题意,在判断第i朵花时,在[1,i-1]中符合转移条件中最大值必定在最后,符合局部单调性。故我们只需要关注第i-1朵花的状态。f[i][0]和f[i][1]的含义同上,根据h[i]>h[i-1],hh[i]<h[i-1],要么选,要么不选;对于h[i]=h[i-1],必定不选,O(n)就行了,方程:f[i][0]=max(f[i-1][1]+1,f[i-1][0]),f[i][1]=f[i-1][1](h[i]>h[i-1]);f[i][0]=f[i-1][0],f[i][[1]=f[i-1][1](h[i]=h[i-1]);f[i][0]=f[i-1][0],f[i][1]=max(f[i-1][0]+1,f[i-1][1])(h[i]<h[i-1])。

 

11、NOIP2014 飞翔的小鸟

        当时什么都没有想就直接打了个50分的暴力,其实现在看来还是比较容易看出这是一道动态规划。同样裸DP只有70分,需要单调队列优化。待补充。

 

12、小结

       现在对动态规划的题目做一个总结。NOIP的难度大概就是上面所述了,我也在慢慢脑补着动规的大洞,动态规划许多题目做题的方式存在共同点:首先我们要根据题目的修改操作来设计状态,要求状态一定是没有重叠单独存在的(就是说对于动态规划数组中的每一个元素都只对应着一个状态),同时我们能够通过线性转移来得到答案。

       有些状态很明显,很容易设计,比如“合唱队形”,“乌龟棋”,“能量项链“感觉能够目测大概的模型;但是有些题目模模糊糊,一定要找到最合适的状态。比如“矩阵取数”,原本我设的是f[i][j]为当前共取走i个,其中从左边取走j个时的最大分数。这符合状态的独立性,但是不方便转移(也有可能是我当时没有想到)。

       然而在一切都大功告成之后,一定要验证一下复杂度是否正确。虽然大多数题目是提供部分分的,但是如果动态规划能够做到AC的话,当然是要去考虑优化的,如“花匠”。类似于最长上升子序列就存在O(n^2)和O(n log n)两种做法。优化的话,可以直接用朴素的方式;同样斜率优化,单调队列也是必要的。

[题解+总结]动态规划大合集II

标签:

原文地址:http://www.cnblogs.com/jinkun113/p/4709397.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!