【bzoj3944】Sum

题意：略

神奇的杜教筛。。。

顺便当做学习好了。。

对于所求的两个函数，均有g=1

然后举一个莫比乌斯函数的例子

欧拉函数同理即可

注意long long和记录已算出的值

#include <bits/stdc++.h>
#define gc getchar()
#define ll long long
#define N 2000000
using namespace std;
ll cnt,T,n,m,pri[N],pd[N],mu[N],phi[N],p[N],P[N],vis[N];
ll read()
{
	ll x=1;
	char ch;
	while (ch=gc,ch<'0'||ch>'9') if (ch=='-') x=-1;
	ll s=ch-48;
	while (ch=gc,ch>='0'&&ch<='9') s=s*10+ch-48;
	return s*x;
}
ll get_phi(ll x)
{
	return (x<N)?phi[x]:P[m/x];
}
ll get_mu(ll x)
{
	return (x<N)?mu[x]:p[m/x];
}
void get(ll n)
{
	ll t=m/n;
	if (n<N||vis[t]) return;
	vis[t]=1;
	p[t]=1,P[t]=(ll)n*(n+1)/2;
	for (ll i=2,j;i<=n;i=j+1)
	{
		j=n/(n/i),get(n/i);
		p[t]-=(ll)(j-i+1)*get_mu(n/i);
		P[t]-=(ll)(j-i+1)*get_phi(n/i);
	}
}
int main()
{
	pd[phi[1]=mu[1]=1]=1;
	for (int i=2;i<N;i++)
	{
		if (!pd[i]) pri[++cnt]=i,mu[i]=-1,phi[i]=i-1;
		for (int j=1;j<=cnt&&pri[j]*i<N;j++)
		{
			pd[pri[j]*i]=1;
			if (i%pri[j]==0)
			{
				phi[pri[j]*i]=phi[i]*pri[j];
				mu[pri[j]*i]=0;
				break;
			}
			mu[pri[j]*i]=-mu[i];
			phi[pri[j]*i]=phi[i]*(pri[j]-1);
		}
	}
	for (int i=1;i<N;i++) mu[i]+=mu[i-1],phi[i]+=phi[i-1];
	T=read();
	while (T--)
	{
		memset(vis,0,sizeof(vis));
		m=n=read();
		if (n<N) printf("%lld %lld\n",phi[n],mu[n]);
		else get(n),printf("%lld %lld\n",P[1],p[1]);
	}
	return 0;
}