整数划分算法原理与实现

本文为原创,如需转载,请注明作者和出处,谢谢!

    整数划分问题是将一个正整数n拆成一组数连加并等于n的形式,且这组数中的最大加数不大于n。
    如6的整数划分为
   
    6
    5 + 1
    4 + 2, 4 + 1 + 1
    3 + 3, 3 + 2 + 1, 3 + 1 + 1 + 1
    2 + 2 + 2, 2 + 2 + 1 + 1, 2 + 1 + 1 + 1 + 1
    1 + 1 + 1 + 1 + 1 + 1
   
    共11种。下面介绍一种通过递归方法得到一个正整数的划分数。
   
    递归函数的声明为 int split(int n, int m);其中n为要划分的正整数,m是划分中的最大加数(当m > n时,最大加数为n),
    1 当n = 1或m = 1时,split的值为1,可根据上例看出,只有一个划分1 或 1 + 1 + 1 + 1 + 1 + 1
    可用程序表示为if(n == 1 || m == 1) return 1;
   
    2 下面看一看m 和 n的关系。它们有三种关系
    (1) m > n
    在整数划分中实际上最大加数不能大于n,因此在这种情况可以等价为split(n, n);
    可用程序表示为if(m > n) return split(n, n);   
    (2) m = n
    这种情况可用递归表示为split(n, m - 1) + 1,从以上例子中可以看出,就是最大加
    数为6和小于6的划分之和
    用程序表示为if(m == n) return (split(n, m - 1) + 1);
    (3) m < n
    这是最一般的情况,在划分的大多数时都是这种情况。
    从上例可以看出,设m = 4,那split(6, 4)的值是最大加数小于4划分数和整数2的划分数的和。
    因此,split(n, m)可表示为split(n, m - 1) + split(n - m, m)
   
    根据以上描述,可得源程序如下:
  
#include <stdio.h>

   
int split(int n, int m)
   {
      
if(n < 1 || m < 1return 0;
      
if(n == 1 || m == 1return 1;
      
if(n < m) return split(n, n);
      
if(n == m) return (split(n, m - 1+ 1);
      
if(n > m) return (split(n, m - 1+ split((n - m), m));
  }

int main()
{
     printf(
"12的划分数: %d", split(1212));
    
return 0;
}

将正整数划分成连续的正整数之和
如15可以划分成4种连续整数相加的形式:
15
7 8
4 5 6
1 2 3 4 5

    首先考虑一般的形式,设n为被划分的正整数,x为划分后最小的整数,如果n有一种划分,那么
结果就是x,如果有两种划分,就是x和x x + 1, 如果有m种划分,就是 x 、x x + 1 、 x x + 1 x + 2 、... 、x x + 1 x + 2 ... x + m - 1
将每一个结果相加得到一个公式(i * x + i * (i - 1) / 2) = n,i为当前划分后相加的正整数个数。
满足条件的划分就是使x为正整数的所有情况。
如上例,当i = 1时,即划分成一个正整数时,x = 15, 当i = 2时, x = 7。
当x = 3时,x = 4, 当x = 4时,4/9,不是正整数,因此,15不可能划分成4个正整数相加。
当x = 5时,x = 1。

    这里还有一个问题,这个i的最大值是多少?不过有一点可以肯定,它一定比n小。我们可以做一个假设,
假设n可以拆成最小值为1的划分,如上例中的1 2 3 4 5。这是n的最大数目的划分。如果不满足这个假设,
那么 i 一定比这个划分中的正整数个数小。因此可以得到这样一个公式i * (i + 1) / 2 <= n,即当i满足
这个公式时n才可能被划分。

综合上述,源程序如下
int split1(int n)
{
    
int i, j, m = 0, x, t1, t2;
   
// 在这里i + 1之所以变为i - 1,是因为i * (i - 1) / 2这个式子在下面多次用到,
  
// 为了避免重复计算,因此将这个值计算完后保存在t1中。并且将<= 号变为了<号。
    for(i = 1; (t1 = i * (i - 1/ 2< n; i++
    {
        t2 
= (n - t1);
        x 
=  t2 / i;
        
if(x <= 0break;
        
if((n - t1) % i == 0)
        {
            printf(
"%d ", x);
            
for(j = 1; j < i; j++)
                printf(
"%d ", x + j);
            printf(
"\n");
            m
++;
        }
    }
    
return m;
}

posted on 2008-05-11 15:39 银河使者 阅读(1322) 评论(5)  编辑 收藏 所属分类: algorithm 原创

评论

#1楼  2008-05-11 16:09 jfan [未注册用户]

呵呵,软考里就有这样的算法.   回复  引用    

#2楼  2008-05-11 20:53 石桥头 [未注册用户]

以前读书时候考数据结构和算法都不及格   回复  引用    

#3楼 [楼主] 2008-05-11 21:02 银河使者      

哈哈,算法是程序设计的基础,锻练思维能力。   回复  引用  查看    

#4楼  2008-05-12 04:36 hoodlum1980      

我希望楼主能严肃对待“原创”一词,你确定这也是你的原创吗?难道我们以前看到的算法是现在才被你发现的吗。   回复  引用  查看    

#5楼 [楼主] 2008-05-12 09:02 银河使者      

基本的算法思路都是一样的,不过程序和文章是我自己写的。其他很多原创的文章,谁能保证自己的文章的思路是世界第一个。恐怕没人能保证。所谓“原创”就是指程序和文章都是通过自己的理解写出的,而不是直接copy and paste。 至于某些思想和解释同书上或网上的内容有类似,也是正常的。

如果“原创”的文章的思想非常要以前从来没人提出过,恐怕“原创”文章会非常少了。哈哈!!   回复  引用  查看    


标题  
姓名  
主页
Email (只有博主才能看到) 
验证码 *  看不清,换一张 [登录][注册]
内容(请不要发表任何与政治相关的内容)  
  登录  使用高级评论  新用户注册  返回页首  恢复上次提交      
该文被作者在 2008-05-22 17:10 编辑过
 
 


<2008年5月>
27282930123
45678910
11121314151617
18192021222324
25262728293031
1234567

导航

统计

与我联系

搜索

 

常用链接

留言簿(4)

我参与的团队

我的标签

随笔分类(108)

随笔档案(52)

相册

专家Blog

积分与排名

最新评论

阅读排行榜

评论排行榜

60天内阅读排行