hdu,1028,整数拆分的理解
#include"iostream"
using namespace std;
int main()
{
int n,i,j,k;
int c[122],temp[122];
//c[] 数组用于储存当前多项式各项系数
//temp[]数组用于暂时储存在运算时的两多项式相加的系数和
while(cin>>n&&n!=0)
{
for(i=0;i<122;i++) //系数初始化,当前c[]所指的多项式是第一个多项式
{c[i]=1; temp[i]=0;}
for(i=2;i<=n;i++) //运算次数的控制
{
for(j=0;j<=n;j++)
for(k=0;k+j<=n;k+=i)
temp[k+j]+=c[j]; //系数相加
for(j=0;j<=n;j++) //新的多项式系数和temp[]初始化
{
c[j]=temp[j];
temp[j]=0;
}
}
cout<<c[n]<<endl;
}
return 0;
}
看了一下午加一晚上的母函数还有代码,原谅我笨,总于明白了一点,所以就和大家分享一下。
源问题是杭电ACM的1028题,大概意思是求 整数 n 的拆分数,大概算法思路有两种,一种是多项式乘法,一种是递推,我起先也是从递推考虑,后来才发现我找不到递推规律,因为它是一个二维数组递推,record[n-m][m]+record[n][m-1] (N>m)(原谅我找不到他的规律)所以就在多项式上着手
先给出公式(这是有好处的):G(n)=(1+x+x^2+....x^i)*(1+x^2+x^4+....x^i)(1+x^3...+x^i).....(1+x^i)( i 是指最接近 n且满足规律且小于等于n )
至于式子,自己理解,大概就是排列组合,这个还是百度一下 母函数 的基础概念比较好(宝宝是数学渣渣)
对于算法的实现:我们可以这样想G1*G2*G3*G4...Gi,可以先算出G1*G2,得到Gh,再让Gh*G3,所以要有一个大的循环控制运算次数
对于多项式的次数可以用数组下标表示,而系数用数组值表示,储存在c[]数组中,而与c[]相乘的多项式可以用第三次循环表示for(k=0;k+j<=n;k+=i){},temp[i]用于暂时存储相乘后次数为 i 的项的系数,k+j 就是表示两个多项式其中的各自的某一项相乘的次数。每次计算完成后可以由temp[]转给c[],temp[]清零后接着与下一个多项项相乘
对于多项式建议大家好好研究研究在排列组合和一些比较麻烦的问题上相当的重要,泰勒公式没学好,宝宝心里苦呀!数学很重要,真的!对于多项式和母函数之类的数学问题等我补好了缺口再与大家分享。
hdu,1028,整数拆分的理解的更多相关文章
- Ignatius and the Princess III HDU - 1028 || 整数拆分,母函数
Ignatius and the Princess III HDU - 1028 整数划分问题 假的dp(复杂度不对) #include<cstdio> #include<cstri ...
- HDU 1028 整数拆分 HDU 2082 找单词 母函数
生成函数(母函数) 母函数又称生成函数.定义是给出序列:a0,a1,a2,...ak,...an, 那么函数G(x)=a0+a1*x+a2*x2+....+ak*xk +...+an* xn 称为序 ...
- HDU 1028 整数拆分问题 Ignatius and the Princess III
Ignatius and the Princess III Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K ...
- HDU 1028(数字拆分 分治)
题意是求所给的数能够被拆分成的不同组合数目. 方法有三种: 一.完全背包. 限制条件:所用数字不大于 n. 目标:求分解种数(组合出 n 的方法数). 令 dp[ i ][ j ] = x 表示 用前 ...
- hdu 1028 整数划分 (母函数)
假如输入44 = 4;4 = 3 + 1;4 = 2 + 2;4 = 2 + 1 + 1;4 = 1 + 1 + 1 + 1;一共5种 假如输入3 用母函数的方法就是写成(1+X+X2+X3)(1+X ...
- HDU 1028 Ignatius and the Princess III(母函数整数拆分)
链接:传送门 题意:一个数n有多少种拆分方法 思路:典型母函数在整数拆分上的应用 /********************************************************** ...
- hdu acm 1028 数字拆分Ignatius and the Princess III
Ignatius and the Princess III Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 65536/32768 K ...
- HDU 4651 Partition(整数拆分)
题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=4651 题意:给出n.求其整数拆分的方案数. i64 f[N]; void init(){ f[0 ...
- hdu 1028 & hdu 1398 —— 整数划分(生成函数)
题目:http://acm.hdu.edu.cn/showproblem.php?pid=1028 整数划分,每个数可以用无限次: 所以构造 f(x) = (1+x+x2+x3+...)(1+x2+x ...
随机推荐
- 奇妙的滚动css+html
<html lang="zh-CN"> <head> <meta charset="UTF-8"> <title> ...
- yum 源本地化 (two)
之前写过一个yum源本地化的文章. 后来发现那个方法有些缺陷, yum install --downloadonly --downloaddir=/tmp/atomicdownload memcach ...
- LeetCode 463. Island Perimeter (岛的周长)
You are given a map in form of a two-dimensional integer grid where 1 represents land and 0 represen ...
- [DB][MySql]关于取得自增字段的值、及@@IDENTITY 与并发性问题
对于刚从Oracle转向MySql的人都会为,MySql中没有Oracle里的Sequence而感到困惑.MySql中没有了Sequence,那么MySql的主键用什么方式来实现最好呢? 主要有下面几 ...
- User32.dll详细介绍
RegisterServiceProcess(ProcessID:Long,Type:Long) 该函数存在于Kernal32.dll中. Process指向进程的ID,Type表示是否向系统注册该进 ...
- Linux I2C设备驱动编写(一)【转】
本文转载自:http://blog.csdn.net/airk000/article/details/21345457 在Linux驱动中I2C系统中主要包含以下几个成员: I2C adapter 即 ...
- timeout in asp.net
Forms authentication timeout vs sessionState timeout They are different things. The Forms Authentica ...
- [Codeplus 2017] 晨跑
[题目链接] https://www.lydsy.com/JudgeOnline/problem.php?id=5105 [算法] 答案为三个数的最小公倍数 [代码] #include<bits ...
- 56. EditorGridPanel和渲染器renderer的使用
转自:https://blog.csdn.net/chenkangwan/article/details/4748716?utm_source=blogxgwz6 1. 在EditorGriddPan ...
- jsp jquery js 的基本路径获取
引子:js中需要当前页面的基础路径,获取不到request,可以通过如下方法来解决! 1.jsp基础路径,在jsp头部加上,获取基础路径http://localhost:8080/project/ ...