P5176 公约数 [莫比乌斯反演][线性筛积性函数]

传送门

来自zxy学长, 强烈推荐哦 

 

 

 

 

 

关于最后一个的证明:

首先p为n的最小质因子, 那么np一定含有p因子的平方项, 如果保留的话对mu是没有贡献的

也就是说对f(n)有贡献的d, 只有d*p对 f(n*p)有贡献


关于线性筛函数, 考虑

等于1

为质数  ->  代入两个端点

不包含最小质因子  ->  积性函数的性质

包含最小质因子  ->  讨论mu的贡献...


#include<bits/stdc++.h>
#define N 20000050
#define Mod 1000000007
#define LL long long
using namespace std;
int prim[N], tot;
bool isp[N]; LL g[N];
int T;
void prework(){
	g[1] = 1;
	for(int i=2;i<=N-50;i++){
		if(!isp[i]) prim[++tot] = i, g[i] = (LL)i * (LL)i % Mod - 1;
		for(int j=1;j<=tot;j++){
			if(prim[j] * i > N - 50) break;
			isp[prim[j] * i] = 1;
			if(i % prim[j] == 0){
				g[i * prim[j]] = g[i] * (LL)prim[j] % Mod * (LL)prim[j] % Mod;
				break;
			}
			g[i * prim[j]] = g[i] * g[prim[j]] % Mod;
		}
	}
	for(int i=2;i<=N-50;i++) g[i] = (g[i] + g[i-1]) % Mod;
}
LL Solve(int n,int m){
	if(n>m) swap(n,m);
	LL ans = 0;
	for(int l=1,r;l<=n;l=r+1){
		int v1 = n/l, v2 = m/l;
		r = min(n/v1, m/v2);
		ans += (LL)v1 * (LL)v2 % Mod * (g[r] - g[l-1]) % Mod;
		ans = (ans % Mod + Mod) % Mod;
	} return ans;
}
int main(){
	prework(); scanf("%d",&T);
	while(T--){
		int n,m,p; scanf("%d%d%d",&n,&m,&p);
		LL ans = (LL)p * Solve(n,m) % Mod;
		ans = (ans + ((LL)n * Solve(m,p)) % Mod) % Mod;
		ans = (ans + ((LL)m * Solve(n,p)) % Mod) % Mod;
		printf("%lld\n",ans);
	} return 0;
}

 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

FSYo

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值