由于这道题比较鬼畜,而且公式巨难打出,所以我粘了两页PoPoQQQ的PPT orz,并对其中一些部分解释一下
最下面的公式就是换了一种枚举的方式,其中u括号里的东西就是把d变成T/p,其中p|T的条件只要带回到原表达式中即可理解。
这里注意这里的每个前缀和代表的意义是对于每一个T满足条件的Σu之和。
下面附上我的代码
#include<cstdio>
#include<cstdlib>
#include<cmath>
#include<ctime>
#include<cstring>
#include<string>
#include<iostream>
#include<iomanip>
#include<algorithm>
using namespace std;
#define M 10001000
bool pd[M];
int u[M];
long long zs[1000000];
int top=0;
int sum[M];
void shai()
{
u[1]=1;
for(int i=2;i<M;i++)
{
if(!pd[i])
{
zs[++top]=i,
u[i]=-1;
}
for(int j=1;(long long)zs[j]*i<M;j++)
{
pd[zs[j]*i]=true;
if(i%zs[j]==0) break;
u[i*zs[j]]=-u[i];
}
}
for(int i=1;i<=top;i++)
for(int j=zs[i];j<M;j+=zs[i])
sum[j]+=u[j/zs[i]];
for(int i=1;i<M;i++) sum[i]+=sum[i-1];
}
long long get_ans(int n,int m)
{
long long re=0;
int last;
for(int i=1;i<=n && i<=m;i=last+1)
{
last=min(n/(n/i),m/(m/i));
re+=(long long)(sum[last]-sum[i-1])*(n/i)*(m/i);
}
return re;
}
int main()
{
int T;
scanf("%d",&T);
shai();
while(T--)
{
int n,m;
scanf("%d%d",&n,&m);
printf("%lld\n",get_ans(n,m));
}
return 0;
}