888. 求组合数 IV (高精乘,预处理阶乘中质数的次数)
作者:
闪回
,
2024-03-17 21:09:17
,
所有人可见
,
阅读 11
#include<bits/stdc++.h>
using namespace std;
const int N = 5010;
int primes[N],cnt;
bool st[N];
int sum[N];
void get_primes(int n)
{
for(int i = 2;i<=n;i++)
{
if(!st[i])primes[cnt++] = i;
for(int j = 0;primes[j] * i<=n;j++)
{
st[primes[j] * i] = true;
if(i%primes[j] == 0)break;
}
}
}
int get(int n,int p)
{
int res = 0;
while(n)
{
res += n/p;
n/=p;
}
return res;
}
vector<int> mul(vector<int> &A,int b)
{
vector<int> C;
int t = 0;
for(int i = 0;i<A.size() || t;i++)
{
if(i<A.size())t += A[i] * b;
C.push_back(t % 10);
t /= 10;
}
while(C.back() == 0 && C.size() > 1)C.pop_back();
return C;
}
int main()
{
int a,b;
cin>>a>>b;
get_primes(a);
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> ans;
ans.push_back(1);
for(int i = 0;i<cnt;i++)
{
for(int j = 0;j<sum[i];j++)
{
ans = mul(ans,primes[i]);
}
}
for(int i = ans.size()-1;i>=0;i--)cout<<ans[i];
return 0;
}