bzoj 3994 约数个数和 —— 反演+数论分块

题目：https://www.lydsy.com/JudgeOnline/problem.php?id=3994

推导过程和这里一样：https://www.cnblogs.com/MashiroSky/p/6365020.html

不用取模但注意开 long long 。

代码如下：

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
int const xn=5e4+;
int cnt,pri[xn];
ll f[xn],mu[xn];
bool vis[xn];
int rd()
{
  int ret=,f=; char ch=getchar();
  while(ch<''||ch>''){if(ch=='-')f=; ch=getchar();}
  while(ch>=''&&ch<='')ret=ret*+ch-'',ch=getchar();
  return f?ret:-ret;
}
void init()
{
  int mx=5e4; mu[]=;
  for(int i=;i<=mx;i++)
    {
      if(!vis[i])pri[++cnt]=i,mu[i]=-;
      for(int j=;j<=cnt&&(ll)i*pri[j]<=mx;j++)
    {
      vis[i*pri[j]]=;
      if(i%pri[j])mu[i*pri[j]]=-mu[i];
      else break;
    }
    }
  for(int i=;i<=mx;i++)mu[i]+=mu[i-];
  for(int n=;n<=mx;n++)
    for(int i=,j;i<=n;i=j+)
      {
    j=n/(n/i);
    f[n]+=(ll)(n/i)*(j-i+);//
      }
}
int main()
{
  int T=rd(); init();
  while(T--)
    {
      int n=rd(),m=rd(); int mn=min(n,m); ll ans=;//
      for(int i=,j;i<=mn;i=j+)
    {
      j=min(n/(n/i),m/(m/i));
      ans+=(mu[j]-mu[i-])*f[n/i]*f[m/i];
    }
      printf("%lld\n",ans);
    }
  return ;
}