一共T+1行
第1行为数据组数T(T<=10)
第2~T+1行每行一个非负整数N,代表一组询问
标签:last -- pre turn scan log images ring blank
题解:
粘自http://blog.csdn.net/skywalkert/article/details/50500009
求莫比乌斯函数的前缀和类似,从开始推就好了
#include <cstdio> #include <cstring> #include <iostream> #include <map> #include <utility> #define MP(A,B) make_pair(A,B) using namespace std; const int m=3000000; typedef long long ll; int n,num; ll phi[m+10],mu[m+10],sp[m+10],sm[m+10],pri[m+10]; bool np[m+10]; typedef pair<ll,ll> pll; map<ll,pll> mp; pll dfs(ll x) { if(x<=m) return MP(sp[x],sm[x]); if(mp.find(x)!=mp.end()) return MP(mp[x].first,mp[x].second); ll rp=x*(x+1)>>1,rm=1,i,last; for(i=2;i<=x;i=last+1) { last=x/(x/i); pll tmp=dfs(x/i); rp-=tmp.first*(last-i+1); rm-=tmp.second*(last-i+1); } mp[x]=MP(rp,rm); return MP(rp,rm); } int main() { int T,i,j; scanf("%d",&T); phi[1]=sp[1]=mu[1]=sm[1]=1; for(i=2;i<=m;i++) { if(!np[i]) pri[++num]=i,phi[i]=i-1,mu[i]=-1; sp[i]=sp[i-1]+phi[i],sm[i]=sm[i-1]+mu[i]; for(j=1;j<=num&&i*pri[j]<=m;j++) { np[i*pri[j]]=1; if(i%pri[j]==0) { phi[i*pri[j]]=phi[i]*pri[j]; mu[i*pri[j]]=0; break; } mu[i*pri[j]]=-mu[i]; phi[i*pri[j]]=phi[i]*(pri[j]-1); } } while(T--) { ll a; scanf("%lld",&a); pll tmp=dfs(a); printf("%lld %lld\n",tmp.first,tmp.second); } return 0; }
标签:last -- pre turn scan log images ring blank
原文地址:http://www.cnblogs.com/CQzhangyu/p/6999122.html