原题链接: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 g c d 变为枚举满足 μ≠0 μ ≠ 0 的 gcd g c d ,化简出来的式子与 Crash C r a s h 的数字表格无异:
其中, p p 表示满足的数,因为 p p 的特殊限制,无法像之前那样线筛,只能暴力筛出。
代码
由于模数为 230 2 3 0 可以直接自然溢出。
#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();}