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 ...
随机推荐
- CODEVS——T 2956 排队问题
http://codevs.cn/problem/2956/ 时间限制: 1 s 空间限制: 32000 KB 题目等级 : 黄金 Gold 题解 查看运行结果 题目描述 Descri ...
- Mybatis错误——Could not find parameter map java.util.Map
错误信息 org.apache.ibatis.builder.IncompleteElementException: Could not find parameter map java.util.Ma ...
- 1. 少了一个PermMissingElem Find the missing element in a given permutation.
少了一个: package com.code; import java.util.Arrays; public class Test03_2 { public static int solution( ...
- GNS3 使用SecureCRT
"D:\SecureCRT\SecureCRT.EXE" /SCRIPT D:\SecureCRT\DyRouter.vbs /ARG %d /T /TELNET %h %p st ...
- Java课程设计——人事管理系统
主界面代码: package PersonSystem; import java.awt.*; import java.awt.event.*; import javax.swing.*; impor ...
- android:descendantFocusability
开发中很常见的一个问题,项目中的GridView不仅仅是简单的文字和图片,常常需要自己定义GridView,自己的Adapter去继承BaseAdapter,在adapter中按照需求进行编写,问题就 ...
- Bing Maps进阶系列一:初识Bing Maps地图服务
Bing Maps进阶系列一:初识Bing Maps地图服务 Bing Maps提供了一组WCF的地图服务,使用这些服务我们可以方便的在自己的应用系统里实现地理位置搜索等相关功能.他们分别是地理编码服 ...
- POJ2492 A Bug's Life 判断二分图
给一个无向图,判断是否为二分图 是否有奇环.用染色法解决,也可以用并查集解决,时间复杂度是相同的,但是并查集用的内存少. 这里给出Bfs染色的代码 #include<iostream> # ...
- Android之NDK开发(转载)
http://www.cnblogs.com/devinzhang/archive/2012/02/29/2373729.html 一.NDK产生的背景 Android平台从诞生起,就已经支持C.C+ ...
- Akka源码分析-Event Bus
akka中的EventBus其实是不常用,也最容易被忽略的一个组件. 但如果你深入Cluster的实现就会发现,这个东西其实还挺有用的,而且它是ActorSystem系统中所有事件消息的一个横切面,通 ...