【数论】莫比乌斯反演巩固1

最新推荐文章于 2024-05-20 15:29:54 发布

conti123

最新推荐文章于 2024-05-20 15:29:54 发布

阅读量1k

点赞数 22

分类专栏：数论文章标签： c++ 洛谷笔记数论

本文链接：https://blog.csdn.net/conti123/article/details/136449739

版权

数论专栏收录该内容

4 篇文章 0 订阅

订阅专栏

文章目录

前言

这里需要对莫反有一些基础。
不会的可以点这里

洛谷P3455

题意：
给定 $a, b, d$
求 $\sum_{x=1}^a\sum_{y=1}^b[\gcd(x,y)=d]$
$1\leq n\leq5\times 10^4$ ， $\leq d \leq a,b \leq 5 \times 10^4$
按照套路：
$\sum_{x=1}^a\sum_{y=1}^b[\gcd(x,y)=d]$
$=\sum_{x=1}^{\frac{a}{d}}\sum_{y=1}^{\frac{b}{d}}[\gcd(x,y)=1]$
$=\sum_{x=1}^{\frac{a}{d}}\sum_{y=1}^{\frac{b}{d}}\sum_{k|\gcd(x,y)}\mu(k)$
$=\sum_{x=1}^{\frac{a}{d}}\sum_{y=1}^{\frac{b}{d}}\sum_{k=1}\mu(k)[k|\gcd(x,y)]$
$=\sum_{k=1}^{\min(\frac{a}{d},\frac{b}{d})}\mu(k)\sum_{x=1}^{\frac{a}{d}}\sum_{y=1}^{\frac{b}{d}}[k|\gcd(x,y)]$
$=\sum_{k=1}^{\min(\frac{a}{d},\frac{b}{d})}\mu(k)\left\lfloor\frac{\frac{a}{d}}{k}\right\rfloor\left\lfloor\frac{\frac{b}{d}}{k}\right\rfloor$

数论分块解决即可。
Code：

#include<bits/stdc++.h>
#define int long long
#define IOS ios::sync_with_stdio(false),cin.tie(NULL),cout.tie(NULL)
using namespace std;
const int N=2e6+1,M=2e6+1;
int t,n,m,k,ans=0;
struct fy{
	int prv[N],cnt,mu[N],sum[N];
	bool pr[N];
	void ola(int x){
		pr[1]=mu[1]=1;
		for(int i=2;i<=x;i++){
			if(!pr[i])
				prv[++cnt]=i,mu[i]=-1;
			for(int j=1;j<=cnt&&i*prv[j]<=x;j++){
				int u=i*prv[j];
				pr[u]=1;
				if(i%prv[j]==0){
					mu[u]=0;
					break;
				}
				else{
					mu[u]=-mu[i];
				}
			}
		}
		for(int i=1;i<=x;i++)
			sum[i]=sum[i-1]+mu[i];
			
	}
}A;
signed main(){
	IOS;
	A.ola(N-1);
	cin>>t;
	while(t--){
		cin>>n>>m>>k;
		n/=k,m/=k;
		ans=0;
		for(int l=1,r;l<=min(n,m);l=r+1){
			r=min(n/(n/l),m/(m/l));
			ans+=(A.sum[r]-A.sum[l-1])*(n/l)*(m/l);
		}
		cout<<ans<<"\n";
	}
	return 0;
}

洛谷P1829

题意：给定 $n, m$ ，求 $\sum_{i=1}^n\sum_{j=1}^m lcm(i,j)$
$1\leq n,m\leq 10^7$

