【bzoj3944】Sum

来源：互联网发布：淘宝联盟优惠券链接编辑：程序博客网时间：2024/06/04 00:47

题意：略

神奇的杜教筛。。。

顺便当做学习好了。。

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

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

欧拉函数同理即可

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

#include <bits/stdc++.h>#define gc getchar()#define ll long long#define N 2000000using 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;}

0 0