吉首 - 组合数(求组合数因子个数)

题目链接:http://120.78.162.102/problem.php?id=6240
时间限制: 1 Sec  内存限制: 128 MB

题目描述

求组合数C(N,M)，以及C(N,M)因子个数。

输入

N和M，其中0<=M<=N<=50,以EOF结束。

输出

该组合数结果。

样例输入

3 2
4 2

样例输出

3 2
6 4

解题思路

先利用杨辉三角求出组合数，然后就是求出因子数了;
求因子数:素数分解的唯一性，一个数可以被分解成若干素数相乘 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循环从1到n枚举i统计。

#include <math.h>
#include <stdio.h>
#include <string.h>
long long a[60][60];
const int maxn = 60;
int vis[maxn], dp[maxn][maxn];
void C()
{
    memset(a, 0, sizeof(a));
    a[0][0] = 1;
    for (int i = 1; i <= 55; i++)
        a[1][i] = i,a[0][i] = 1;
    for (int r = 2 ; r <= 55 ; r++)
        for (int c = r; c <= 55; c++)
            a[r][c] = a[r][c - 1] + a[r - 1][c - 1];
}
void sieve(int n)
{
    int m = sqrt(n + 0.5);
    memset(vis, 0, sizeof(vis));
    vis[0] = vis[1] = 1;
    for (int i = 2; i <= m; i++)
        if (!vis[i])
            for (int j = i * i; j <= n; j += i)
                vis[j] = 1;
}
int get(int n, int m)
{
    int sum = 0;
    while (n)
    {
        sum += n / m;
        n /= m;
    }
    return sum;
}
void pre(int n)
{
    for (int i = 2; i <= n; i++)
        for (int j = 2; j <= i; j++)
            dp[i][j] = get(i, j);
}
int MM(int n, int m)
{
    pre(55);
    sieve(55);
    long long ans = 1;
    for (int i = 2; i <= n; i++)
        if (!vis[i])
            ans *= dp[n][i] - dp[m][i] - dp[n - m][i] + 1;
    return ans;
}
int main()
{
    C();
    int n, m;
    while (~scanf("%d%d", &n, &m))
        printf("%lld %lld\n", a[m][n], MM(n, m));
    return 0;
}