Luogu 1313（组合数学）（NOIP2011）

最新推荐文章于 2018-05-23 17:19:58 发布

嘉伟森的猫

最新推荐文章于 2018-05-23 17:19:58 发布

阅读量278

点赞数

分类专栏： OI-数学文章标签：组合数学

本文链接：https://blog.csdn.net/KGV093/article/details/77586865

版权

OI-数学专栏收录该内容

38 篇文章 0 订阅

订阅专栏

传送门
NOIP 2011 D2T1（快速幂，组合数，费马小定理，逆元）
题解就在代码之前，发现第一行所述性质就很容易搞定了。

/*
    在多项式(by+ax)^k中，某一项(x^n)*(y^m)中m=k-n恒成立
    ans=C(k,n)*(a^n)*(b^(k-n))
*/
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const ll MOD=1e4+7;
ll a,b,k,n,m;
ll fac[1002]={1,1};
ll fpow(ll x,ll y,ll p) {
    ll ret=1;
    while (y) {
        if (y&1) ret=ret*x%p;
        y>>=1,x=x*x%p;
    }
    return ret;
}
int main() {
    scanf("%lld%lld%lld%lld%lld",&a,&b,&k,&n,&m);
    for (int i=2;i<=1000;++i) fac[i]=fac[i-1]*(ll)i%MOD;
    ll comb=fac[k]*fpow(fac[k-n]*fac[n]%MOD,MOD-2,MOD)%MOD;
    printf("%lld\n",comb*fpow(a,n,MOD)%MOD*fpow(b,k-n,MOD)%MOD);
}