【BZOJ3994】[SDOI2015]約數個數和

Description

設d(x)為x的約數個數，給定N、M，求 ?

Input

輸入文件包含多組測試數據。

第一行，一個整數T，表示測試數據的組數。 接下來的T行，每行兩個整數N、M。

Output

?T行，每行一個整數，表示你所求的答案。

Sample Input

2
7 4
5 6

Sample Output

110
121

HINT

?1<=N, M<=50000

1<=T<=50000

題解：依舊是這個結論

但由于這次是多組詢問，我們先O(nsqrt(n))預處理，然后就能O(sqrt(n))回答詢問了。

#include <cstdio> #include <cstring> #include <iostream> using namespace std; typedef long long ll; const int m=50000; int sm[m+10],mu[m+10],pri[m+10]; bool np[m+10]; int T,num; ll x,y; ll f[m+10]; int main() {ll i,j,last,ans,x,y;sm[1]=mu[1]=1;for(i=2;i<=m;i++){if(!np[i]) pri[++num]=i,mu[i]=-1;sm[i]=sm[i-1]+mu[i];for(j=1;j<=num&&i*pri[j]<=m;j++){np[i*pri[j]]=1;if(i%pri[j]==0){mu[i*pri[j]]=0;break;}mu[i*pri[j]]=-mu[i];}}for(i=1;i<=m;i++) for(j=1;j<=i;j=last+1) last=i/(i/j),f[i]+=(last-j+1)*(i/j);scanf("%d",&T);while(T--){scanf("%lld%lld",&x,&y),ans=0;if(x<y) swap(x,y);for(i=1;i<=y;i=last+1){last=min(x/(x/i),y/(y/i));ans+=(sm[last]-sm[i-1])*f[x/i]*f[y/i];}printf("%lld\n",ans);}return 0; }

轉載于:https://www.cnblogs.com/CQzhangyu/p/6999507.html

總結

以上是生活随笔為你收集整理的【BZOJ3994】[SDOI2015]约数个数和 莫比乌斯反演的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。