bzoj 3944: Sum

传说中的杜教筛，，，

其实就是用狄利克雷卷构造出一个好求的东西，然后表示出要求的东西（好简单啊看起来2333）

虚的不行，这些东西肯定2天之后就什么都不知道了2333

#include<bits/stdc++.h>
#define LL long long
using namespace std;


int n,m,cnt,prime[10000000+5]; 
LL phi[10000000+5];
int mo[10000000+5]; 

void mobius()
{
    phi[1]=mo[1]=1; m=10000000;
    for (int i=2; i<=m; i++)
    {
        if (!phi[i]) phi[i]=i-1,mo[i]=-1,prime[++cnt]=i;
        for (int j=1; j<=cnt && i*prime[j]<=m; j++)
        {
            if (i%prime[j]) phi[i*prime[j]]=phi[i]*(prime[j]-1),mo[i*prime[j]]=-mo[i];
            else {phi[i*prime[j]]=phi[i]*prime[j],mo[i*prime[j]]=0; break;}
        }
    }
    for (int i=1; i<=m; i++) phi[i]+=phi[i-1];
    for (int i=1; i<=m; i++) mo[i]+=mo[i-1];
}

map<int , LL > p,q;
map<int , bool > vis;
LL get_phi(int x){return (x<=m)?phi[x]:p[x];}
LL get_mo(int x){return (x<=m)?mo[x]:q[x];}
void solve(int x)
{
    if (x<=m) return; 
    int j=1,i;
    if (vis[x]) return; vis[x]=1;
    p[x]=((LL)x+1)*x>>1; q[x]=1;
    while (j<x)
    {
        i=j+1; j=x/(x/i); solve(x/i);
        p[x]-=get_phi(x/i)*(j-i+1); q[x]-=get_mo(x/i)*(j-i+1);
    }
}
int main()
{
    int T;
    scanf("%d",&T); 
    mobius();
    while (T--)
    {
        scanf("%d",&n); vis.clear();
        if (n<=m) printf("%lld %lld
",phi[n],mo[n]);
        else{
            solve((LL)n); printf("%lld %lld
",p[n],q[n]);
        }
    }
    return 0;
}