题意:略
神奇的杜教筛。。。
顺便当做学习好了。。
对于所求的两个函数,均有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;
}