题目链接: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;
}
