acm-(数论、莫比乌斯反演、卡空间)LuoGu P5221 Product

最新推荐文章于 2020-10-20 15:59:00 发布

&*^*&

最新推荐文章于 2020-10-20 15:59:00 发布

阅读量179

点赞数

文章标签： acm竞赛算法

本文链接：https://blog.csdn.net/qq_30361651/article/details/109138713

版权

数论同时被 2 个专栏收录

11 篇文章 0 订阅

订阅专栏

莫比乌斯反演

9 篇文章 0 订阅

订阅专栏

传送门
由于 $\mathbf{\prod_{i=1}^n\prod_{j=1}^n\frac{lcm(i,j)}{gcd(i,j)}=\prod_{i=1}^n\prod_{j=1}^n\frac{ij}{gcd(i,j)^2}=\frac{\prod_{i=1}^n\prod_{j=1}^nij}{\prod_{i=1}^n\prod_{j=1}^ngcd(i,j)^2}}$ ，因此我们可以分成上下两部分来求，对于分子我们有 $\mathbf{\prod_{i=1}^n\prod_{j=1}^nij=\prod_{i=1}^ni^n\prod_{j=1}^nj=\prod_{i=1}^ni^nn!=(n!)^n\prod_{i=1}^ni^n=(n!)^n(n!)^n=(n!)^{2n}}$ 。
然后重头戏是下部分，类似于普通的莫比乌斯反演，我们更改枚举的对象，设 $\mathbf{d=gcd(i,j)}$ ，我们枚举这个 $\mathbf{d}$ 从1~n，于是有 $\mathbf{\prod_{i=1}^nd^{k}}$ ，这个k使我们要确定的，它代表着在所有的 $\mathbf{(i,j)}$ 对中有多少对满足 $\mathbf{gcd(i,j)=d}$ ，也就是说 $\mathbf{k=\sum_{i=1}^n\sum_{j=1}^n[gcd(i,j)=d]}$ ，然后根据结论基础数论学习笔记(下)第十部分6.18我们能立刻得到 $\mathbf{k=2sum\varphi(\lfloor \frac nd\rfloor )-1}$ ， $\mathbf{sum\varphi(n)}$ 代表 $\mathbf{\varphi}$ 的1~n的前缀和。
然后剩下的工作就简单了，我们要计算的分母其实就是这么一个式子：
$\mathbf{\prod_{d=1}^nd^{2sum\varphi(\lfloor \frac nd\rfloor)-1}}$ ,
显然我们可以对这个式子的 $\mathbf{2sum\varphi(\lfloor \frac nd\rfloor)-1}$ 这个部分进行分块处理，然后由于这题卡空间，暴力算分块过程中l到r的阶乘即可，线性筛复杂度为 $\mathbf{O(n)}$ ，分块由于要暴力算阶乘复杂度仍为 $\mathbf{O(n)}$ ，故总复杂度为 $\mathbf{O(n)}$ ，然后要注意空间，本题空间卡得很紧，线性筛的时候标记数组最好开成bitset。

const int mod = 104857601;
inline void add(int &x,int y){x+=y,x>=(mod-1)&&(x-=mod-1);}
int n,prim[maxn],tot,phi[maxn];
void init(int n){
	bitset<maxn>flag;
	phi[1]=1;
	FOR(i,2,n+1){
		if(!flag[i])prim[tot++]=i,phi[i]=i-1;
		for(register int j =0;j<tot && prim[j]*i<=n;++j){
			flag[i*prim[j]]=1;
			if(i%prim[j]==0){
				phi[i*prim[j]]=phi[i]*prim[j];
				break;
			}
			phi[i*prim[j]]=phi[i]*(prim[j]-1);
		}
	}
	FOR(i,1,n+1)add(phi[i],phi[i-1]);
}
int main(){
	rd(&n);
	init(n);
	register int i=1,j,c,a=1,b=1;
	while(i<=n){
		j=n/(n/i),c=1;
		FOR(k,i,j+1)c=1ll*c*k%mod,a=1ll*a*k%mod;
		b=1ll*b*qpow(c,4*phi[n/i]-2,mod)%mod;
		i=j+1;
	}
	a=qpow(a,2*n,mod);
	b=1ll*a*qpow(b,mod-2,mod)%mod;
	wrn(b);
}