话说今天百度面试,可能是由于我表现的不太好,面试官显得有点不耐烦,说话的语气也很具有嘲讽的意思,搞得我有点不爽。Whatever,面试中有问到整数划分问题,回答这个问题过程中被面试官搞的不胜其烦,最后也给出了其动态规划的算法,但是显然,醉翁之意不在动态规划而在于生成函数(generating function)。下面开始吧:

 

参考:http://www.skymoon.biz/?p=192  (问题定义以及动态规划)

         http://www.artofproblemsolving.com/Resources/Papers/LaurendiPartitions.pdf (生成函数)

         http://www.cnblogs.com/hoodlum1980/archive/2008/10/11/1308493.html(部分代码实现)

 

本系列有两篇文章,会分为下面三个部分,首先我会简单介绍一下动态规划的解法,然后会介绍生成函数的解法,第三部分(在下一篇文章中说明)做一下生成函数的扩展。要求读者具备基本的动态规划常识以及生成函数概念。

 

文章最后有我最喜欢的东西。

 

先说明一下问题,什么是整数划分?

  • n=m1+m2+...+mi; (其中mi为正整数,并且1 <= mi <= n),则{m1,m2,...,mi}为n的一个划分。
  • 如果{m1,m2,...,mi}中的最大值不超过m,即max(m1,m2,...,mi)<=m,则称它属于n的一个m划分。这里我们记n的m划分的个数为f(n,m);
  • 举个例子,当n=5时我们可以获得以下这几种划分(注意,例子中m>=5)

5 = 5
   = 4 + 1
   = 3 + 2
   = 3 + 1 + 1
   = 2 + 2 + 1
   = 2 + 1 + 1 + 1
   = 1 + 1 + 1 + 1 + 1

 

一、 动态规划解法

根据n和m的关系,考虑以下几种情况:
1. 当n=1时,不论m的值为多少(m>0),只有一种划分即{1};
2. 当m=1时,不论n的值为多少,只有一种划分即n个1,{1,1,1,...,1};
3. 当n=m时,根据划分中是否包含n,可以分为两种情况:
    (1) 划分中包含n的情况,只有一个即{n};
    (2) 划分中不包含n的情况,这时划分中最大的数字也一定比n小,即n的所有(n-1)划分。因此 f(n,n) =1 + f(n,n-1);
4. 当n<m时,由于划分中不可能出现负数,因此就相当于f(n,n);
5. 但n>m时,根据划分中是否包含最大值m,可以分为两种情况:
    (1) 划分中包含m的情况,即{m, {x1,x2,...xi}}, 其中{x1,x2,... xi} 的和为n-m,可能再次出现m,因此是(n-m)的m划分,因此这种划分个数为f(n-m, m);
    (2) 划分中不包含m的情况,则划分中所有值都比m小,即n的(m-1)划分,个数为f(n,m-1);因此 f(n, m) = f(n-m, m)+f(n,m-1);

 

综合以上情况,我们可以看出,上面的结论具有递归定义特征,其中(1)和(2)属于回归条件,(3)和(4)属于特殊情况,将会转换为情况(5)。而情况(5)为通用情况,属于递推的方法,其本质主要是通过减小m以达到回归条件,从而解决问题。其递推表达式如下:

  • f(n, m)= 1; (n=1 or m=1)
  • f(n, m)=f(n, n); (n<m)
  • 1+ f(n, m-1); (n=m)
  • f(n-m,m)+f(n,m-1); (n>m)

 

据此我们获得了动态规划的代码:

   1:  #include<iostream>
   2:   
   3:  using namespacestd;
   4:   
   5:  int equationCount(intn,intm)
   6:  {
   7:      if(n==1||m==1)
   8:          return 1;
   9:      else if(n<m)
  10:          return equationCount(n,n);
  11:      else if(n==m)
  12:          return 1+equationCount(n,n-1);
  13:      else
  14:          return equationCount(n,m-1)+equationCount(n-m,m);
  15:  }
  16:   
  17:  int main(void)
  18:  {
  19:      in tn;
  20:      while(scanf("%d",&n)!=EOF&&(n>=1&&n<=120))
  21:      {
  22:          printf("%d\n",equationCount(n,n));
  23:      }
  24:      return 0;
  25:  }

