[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_{2ile n}varphi(2i)$

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

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

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

对于$sum_{2ile n}varphi(2i)$，对$i$的奇偶性分类讨论，即$sum_{2ile n}varphi(2i)=sum_{4ile n}varphi(2i)+sum_{4i+2le 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}}log_{2}n)$，可以通过

#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
",ans);
    }
}