快速幂+组合数

来源:互联网 发布:新奔腾计价软件 编辑:程序博客网 时间:2024/06/05 07:47

此题以NOIP 2011 day2 t1为例
主要是记录模板,就不多加解释

#include<iostream>#include<cstring>#include<algorithm>#include<cstdio>#include<cmath>using namespace std;int n,m,k,a,b;int C[1005][1005]; inline int mi(int a,int b){    long long w=a,s=1;    for(;b;b>>=1){        if(b&1) s=(s*w)%10007;        w=(w*w)%10007;    }    return s;}int main(){    scanf("%d%d%d%d%d",&a,&b,&k,&n,&m);    C[0][0]=1;    for(int i=1;i<=k;i++) C[i][0]=1,C[i][i]=1;    for(int i=1;i<=k;i++)     for(int j=1;j<i;j++)      C[i][j]=(C[i-1][j]+C[i-1][j-1])%10007;    printf("%d",((C[k][m]*mi(a,n)%10007)*mi(b,m))%10007);    return 0;}
原创粉丝点击