组合问题(高精度)

http://acm.hdu.edu.cn/showproblem.php?pid=1261

字串数

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


Problem Description
一个A和两个B一共可以组成三种字符串:"ABB","BAB","BBA".
给定若干字母和它们相应的个数,计算一共可以组成多少个不同的字符串.
 
Input
每组测试数据分两行,第一行为n(1<=n<=26),表示不同字母的个数,第二行为n个数A1,A2,...,An(1<=Ai<=12),表示每种字母的个数.测试数据以n=0为结束.
 
Output
对于每一组测试数据,输出一个m,表示一共有多少种字符串.
 
Sample Input
2 1 2 3 2 2 2 0
 
Sample Output
3 90
 
Source
//#include <bits/stdc++.h>
#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <iostream>
#include <string>
#include <stdio.h>
#include <queue>
#include <stack>
#include <map>
#include <set>
#include <string.h>
#include <vector>
#define ME(x , y) memset(x , y , sizeof(x))
#define SF(n) scanf("%d" , &n)
#define rep(i , n) for(int i = 0 ; i < n ; i ++)
#define INF  0x3f3f3f3f
#define mod 1000000007
#define PI acos(-1)
using namespace std;
typedef long long ll ;
int a[29];
int num[159];


int main()
{
    int n ;
    while(~scanf("%d" , &n) && n)
    {
        int sum = 0 ;
        for(int i = 0 ; i < n ; i++)
        {
            scanf("%d" , &a[i]);
            sum += a[i];
        }

        memset(num , 0 , sizeof(num));
        num[0] = 1 ;
        int len = 1 ;
        for(int i = 2 ; i <= sum ; i++)
        {
            int carry = 0 ;
            for(int j = 0 ; j < len ; j++)
            {
                num[j] *= i;
                num[j] += carry;
                if(num[j] >= 10000)
                {
                    carry = num[j] / 10000 ;
                    num[j] %= 10000 ;
                }
                else
                    carry = 0 ;
            }
            if(carry) num[len] = carry ;
            while(num[len]) len++;
        }
        for(int i = 0 ; i < n ; i++)
        {
            for(int j = 2 ; j <= a[i] ; j++)
            {
                int t = 0 ;
                int carry = 0 ;
                for(int k = len - 1 ; k >= 0 ; k--)
                {
                    carry = t ;
                    t = (num[k] + t * 10000) % j ;
                    num[k] = (num[k] + carry * 10000) / j;
                }
                while(num[len-1] == 0) len--;
            }
        }
        int i ;
        for(i = len ; i >= 0 ; i--)
        {
            if(num[i])
            {
                printf("%d" , num[i]);
                break ;
            }
        }
        for(i-- ; i >= 0 ; i--)
        {
            printf("%04d" , num[i]);
        }
        cout << endl ;
    }


    return 0;
}

猜你喜欢

转载自www.cnblogs.com/nonames/p/11723437.html