排列组合 + 指数型母函数

排列组合

Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)
Total Submission(s): 1569    Accepted Submission(s): 659

Problem Description

有n种物品，并且知道每种物品的数量。要求从中选出m件物品的排列数。例如有两种物品A,B，并且数量都是1，从中选2件物品，则排列有"AB","BA"两种。

 

Input

每组输入数据有两行，第一行是二个数n,m(1<=m,n<=10)，表示物品数，第二行有n个数，分别表示这n件物品的数量。

 

Output

对应每组数据输出排列数。(任何运算不会超出2^31的范围)

 

Sample Input

  
  

   
   2 2
1 1
  
  

 

Sample Output

  
  

   
   2
  
  
  
  

   
    
  
  
  
  

   
   指数型母函数的详解  ： 
   
   http://blog.csdn.net/wenlovingliu/article/details/8966293
  
  
  
  

   
    
  
  
  
  

   
   
#include<stdio.h>
#include<string.h>
double fac[]={1,1,2,6,24,120,720,5040,40320,362880,3628800};      // 0~10的阶乘
int a[11];
double b[11],c[11];
int main()
{
	int  m,n,i,len,j,k;
	while(~scanf("%d%d",&n,&m))
	{
		for(i=1;i<=n;i++) scanf("%d",&a[i]);
		for(i=0;i<11;i++) b[i]=c[i]=0;
		for(i=0;i<=a[1];i++) c[i]=1.00/fac[i];  //和普通母函数的区别                     
		for(i=2;i<=n;i++)
		{
			for(j=0;j<=m;j++)
			{
				for(k=0;j+k<=m&&k<=a[i];k++)
					b[j+k]+=c[j]/fac[k];                             
			}
			for(j=0;j<=m;j++)
			{
				c[j]=b[j];
				b[j]=0;
			}
		}
		printf("%.0lf\n",c[m]*fac[m]);  
		
		//语句若写成printf("%d\n",int(c[m]*fac[m]));则会WA，double转换成int类型有精度损失，而只取double的整数部分则会AC
	}
	return 0;
}
  
  
  
  

   
   
#include<iostream>
#include<iomanip>
using namespace std ;
#define MAX 15
double c1[MAX] ,c2[MAX] ,f[10];
void fun() //求1-10之间的阶乘
{
	int i , j ;
	f[0] = 1 ;
	for( i = 1 ; i <10 ; i++)
	{
		int sum = 1;
		for( j = 1 ; j <= i ; j++)
			sum *=j;
		f[i] = sum ;
	}
	return ;
}
int main(void)
{
	fun();
	int n , m ;
	while(cin >> n >>m)
	{
		int i,j,k,a[MAX];
		for( i = 0 ; i < n ; i++)
			cin >> a[i] ;
		for( i = 0 ; i < MAX ; i++)
			c1[i] = c2[i] = 0 ;
		for( i = 0 ; i <= a[0] ; i ++)
			c1[i] = 1.00/f[i] ;
		for( i = 1 ; i< n ; i ++)
		{
			for( j = 0 ; j <= m ; j++)
				for( k = 0 ;j+k <= m &&  k<= a[i] ; k++)
					c2[j + k] += c1[j] / f[k] ;
			for( j = 0 ; j<= m ; j ++)
			{
				c1[j] = c2[j] ;
				c2[j] = 0 ;
			}
				
		}
		cout<<fixed<<setprecision(0)<<c1[m]*f[m]<<endl;
	}
	return 0;
}