$\sum_{i=1}^n\sum_{j=1}^m lcm(i,j)$
$=\sum_{i=1}^n\sum_{j=1}^m\frac{ij}{\gcd(i,j)}$
令 $o=\min(n,m)$
$=\sum_{d=1}^o\sum_{i=1}^n\sum_{j=1}^m[\gcd(i,j)=d]\frac{ij}{d}$
令 $a=\left\lfloor\frac{n}{d}\right\rfloor,b=\left\lfloor\frac{m}{d}\right\rfloor,c=\min(a,b)$
$=\sum_{d=1}^o d\sum_{i=1}^a\sum_{j=1}^b[\gcd(i,j)=1]ij$
$=\sum_{d=1}^o d\sum_{i=1}^a\sum_{j=1}^bij\sum_{k|\gcd(i,j)}\mu(k)$
$=\sum_{d=1}^o d\sum_{k=1}^c\mu(k)\sum_{i=1}^a\sum_{j=1}^bij[k|\gcd(i,j)]$
$=\sum_{d=1}^o d\sum_{k=1}^c\mu(k)\sum_{i=1}^{\left\lfloor\frac{a}{k}\right\rfloor}\sum_{j=1}^{\left\lfloor\frac{b}{k}\right\rfloor}ik\cdot jk$
$=\sum_{d=1}^o d\sum_{k=1}^c\mu(k)\sum_{i=1}^{\left\lfloor\frac{a}{k}\right\rfloor}\sum_{j=1}^{\left\lfloor\frac{b}{k}\right\rfloor}ijk^2$
$=\sum_{d=1}^o d\sum_{k=1}^c\mu(k)k^2\sum_{i=1}^{\left\lfloor\frac{a}{k}\right\rfloor}i\sum_{j=1}^{\left\lfloor\frac{b}{k}\right\rfloor}j$
令 $sum(n)=\sum_{i=1}^ni,T=kd$
$=\sum_{d=1}^o d\sum_{k=1}^c\mu(k)k^2sum(\left\lfloor\frac{a}{k}\right\rfloor)sum(\left\lfloor\frac{b}{k}\right\rfloor)$
$=\sum_{d=1}^o d\sum_{k=1}^{\min(\left\lfloor\frac{n}{d}\right\rfloor,\left\lfloor\frac{m}{d}\right\rfloor)}\mu(k)k^2sum(\left\lfloor\frac{n}{kd}\right\rfloor)sum(\left\lfloor\frac{m}{kd}\right\rfloor)$
$=\sum_{T=1}^osum(\frac{n}{T})sum(\frac{m}{T})T\sum_{d|T}d\mu(d)$
令 $g(n)=\sum_{d|n}d\mu(d)$
$=\sum_{T=1}^osum(\frac{n}{T})sum(\frac{m}{T})Tg(T)$

考虑线性筛出 $g (n)$
还可以用数论分块优化，但也可以不用。
Code：

#include<bits/stdc++.h>
#define int long long
#define IOS ios::sync_with_stdio(false),cin.tie(NULL),cout.tie(NULL)
using namespace std;
const int N=1e7+1,mod=20101009;
int t,n,m,k;
struct fy{
	int prv[N],cnt,g[N],s[N];
	bool pr[N];
	void ola(int x){
		pr[1]=g[1]=1;
		for(int i=2;i<=x;i++){
			if(!pr[i])
				prv[++cnt]=i,g[i]=(1-i+mod)%mod;
			for(int j=1;j<=cnt&&i*prv[j]<=x;j++){
				int u=i*prv[j];
				pr[u]=1;
				if(i%prv[j]==0){
					g[u]=g[i];
					break;
				}
				else{
					g[u]=g[i]*g[prv[j]]%mod;
				}
			}
		}
	}
	void getsum(int x){
		for(int i=1;i<=x;i++)
			s[i]=(s[i-1]+i)%mod;
	}
}A;
signed main(){
	IOS;
	A.ola(N-1);
	A.getsum(N-1);
	cin>>n>>m;
	int ans=0;
	for(int i=1;i<=min(n,m);i++){
		ans+=A.s[n/i]*A.s[m/i]%mod*i%mod*A.g[i]%mod;
		ans%=mod;
	}
	cout<<ans;
	return 0;
}
#include<bits/stdc++.h>
#define int long long
#define IOS ios::sync_with_stdio(false),cin.tie(NULL),cout.tie(NULL)
using namespace std;
const int N=1e7+1,mod=20101009;
int t,n,m,k;
struct fy{
	int prv[N],cnt,g[N],s[N];
	bool pr[N];
	void ola(int x){
		pr[1]=g[1]=1;
		for(int i=2;i<=x;i++){
			if(!pr[i])
				prv[++cnt]=i,g[i]=(1-i+mod)%mod;
			for(int j=1;j<=cnt&&i*prv[j]<=x;j++){
				int u=i*prv[j];
				pr[u]=1;
				if(i%prv[j]==0){
					g[u]=g[i];
					break;
				}
				else{
					g[u]=g[i]*g[prv[j]]%mod;
				}
			}
		}
	}
	void getsum(int x){
		for(int i=1;i<=x;i++)
			s[i]=(s[i-1]+i)%mod,g[i]=(i*g[i]+g[i-1])%mod;
	}
}A;
signed main(){
	IOS;
	A.ola(N-1);
	A.getsum(N-1);
	cin>>n>>m;
	int ans=0;
	for(int l=1,r;l<=min(n,m);l=r+1){
		r=min(n/(n/l),m/(m/l));
		ans+=A.s[n/l]*A.s[m/l]%mod*((A.g[r]-A.g[l-1]+mod)%mod)%mod;
		ans%=mod;
	}
	cout<<ans;
	return 0;
}

