bzoj 3944 sum 杜教筛

题意，求1-n的欧拉函数和莫比乌斯函数之和。
裸的杜教筛，比较简单，直接做就好。但是存储的时候用了两个模拟链表，虽然空间没炸，但是很明显空间较大的时候会很大程度上影响时间，所以我理所当然的T了。
不得已用了一发Po姐的方法，用map，就会好一点，但是时间玄学，交了5遍过了3遍。。
真是玄学。。

#include<cstdio>
#include<cstdlib>
#include<algorithm>
#include<cstring>
#define fo(i,a,b) for(int i=a;i<=b;i++)
#define fd(i,a,b) for(int i=a;i>b=;i--)
using namespace std;
typedef long long ll;
ll n,m;
const int M=5000000;
const int mo=857777;
ll mu[M],phi[M],bz[M],a[M],tot=0;
ll next[mo],go[mo],head[mo],v[mo],cnt=0;
ll next1[mo],go1[mo],head1[mo],v1[mo];
inline void add(ll x,ll y,ll z)
{
    go[++tot]=y;
    v[tot]=z;
    next[tot]=head[x];
    head[x]=tot;
}
inline void add1(ll x,ll y,ll z)
{
    go1[++cnt]=y;
    v[cnt]=z;
    next[cnt]=head[x];
    head[x]=cnt;
}
inline ll calcphi(ll x)
{
    if (x<=M)return phi[x];
    ll k=x%mo;
    ll ans=n*(n+1)>>1;
    ll i=head[k];
    while (i)
    {
        if (go[i]==x)return v[i];
        i=next[i];
    }
    for(ll l=2,r;l<=x;l=r+1)
    {
        r=x/(x/l);
        ans-=(r-l+1)*calcphi(x/l);
    }
    add(k,x,ans);
    return ans;
}
inline ll calcmu(ll x)
{
    if (x<=M)return mu[x];
    ll ans=1,k=x%mo;
    ll i=head[k];
    while (i)
    {
        if (go[i]==x)return v[i];
        i=next[i];
    }
    for(ll l=2,r;l<=x;l=r+1)
    {
        r=x/(x/l);
        ans-=(r-l+1)*calcmu(x/l); 
    }
    add1(k,x,ans);
    return ans;
}
inline void pre()
{
    phi[1]=mu[1]=1;
    tot=0;
    fo(i,2,M)
    {
        if (!bz[i])a[++tot]=i,mu[i]=-1,phi[i]=i-1;
        fo(j,1,tot)
        {
            if (i*a[j]>M)break;
            bz[i*a[j]]=1;
            if (!(i%a[j]))
            {
                mu[i*a[j]]=0;
                phi[i*a[j]]=phi[i]*a[j];
                break;
            }
            phi[i*a[j]]=phi[i]*(a[j]-1);
            mu[i*a[j]]=-mu[i];
        }
    }
    fo(i,1,M)
    {
        phi[i]+=phi[i-1];
        mu[i]+=mu[i-1];
    }
}
int main()
{
    int cas,n;
    pre();
    scanf("%d",&cas);
    fo(fuck,1,cas)
    {
        scanf("%d",&n);
        tot=0;
        cnt=0;
        //memset(head,0,sizeof(head));
        //memset(head1,0,sizeof(head1));
        printf("%lld %lld\n",calcphi(n),calcmu(n));
    }
    return 0;
}