Luogu2257 YY的GCD

原题链接：https://www.luogu.org/problemnew/show/P2257

YY的GCD

题目描述

神犇YY虐完数论后给傻×kAc出了一题

给定N, M,求1<=x<=N, 1<=y<=M且gcd(x, y)为质数的(x, y)有多少对

kAc这种傻×必然不会了，于是向你来请教……

多组输入

输入输出格式

输入格式：

第一行一个整数T 表述数据组数

接下来T行，每行两个正整数，表示N, M

输出格式：

T行，每行一个整数表示第i组数据的结果

输入输出样例

输入样例#1：

2
10 10
100 100

输出样例#1：

30
2791

说明

T = 10000

N, M <= 10000000

题解

日常推式子。。。

设$p$为素数，$P$表示素数集合：

$$ans=\sum_P\sum_{i=1}^{n}\sum_{j=1}^{m}[gcd(i,j)=p]$$

日常提出$p$：

$$ans=\sum_P\sum_{i=1}^{\lfloor\frac{n}{p}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{p}\rfloor}[gcd(i,j)=1]$$

日常替换$[gcd(i,j)=1]$：

$$ans=\sum_P\sum_{i=1}^{\lfloor\frac{n}{p}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{p}\rfloor}\sum_{d|gcd(i,j)}\mu(d)$$

日常枚举$d$：

$$ans=\sum_{d|gcd(i,j)}\mu(d)\sum_P\sum_{i=1}^{\lfloor\frac{n}{p}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{p}\rfloor}[d|gcd(i,j)]$$

日常枚举$d$的倍数：

$$ans=\sum_P\sum_{d=1}^{min(\lfloor\frac{n}{p}\rfloor,\lfloor\frac{m}{p}\rfloor)}\mu(d)\sum_{i=1}^{\lfloor\frac{n}{dp}\rfloor}\sum_{j=1}^{\lfloor\frac{m}{dp}\rfloor}$$

$$ans=\sum_P\sum_{d=1}^{min(\lfloor\frac{n}{p}\rfloor,\lfloor\frac{m}{p}\rfloor)}\mu(d)\lfloor\frac{n}{dp}\rfloor\lfloor\frac{m}{dp}\rfloor$$

设$T=dp$，代入：

$$ans=\sum_P\sum_{d=1}^{min(\lfloor\frac{n}{p}\rfloor,\lfloor\frac{m}{p}\rfloor)}\mu(\frac{T}{p})\lfloor\frac{n}{T}\rfloor\lfloor\frac{m}{T}\rfloor$$

改为枚举$T$：

$$ans=\sum_{T=1}^{min(n,m)}\sum_{p|T}\mu(\frac{T}{p})\lfloor\frac{n}{T}\rfloor\lfloor\frac{m}{T}\rfloor$$

设$f(T)=\sum_{p|T}\mu(\frac{T}{p}),T=i\times p[j]$，考虑线筛：

1.当$T$本身为素数时，显然有$f(T)=\mu(1)=1$；

2.当$T$拥有多个最小质因子（即$i\ mod\ p[j]=0$）时：

(1)当$i$本身无多个相同质因子时，那么当且仅当我们枚举的素数等于$p[j]$时，$\mu(\frac{T}{p})$的值不为$0$，此时$p=p[j],T=i\times p[j]$，所以有$\mu(\frac{T}{p})=\mu(i)$；

(2)当$i$本身就是拥有多个因子的数时，无论我们怎么枚举，$\mu(\frac{T}{p})$的值都为$0$，此时，仍然有$\mu(\frac{T}{p})=\mu(i)$；

3.当$T$的最小质因子只有一个时：

此时，$T$就比$i$多了一个质因数$p[j]$，因为有：

$$f(i)=\sum_{p|i}\mu(\frac{i}{p})\qquad f(T)=\sum_{p|T}\mu(\frac{i\times p[j]}{p})$$

因为$i$中没有质因数$p[j]$，所以有$\mu(\frac{i\times p[j]}{p})=-\mu(\frac{i}{p})$，在此基础上，$T$还多了一个$\mu(\frac{i\times p[j]}{p[j]})=\mu(i)$。可得，$f(T)=-f(i)+\mu(i)$；

综上，我们得到线筛方程：

$$f(T)=\left\{ \begin{array}\\ \mu(1)&&{T\in P}\\ \mu(i)&&{i\ mod\ p[j]=0}\\ -f(i)+\mu(i)&&{i\ mod\ p[j]\ne 0} \end{array} \right.$$

再求一下$f(T)$的前缀和，利用下底分块，我们就可以$O(\sqrt{n})$解决询问了。

代码

#include<bits/stdc++.h>
#define R register int
using namespace std;
const int M=1e7+5,N=1e7;
int p[M],miu[M],f[M],T;
bool check[M];
void get()
{
    R i,j,t;
    check[1]=miu[1]=1;
    for(i=2;i<=N;++i)
    {
        if(!check[i])f[i]=1,miu[i]=-1,p[++p[0]]=i;
        for(j=1;j<=p[0];++j)
        {
            t=i*p[j];if(t>N)break;
            check[t]=1;
            if(i%p[j]==0){f[t]=miu[i];miu[t]=0;break;}
            miu[t]=-miu[i],f[t]=miu[i]-f[i];
        }
        f[i]+=f[i-1];
    }
}
long long ans;
void query(int n,int m)
{
    R l,r;ans=0;
    if(n>m)swap(n,m);
    for(l=1;l<=n;l=r+1)r=min(n/(n/l),m/(m/l)),ans+=1ll*(f[r]-f[l-1])*(n/l)*(m/l);
}
void in(){get();scanf("%d",&T);}
void ac()
{
    R i,a,b;
    for(i=1;i<=T;++i)
    scanf("%d%d",&a,&b),query(a,b),printf("%lld\n",ans);
}
int main()
{
    in();ac();
    return 0;
}