bzoj3944: Sum//杜教筛

bzoj3944: Sum
前几天的blog，没写完，现在补上吧…

---

题意

求莫比乌斯函数和欧拉函数的前缀和。N<=2^31-1。T(<=10)组询问。

---

题解

杜教筛。
然而这篇博客不是写题解的。这主要是因为
1.网上有很多详细题解。
2.不想打公式。
3.表达能力拙计。

---

注意事项

今天真是涨姿势…居然还有这种挂题方法…
首先CE的时候显示 File size limit exceeded。问了一圈，不解其意…
改了几个地方，把N调小了一点，莫名过了编译，但是RE了？
查错，发现l会变成负的，于是从

long long l l=(long long)n*(n+1)>>1;

改成了

long long l=(long long)n*(long long)(n+1)>>1;

诶怎么还是负的？想了想发现自己弱智……
再交，又CE，还是File size limit exceeded（FSLE??），二分了几发试了一下内存，发现好像是MLE①的另一种表达方式……真是十分奇怪的体验

UPD (2018.3.28)：昨天又CE了一道题后…发现File size limit exceeded是因为数组用花括号初始化会占内存 QAQ

①：bz上MLE的表达方法还有TLE…大概

---

代码

#include<cstdio>
#include<map>
#define N 5125000
using namespace std;
typedef long long ll;
typedef pair<ll,int> pr;
pr b;
ll f[N]={0,1};
int T,a,p[N],g[N]={0,1},tot;
bool pri[N];
map<int,pr>m;
map<int,pr>::iterator s;
pr G(int n)
{
    if(n<N)return make_pair(f[n],g[n]);
    if((s=m.find(n))!=m.end())
    return s->second;
    ll l=n*((ll)n+1)>>1,t;
    int r=1;
    for(ll i=2;i<n;i=t)
    t=n/(n/i)+1,b=G(n/i),
    l-=(t-i)*b.first,r-=(t-i)*b.second;
    return m[n]=make_pair(l,r);
}
int main()
{
    for(int i=2;i<N;i++)
    {
        if(!pri[i])p[tot++]=i,g[i]=-1,f[i]=i-1;
        for(int j=0;i*p[j]<N;j++)
        {
            pri[a=i*p[j]]=1,f[a]=f[i]*p[j];
            if(!(i%p[j]))break;
            g[i*p[j]]=-g[i],f[a]-=f[i];
        }
        g[i]+=g[i-1],f[i]+=f[i-1];
    }
    scanf("%d",&T);
    while(T--)
    scanf("%d",&a),b=G(a),
    printf("%lld %d\n",b.first,b.second);
}