莫比乌斯反演

学习课件：popoqqq课件，可以百度下载

Gcd

Bzoj2818:
题解来源：https://www.cnblogs.com/Emiya-wjk/p/9991997.html
题意：给定整数N(1<=N<=10000000)，求1<=x,y<=N,gcd(x,y)为素数的对数
题解：

#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdio>
#include<iomanip>
#include<cstdlib>
#define MAXN 0x7fffffff
typedef long long LL;
const int N=10000005;
using namespace std;
inline int Getint(){register int x=0,f=1;register char ch=getchar();while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}while(isdigit(ch)){x=x*10+ch-'0';ch=getchar();}return x*f;}
int mu[N],prime[N]; 
bool vis[N];

int main(){
    int n=Getint();
    mu[1]=1;
    for(int i=2;i<n;i++){
        if(!vis[i])prime[++prime[0]]=i,mu[i]=-1;
        for(int j=1;j<=prime[0]&&i*prime[j]<n;j++){
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)break;
            mu[i*prime[j]]=-mu[i];
        }
    }
    LL ans=0;
    for(int i=1;i<=prime[0];i++)
        for(int j=1;j*prime[i]<=n;j++)
            ans+=(LL)mu[j]*(n/j/prime[i])*(n/j/prime[i]);
    cout<<ans;
    return 0;
}

能量采集

noi 2010
题解来源：https://www.cnblogs.com/emiya-wjk/p/9996805.html

#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdio>
#include<iomanip>
#include<cstdlib>
#define MAXN 0x7fffffff
typedef long long LL;
const int N=100005;
using namespace std;
inline int Getint(){register int x=0,f=1;register char ch=getchar();while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}while(isdigit(ch)){x=x*10+ch-'0';ch=getchar();}return x*f;}
int mu[N],prime[N];
bool vis[N];
int main(){
    int n=Getint(),m=Getint();if(n>m)swap(n,m); 
    mu[1]=1;
    for(int i=2;i<=n;i++){
        if(!vis[i])prime[++prime[0]]=i,mu[i]=-1;
        for(int j=1;j<=prime[0]&&1ll*i*prime[j]<=n;j++){
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)break;
            mu[i*prime[j]]=-mu[i];
        }
    }
    LL ans=0;
    for(int d=1;d<=n;d++){
        LL ret=0;
        for(int x=1;x*d<=n;x++)
            ret+=1ll*mu[x]*(n/x/d)*(m/x/d);
        ans+=ret*d;
    }
    cout<<2*ans-1ll*n*m;
    return 0;
}

Problem b

Bzoj 2301
题解来源：https://blog.csdn.net/outer_form/article/details/50590197
题意：n次询问，每次询问给出a,b,c,d,k，求在a<=x<=b,c<=y<=d，满足gcd(x,y)是k的(x,y)的对数。
令g(n,m,k)表示在1<=x<=n,1<=y<=m，满足gcd(x,y)是k的(x,y)的对数。 那么由容斥原理可得
$ans=g(c,d,k)--g(a-1,d,k)-g(b,c-1,k)+g(a-1,c-1,k)$
其中$g(n,m,k)=g(n/k,m/k,1)$，我们根据gcd那道，
令f[i]为$gcd(x,y)==i$的对数，$g[i]$为$gcd(x,y)==x\ast i$的对数，$x\in Z$，$g[i]=(n/i)\ast (m/i)$，$g[x]={\sum }_{x|d}^{n}f[d]$，反演后
$f[x]={\sum }_{x|d}^{n}u(d/x)\ast g[d]={\sum }_{x|d}^{n}u(d/x)\ast \lfloor n/x\rfloor \ast \lfloor m/x\rfloor$
当x=1时，$f[1]={\sum }_{d=1}^{min(n,m)}u(d)\lfloor n/d\rfloor \lfloor m/d\rfloor$，由于$\lfloor n/x\rfloor$取值最多只有$2\sqrt{n}$个，我们算出$u$的前缀和sum，然后只需要在$O(2(\sqrt{n}+\sqrt{m}))$的时间回答每次查询即可。