.csharpcode, .csharpcode pre
{
font-size: small;
color: black;
font-family: consolas, "Courier New", courier, monospace;
background-color: #ffffff;
/*white-space: pre;*/
}
.csharpcode pre { margin: 0em; }
.csharpcode .rem { color: #008000; }
.csharpcode .kwrd { color: #0000ff; }
.csharpcode .str { color: #006080; }
.csharpcode .op { color: #0000c0; }
.csharpcode .preproc { color: #cc6633; }
.csharpcode .asp { background-color: #ffff00; }
.csharpcode .html { color: #800000; }
.csharpcode .attr { color: #ff0000; }
.csharpcode .alt
{
background-color: #f4f4f4;
width: 100%;
margin: 0em;
}
.csharpcode .lnum { color: #606060; }

 

几个变种:

(一)要求1,2,3,4..,m中每个数只允许使用一次的时?

此时我们需要调整我们的状态转换公式。

f(n-m,m)+f(n,m-1); (n>m) 应该更改为:f(n-m,m-1)+f(n,m-1); (n>m)

为什呢?因为每个数最多使用一次,f(n-m,m-1)表示我们取了数m,f(n,m-1)表示我们没取,但是无论取不取数m我们以后都不会再次取数m了。

当然喽,我们还需要调整边界状态:当m=1时,f(n,m)=1;当n=1而m>1时,f(n,m)=0。

其他不变!

 

(二)要求只能取1,2,3,4,..,m中的奇数?(默认m为奇数,如果不是则m=m-1)

这个呢,我们首先需要调整边界状态:当m=1时,f(n,m)=1;当n=1而m>1时,f(n,m)=0

其次,我们需要调整状态转换公式:

f(n-m,m)+f(n,m-1); (n>m) 应该更改为:f(n-m,m)+f(n,m-2); (n>m)

这是因为我们不能取偶数,故而当m为奇数的时候,m-1为偶数(只能被选择0次),f(n,m-1)=f(n,m-2);

(三)要求我们所取的 (n=m1+m2+...+mi )中  m1 m2 ... mi连续,比如5=1+4就不符合要求了。

这个的话,需要做一下转换,留待下一篇文章说明

注意,一般而言动态规划算法是用非递归从下往上计算的,上述代码采用递归形式只是为了便于理解,真正实现的话最好采用非递归形式。

 

二、 生成函数解法

先简单说明一下生成函数吧,下面是一个生成函数,

xk的系数ak代表了可获得数字k的组合数。

 

那么回到我们的问题中,我们怎么用生成函数去解决呢?

类似的,我们可以计算生成函数:

计算出来的xk的系数ak就是我们需要的划分数了。我们来详细说明一下为什么是这样。考虑一下x3,我们可以在第一个括号里面选择x、在第二个括号里面选择x2、其他括号选择1,这种方式来得到;同样的,我们也可以通过在第三个括号里面选择x3,其他扩号里面选择1,这种方式来获得。那么这跟生成函数有什么关系呢?

我们来做一个说明,第i个括号(1+xi+x2i+x3i · · ·)选择的元素代表了数字i在我们最终的划分中出现的次数,具体而言,如果我们在第i个括号中选择了元素 xc(i) * i  则表示数字 i 在我们最终的划分中出现了c(i)次。如果我们把最终从每个括号里面选择出来的元素相乘 x1 * c(1) · x2 * c(2) · x3 * c(3) · · · = x c(1) + 2 * c(2) + 3 * c(3)···.那么,xn的系数就是我们可以获得多少种不同的方式使得c(1) + 2 * c(2) + 3 * c(3)··· = n,也就是n的划分数(其中,c(i)代表了在一次划分中数字 i 的出现次数)。比如说25=6+4+4+3+2+2+2+1+1,用上式表示就是25=1(2)+2(3)+3(1)+4(2)+5(0)+6(1),也就是在划分中有两个1,三个2,一个3,两个4,0个5以及一个6。

假设x<1,那么我们可以将上面的生成函数表示为:

这将给我们第三部分的扩展带来方便。

考虑到我们能选择的最大的数 i 是 m ,所以真正计算的时候我们需要对上面的生成函数式子做一下修改:

G(x) = (1+x+x^2+x^3+...+x^n) (1+x^2+x^4+...) (1+x^3+x^6+...) ... (1+x^m)

       = g(x,1) g(x,2) g(x,3) ... g(x, n)

 

代码实现如下:

   1:  #define N 130
   2:  unsigned long a[N];/*多项式a的系数数组*/
   3:  unsigned long b[N];/*多项式b的系数数组*/
   4:  unsigned long c[N];/*存储多项式a*b的结果*/
   5:   
   6:  /*两个多项式进行乘法,系数分别在a和b中,结果保存到c ,项最大次数到N */
   7:  /*注意这里我们只需要计算到前N项就够了。*/
   8:  void Poly()
   9:  {
  10:      int i,j;
  11:      memset(c,0,sizeof(c));
  12:      for(i=0; i<N; i++)
  13:              for(j=0; j<N-i; j++) /*y<N-i: 确保i+j不会越界*/
  14:                    c[i+j] += a[i]*b[j];
  15:  }
  16:   
  17:  /*计算出前N项系数!即g(x,1) g(x,2)... g(x,n)的展开结果*/
  18:  void Init()
  19:  {
  20:      int i,k;
  21:      memset(a,0,sizeof(a));
  22:      memset(c,0,sizeof(c));
  23:      for(i=0;i<N;i++) a[i]=1; /*第一个多项式:g(x, 1) = x^0 + x^1 + x^2 + x^3 +  */
  24:      for(k=2;k<N;k++)
  25:      {
  26:          memset(b,0,sizeof(b));
  27:          for(i=0;i<N;i+=k) b[i]=1;/*第k个多项式:g(x, k) = x^0 + x^(k) + x^(2k) + x^(3k) +  */
  28:          Poly(); /* 多项式乘法:c= a*b */
  29:          memcpy(a,c,sizeof(c)); /*把相乘的结果从c复制到a中:c=a; */
  30:      }
  31:  }

.csharpcode, .csharpcode pre
{
font-size: small;
color: black;
font-family: consolas, "Courier New", courier, monospace;
background-color: #ffffff;
/*white-space: pre;*/
}
.csharpcode pre { margin: 0em; }
.csharpcode .rem { color: #008000; }
.csharpcode .kwrd { color: #0000ff; }
.csharpcode .str { color: #006080; }
.csharpcode .op { color: #0000c0; }
.csharpcode .preproc { color: #cc6633; }
.csharpcode .asp { background-color: #ffff00; }
.csharpcode .html { color: #800000; }
.csharpcode .attr { color: #ff0000; }
.csharpcode .alt
{
background-color: #f4f4f4;
width: 100%;
margin: 0em;
}
.csharpcode .lnum { color: #606060; }

 

好了,我们来看几个变种:

(一)要求1,2,3,4..,m中每个数只允许使用一次的时?

这个简单,将生成函数式子改为

G(x) = (1+x) (1+x^2) (1+x^3) ... (1+x^m) = g(x,1) g(x,2) g(x,3) ... g(x, n)就可以了

(二)要求只能取1,2,3,4,..,m中的奇数?(默认m为奇数,如果不是则m=m-1)

这个也简单,调整生成函数的式子为:

G(x) = (1+x+x^2+x^3+...+x^n)(1+x^3+x^6+...) ... (1+x^m)

       = g(x,1) g(x,2) g(x,3) ... g(x, n)

(三)要求我们所取的 (n=m1+m2+...+mi )中  m1 m2 ... mi连续,比如5=1+4就不符合要求了。

这个呢,我们需要做一点转换

 

 

最后我们来证明一下变种一和变种二是一样的:

整数划分 Integer Partition(一)的更多相关文章

  1. 整数划分 Integer Partition(二)

    本文是整数划分的第二节,主要介绍整数划分的一些性质. 一 先来弥补一下上一篇文章的遗留问题:要求我们所取的 (n=m1+m2+...+mi )中  m1 m2 ... mi连续,比如5=1+4就不符合 ...

  2. 51nod p1201 整数划分

    1201 整数划分 基准时间限制:1 秒 空间限制:131072 KB 分值: 80 难度:5级算法题 将N分为若干个不同整数的和,有多少种不同的划分方式,例如:n = 6,{6} {1,5} {2, ...

  3. 2014北大研究生推免机试(校内)-复杂的整数划分(DP进阶)

    这是一道典型的整数划分题目,适合正在研究动态规划的同学练练手,但是和上一个随笔一样,我是在Coursera中评测通过的,没有找到适合的OJ有这一道题(找到的ACMer拜托告诉一声~),这道题考察得较全 ...

  4. 整数划分 (区间DP)

    整数划分(四) 时间限制:1000 ms  |  内存限制:65535 KB 难度:3   描述 暑假来了,hrdv 又要留学校在参加ACM集训了,集训的生活非常Happy(ps:你懂得),可是他最近 ...

  5. nyoj 90 整数划分

    点击打开链接 整数划分 时间限制:3000 ms  |  内存限制:65535 KB 难度:3 描述 将正整数n表示成一系列正整数之和:n=n1+n2+-+nk,  其中n1≥n2≥-≥nk≥1,k≥ ...

  6. 51nod1201 整数划分

    01背包显然超时.然后就是一道神dp了.dp[i][j]表示j个数组成i的方案数.O(nsqrt(n)) #include<cstdio> #include<cstring> ...

  7. NYOJ-571 整数划分(三)

    此题是个非常经典的题目,这个题目包含了整数划分(一)和整数划分(二)的所有情形,而且还增加了其它的情形,主要是用递归或者说是递推式来解,只要找到了递推式剩下的任务就是找边界条件了,我觉得边界也是非常重 ...

  8. BZOJ1263: [SCOI2006]整数划分

    1263: [SCOI2006]整数划分 Time Limit: 1 Sec  Memory Limit: 162 MBSubmit: 677  Solved: 332[Submit][Status] ...

  9. BZOJ 1263: [SCOI2006]整数划分( 高精度 )

    yy一下发现好像越小越好...分解成3*3*3*3……这种形式是最好的...然后就是高精度了 ----------------------------------------------------- ...

随机推荐

  1. 3.Windows Server 2012 R2数据库部署

    很多人竟然不会安装数据库....好吧,来个图文教程,其实和windows里面一样安装,和安装2008一样的 先安装3.5:http://www.cnblogs.com/dunitian/p/53487 ...

  2. ButterKnife Zelezny从配置到使用

    插件介绍:ButterKnife是一个专注于Android系统的View注入框架,可以减少大量的findViewById以及setOnClickListener代码,可视化一键生成.又一神器,完美告别 ...

  3. 八款常见的Android游戏引擎

    原文地址:http://bbs.csdn.net/topics/380203732 1.Angle  Angle是一款专为Android平台设计的,敏捷且适合快速开发的2D游戏引擎,基于OpenGL  ...

  4. MongoDB的基础知识

    本人只是软件开发的一个菜鸟,在学习MongoDB,总结了一点自己学习的知识,监督自己学习. 如果文章中有不足的地方,还请大神指点迷津,纠正改错,谢谢. 一.MongoDB简介 MongoDB是一个基于 ...

  5. CF#335 Freelancer&#39;s Dreams

    Freelancer's Dreams time limit per test 2 seconds memory limit per test 256 megabytes input standard ...

  6. 编程之美--2. Trie树 (Trie图)

    #1014 : Trie树 时间限制:10000ms 单点时限:1000ms 内存限制:256MB 描述 小Hi和小Ho是一对好朋友,出生在信息化社会的他们对编程产生了莫大的兴趣,他们约定好互相帮助, ...

  7. Eclipse tomcat插件禁用热加载

    Eclipse中的tomcat插件默认是开户了热加载,只要是修改了java文件一保存,tomcat自动编译.加载.发布,很吃内存. 关闭方法: 打开eclipse,找到server项: 双击打开,修改 ...

  8. “康园圈--互联网+校园平台“项目之Sprint3

    sprint3任务列表 1.部署项目内测版上线,并绑定相应域名. 2.添加测试用的模拟数据. 3.写已实现功能及效果汇报文档. 4.写项目总结博客. 5.记录并分析用户反馈. 看板动态 第一次 第二次 ...

  9. Android学习一(入门)

    一.Android 系统简介 1.1.1G-4G 1G:模拟制式手机,1995年问世的第一代模拟制式手机,只能进行语音通话, 2G:手机使用GSM,CDMA(9K/s),增加了接收数据的功能 2.5G ...

  10. OAF中的面包屑(breadcrumbs)始无法显示(转)

    原文地址:OAF中的面包屑(breadcrumbs)始无法显示 OAF中面包屑是一种重要的导航工具.可以提示用户抵达当前页面的路径,也可以方便的切换到之前的节点. 开始做面包屑开发的时候发现面包屑总是 ...