洛谷P4449

题意：给定 $n, m, k$ ，求 $\sum_{i=1}^n\sum_{j=1}^m\gcd(i,j)^k$
$\leq T \leq 2 \times 10^3,1\leq n, m, k \leq 5\times10^6$
$\sum_{i=1}^n\sum_{j=1}^m\gcd(i,j)^k$
$a, b, c, o$ 的含义同上文。
$=\sum_{d=1}^od^k\sum_{i=1}^n\sum_{j=1}^m[\gcd(i,j)=d]$
$=\sum_{d=1}^od^k\sum_{p=1}^c\mu(p)\left\lfloor\frac{a}{p}\right\rfloor\left\lfloor\frac{b}{p}\right\rfloor$
令 $T = d p$
$=\sum_{d=1}^od^k\sum_{p=1}^c\mu(p)\left\lfloor\frac{n}{T}\right\rfloor\left\lfloor\frac{n}{T}\right\rfloor$
$=\sum_{T=1}^o\left\lfloor\frac{n}{T}\right\rfloor\left\lfloor\frac{n}{T}\right\rfloor\sum_{d|T}d^k\mu(\frac{T}{d})$

就跟上一题一样了，线性筛+数论分块。
Code：

#include<bits/stdc++.h>
#define int long long
#define IOS ios::sync_with_stdio(false),cin.tie(NULL),cout.tie(NULL)
using namespace std;
const int N=5e6+1,mod=1e9+7;
int t,n,m,k;
struct fy{
	int prv[N],cnt,g[N],s[N];
	bool pr[N];
	inline int qmi(int x,int y){
		int res=1;
		while(y>0){
			if(y&1)
				res=res*x%mod;
			x=x*x%mod,y>>=1;
		}
		return res;
	}
	void ola(int x){
		pr[1]=g[1]=1;
		for(int i=2;i<=x;i++){
			if(!pr[i])
				prv[++cnt]=i,g[i]=(qmi(i,k)-1+mod)%mod;
			for(int j=1;j<=cnt&&i*prv[j]<=x;j++){
				int u=i*prv[j];
				pr[u]=1;
				if(i%prv[j]==0){
					g[u]=g[i]*qmi(prv[j],k)%mod;
					break;
				}
				else{
					g[u]=g[i]*g[prv[j]]%mod;
				}
			}
		}
	}
	void getsum(int x){
		for(int i=1;i<=x;i++)
			g[i]=(g[i]+g[i-1])%mod;
	}
}A;
signed main(){
	IOS;
	cin>>t>>k;
	A.ola(N-1);
	A.getsum(N-1);
	while(t--){
		cin>>n>>m;
		int ans=0;
		for(int l=1,r;l<=min(n,m);l=r+1){
			r=min(n/(n/l),m/(m/l));
			ans+=(n/l)*(m/l)%mod*((A.g[r]-A.g[l-1]+mod)%mod)%mod;
			ans%=mod;
		}
		cout<<ans<<"\n";
	}
	return 0;
}