#include<iostream>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<cstdio>
#include<iomanip>
#include<cstdlib>
#define MAXN 0x7fffffff
typedef long long LL;
const int N=50005;
using namespace std;
inline int Getint(){register int x=0,f=1;register char ch=getchar();while(!isdigit(ch)){if(ch=='-')f=-1;ch=getchar();}while(isdigit(ch)){x=x*10+ch-'0';ch=getchar();}return x*f;}
int mu[N],prime[N],sum[N];
bool vis[N];
int Cal(int n,int m)
{
	if(n>m) swap(n,m);
	int last,res=0;
	for(int i=1;i<=n;i=last+1){
		last=min(n/(n/i),m/(m/i));
		res+=(n/i)*(m/i)*(sum[last]-sum[i-1]);
	}
	return res;
}
int n,a,b,c,d,k;
int main(){
    n=50000;
    sum[1]=mu[1]=1;sum[0]=0;
    for(int i=2;i<=n;i++){
        if(!vis[i])prime[++prime[0]]=i,mu[i]=-1;
        for(int j=1;j<=prime[0]&&1ll*i*prime[j]<=n;j++){
            vis[i*prime[j]]=1;
            if(i%prime[j]==0)break;
            mu[i*prime[j]]=-mu[i];
        }
        sum[i]=mu[i]+sum[i-1];
    }
    n=Getint();
    for(int i=1;i<=n;i++){
    	a=Getint();b=Getint();c=Getint();d=Getint();
    	k=Getint();a--;c--;
    	int ans=Cal(b/k,d/k)+Cal(a/k,c/k)-Cal(b/k,c/k)-Cal(a/k,d/k);
    	cout<<ans<<endl;
	}
    return 0;
}

YY的GCD

Bzoj 2820
题解来源：http://hzwer.com/6142.html


可以线筛预处理F,按照素数粗略个数n/logn以及调和级数求和nlogn来看暴力的复杂度接近On，处理完F以后就是喜闻乐见的下底函数分块。
莫得提交，代码正确性待验证。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=50005;

inline int read()
{
	int res=0;
	char ch;
	while((ch=getchar())<'0'||ch>'9');
	while(ch>='0'&&ch<='9') res=(res<<3)+(res<<1)+ch-'0',ch=getchar();
	return res; 
}
int mu[maxn],prime[maxn];
int n,m;
bool vis[maxn];
ll f[maxn];
void getphi()
{
	memset(vis,0,sizeof(vis));
	memset(mu,0,sizeof(mu));
	mu[1]=1;
	for(int i=2;i<maxn;i++){
		if(!vis[i]) prime[++prime[0]]=i,mu[i]=-1;
		for(int j=1;j<=prime[0]&&1LL*i*prime[j]<maxn;j++){
			vis[i*prime[j]]=1;
			if(i%prime[j]==0) break;
			mu[i*prime[j]]=-mu[i];
		}
	}
	memset(f,0,sizeof(f));
	for(int i=1;i<=prime[0];i++){
		int p=prime[i];
		for(int j=1;j*p<maxn;j++) f[j*p]+=mu[j];
	}
	for(int i=1;i<maxn;i++) f[i]+=f[i-1];
}
int main(){
	getphi();
    int t=read();
    while(t--){
    	ll ans=0;
    	n=read();m=read();
    	if(n>m) swap(n,m);
    	for(int i=1,j;i<=n;i=j+1){
    		j=min(n/(n/i),m/(m/i));
    		ans+=(f[j]-f[i-1])*(n/i)*(m/i);
		}
		printf("%lld\n",ans);
	}
    return 0;
}

数表

Bzoj 3529
参考代码：https://www.cnblogs.com/gzez181027/p/bzoj3529.html
题解：



PS:求一个数x的约数和，令$x=p_1^{a_1}\ast p_2^{a_2}\ast ...\ast p_n^{a_n}$，那么x的约数和为
$sum=(p_1^0+p_1^1+..+p_1^{a_1})\ast (p_2^0+p_2^1+..+p_2^{a_2})\ast ...\ast (p_n^0+..+p_n^{a_n})=(\frac{p_1^{a_1+1}-1}{p_1-1})\ast (\frac{p_2^{a_2+1}-1}{p_2-1})\ast ...\ast (\frac{p_n^{a_n+1}-1}{p_n-1})$
求数的约数和除了用上面提到的方法，这题也可以用朴素的暴力方法，有线性的方法，但弱鸡看不懂。。

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int maxn=100001;
const int maxm=20010;
inline int read()
{
	int res=0;
	char ch;
	while((ch=getchar())<'0'||ch>'9');
	while(ch>='0'&&ch<='9') res=(res<<3)+(res<<1)+ch-'0',ch=getchar();
	return res; 
}
int mu[maxn],prime[maxn];
int n,m,a;
bool vis[maxn];
struct node{
	int id,val;
}f[maxn];
bool cmp(const node &x,const node &y)
{
	return x.val<y.val;
}
struct Query{
	int id,n,m,a;
}q[maxm];
int ans[maxm];
bool cmp2(const Query &x,const Query &y)
{
	return x.a<y.a;
}
void getphi()
{
	memset(vis,0,sizeof(vis));
	memset(mu,0,sizeof(mu));
	mu[1]=1;
	for(int i=2;i<maxn;i++){
		if(!vis[i]) prime[++prime[0]]=i,mu[i]=-1;
		for(int j=1;j<=prime[0]&&1LL*i*prime[j]<maxn;j++){
			vis[i*prime[j]]=1;
			if(i%prime[j]==0) break;
			mu[i*prime[j]]=-mu[i];
		}
	}
	for(int i=1;i<maxn;i++) f[i].id=i,f[i].val=0;
	for(int i=1;i<maxn;i++)
		for(int j=i;j<maxn;j+=i)
			f[j].val+=i; 
	sort(f+1,f+maxn,cmp);
}
int s[maxn];
int lowbit(int x)
{
	return x&(-x);
}
void Add(int x,int val)
{
	for(;x<maxn;x+=lowbit(x))
		s[x]+=val;
	return;
}
int query(int x)
{
	int res=0;
	for(;x;x-=lowbit(x)) res+=s[x];
	return res;
}
int solve(int n,int m)
{
	int res=0;
	if(n>m) swap(n,m);
	for(int i=1,j;i<=n;i=j+1){
		j=min(n/(n/i),m/(m/i));
		res+=(query(j)-query(i-1))*(n/i)*(m/i);
	}
	return res&0x7fffffff;//这里模为（1LL<<31）等价下面 减法取模
//	return (res+(1LL<<31))%(1LL<<31);
}
int main()
{
	getphi();
	scanf("%d",&m);
	for(int i=1;i<=m;i++){
		scanf("%d%d%d",&q[i].n,&q[i].m,&q[i].a);
		q[i].id=i;
	}
	sort(q+1,q+m+1,cmp2);
	for(int i=1,j=1;i<=m;i++){
		for(;j<maxn&&f[j].val<=q[i].a;j++)
			for(int k=f[j].id;k<maxn;k+=f[j].id)
				Add(k,f[j].val*mu[k/f[j].id]);
		ans[q[i].id]=solve(q[i].n,q[i].m);
	}
	for(int i=1;i<=m;i++) printf("%d\n",ans[i]);
	return 0;
}

