动态规划之求组合数

时间:2014.05.29

地点:基地

--------------------------------------------------------

一、关于动态规划

如果问题是由交叠的子问题构成,则可用动态规划的方式求解。我们在将一个大问题划分为子问题的的过程中,如果递推关系中包含的子问题和大问题具有相同的形式,但由于子问题的交叠性质,我们用递归解决的代价往往很大,这时可考虑动态规划,即对每个交叠的子问题只求解一次,并把结果存储在记录表中,最后得出原始问题的解。直观上,好像动态规划是采取空间换时间的权衡技术,但经过改进,可以避免使用额外的空间,一个典型的应用就是求斐波那契数列的第n项和组合数。这两个应用都属于经典的从底到上动态规划的应用,即需要解出给定问题的所有较小子问题。

--------------------------------------------------------

二、求组合数

C(n,k)=C(n-1,k)+C(n-1,k-1)                              ————————F(n)=F(n-1)+F(n-2)也类似

像这样,对于给定问题可分解为子问题,而子问题又存在交叠,适合于动态规划,更具体的这种形式适合采取自底向上的动态规划策略。比如这里:

我们先从C(0,0)算起,然后C(1,0),C(1,1),C(2,0),C(2,1),C(2,2)......C(n,k)

这样还可以为之建立一张表,我们对数据进行动态更新,即没一次迭代,我们都去更具公式计算出合适的值。

计算组合的公式在程序上可实现为:正上方值+左上方值

为了节省存储空间,对于计算C(n,k)实际上我们只需要一个k+1维的数组即可,通过对该数组的n+1次维护即可得出目标值,当然注意到0号元素总是1可以对问题略加简化:

--------------------------------------------------------

三、动态求解组合数的伪代码

Binomial(n,k)
  for i=[1..n]
      left_top=1         //计算每一行最开始时,left_top初始化为1,因为C(n,0)总是为1,所以每行计算第二个值时,它左上方元素值总是为1
       for j=[1...k]     //从1号元素算起,0号元素总是为1
          if   j=k
               C[j]=1    //计算对角线,即C(j,j)
               break     //结束当前行(轮)计算
           else
               temp=C[j]             //必须保存旧值,因为当前值将被更新,然而我们计算下一个值时,该位置的旧值是下一个元素的左上方值
               C[J]+=left_top        //更新当前值=正上方值(旧值)+左上方值
               left_top=temp         //更新左上方值用于下次计算
  return C[k]

--------------------------------------------------------

四、源码实现

int Binomial(size_t n, size_t k)
{
	//Precontion:
	//Postcontion:
	//Facilities used:
	assert(k>=0&&k <= n);
	int* arr_ptr = new int[k+1];
	arr_ptr[0] = 1;
	int temp, left_top;
	for (size_t i = 1; i < n+1; ++i)
	{
		left_top = 1;
		for (size_t j = 1; j < k+1; ++j)
		{
			if (j == i)
			{
				arr_ptr[j] = 1;
				break;
			}
			else
			{
				temp = arr_ptr[j];
				arr_ptr[j] += left_top;
				left_top = temp;
			}
		}
	}
	temp=arr_ptr[k];
	delete[] arr_ptr;
	return temp;
}

动态规划之求组合数,布布扣,bubuko.com

时间: 2024-10-09 22:25:07

动态规划之求组合数的相关文章

求组合数

组合数的计算虽说简单但也不乏有些陷阱,这主要是因为语言中的数据类型在表示范围上是有限的.更何况还有中间结果溢出的现象,所以千万要小心. 输入 求组合数的数据都是成对(M与N)出现的,每对整数M和N满足0<m, n≤20,以EOF结束. 输出 输出该组合数.每个组合数换行. 样例输入 5 2 18 13 样例输出 10 8568 代码 #include<stdio.h> int main(){int isum=1;int m,n,k;while(scanf("%d%d"

POJ 2992 Divisors 求组合数因子个数

题目来源:POJ 2992 Divisors 题意:... 思路:素数分解的唯一性 一个数可以被分解成若干素数相乘 p1^x1*p2^x2*...*pn^xn 根据乘法原理 因子数为 (x1+1)*(x2+1)*...*(xn+1) 不能直接求出组合数 会溢出 也不能把每个乘的数分解因子 这样会超时 C(N,M)=N!/(M!*(N-M)!) 另dp[i][j] 代表为i的阶乘中j因子的个数(j是素数) 那么i素数的个数为dp[n][i]-dp[m][i]-dp[n-m][i] 最后for循环从

[2011山东ACM省赛] Binomial Coeffcients(求组合数)

Binomial Coeffcients nid=24#time" style="padding-bottom:0px; margin:0px; padding-left:0px; padding-right:0px; color:rgb(83,113,197); text-decoration:none; padding-top:0px"> Time Limit: 1000ms   Memory limit: 65536K  有疑问?点这里^_^ 题目描写叙述 输入

hdu 2519 求组合数

求组合数 如果求C5 3 就是5*4*3/3*2*1 也就是(5/3)*(4/2)*(3/1) Sample Input5 //T3 2 //C3 25 34 43 68 0 Sample Output310101 1 # include <iostream> 2 # include <cstdio> 3 # include <cstring> 4 # include <algorithm> 5 # include <cmath> 6 # def

hdu 4927 java求组合数(大数)

import java.util.Scanner; import java.math.BigInteger; public class Main { private static int [] a = new int[3005]; private static int T; private static int n; public static void solve() { BigInteger rs = BigInteger.ZERO; BigInteger temp = BigInteger

求组合数小结

今天学了一天数学,觉得自己都要转竞了23333 题目链接https://vjudge.net/contest/282927#problem/E 这里说一说求组合数的方法吧 其实就是求阶乘及其逆元的方法: 规定mod为模数,n为数据规模 1.mod为素数 费马小定理:nlogn 线性求逆元(n较小) 2.mod不为素数 欧拉筛出两个阶乘的素数,再计算出每个素数的次幂,最后快速幂乘起来即可 详见 n = read(),p = read(); ans = 1; for(int i = 2;i <= 2

求组合数C(m,n)的多种计算方法

https://ac.nowcoder.com/discuss/187813?type=101&order=0&pos=1&page=0 https://blog.csdn.net/shadandeajian/article/details/82084087 1.简单法---适合n,m很小 #include<bits/stdc++.h> using namespace std; const int MAXN = 1000; int C[MAXN+1][MAXN+1];

常见的求组合数的板子

1 LL C[3010][3010]; 2 3 void init() { 4 C[0][0] = 1; 5 for(int i = 1; i < 3010; i++) { 6 C[i][0] = 1; 7 for(int j = 1; j <= i; j++) { 8 C[i][j] = (C[i - 1][j] + C[i - 1][j - 1]) % MOD; 9 } 10 } 11 } 费马小定理加快速幂进行优化求组合数 LL da[MAXN];//G++ long long void

[组合数]求组合数的几种方法总结

求C(n,m)%mod的方法总结 1.当n,m都很小的时候可以利用杨辉三角直接求. C(n,m)=C(n-1,m)+C(n-1,m-1): 2.利用乘法逆元. 乘法逆元:(a/b)%mod=a*(b^(mod-2)) mod为素数. 逆元可以利用扩展欧几里德或欧拉函数求得: 1).扩展欧几里德:b*x+p*y=1 有解,x就是所求 2).费马小定理:b^(p-1)=1(mod p),故b*b^(p-2)=1(mod p),所以x=b^(p-2) 1. n!/(m!*(n-m)! =x%p ,先对