当前位置:小鱼儿玄机二站 > 计算机网络 > bzoj2982,02786872982

bzoj2982,02786872982

文章作者:计算机网络 上传时间:2019-11-12

bzoj2982,02786872982

Lucas定理裸题。。

Lucas定理:C(n,m)=C(n%p,m%p)*C(n/p,m/p)%p

预管理出阶乘、逆元的阶乘就足以了。

 

代码:

图片 1

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5 #include<cmath>
 6 using namespace std;
 7 #define p 10007
 8 int i,j,k,n,m,inv[p],fac[p],fac2[p],t;
 9 inline void Init(){
10     for(fac[0]=fac[1]=1,i=2;i<p;i++)fac[i]=fac[i-1]*i%p;
11     for(inv[0]=inv[1]=fac2[0]=fac2[1]=1,i=2;i<p;i++)inv[i]=inv[p%i]*(p-p/i)%p,fac2[i]=fac2[i-1]*inv[i]%p;
12 }
13 inline int C(int n,int m){
14     if(n<m)return 0;
15     return fac[n]*fac2[m]%p*fac2[n-m]%p;
16 }
17 inline int Get(int n,int m){
18     if(m==0)return 1;
19     return C(n%p,m%p)*Get(n/p,m/p)%p;
20 }
21 int main(){
22     scanf("%d",&t);
23     Init();
24     while(t--){
25         scanf("%d%d",&n,&m);
26         printf("%dn",Get(n,m));
27     }
28     return 0;
29 }

bzoj2982

 

Lucas定理裸题。。 Lucas定理:C(n,m)=C(n%p,m%p)*C(n/p,m/p)%p 预管理出阶乘、逆元的阶乘就能够了。 代码: 1 #includeiostream 2 #includ...

Lucas定理裸题。。

Lucas定理:C(n,m)=C(n%p,m%p)*C(n/p,m/p)%p

预管理出阶乘、逆元的阶乘就足以了。

 

代码:

图片 2图片 3

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5 #include<cmath>
 6 using namespace std;
 7 #define p 10007
 8 int i,j,k,n,m,inv[p],fac[p],fac2[p],t;
 9 inline void Init(){
10     for(fac[0]=fac[1]=1,i=2;i<p;i++)fac[i]=fac[i-1]*i%p;
11     for(inv[0]=inv[1]=fac2[0]=fac2[1]=1,i=2;i<p;i++)inv[i]=inv[p%i]*(p-p/i)%p,fac2[i]=fac2[i-1]*inv[i]%p;
12 }
13 inline int C(int n,int m){
14     if(n<m)return 0;
15     return fac[n]*fac2[m]%p*fac2[n-m]%p;
16 }
17 inline int Get(int n,int m){
18     if(m==0)return 1;
19     return C(n%p,m%p)*Get(n/p,m/p)%p;
20 }
21 int main(){
22     scanf("%d",&t);
23     Init();
24     while(t--){
25         scanf("%d%d",&n,&m);
26         printf("%dn",Get(n,m));
27     }
28     return 0;
29 }

bzoj2982

 

本文由小鱼儿玄机二站发布于计算机网络,转载请注明出处:bzoj2982,02786872982

关键词: