【BZOJ2693】jzptab(莫比乌斯反演)
题面
讨厌权限题,只能跑到别的OJ上交
和这题是一样的
多组数据
求$$sum_{i=1}^nsum_{j=1}^mlcm(i,j)$$
题解
前面的部分直接看上面的那个链接
[ans=sum_{d=1}^ndsum_{i=1}^{n/d}i^2S(frac{n}{id})S(frac{m}{id})mu(i)
]
其中(S(x)=1+2+...x=frac{x(x+1)}{2})
令(T=id)
[ans=sum_{T=1}^nS(frac{n}{T})S(frac{m}{T})sum_{d|T}d^2frac{T}{d}mu(d)
]
设(f(x)=x^2mu(x)),(g(x)=frac{T}{x})
很显然,这两个都是积性函数
所以,后面的东西也是积性函数
直接线性筛出来就行了
#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<cmath>
#include<algorithm>
#include<set>
#include<map>
#include<vector>
#include<queue>
using namespace std;
#define MOD 100000009
#define MAX 10000000
inline int read()
{
int x=0,t=1;char ch=getchar();
while((ch<'0'||ch>'9')&&ch!='-')ch=getchar();
if(ch=='-')t=-1,ch=getchar();
while(ch<='9'&&ch>='0')x=x*10+ch-48,ch=getchar();
return x*t;
}
bool zs[MAX+10];
int sum[MAX+10],pri[MAX],tot;
void pre()
{
zs[1]=true;sum[1]=1;
for(int i=2;i<=MAX;++i)
{
if(!zs[i])pri[++tot]=i,sum[i]=(i-1ll*i*i%MOD+MOD)%MOD;
for(int j=1;j<=tot&&i*pri[j]<=MAX;++j)
{
zs[i*pri[j]]=true;
if(i%pri[j]==0){sum[i*pri[j]]=1ll*sum[i]*pri[j]%MOD;break;}
else sum[i*pri[j]]=1ll*sum[i]*sum[pri[j]]%MOD;
}
}
for(int i=1;i<=MAX;++i)sum[i]=(sum[i-1]+sum[i])%MOD;
}
int n,m;
int main()
{
pre();
int T=read();
while(T--)
{
n=read();m=read();
if(n>m)swap(n,m);
int i=1,j;
long long ans=0;
while(i<=n)
{
j=min(n/(n/i),m/(m/i));
int tt=(1ll*(1+n/i)*(n/i)/2%MOD)*(1ll*(1+m/i)*(m/i)/2%MOD)%MOD;
ans+=1ll*(sum[j]-sum[i-1]+MOD)%MOD*tt%MOD;
ans%=MOD;
i=j+1;
}
printf("%lld
",(ans+MOD)%MOD);
}
return 0;
}