标签:
Description
Input
Output
Sample Input
Sample Output
可以轻易推出公式 :(n1+n2+n3+...nn)!/(n1!*n2!*...*nn!);
因为15!还在long long的范围之内,可以先定义一个数组f[15]保存1~15的阶乘,接着就是将(n1+n2+n3+...nn)!计算出来并存到数组内,接着就是大数除法了(相当于一个大数除一个小数)。
#include <iostream> #include <algorithm> #include <cstring> #include <cstdio> #include <cctype> #include <cstdlib> #include <stack> #include <cmath> #include <string> #include <queue> using namespace std; #define SIZE 30 typedef long long ll ; int d[SIZE] ; int ans[1000] , f[15]; void multiply(int c){ ans[0] = ans[1] = 1 ; for(int i = 2 ; i <= c ; ++i){ int r = 0 ; for(int j = 1 ; j <= ans[0] ; ++j){ ans[j] *= i ; ans[j] += r ; r = ans[j]/10 ; ans[j] %= 10 ; } if(r != 0){ while(r){ ans[ans[0]+1] += r%10 ; ans[0] = ans[0]+1 ; r /= 10 ; } } } } void divide(int n){ for(int i = 0 ; i < n ; ++i){ if(d[i] == 1) continue ; ll r = 0 ; for(int j = ans[0] ; j > 0 ; --j){ r = r*10 + ans[j] ; ans[j] = (int)(r/f[d[i]]) ; r %= f[d[i]] ; } int j = ans[0] ; while(!ans[j--]) ; ans[0] = j+1 ; } } int main(){ int n ; f[0] = f[1] = 1 ; for(int i = 2 ; i < 15 ; ++i) f[i] = f[i-1]*i ; while(scanf("%d",&n) && n){ int c = 0; memset(ans,0,sizeof(ans)) ; for(int i = 0 ; i < n ; ++i){ scanf("%d",&d[i]) ; c += d[i] ; } multiply(c) ; divide(n) ; for(int i = ans[0] ; i > 0 ; --i) printf("%d",ans[i]) ; puts("") ; } return 0 ; }
标签:
原文地址:http://www.cnblogs.com/Asimple/p/5698995.html