[SDOI2010]古代猪文

最新推荐文章于 2020-10-19 20:39:51 发布

qq_43520313

最新推荐文章于 2020-10-19 20:39:51 发布

阅读量173

点赞数

分类专栏：组合数学初等数论

本文链接：https://blog.csdn.net/qq_43520313/article/details/108232124

版权

组合数学同时被 2 个专栏收录

34 篇文章 1 订阅

订阅专栏

初等数论

12 篇文章 0 订阅

订阅专栏

题目：
https://ac.nowcoder.com/acm/problem/20333

求

$G^{\sum_{d|n}{n\choose i}}\ mod\;999911659$

思路：

$\begin{aligned} &G^{\sum_{d|n}{n\choose i}}\ mod\;p\\ =&G^{(\sum_{d|n}{n\choose i})mod\;p-1}\ mod\;p\\ \end{aligned}$
$p - 1 = 2 \times 3 \times 4697 \times 35617$
中国剩余定理 $+ L u c a s$ 求
$\sum_{d|n}{n\choose i}\equiv x\ mod\;p-1\Longleftrightarrow \begin{cases} \sum_{d|n}{n\;mod\;2\choose i\;mod\;2}\equiv x\ mod\;2\\ \sum_{d|n}{n\;mod\;3\choose i\;mod\;3}\equiv x\ mod\;3\\ \sum_{d|n}{n\;mod\;4697\choose i\;mod\;4697}\equiv x\ mod\;4697\\ \sum_{d|n}{n\;mod\;35617\choose i\;mod\;35617}\equiv x\ mod\;35617\\ \end{cases}$

最后再用一下快速幂即可。

注意：
当 $G = p$ 时，输出 $0$ 。

#include <bits/stdc++.h>
#define LL long long
using namespace std;
const int mod=999911658;
LL n,G,farc[50010],a[5],b[5]={0,2,3,4679,35617},val;

LL fast_pow(LL a,LL b,LL p)//快速幂
{
	LL ret=1;
	for(;b;b>>=1,a=a*a%p)
		ret=ret*(b&1?a:1)%p;
	return ret;
}

void init(LL p)//预处理
{
	farc[0]=1;
	for(LL i=1;i<=p;i++)
		farc[i]=farc[i-1]*i%p;
}

LL C(LL n,LL m,LL p)//组合数
{
	if(n<m) return 0;
	return farc[n]*fast_pow(farc[m],p-2,p)%p*fast_pow(farc[n-m],p-2,p)%p;
}

LL Lucas(LL n,LL m,LL p)//Lucas定理
{
	if(n<m) return 0;if(!n) return 1;
	return Lucas(n/p,m/p,p)*C(n%p,m%p,p)%p;
}

void CRT()//中国剩余定理
{
	for(LL i=1;i<=4;i++)
		val=(val+a[i]*(mod/b[i])%mod*fast_pow(mod/b[i],b[i]-2,b[i]))%mod;
}

int main()
{
	scanf("%lld%lld",&n,&G);
	if(G%(mod+1)==0){
		printf("0\n");
		return 0;
	}//特判
	for(LL k=1;k<=4;k++){
		init(b[k]);
		for(LL i=1;i*i<=n;i++){
			if(n%i==0){
				a[k]=(a[k]+Lucas(n,i,b[k]))%b[k];
				if(i*i!=n){
					a[k]=(a[k]+Lucas(n,n/i,b[k]))%b[k];
				}
			}
		}
	}//逐一枚举n的约数
	CRT();
	printf("%lld\n",fast_pow(G,val,mod+1));//注意mod要+1
	return 0;
}