HDU 5663 （莫比乌斯反演）

分析：

这道题就是一道公式推导题，我再来推一遍公式：

设$f(x)$函数：

$$f(x) = \begin{cases} 1, & \text{x是完全平方数} \\ 0, & \text{否则} \end{cases}$$

则$ANS={n}\times{m}-\sum_{i=1}^{n}\sum_{j=1}^{m}f(gcd(i,j))$这个时候只需要维护后面一块就可以了

$$\begin{align} TEST&= \sum_{i=1}^{n}\sum_{j=1}^{m}f(gcd(i,j))\\ &=\sum_{i=1}^{n}\sum_{j=1}^{m}\sum_{d=gcd(i,j)}f(d)\\ &=\sum_{d=1}^{min(n,m)}f(d)\sum_{i=1}^{n}\sum_{j=1}^{m}(gcd(i,j)=d)\\ &=\sum_{d=1}^{min(n,m)}f(d)\sum_{i=1}^{\lfloor{{n}\over{d}}\rfloor}\sum_{j=1}^{\lfloor{{m}\over{d}}\rfloor}(gcd(i,j)=1)\\ &=\sum_{d=1}^{min(n,m)}f(d)\sum_{i=1}^{\lfloor{{n}\over{d}}\rfloor}\sum_{j=1}^{\lfloor{{m}\over{d}}\rfloor}\sum_{t|gcd(i,j)}\mu(t)\\ &=\sum_{d=1}^{min(n,m)}f(d)\sum_{t=1}^{min({\lfloor{{n}\over{d}}\rfloor},{\lfloor{{m}\over{d}}\rfloor})}\mu(t)\times{\lfloor{{n}\over{dt}}\rfloor}\times{\lfloor{{m}\over{dt}}\rfloor}，&\text{使$g=dt$}\\ &=\sum_{g=1}^{min(n,m)}{\lfloor{{n}\over{g}}\rfloor}\times{\lfloor{{m}\over{g}}\rfloor}\sum_{d|g}\mu(d)\times{f({{g}\over{d}})} \end{align}$$
这个时候就搞定了，后面一块 $\sum_{g}\sum_{d|g}\mu(d)\times{f({{g}\over{d}})}$维护这个式子和，复杂度是 $o(nlogn)$但是还可降低： $\sum_{g}\sum_{x^2|g}\mu(x^2)\times{f({{g}\over{x^2}})}\implies\sum_{x}\sum_{x^2|g}\mu(d)$，这个时候复杂度就是 $o(n)$（复杂度计算就是平方的倒数和，泰勒展开）。
然后维护 $\sum_{g}\sum_{d|g}\mu(d)\times{f({{g}\over{d}})}$的前缀和，就可以分块做了，分块的复杂度 $o(\sqrt{n})$

总的复杂度是$o(n+T\sqrt{n})$

代码：

#include <cstdio>
#include <cstdlib>
#include <cmath>
#include <cstring>
#include <algorithm>
#include <set>
#include <map>
#include <queue>
#include <vector>
#include <string>

using namespace std;

typedef long long LL;
typedef vector <int>    VI;
typedef pair <int,int>  PII;
#define FOR(i,x,y)  for(int i = x;i < y;++ i)
#define IFOR(i,x,y) for(int i = x;i > y;-- i)
#define pb  push_back
#define mp  make_pair
#define fi  first
#define se  second 

const int maxn = 10000010;

int mu[maxn],prime[maxn];
bool check[maxn];

void Mobius(){
    memset(check,false,sizeof(check));
    prime[0] = 0;
    mu[1] = 1;
    FOR(i,2,maxn){
        if(!check[i])   {mu[i] = -1;prime[++prime[0]] = i;}
        FOR(j,1,prime[0]+1){
            if(i*prime[j] > maxn)  break;
            check[i*prime[j]] = true;
            if(i% prime[j] == 0)    {mu[i*prime[j]] = 0;break;}
            else    {mu[i*prime[j]] = -mu[i];}
        }
    }
}

LL sum[maxn];
int n,m;

LL calc(int x,int y){
    LL ans = 0;
    int i = 1;
    while(i <= x){
        int p = x/i,q = y/i;
        int j = min(x/p,y/q);
        ans += (LL)p*q*(sum[j]-sum[i-1]);
        i = j+1;
    }
    return ans;
}

int main(){
    Mobius();
    memset(sum,0,sizeof(sum));
    FOR(i,1,maxn){
        int k = i*i;
        if(k >= maxn)   break;
        for(int j = k;j < maxn;j += k)  sum[j] += mu[j/k];
    }
    FOR(i,1,maxn)   sum[i] += sum[i-1];
    int T;  scanf("%d",&T);
    while(T--){
        scanf("%d%d",&n,&m);
        if(n < m)   swap(n,m);
        printf("%I64d\n",(LL)n*m-calc(m,n));
    }
    return 0;
}