rqnoj[NOIP2011提高组]计算系数

链接:http://www.rqnoj.cn/problem/659

题意:中文题。。。

分析:不能再裸的组合数了。。练下逆元。

代码:

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<math.h>
#include<cstdio>
#include<vector>
#include<string>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=410;
const int MAX=151;
const int MOD1=100000007;
const int MOD2=100000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=10007;
const ll INF=10000000010;
typedef unsigned long long ull;
int inv[1010];
int main()
{
    int a,b,k,n,m,i;
    ll ans=1;
    scanf("%d%d%d%d%d", &a, &b, &k, &n, &m);
    for (i=1;i<=n;i++) ans=(ans*a)%MOD;
    for (i=1;i<=m;i++) ans=(ans*b)%MOD;
    inv[1]=1;
    for (i=2;i<=n;i++) inv[i]=(MOD-MOD/i)*inv[MOD%i]%MOD;
    for (i=k-n+1;i<=k;i++) ans=(ans*i)%MOD;
    for (i=1;i<=n;i++) ans=(ans*inv[i])%MOD;
    printf("%lld\n", (ans+MOD)%MOD);
    return 0;
}

/*
1 1 3 1 2
*/


你可能感兴趣的:(rqnoj[NOIP2011提高组]计算系数)