HDU5608 function

设
$G(i)=i^2−3i+2$
有莫比乌斯反演可得出$F(i)$关于$G(i)$的关系
带入得
$∑_{i=1}^nF(i)=∑_{d=1}^nu(d)∑_{i=1}^{n/d}G(i)$
后面可以O(1)
前面就是经典的杜教筛了

#include<cstdio>
#include<iostream>
#include<cstring>
#include<cstdlib>
#include<algorithm>
#include<map>
using namespace std;
#define ll long long
char c;
inline void read(ll &a)
{
    a=0;do c=getchar();while(c<'0'||c>'9');
    while(c<='9'&&c>='0')a=(a<<3)+(a<<1)+c-'0',c=getchar();
}

const
    ll Mod=(int)1e9+7,MAXN=(int)3e6+1;
bool check[MAXN];
int mu[MAXN],mu2[MAXN],mu3[MAXN],prime[MAXN],tot;

inline int add(int x,int y)
{
    int res=x+y;
    if(res>=Mod)
        res-=Mod;
    if(res<0)res+=Mod;
    return res;
}
inline int mul(int x,int y){return x*1ll*y%Mod;}
inline 
int A(ll x)
{
    if(x&1)
        return (x%Mod)*((x+1>>1)%Mod)%Mod;
    return ((x>>1)%Mod)*((x+1)%Mod)%Mod;
}
inline 
int cn(ll x,ll y){return A(y)-A(x-1);}
map<ll ,int >M;
int Mu(ll x)
{
    if(x<MAXN)return mu[x];
    if(M[x])return M[x];
    ll i=2,j;
    int res=0;
    while(i<=x)
    {
        j=x/(x/i);
        res=add(res,mul(add(j,-i+1),Mu(x/i))%Mod);
        i=j+1;
    }
    return M[x]=add(1,Mod-res);
}

map<ll ,int >M2;
int Mu2(ll x)
{
    if(x<MAXN)return mu2[x];
    if(M2[x])return M2[x];
    ll i=2,j;
    int res=0;
    while(i<=x)
    {
        j=x/(x/i);
        res=add(res,mul(cn(i,j),Mu2(x/i))%Mod);
        i=j+1;
    }
    return M2[x]=add(1,Mod-res);
}
const
    int rev=166666668;
int A3(int x)
{
    return mul(mul(mul(add(mul(x,2),1),add(x,1)),x),rev);
}

int cn3(int x,int y)
{
    return add(A3(y),-A3(add(x,-1)));
}
map<ll ,int >M3;
int Mu3(ll x)
{
    if(x<MAXN)return mu3[x];
    if(M3[x])return M3[x];
    ll i=2,j;
    int res=0;
    while(i<=x)
    {
        j=x/(x/i);
        res=add(res,mul(cn3(i,j),Mu3(x/i))%Mod);
        i=j+1;
    }
    return M3[x]=add(1,Mod-res);
}

int G(ll x)
{
    x%=Mod;
    return add(add(cn3(1,x),-mul(cn(1,x),3)),mul(2,x));
}

ll S(ll x)
{
    ll i=1,j,res=0;
    while(i<=x)
    {
        j=x/(x/i);
        int cc=x/i%Mod;
        res=add(res,mul(G(cc),add(Mu(j),-Mu(i-1))) );
        i=j+1;
    }
    return res;
}

void out(int x)
{
    if(x<0)putchar('-'),x=-x;
    if(x>9)out(x/10);
    putchar('0'+x%10);
}

int main()
{
    mu2[1]=mu3[1]=mu[1]=1;
    for(ll i=2;i<MAXN;i++)
    {
        if(!check[i])prime[++tot]=i,mu[i]=-1;
        for(ll j=1,k;j<=tot&&(k=prime[j]*i)<MAXN;j++)
        {
            check[k]=true;
            if(i%prime[j]==0)
            {
                mu[k]=0;
                break;
            }
            mu[k]=-mu[i];
        }
        mu2[i]=mul(mu[i],i);
        mu3[i]=mul(mu2[i],i);
        mu[i]=add(mu[i],mu[i-1]);
        mu2[i]=add(mu2[i],mu2[i-1]);
        mu3[i]=add(mu3[i],mu3[i-1]);
    }
    ll t;
    read(t);
    while(t--)
    {
        ll n;
        read(n);
        out(S(n));
        puts("");
    }
    return 0;
}