Codeforces Round #232 (Div. 1) 解题报告

遇到这种都是数学题的，非常合理的爆了一个零。然后A题的原因是没有想到用map计数，还有数组开小了。。。

Problem A On Number of Decompositions into Multipliers

题意：给一个n个元素的数组ai求所有元素的乘积分解为n个数的情况总数。

思路：一次拆分数组元素，由于又大素数存在所以bixuyongmap标记然后推出公式这个直接看程序吧，不好打。

代码如下：

/**************************************************
 * Author     : xiaohao Z
 * Blog     : http://www.cnblogs.com/shu-xiaohao/
 * Last modified : 2014-02-26 23:26
 * Filename     :
 * Description     :
 * ************************************************/

#include <iostream>
#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <algorithm>
#include <queue>
#include <stack>
#include <vector>
#include <set>
#include <map>
#define MP(a, b) make_pair(a, b)
#define PB(a) push_back(a)

using namespace std;
typedef long long ll;
typedef pair<int, int> pii;
typedef pair<unsigned int,unsigned int> puu;
typedef pair<int, double> pid;
typedef pair<ll, int> pli;
typedef pair<int, ll> pil;

const int INF = 0x3f3f3f3f;
const double eps = 1E-6;
const int MOD = 1000000007;
ll C[20010][510];
map<ll, ll> mp;

const int MAXN=1000000;
int prime[MAXN+1];
void getPrime()
{
    memset(prime,0,sizeof(prime));
    for(int i=2; i<=MAXN; i++)
    {
        if(!prime[i])prime[++prime[0]]=i;
        for(int j=1; j<=prime[0]&&prime[j]<=MAXN/i; j++)
        {
            prime[prime[j]*i]=1;
            if(i%prime[j]==0) break;
        }
    }
}

void make(int n){
    int temp = sqrt(n);
    for(int i=1; prime[i] <= temp; i++){
        while(n%prime[i]==0){
            mp[prime[i]]++;
            n/=prime[i];
        }
    }
    if(n>1) mp[n] ++;
}

void init(){
    memset(C, 0, sizeof C);
    C[0][0] = 1;
    for(int i=0; i<20010; i++){
        C[i][0] = 1;
        if(i < 510)C[i][i] = 1;
        for(int j=1; j<min(510, i); j++){
            C[i][j] = (C[i-1][j] + C[i-1][j-1]) % MOD;
        }
    }
}

int main()
{
//    freopen("in.txt", "r", stdin);

    int n, tnum;
    while(scanf("%d", &n)!=EOF){
        init();getPrime();mp.clear();
        for(int i=0; i<n; i++){
            scanf("%d", &tnum);
            make(tnum);
        }
        ll ans = 1, tans;
        for(map<ll, ll>::iterator it=mp.begin(); it!=mp.end();++it){
            tans = 0;
            for(int j=0; j<n; j++) {
                tans = (tans+(C[n][j+1]*C[it->second-1][j])%MOD)%MOD;
            }
            ans = ans * tans;
            ans%=MOD;
        }
        printf("%I64d
", ans);
    }
    return 0;
}