阶乘分解——质因数分解+组合数的计算

首先介绍算数基本定理

$\mathbf{m}={\mathbf{P}}_1^{{\mathbf{\alpha }}_1}{\mathbf{P}}_2^{{\mathbf{\alpha }}_2}{\mathbf{P}}_3^{{\mathbf{\alpha }}_3}{\mathbf{P}}_4^{{\mathbf{\alpha }}_4}.....{\mathbf{P}}_{\mathbf{k}}^{{\mathbf{\alpha }}_{\mathbf{k}}}$

1. 其中$P_1,P_2......P_k$都是质数，即为任何一个数都可以表示成为质数的组合乘积

接下来介绍 $n！$的质数分解

• 首先我们以一个数字为例子 $12!$
  $$12!=12\times 11\times 10\times 9\times 8\times 7\times 6\times 5\times 4\times 3\times 2\times 1$$

• 观察他的质因数为$$2,3,5,11,7$$

• 我们在观察一下这些质因数可能得倍数：
  $$2:2,4,6,8,10,12$$$$3:3,6,9,12$$$$5:5,10$$$$11:11$$$$7:7$$

• 根据我们发掘出来的质因数$12!$可以写为
  $$12!=(3^1\times 2^2)\times 11^1\times (5^1\times 2^1)\times (3^2)\times (1\times 2^3)\times 7^1\times (3\times 2^1)\times 5\times 2^2\times 3^1\times 2^1\times 1$$

• 咱们在进行一些整理,看看2的不同次方的个数；
  个数为(以下除法均为向下取整)：$$\frac{12}{2^3}个+\frac{12}{2^2}个+\frac{12}{2^1}个$$
  $$一共为1+3+6=10个$$

• 上面发现质数$$P的个数=\lfloor \frac{a}{p^1}\rfloor +\lfloor \frac{a}{p^2}\rfloor +\lfloor \frac{a}{p^3}\rfloor +.......\lfloor \frac{a}{p^k}\rfloor$$

• 所以我们就可以算出所有的质因子个数然后再根据基本算数定理就可以求出$a!$

接下来咱们开始利用这个性质来做组合数的题

选自AcWing 888. 求组合数 IV

• 输入$a,b$ 求$C_{a}b$的数值

基本思路：

1.因为$C_{a}b=\frac{a!}{b!\times (a-b)!}$

2.根据 基本算数定理 我们把$C_{a}b=P_1^{a_1}\times P_2^{a_2}\times .......\times P_k^{ak}$

3.除法无非就是把对应的素数个数减去即为：
$$P_i=\sum _{j=1}^n\frac{a}{p_i^j}-\sum _{j=1}^n\frac{a-b}{p_i^j}-\sum _{j=1}^n\frac{b}{p_i^j}(公式如有错误请指正)$$
4.再把所有的$P_i$相乘就是最终的答案；
5.上代码

#include<iostream>
#include<algorithm>
#include<vector>

using namespace std;

const int N=5010;

int primes[N],cnt;
bool st[N];
int sum[N];//存放素数的个数；

void get_primes(int x)
{
    for(int i=2;i<=x;i++){
        
        if(!st[i])primes[cnt++]=i;
        
        for(int j=0;primes[j]<=x/i;j++)
        {
            st[primes[j]*i]=true;
            if(i%primes[j]==0)break;
        }
        
    }
}

int get(int x,int p)
{
    int res=0;
   while(x)
   {
       res+=x/p;
       x/=p;
   }
   return res;
}

vector<int> mul(vector<int>a,int b)
{
    vector<int>c;
    int t=0;
    for(int i=0;i<(int)a.size();i++)
    {
        t+=a[i]*b;
        c.push_back(t%10);
        t/=10;
    }
    while(t)
    {
        c.push_back(t%10);
        t/=10;
    }
    return c;
}
int main()
{
    int a,b;
    cin>>a>>b;
    get_primes(a);//把可能出现的次数存入primes数组里；
    for(int i=0;i<cnt;i++)
    {
        int p=primes[i];
        sum[i]=get(a,p)-get(a-b,p)-get(b,p);
    }//获取每个次数的个数；
    //下面用高精度来算乘法；
    vector<int>end;
    end.push_back(1);
    for(int i=0;i<cnt;i++)
    {
        for(int j=0;j<sum[i];j++)
        {
            end=mul(end,primes[i]);
        }
    }
    for(int i=(int)end.size()-1;i>=0;i--)cout<<end[i];
}