【51Nod1227】平均最小公倍数-杜教筛

测试地址：平均最小公倍数
做法：这一题需要用到杜教筛。这一题推式子的过程比较经典，是做杜教筛更难题的基础。

首先定义几个接下来要用到的符号和函数：
幂函数：$id^k(n)=n^k$（完全积性）
$lcm(i,j)$和$\gcd(i,j)$：$i$和$j$的最小公倍数和最大公约数（这个不懂就……）
函数$f$和$g$的狄利克雷卷积：$f*g$（$(f*g)(n)=\sum_{d|n}f(d)g(n/d)$）
函数$f$和$g$的普通乘法：$f\cdot g$（$(f\cdot g)(n)=f(n)g(n)$）
$[S]$：若式子$S$为真则$[S]=1$，否则$[S]=0$。

我们用公式表达要求的函数：$A(n)=\frac{1}{n}\sum_{i=1}^{n}lcm(n,i)$，求$F(n)=\sum_{i=1}^{n}A(i)$。$A$不是积性函数，所以我们要想办法将式子化成积性函数前缀和的形式。
我们知道$lcm(i,j)\times\gcd(i,j)=ij$，所以$A(n)=\frac{1}{n}\sum_{i=1}^{n}\frac{ni}{\gcd(n,i)}=\sum_{i=1}^{n}\frac{i}{\gcd(n,i)}$。但是$\frac{i}{\gcd(n,i)}$仍然不是积性的，所以我们进行进一步化简：

$$A(n)=\sum_{i=1}^{n}\frac{i}{\gcd(n,i)}=\sum_{i=1}^{n}\sum_{d|n}[\gcd(n,i)=d]\times\frac{i}{d}\\=\sum_{d|n}\sum_{\frac{i}{d}=1}^{\frac{n}{d}}[\gcd(\frac{n}{d},\frac{i}{d})=1]\times\frac{i}{d}$$
注意到 $\sum_{i/d=1}^{n/d}[\gcd(\frac{n}{d},\frac{i}{d})=1]\times\frac{i}{d}$其实就是问小于等于 $\frac{n}{d}$的数中和它互质的数的和，又因为如果 $i$和$x$互质，那么 $x-i$也和 $x$互质（这个很显然，随便证证就可以了），推测出小于等于$x$的与 $x$互质的数是成对出现的，而且每一对的贡献都是$x$，所以小于等于 $x$的与$x$互质的数的和为 $x\times\varphi(x)\times\frac{1}{2}$，特别地，小于等于 $1$的与$1$互质的数的和为 $1$（若用上面公式算将算出$\frac{1}{2}$），代入上面的式子中得：
$$A(n)=\frac{1}{2}(1+\sum_{d|n}d\times\varphi(d))$$
我们发现这个式子已经有点积性函数的意思了，我们再用这个式子来算 $F(n)$：
$$F(n)=\sum_{i=1}^{n}A(i)=\frac{1}{2}(n+\sum_{i=1}^{n}\sum_{d|i}d\times\varphi(d))=\frac{1}{2}(n+\sum_{i=1}^{n}\sum_{d=1}^{\lfloor\frac{n}{i}\rfloor}d\times\varphi(d))$$
那么 $G(n)=\sum_{i=1}^{n}i\times\varphi(i)$就是一个积性函数前缀和了。
要求函数 $id\cdot\varphi$的前缀和，考虑卷上一个好算前缀和的函数，从而得到另一个好算前缀和的函数。显然 $(id\cdot\varphi)*id=id^2$，因为 $\sum_{d|n}d\times\varphi(d)\times\frac{n}{d}=n\sum_{d|n}\varphi(d)=n^2$，那么：
$$\frac{n(n+1)(2n+1)}{6}=\sum_{i=1}^{n}i^2=\sum_{i=1}^{n}\sum_{d|i}d\times\varphi(d)\times\frac{i}{d}=\sum_{\frac{i}{d}=1}^{n}\frac{i}{d}\sum_{d=1}^{\Big\lfloor\frac{n}{\frac{i}{d}}\Big\rfloor}d\times\varphi(d)\\=\sum_{i=1}^{n}i\times G(\lfloor\frac{n}{i}\rfloor)$$
看不懂上述式子的同学可以尝试从约数、倍数和贡献的方面来理解一下……
将 $G(n)$一项提出来，得：
$$G(n)=\frac{n(n+1)(2n+1)}{6}-\sum_{i=2}^{n}i\times G(\lfloor\frac{n}{i}\rfloor)$$
再结合这个式子：
$$F(n)=\frac{1}{2}(n+\sum_{i=1}^{n}G(\lfloor\frac{n}{i}\rfloor))$$
这样就可以使用杜教筛来计算了，使用哈希表判重，再预处理前 $n^{\frac{2}{3}}$个 $G$值，就可以达到$O(n^{\frac{2}{3}})$的时间复杂度了。注意因为答案要模 $10^9+7$，所以不能直接将结果乘上 $\frac{1}{2}$，而应该乘上 $2$关于模$10^9+7$的逆元： $5\times10^8+4$。
以下是本人代码：

#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <iostream>
#include <algorithm>
#define hashsize 5000000
#define limit 1000000
#define mod 1000000007
#define ll long long
using namespace std;
ll a,b,phi[limit+5],sum[limit+5],h[hashsize+5]={0},f[hashsize+5];
bool prime[limit+5]={0};

int hash(ll x)
{
  ll s=x%hashsize;
  while(h[s]&&h[s]!=x) s=(s+1)%hashsize;
  return s;
}

void init()
{
  for(int i=1;i<=limit;i++) phi[i]=i;
  for(ll i=2;i<=limit;i++)
    if (!prime[i])
    {
      for(ll j=1;i*j<=limit;j++)
      {
        prime[i*j]=1;
        phi[i*j]=phi[i*j]/i*(i-1);
      }
    }
  sum[0]=0;
  for(int i=1;i<=limit;i++)
    sum[i]=(sum[i-1]+i*phi[i])%mod;
}

ll gcd(ll a,ll b)
{
  return (b==0)?a:gcd(b,a%b);
}

ll sum1(ll a,ll b)
{
  return (b*(b+1)/2-a*(a-1)/2)%mod;
}

ll sum2(ll x)
{
  ll a=x,b=x+1,c=2*x+1,s=6,g;
  g=gcd(a,s),a/=g,s/=g;
  g=gcd(b,s),b/=g,s/=g;
  g=gcd(c,s),c/=g,s/=g;
  return (((a*b)%mod)*c)%mod;
}

ll G(ll x)
{
  if (x<=limit) return sum[x];
  int pos=hash(x);
  if (h[pos]==x) return f[pos];
  ll i=2,next,s=0;
  while(i<=x)
  {
    next=x/(x/i);
    s=(s+sum1(i,next)*G(x/i))%mod;
    i=next+1;
  }
  h[pos]=x,f[pos]=(sum2(x)-s+mod)%mod;
  return f[pos];
}

ll F(ll x)
{
  ll i=1,next,s=0;
  while(i<=x)
  {
    next=x/(x/i);
    s=(s+(next-i+1)*G(x/i))%mod;
    i=next+1;
  }
  return (500000004*(s+x))%mod;
}

int main()
{
  init();
  scanf("%lld%lld",&a,&b);
  printf("%lld",(F(b)-F(a-1)+mod)%mod);

  return 0;
}