リンクのタイトルへ
の質問の意味:C [A] [B]の値を求めて、A、Bを入力します。
結果は、高精度計算を使用する大規模な、必要であってもよいことに留意されたいです。
入力フォーマット
二つの整数aとbを備えた共通線、。
出力フォーマット
共通線、出力C [A] [B]値。
データ範囲
1≤b≤a≤5000
入力サンプル:
53
のサンプル出力:
10の
アイデア:品質係数の分解:
C [A] [B] ^ P1 = M1 * M2 * P2 ^ ... ^ * PK MK;ここP1 。PKが素数である1〜A〜中、M1〜MKことで回数で、素数対応するそれぞれの数である-の数- - bが数である(b)は、
コード実装:
#include<iostream>
#include<vector>
#include<algorithm>
using namespace std;
const int N = 5010;
int primes[N], cnt;
int sum[N];
bool st[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] <= n / i; j ++){
st[primes[j] * i] = 1;
if(i % primes[j] == 0) break;
}
}
}
//求出n!里面包含p的个数
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(); 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);
//枚举一下每一个质数
for(int i = 0; i < cnt; i ++){
int p = primes[i];
sum[i] = get(a, p) - get(b, p) - get(a - b, p);
}
vector<int> res;
res.push_back(1);
for(int i = 0; i < cnt; i ++)
for(int j = 0; j < sum[i]; j ++)
res = mul(res, primes[i]);
for(int i = res.size() - 1; i >= 0; i --)
printf("%d", res[i]);
cout << endl;
return 0;
}