Crash的数字表格

Bzoj 2154
题意：求${\sum }_{i=1}^n{\sum }_{j=1}^{m}lcm(i,j),n,m<=1e7$
题解&代码：https://www.cnblogs.com/iwtwiioi/p/4268926.html
题解：
下面设n<m
${\sum }_{i=1}^n{\sum }_{j=1}^{m}lcm(i,j)={\sum }_{i=1}^n{\sum }_{j=1}^m\frac{ij}{(i,j)}$
然后枚举$d=(i,j)$，那么确定i,j怎么取，显然我们只只需要先除去i和j的d，也就是$(i/d,j/d)=1$就行了，那么设$F(x,y)={\sum }_{i-1}^x{\sum }_{j=1}^{y}ij[(i,j)==1]$，那么原式变成${\sum }_{d=1}^n\frac{d^{2}F(\lfloor \frac{n}{d}\rfloor ,\lfloor \frac{m}{d}\rfloor )}{d}$
考虑$F(x,y)$
$$F(x,y)=\sum _{i-1}^x\sum _{j=1}^{y}ij[(i,j)==1]$$
$$=\sum _{i=1}^x\sum _{j=1}^{y}ij\sum _{d|(i,j)}u(d)$$
$$=\sum _{d=1}^{x}u(d)\sum _{d|i}^{x}i\sum _{d|j}^{y}j$$
$$=\sum _{d=1}^{x}u(d)d^2\sum _{i=1}^{\lfloor \frac{x}{d}\rfloor }i\sum _{j=1}^{\lfloor \frac{y}{d}\rfloor }j$$
$$=\sum _{d=1}^{x}u(d)d^2\frac{\lfloor \frac{x}{d}\rfloor (\lfloor \frac{x}{d}\rfloor +1)}{2}\frac{\lfloor \frac{y}{d}\rfloor (\lfloor \frac{y}{d}\rfloor +1)}{2}$$
带回原式得：
$$\sum _{d=1}^{n}F(\lfloor \frac{n}{d}\rfloor ,\lfloor \frac{m}{d}\rfloor )$$
$$=\sum _{d=1}^{n}d\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }u(i)i^2\frac{\lfloor \frac{\lfloor \frac{n}{d}\rfloor }{i}\rfloor (\lfloor \frac{\lfloor \frac{n}{d}\rfloor }{i}\rfloor +1)}{2}\frac{\lfloor \frac{\lfloor \frac{m}{d}\rfloor }{i}\rfloor (\lfloor \frac{\lfloor \frac{m}{d}\rfloor }{i}\rfloor +1)}{2}$$
$$=\sum _{d=1}^{n}d\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }u(i)i^2\frac{\lfloor \frac{n}{di}\rfloor (\lfloor \frac{n}{di}\rfloor +1)}{2}\frac{\lfloor \frac{m}{di}\rfloor (\lfloor \frac{m}{di}\rfloor +1)}{2}$$
现在可以$O(\sqrt{n}\sqrt{n})=O(n)$查询了，可以继续化简，令$T=di$，则$i|T,d=T/i$，换掉指标，得
$$=\sum _{d=1}^{n}d\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }u(i)i^2\frac{\lfloor \frac{n}{di}\rfloor (\lfloor \frac{n}{di}\rfloor +1)}{2}\frac{\lfloor \frac{m}{di}\rfloor (\lfloor \frac{m}{di}\rfloor +1)}{2}$$
$$=\sum _{T=1}^n\frac{\lfloor \frac{n}{T}\rfloor (\lfloor \frac{n}{T}\rfloor +1)}{2}\frac{\lfloor \frac{m}{T}\rfloor (\lfloor \frac{m}{T}\rfloor +1)}{2}\sum _{i|T}\frac{T}{i}u(i)i^2$$
$$=\sum _{T=1}^n\frac{\lfloor \frac{n}{T}\rfloor (\lfloor \frac{n}{T}\rfloor +1)}{2}\frac{\lfloor \frac{m}{T}\rfloor (\lfloor \frac{m}{T}\rfloor +1)}{2}T\sum _{i|T}u(i)i$$
设$g(T)=T{\sum }_{i|T}u(i)i$，$f(T)={\sum }_{i|T}u(i)i$,则$g(T)=T\ast f(T)$，考虑求$f(T)$，在线性筛中，外层为k，内层为$p_y$，求$f(k{p}_y)={\sum }_{i|k{p}_y}u(i)i$
当$p_y|k$时，如果i取的因子不包含新加入的$p_y$，答案为$f(k)$；当i取的因子包含新加入的$p_y$时，此时$p_y$已经$>=2$，所以$u(i)=0$，贡献为0，综上，当$p_y|k$答案为$f(k)$
当$p_y\nmid k$时，当i取的因子不包含新加入的$p_y$，答案为$f(k)$；当i取的因子包含新加入的$p_y$时，指数为1,考虑$i=a{p}_y$
则$$\sum _{i|T}u(i)i=\sum _{a{p}_y|k{p}_y}u(a{p}_y)a{p}_y$$
$$=p_y\sum _{a|k}u(a)u(p_y)a$$
$$=-p_y\sum _{a|k}u(a)a=-p_{y}f(k)$$
综上，当$p_y\nmid k$时，答案为$(1-p_y)f(k)$
预处理求出$f(T)$，进而求出$g[T)$，最后答案为
$$\sum _{T=1}^n\frac{\lfloor \frac{n}{T}\rfloor (\lfloor \frac{n}{T}\rfloor +1)}{2}\frac{\lfloor \frac{m}{T}\rfloor (\lfloor \frac{m}{T}\rfloor +1)}{2}g(T)$$
再预处理$g(T)$的前缀和，分块查询即可

