HDU 6706 huntian oy(杜教筛)

题意

给你n计算函数
$$f(n,a,b)=\sum _{i=1}^n\sum _{j=1}^{i}gcd(i^a-j^a,i^b-j^b)[gcd(i,j)==1]\mathrm{m}\mathrm{o}\mathrm{d}1e^9+7$$

思路

数论中的若干定理，也可以通过打表发现
若$n>m$,$gcd(n,m)=1$,那么$gcd(n^a-m^a,n^b-m^b)=n^{gcd(a,b)}-m^{gcd(a,b)}$
那又有题目中说了$a$与$b$互质，那么原式就转变为
$$\sum _{i=1}^n\sum _{j=1}^{i}i-j[gcd(i,j)==1]$$
我们先只考虑后面这个式子
$$\sum _{j=1}^{i}i-j[gcd(i,j)==1]$$
这个式子的含义就是$1$到$i$中与$i$互质的数与$i$的差的和。我们也知道$1$到$n(n>1)$中与$n$互质的数的和为$\frac{n\phi (n)}{2}$,与$n$互质的数的个数就是$\phi (n)$。那么式子就为
$$\sum _{j=1}^{i}i-j[gcd(i,j)==1]=i\phi (i)-\frac{i\phi (i)}{2}(i>1)$$
原式就等于$(i=1)$答案为$0$没有贡献
$$\sum _{i=2}^{n}i\phi (i)-\frac{i\phi (i)}{2}=\sum _{i=2}^n\frac{i\phi (i)}{2}$$
对这个式子我们可以换成求解
$$\sum _{i=1}^{n}i\phi (i)$$
用狄利克雷卷积求杜教筛的方法，对于原函数${\sum }_{i=1}^{n}i\ast \phi (i)$
令$f(x)=x\phi (x)$，考虑这个函数卷上一个什么函数求前缀和会比较容易呢
已知$(\phi \ast I)(n)=id$，
有一个小技巧就是把卷上一个函数能把前面的系数给消掉变为常数的，那么我们再卷上一个$x$就可以了
令$$g(x)=x$$
那么

$$\sum _{d|n}g(d)f(\frac{n}{d})=\sum _{d|n}d\frac{n}{d}\phi (\frac{n}{d})=n\sum _{d|n}\phi (\frac{n}{d})=n^2$$
所以这两个函数的卷积的前缀和是很好求的
令$S(n)={\sum }_{i=1}^{n}f(i)$有

$$\sum _{i=1}^{n}g\ast f=\sum _{i=1}^n\sum _{d|i}g(d)f(\frac{i}{d})=\sum _{i=1}^n{i}^2$$
$$\sum _{i=1}^n\sum _{d|i}g(d)f(\frac{i}{d})$$
对于这个式子是枚举1到n中每个数的因子的，他和枚举1到n中每个数的倍数是等价的那就有
$$\sum _{i=1}^n\sum _{d|i}g(d)f(\frac{i}{d})=\sum _{d=1}^n\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }g(d)f(i)=\sum _{d=1}^{n}g(d)\sum _{i=1}^{\lfloor \frac{n}{d}\rfloor }f(i)=\sum _{d=1}^{n}g(d)S(\lfloor \frac{n}{d}\rfloor )$$
$$\sum _{d=1}^{n}g(d)S(\lfloor \frac{n}{d}\rfloor )=g(1)S(n)+\sum _{d=2}^{n}g(d)S(\lfloor \frac{n}{d}\rfloor )=\sum _{i=1}^n{i}^2$$
$$g(1)S(n)=\sum _{i=1}^n{i}^2-\sum _{d=2}^{n}g(d)S(\lfloor \frac{n}{d}\rfloor )$$
由于$g(x)=x$，所以
$$S(n)=\frac{n(n+1)(2n+1)}{6}-\sum _{d=2}^{n}dS(\lfloor \frac{n}{d}\rfloor )$$
套杜教筛就可以了
求完$S(n)$后，那么答案就是
$$\frac{S(n)-1}{2}$$
把多出来的一个$1$减掉就是$$\sum _{i=2}^{n}i\phi (i)$$
除以二就是答案了

#include
using namespace std;
const int N=1e6+5;
const int mod=1e9+7;
const int inv6=166666668;
const int inv2=500000004;
long long phi[N];
bool vis[N];
int prime[N];
int cnt;
void init()
{
    phi[1]=1;
    cnt=0;
    for(int i=2; iP;
long long S(long long n)
{
    if(n