整理:整数划分问题

http://blog.csdn.net/majestyhao/article/details/7079314

http://www.kuqin.com/algorithm/20080511/8342.html

整数划分问题是将一个正整数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  <   1 return   0 ;
      
if (n  ==   1   ||  m  ==   1 return   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( 12 12 ));
    
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  <=   0 break ;
        
if ((n  -  t1)  %  i  ==   0 )
        {
            printf(
" %d  " , x);
            
for (j  =   1 ; j  <  i; j ++ )
                printf(
" %d  " , x  +  j);
            printf(
"   " );
            m
++ ;
        }
    }
    
return  m;
}




输出所有划分

[cpp]  view plain cop
  1. /* 
  2. Subject:计算机算法设计与分析 
  3. Title:例2-5 整数划分问题:输出一个整数的所有划分并统计总划分数 
  4. Coder:Hao 
  5. Class:计科0906 
  6. Num:0304090614 
  7. Date:Sept 20th,2011  
  8. */  
  9.   
  10. #include <iostream>  
  11. using namespace std;  
  12.   
  13. //用于打印(输出)的函数  
  14. //result为存储某划分结果的数组,length为此划分所占长度-1(从0开始)  
  15. void display(int *result,int length)         
  16. {     
  17.     for(int i=0;i<length;i++)  
  18.         cout<<result[i]<<" ";  
  19.     cout<<endl;  
  20. }  
  21.   
  22. //主划分函数q(int n,int m,int *result,int length)  
  23. //算法参考了书例2-5并为了打印各划分做了一定修改  
  24. //n为待划分整数,m为最大加数上限,result和length意以同display函数同//名变量定义  
  25. //将q分成五种情况分类讨论,其中包含递归调用  
  26. int q(int n,int m,int *result,int length)  
  27. {  
  28. //当n>=1并且m=1时,q(n,m,result,length)=q(n-1,m,result,length)      
  29.      if(n>=0&&m==1)                                                                                       
  30.     {  
  31. //直至n=0并且m=1时,输出  
  32.         if(n==0)                                                            display(result,length);  
  33.     }  
  34.         else  
  35.         {  
  36.         result[length]=1;  
  37.         q(n-1,m,result,length+1);  
  38.         }  
  39.         return 1;  
  40.     }  
  41.   
  42.   
  43. // 当 n=1并且m>1 时,分解已经完成,进行输出  
  44.     else if(n==1&&m>1)                                 
  45.     {  
  46.         result[length]=n;  
  47.         display(result,length+1);  
  48.         return 1;  
  49.     }  
  50.   
  51. //当n<m时,q(n,m,result,length)=q(n,n,result,length)  
  52.     else if(n<m)                                        
  53.    {  
  54.         return q(n,n,result,length);  
  55.     }  
  56.   
  57. //当n=m时,q(n,m,result,length)=q(n,m-1,result,length)+1(划分数目)  
  58.     else if(n==m)                                           
  59. {  
  60.         result[length]=m;  
  61.         display(result,length+1);  
  62.         return q(n,m-1,result,length)+1;  
  63.     }  
  64.   
  65. //当n>m>1时,//q(n,m,result,length)=q(n-m,m,result,length+1)+q(n,m-1,result,length)  
  66.     else                                                           
  67.     {  
  68.         result[length]=m;  
  69.         return q(n-m,m,result,length+1)+q(n,m-1,result,length);  
  70.     }  
  71. }  
  72.   
  73.   
  74.   void main()  
  75.   {  
  76.       int n;                            //定义待划分整数  
  77.       int result[100]={0},length=0;        //初始化  
  78.       cout<<"please input the integer:";    
  79.       cin>>n;  
  80.       cout<<"整数"<<n<<"的划分个数为"<<q(n,n,result,length)<<endl;  
  81.   }  

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值