（排列组合模板题）南华大学第十五届ACM程序设计竞赛 F题有趣的数字

最新推荐文章于 2019-05-01 10:53:02 发布

irvix

最新推荐文章于 2019-05-01 10:53:02 发布

阅读量232

点赞数

分类专栏：模板题数学

本文链接：https://blog.csdn.net/qq_40526226/article/details/89424012

版权

数学同时被 2 个专栏收录

4 篇文章 0 订阅

订阅专栏

模板题

3 篇文章 0 订阅

订阅专栏

Input

3 2 1 4 1 2 1

Output

144

题解：很明显，就是要求出 C(k, e) * C(k-e, f) * (a^e) * (b^f) * (c^g) 的值，那么我们就需要使用快速幂和组合数这两部分来完成了。

我的代(la)码(ji)：

第17-53行的代码可以直接当模板拿来打表。（使用费马小定理）

#include <cstdio>
#include <cstring>
#include <cstdlib>
#include <iostream>
#include <algorithm>
#include <string>
#include <vector>
#include <set>
using namespace std;

#define LL          long long
#define zero(a)     memset(a,0,sizeof(a))

const int mod = 10007;
typedef long long ll;

const int N = 1e3+5;                        //这里使用了组合数打表的模板，N代表C底数的上限
ll a[N] = {0};                              //打表后的结果
ll Pow(ll a,ll b){                          //记得main函数开始要调用init();
    a%=mod;
    ll ans = 1;
    while(b)
    {
        if(b&1)
            ans = (ans%mod)*(a%mod)%mod;
        a = (a%mod)*(a%mod)%mod;
        b/=2;
    }
    return ans%mod;
}
ll Quk(ll a,ll b){
    a%=mod;
    ll ans = 0;
    while(b)
    {
        if(b&1)
            ans = (ans+a)%mod;
        a = (a+a)%mod;
        b/=2;
    }
    return ans%mod;
}
ll C(ll n,ll m){
    return Quk(Quk(a[n],Pow(a[n-m],mod-2)),Pow(a[m],mod-2))%mod;
}
ll A(ll n,ll m){
    return Quk(a[n],Pow(a[n-m],mod-2))%mod;
}
void init(){
    a[0]=a[1]=1;
    for(ll i=2;i<=N-5;i++)
        a[i]=Quk(a[i-1],i);
}
int main()
{
    init();
    ll a,b,c,k,e,f,g;
    scanf("%lld%lld%lld%lld%lld%lld%lld",&a,&b,&c,&k,&e,&f,&g);

    LL ans = 1;
    ans *= C(k, e);
    ans *= C(k-e, f);

    ans = ((Pow(a,e)%mod) * ans)%mod;
    ans = ((Pow(b,f)%mod) * ans)%mod;
    ans = ((Pow(c,g)%mod) * ans)%mod;
    printf("%lld",ans);

    return 0;
}

大佬的代码：

快速幂的模板+阶乘的模板不错。

另外由组合数公式：直接把两个组合数乘积化简转换成求 k! / (e!*f!*g!) %mod，也就是k! * (那一坨的逆元)，也就是k! * ((e!*f!*g!) ^ (mod-2)) 就得到了比我铁憨憨地直接求写着简单多了，高中数学都丢了。

#include<cstdio>
#include<algorithm>
#include<set>
#define rep(i,s,t) for(register int i=s;i<=t;++i)
using namespace std;
const int N=1e6+11,mod=10007;
typedef long long ll;
int a,b,c,k,e,f,g;
inline int fp(int a,int b){    //快速幂
    int res=1;
    for(;b;b>>=1,a=a*a%mod)
        if(b&1)
            res=res*a%mod;
    return res;
}
inline int fac(int x){         //阶乘
    int ans=1;
    rep(i,1,x)
        ans=ans*i%mod;
    return ans;
}
int main(){
    scanf("%d%d%d%d%d%d%d",&a,&b,&c,&k,&e,&f,&g);
    int ans=1ll*fac(k)*fp(1ll*fac(e)*fac(g)*fac(f)%mod,mod-2)%mod;
    ans=1ll*ans*fp(a%mod,e)*fp(b%mod,f)*fp(c%mod,g)%mod;
    printf("%d\n",ans);
    return 0;
}

总结：这题就算是用来熟悉模板的