link
题目大意：有N个数，每个数都在区间[L,H]之间，请求出所有数的gcd恰好为K的方案数
推式子
首先可以把[L,H]之间的数字gcd恰好为K转化为[(L-1)/K+1,H/K]之间数字gcd恰好为1
然后就可以反演了
下面手误把所有的H都打成了R
\(\sum_{i_1=L}^R\sum_{i_2=L}^R\dots\sum_{i_N=L}^R[\gcd(i_1,i_2,\dots,i_N)=1]\)
\(\sum_{i_1=L}^R\sum_{i_2=L}^R\dots\sum_{i_N=L}^R\sum_{d|i_1,d|i_2,\dots,d|i_N}\mu(d)\)
\(\sum_{d=1}^R\mu(d)\left(\frac{R}d-\frac{L-1}d\right)^N\)
显然可以打数论分块
但是这道题H范围达到10的9次方，果断杜教筛
注意由于我们的d是枚举到R的（因为右边式子的关系不是相乘，而是相减，所以大于L的项也有效）而>=L的数整除L-1的值为0，所以会导致除法爆炸，所以需要特判，详细请见代码
#include
#include
#include