#include<bits/stdc++.h>
using namespace std;
typedef long long ll;
const int mod=20101009;
const int maxn=1e7+10;
 
inline int read()
{
    int res=0;
    char ch;
    while((ch=getchar())<'0'||ch>'9');
    while(ch>='0'&&ch<='9') res=(res<<3)+(res<<1)+ch-'0',ch=getchar();
    return res; 
}
int prime[maxn];
int n,m;
bool vis[maxn];
ll g[maxn];
void getphi()
{
//  memset(vis,0,sizeof(vis));
    g[1]=1;
    for(int i=2;i<maxn;++i){
        if(!vis[i]) prime[++prime[0]]=i,g[i]=1-i;
        for(int j=1;j<=prime[0];j++){//&&1LL*i*prime[j]<maxn 
            int tmp=i*prime[j];if(tmp>=maxn) break;
            vis[tmp]=1;
            if(i%prime[j]==0){
                g[tmp]=g[i];
                break;
            }
            g[tmp]=g[i]*(1-prime[j]);//%mod
        }
    }
//  for(int i=2;i<maxn;++i) g[i]=g[i]*i%mod;
    for(int i=2;i<maxn;++i) g[i]=g[i]*i;
    for(int i=1;i<maxn;++i) g[i]+=g[i-1],g[i]%=mod;
}
 
int main(){
    getphi();
    int t=1;
    while(t--){
        scanf("%d%d",&n,&m);if(n>m) swap(n,m);
        ll ans=0,t1,t2;
        for(int i=1,j;i<=n;i=j+1){
            j=min(n/(n/i),m/(m/i));
            t1=(1LL*(n/i)*(n/i+1)/2)%mod;
            t2=(1LL*(m/i)*(m/i+1)/2)%mod;
            
            ans+=(g[j]-g[i-1]+mod)%mod*(t1*t2%mod)%mod;
            ans%=mod;
        }
        printf("%lld\n",((ans%mod)+mod)%mod);
    }
    return 0;
}

jzptab

Bzoj 2693
和上题相同，多次查询，上述代码可以解决。