BZOJ4659 Lcm

原题链接：https://www.lydsy.com/JudgeOnline/problem.php?id=4659

Lcm

Description

给出A，B，考虑所有满足l<=a<=A，l<=b<=B，且不存在n>1使得n^2同时整除a和b的有序数
对(a，b)，求其lcm(a，b)之和。答案模2^30。

Input

第一行一个整数T表示数据组数。接下来T行每行两个整数A，B表示一组数据。
T ≤ 2000,A,B ≤ 4 × 10^6

Output

对每组数据输出一行一个整数表示答案模2^30的值

Sample Input

5
2 2
4 6
3 4
5 1
23333 33333

Sample Output

7
148
48
15
451085813

题解

此题与BZOJ2154 Crash的数字表格类似，只是需要从单纯的枚举$gcd$变为枚举满足$\mu\ne 0$的$gcd$，化简出来的式子与$Crash$的数字表格无异：

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

其中，$p$表示满足$\mu(p)\ne 0$的数，因为$p$的特殊限制，无法像之前那样线筛，只能暴力$O(nln(n))$筛出。

代码

由于模数为$2^30$可以直接自然溢出。

#include<bits/stdc++.h>
using namespace std;
const int M=4e6+5,mod=1<<30;
int T,miu[M],p[M/3],f[M],n,m;
bool isp[M];
void in(){scanf("%d",&T);}
void pre()
{
    miu[1]=isp[1]=1;
    for(int t,i=2;i<M;++i)
    {
        if(!isp[i])p[++p[0]]=i,miu[i]=-1;
        for(int j=1;j<=p[0];++j)
        {
            t=i*p[j];if(t>=M)break;
            isp[t]=1;
            if(i%p[j]==0){miu[t]=0;break;}
            miu[t]=-miu[i];
        }
    }
    for(int t,i=1;i<M;++i)if(miu[i])for(int j=1;;++j){t=i*j;if(t>=M)break;f[t]+=t*j*miu[j];}
    for(int i=2;i<M;++i)f[i]+=f[i-1];
}
int sum(int x){return x*(x+1)/2;}
int ask()
{
    if(n>m)swap(n,m);int ans=0;
    for(int l=1,r;l<=n;l=r+1)r=min(n/(n/l),m/(m/l)),ans+=sum(n/l)*sum(m/l)*(f[r]-f[l-1]);
    return (ans%mod+mod)%mod;
}
void ac(){pre();while(T--)scanf("%d%d",&n,&m),printf("%d\n",ask());}
int main(){in();ac();}