/* 整数划分问题是算法中的一个经典命题之一,有关这个问题的讲述在讲解到递归时基本都将涉及。 所谓整数划分,是指把一个正整数n写成如下形式: 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=4时,他有5个划分,{4},{3,1},{2,2},{2,1,1},{1,1,1,1}; 注意4=1+3 和 4=3+1被认为是同一个划分。 该问题是求出n的所有划分个数,即f(n, n)。下面我们考虑求f(n,m)的方法; 递归法: 根据n和m的关系,考虑以下几种情况: (1)当n=1时,不论m的值为多少(m>0),只有一种划分即{1}; (2)当m=1时,不论n的值为多少,只有一种划分即n个1,{1,1,1,...,1}; (3)当n=m时,根据划分中是否包含n,可以分为两种情况: (a)划分中包含n的情况,只有一个即{n}; (b)划分中不包含n的情况,这时划分中最大的数字也一定比n小,即n的所有(n-1)划分。 因此 f(n,n) =1 + f(n,n-1); (4)当n<m时,由于划分中不可能出现负数,因此就相当于f(n,n); (5)但n>m时,根据划分中是否包含最大值m,可以分为两种情况: (a)划分中包含m的情况,即{m, {x1,x2,...xi}}, 其中{x1,x2,... xi} 的和为n-m,因此这情况下 为f(n-m,m) (b)划分中不包含m的情况,则划分中所有值都比m小,即n的(m-1)划分,个数为f(n,m-1); 因此 f(n, m) = f(n-m, m)+f(n,m-1); 综上所述: 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) */ #include<iostream> #include<cstdio> using namespace std; int fun(int n,int m) { if(n==1||m==1) return 1; else if(n<m) return fun(n,n); else if(n==m) return 1+fun(n,n-1); else return fun(n,m-1)+fun(n-m,m); }
poj1221 单峰回文序列
所谓回文序列就是从左至右和从右至左读都是一样的,而单峰的意思就是说从左至中间的那个数是递增的
从中间至右是递减的,比如 1221,121满足条件,而2112不满足
题目大意:给你一个整数n,求出他的整数划分序列中满足单峰回文的个数
思路:
对于奇数而言:它满足的划分序列肯定是奇数个数组成的(如果是偶数个,因为又要满足回文那么这些数加起来肯定就是偶数了)
对于偶数:它满足的划分序列可以使奇数个或偶数个组成
对于数n:
1. 由奇数个组成的回文序列....k.....,那么除了中间数K外,其它的数的和肯定是一个偶数
因此n-k为偶数,因此我们只要从1--n枚举k就好了,并且k的左边所有数的和是(n-k)/2就好了,
剩下的我们只要求出fun((n-k/2),k)了
2. 偶数个组成的回文序列....KK.... K<=n/2,左半边的和为n/2,并且最大值小于等于n/2
所以只需要求fun(n/2,n/2)就好了
具体代码如下:
#include<iostream> #include<cstdio> using namespace std; int temp[1000][1000]; int fun(int n,int m) { //优化否则会超时,计算过的就不重新计算了 if(temp[n][m]) return temp[n][m]; if(n==1||m==1) return temp[n][m]=1; else if(n<m) return temp[n][m]=fun(n,n); else if(n==m) return temp[n][m]=1+fun(n,n-1); else return temp[n][m]=fun(n,m-1)+fun(n-m,m); } int main() { freopen("in.txt","r",stdin); int n; while(scanf("%d",&n)!=EOF&&n) { __int64 ans=1;//本身就是一种 for(int i=1;i<n;i++) if((n-i)%2==0) ans+=fun((n-i)/2,i); if(n%2==0) ans+=fun(n/2,n/2); printf("%d %I64d ",n,ans); } }