f(n)表示某一范围内(x,y)=n的数对的数量,F(n)表示某一范围内n|(x,y)的数对的数量
那么直接求f(n)并不是很好求,而F(n)求起来相对无脑一些,我们可以通过对F(n)进行莫比乌斯反演来求得f(n)
线性筛μ模板:
void Prime(int N)
{
mu[1]=1;int cnt=0;
for(int i=2;i<=N;i++)
{
if(!v[i]) p[++cnt]=i,mu[i]=-1;
for(int j=1,k;j<=cnt&&p[j]*i<=N;j++)
{
v[k=p[j]*i]=1;
if(i%p[j]==0) {mu[k]=0;break;}
mu[k]=-mu[i];
}
}
}
基本解题思路:枚举gcd
例题1:
求1<=i<=n, 1<=j<=m, gcd(i,j)=k 的(i,j)对数
该问题等价于求1<=i<=n/k (N), 1<=j<=m/k (M) gcd(i,j)=1 的(i,j)对数
由的取值是级的,可以得出:
优化技巧一: 前缀和+分块优化
long long ans=0;
for(int i=1,j;i<=n;i=j+1)
{
j=min(n/(n/i),m/(m/i));
ans+=1ll*(n/i)*(m/i)*(sum[j]-sum[i-1]);
}
例题2:
求1<=i<=n, 1<=j<=m, gcd(i,j)为质数 的(i,j)对数
优化技巧二: 切换枚举次序
令上式的pd=T, 则有
应用技巧一即可
例题3:
设为的约数个数,给定N,M,T
求
N,M,T<=50000
题目解析见这里:约数个数和 (莫比乌斯反演)
于是算不上技巧但是很神奇的技巧三:
技巧三:替换条件框
例题3 代码
#include<cstdio>
#include<algorithm>
#define maxn 50005
using namespace std;
int T,n,m,p[maxn],mu[maxn],sum[maxn],mx[maxn];//线性筛约数个数
bool v[maxn];
void Prime(int N)
{
mu[1]=sum[1]=1;int cnt=0;
for(int i=2;i<=N;i++)
{
if(!v[i]){
p[++cnt]=i,mu[i]=-1;
sum[i]=mx[i]=2;
}
for(int j=1,k;j<=cnt&&p[j]*i<=N;j++)
{
v[k=p[j]*i]=1;
if(i%p[j]==0){
mu[k]=0;
mx[k]=mx[i]+1;
sum[k]=sum[i]/mx[i]*mx[k];
break;
}
mu[k]=-mu[i];
mx[k]=2;
sum[k]=sum[i]*2;
}
}
for(int i=1;i<=N;i++) mu[i]+=mu[i-1],sum[i]+=sum[i-1];
}
int main()
{
Prime(maxn-5);
scanf("%d",&T);
while(T--)
{
scanf("%d%d",&n,&m);if(n>m) swap(n,m);
long long ans=0;
for(int i=1,j;i<=n;i=j+1)
{
j=min(n/(n/i),m/(m/i));
ans+=1ll*sum[n/i]*sum[m/i]*(mu[j]-mu[i-1]);
}
printf("%lld\n",ans);
}
}