[luogu6860]象棋与马

根据扩欧$(a,b)=1$必须要满足，同时，若$a+b$为偶数则格子的”奇偶性“不变，因此$a+b$必须为奇数

反过来，容易证明满足$(a,b)=1$且$a+b$为奇数则一定可行（构造从$(0,0)$到$(0,1)$的一组解即可）

不妨假设$a$为奇数、$b$为偶数（答案再乘以2），分两类考虑：

1.$a<b$，显然与$b$互素的数必然是奇数，因此即$\sum_{2i\le n}\varphi(2i)$

2.$a>b$，由于$(a,b)=1$等价于$(a,a-b)=1$，因此每一个小于$a$且与$a$互素的奇数与另一个偶数对应，因此即$\sum_{2i+1\le n}\frac{\varphi(2i+1)}{2}$

不妨先把答案的2乘进去，那么即$f_{n}=\sum_{i=1}^{n}\varphi(i)+\sum_{2i\le n}\varphi(2i)$（偶数要算两次）

根据积性或$\varphi$的计算过程，若$i$为奇数，则有$\varphi(2i)=\varphi(i)$，若$i$为偶数，则$\varphi(2i)=2\varphi(i)$

对于$\sum_{2i\le n}\varphi(2i)$，对$i$的奇偶性分类讨论，即$\sum_{2i\le n}\varphi(2i)=\sum_{4i\le n}\varphi(2i)+\sum_{4i+2\le n}\varphi(2i+1)=f_{\frac{n}{2}}$

根据$f_{n}=\sum_{i=1}^{n}\varphi(i)+f_{\frac{n}{2}}$，再用杜教筛优化，时间复杂度为$o(Tn^{\frac{2}{3}})$，可以通过

#include<bits/stdc++.h>
using namespace std;
#define N 40000005
#define ull unsigned long long
int t,vis[N],p[N];
ull n,ans,phi[N];
map<ull,ull>mat;
ull calc(ull n){
    if (n<N-4)return phi[n];
    if (mat[n])return mat[n];
    ull ans;
    if (n&1)ans=(n+1)/2*n;
    else ans=n/2*(n+1);
    for(ull i=2,j;i<=n;i=j+1){
        j=n/(n/i);
        ans-=(j-i+1)*calc(n/i);
    }
    return mat[n]=ans;
}
int main(){
    phi[1]=1;
    for(int i=2;i<N-4;i++){
        if (!vis[i]){
            p[++p[0]]=i;
            phi[i]=i-1;
        }
        for(int j=1;(j<=p[0])&&(i*p[j]<N-4);j++){
            vis[i*p[j]]=1;
            if (i%p[j])phi[i*p[j]]=phi[i]*phi[p[j]];
            else{
                phi[i*p[j]]=phi[i]*p[j];
                break;
            } 
        }
    }
    for(int i=2;i<N-4;i++)phi[i]=phi[i]+phi[i-1];
    scanf("%d",&t);
    while (t--){
        scanf("%llu",&n);
        ans=0;
        while (n>1){
            ans+=calc(n);
            n/=2;
        }
        printf("%llu\n",ans);
    }
}