[BZOJ3944]SUM 杜教筛

题意

---

求$\sum_{i=1}^{n}\mu(i)$和$\sum_{i=1}^{n}\phi(i)$

---

1.
$\sum_{i=1}^{n}\sum_{d|i}\mu(d)=1=\sum_{i=1}^{n}\sum_{j=1}^{\lfloor\dfrac{n}{i}\rfloor}\mu(j)$

$1=\sum_{i=1}^{n}S_\mu(\lfloor\dfrac{n}{i}\rfloor)$

$S_\mu(n)=1-\sum_{i=2}^{n}S_\mu(\lfloor\dfrac{n}{i}\rfloor)$

2.
$\sum_{i=1}^{n}\sum_{d|i}\phi(d)=\sum_{i=1}^{n}i=\sum_{i=1}^{n}\sum_{j=1}^{\lfloor\dfrac{n}{i}\rfloor}\phi(j)$

$\dfrac{n*(n+1)}{2}=\sum_{i=1}^{n}S_\phi(\lfloor\dfrac{n}{i}\rfloor)$

$S_\phi(n)=\dfrac{n*(n+1)}{2}-\sum_{i=2}^{n}S_\phi(\lfloor\dfrac{n}{i}\rfloor)$

套杜教筛。
题目卡常数有点坑，发现unsigned int要比int快很多。。

#include <cstdio>
#include <map>
#include <iostream>
#include <cstring>
#include <string>
#include <algorithm>
#define N 6500010
#define M 100011

using namespace std;

typedef long long ll;
typedef unsigned int uint;

int t,n,p[N],vp[M],vm[M];
ll phi[N],mu[N],rp[M],rm[M];

inline char C(){
    static char buf[100000],*p1=buf,*p2=buf;
    if(p1==p2){
        p2=(p1=buf)+fread(buf,1,100000,stdin);
        if(p1==p2)return EOF;
    }
    return *p1++;
}

inline void reaD(int &x){
    char Ch=C();x=0;
    for(;Ch>'9'||Ch<'0';Ch=C());
    for(;Ch>='0'&&Ch<='9';x=x*10+Ch-'0',Ch=C());
}

inline void first(const int x){
    mu[1]=phi[1]=1;
    ll i,j;
    for(i=2;i<=x;i++){
        if(!p[i]) p[++p[0]]=i,mu[i]=-1,phi[i]=i-1;
        for(j=1;i*p[j]<=x;j++)
            if(p[p[j]*i]=1,i%p[j]) mu[i*p[j]]=-mu[i],phi[i*p[j]]=phi[i]*(p[j]-1);
            else {mu[i*p[j]]=0,phi[i*p[j]]=phi[i]*p[j];break;}
    }
    for(i=1;i<=x;i++)phi[i]+=phi[i-1],mu[i]+=mu[i-1];
}

ll calcphi(uint x){
    if(x<=N-10) return phi[x];
    int nw=n/x;
    if(vp[nw]) return rp[nw];
    ll ret=1ll*x*(x+1)>>1;register uint i,j;
    for(i=2,j;i<=x;i=j+1){
        j=x/(x/i);
        ret-=(j-i+1)*calcphi(x/i);
    }
    vp[nw]=1;rp[nw]=ret;
    return ret;
}

ll calcmu(uint x){
    if(x<=N-10) return mu[x];
    ll ret=1;register uint i,j;
    int nw=n/x;
    if(vm[nw]) return rm[nw];
    for(i=2,j;i<=x;i=j+1){
        j=x/(x/i);
        ret-=(j-i+1)*calcmu(x/i);
    }
    vm[nw]=1;rm[nw]=ret;
    return ret;
}

int main(){
    #ifndef ONLINE_JUDGE
    freopen("1.in","r",stdin);
    freopen("1.out","w",stdout);
    #endif
    first(N-10);
    for(reaD(t);t;t--){
        memset(vm,0,sizeof(vm));memset(vp,0,sizeof(vp));
        reaD(n);
        ll Ans1=calcphi(n),Ans2=calcmu(n);;
        printf("%lld %lld\n",Ans1,Ans2);
    }return 0;
}