[hdu6588]Function

令$m=\lfloor \sqrt[3]{n} \rfloor-1$ 
   $\sum_{i=1}^{n}gcd(floor(\sqrt[3]{i}),i)$
=$\sum_{i=1}^{m}\sum\limits_{j=i^{3}}^{(i+1)^{3}-1}gcd(i,j)+\sum\limits_{i=(m+1)^{3}}^{n}gcd(m+1,i)$
考虑该式的前缀和，即
   $\sum_{i=1}^{n}gcd(m,i)$(上式的m和n和之前的无关系)
=$\sum_{d|m}d\sum_{i=1}^{n/d}\varepsilon(gcd(i,m/d))$
=$\sum_{t|m}\mu(t)\sum_{dt|m}n/dt\cdot d$
=$\sum_{T|m}n/T \sum_{d|T}d\cdot \mu(T/d)$
=$\sum_{T|m}n/T\cdot \varphi(T)$
对$\varphi$线性筛，原式后半部分可以用o(m)，考虑前半部分
=$\sum_{i=1}^{m}\sum_{T|i}\varphi(T)\cdot (((i+1)^{3}-1)/T-(i^3-1)/T)$
=$\sum_{T=1}^{m}\varphi(T)\sum_{i=1}^{m/T}((iT+1)^{3}-1)/T-(iT^3-1)/T$
=$\sum_{T=1}^{m}\varphi(T)\sum_{i=1}^{m/T}3Ti^{2}+3i+1$
这个就也可以o(m)计算了(后面的sigma可以预处理，也可以套公式)，总时间复杂度即o(m)

#include<bits/stdc++.h>
using namespace std;
#define N 10000005
#define mod 998244353
#define ll __int128
int T,m,ans,s1[N],s2[N],vis[N],phi[N],p[N];
char s[101];
int main(){
    s1[1]=3;
    s2[1]=4;
    phi[1]=1;
    for(int i=2;i<N-4;i++){
        s1[i]=(s1[i-1]+3LL*i*i)%mod;
        s2[i]=(s2[i-1]+3LL*i+1)%mod;
        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]*(p[j]-1); 
            else{
                phi[i*p[j]]=phi[i]*p[j];
                break;
            }
        }
    }
    scanf("%d",&T);
    while (T--){
        scanf("%s",s);
        ll n=0;
        for(int i=0;s[i];i++)n=n*10+(s[i]-'0');
        for(m=1;(ll)m*m*m<=n;m++);
        m-=2;
        int ans=0;
        for(int i=1;i<=m+1;i++)
            if ((m+1)%i==0)ans=(ans+(n/i-((ll)(m+1)*(m+1)*(m+1)-1)/i)%mod*phi[i])%mod;
        for(int i=1;i<=m;i++)ans=(ans+1LL*phi[i]*(1LL*i*s1[m/i]+s2[m/i]+mod))%mod;
        printf("%d\n",ans);
    }
}