数学:整数拆分

使用五边形数定理进行整数拆分

先介绍一个叫做分拆数的概念:

将一个数用一个或多个正整数的无序和来表示

例如4的分拆有5种:4 , 3+1 , 2+2 , 2+1+1 , 1+1+1+1

还有一个概念叫做限制分拆,就是给分拆加上限制条件,下面给出一些常见的结论

n的分拆数中最大部分为m的个数=把n分拆成m部分的个数

n的分拆数中每一部分都小于等于m的个数=把n分成m份或更小

的分拆数中每部分的数都相等的个数=n 的因子个数
Eg. 6=2+2+2, 6=3+3,6=1+1+1+1+1+1

的分拆数中每部分都是1或2(或者把n分拆成1或2部分)的个数=floor(n/2+1);
Eg. 6=1+1+1+1+1+1, 6=1+1+1+1+2, 6=1+1+2+2, 6=2+2+2
floor()是向下取整

n的分拆数中每部分都是1或2或3(或者把n分拆成1或2或3部分)的个数=(n+3)^2/12;

然后是要点总结部分,转自https://blog.csdn.net/whai362/article/details/42530243

至于问题的变式,等介绍完母函数(生成函数)之后再进行介绍

HDU4651

 1 #include<cstdio>
 2 using namespace std;
 3 const int MOD=1e9+7;
 4 int n;
 5 int dp[100005];
 6 void init()
 7 {
 8     dp[0]=1;
 9     for(int i=1;i<=100000;i++)
10     {
11         for(int j=1,r=1;i-(3*j*j-j)/2>=0;j++,r*=-1)
12         {
13             dp[i]+=dp[i-(3*j*j-j)/2]*r;
14             dp[i]%=MOD;
15             dp[i]=(dp[i]+MOD)%MOD;
16             if(i-(3*j*j+j)/2>=0)
17             {
18                 dp[i]+=dp[i-(3*j*j+j)/2]*r;
19                 dp[i]%=MOD;
20                 dp[i]=(dp[i]+MOD)%MOD;
21             }
22         }
23     }
24 }
25 int main()
26 {
27     init();
28     int T;
29     scanf("%d",&T);
30     while(T--)
31     {
32         scanf("%d",&n);
33         printf("%d\n",dp[n]);    
34     }
35     return 0;
36 }

 

posted @ 2018-08-17 10:08  静听风吟。  阅读(1027)  评论(0编辑  收藏  举报