bzoj3944 Sum 杜教筛

模板题，注意杜教筛化式子的技巧：

1、构造求和函数

2、利用积性函数约数性质

3、做差算出单项

4、约数倍数转化贡献

5、分块求商

码：

#include<iostream>
#include<cstdio>
#include<map>
#include<cmath>
using namespace std;
#define M 3000000
#define ll long long 
int tot,i,j,T,su[M+5];
ll ans,phe[M+5],n[999],N,mhe[M+5],mu[M+5],phi[M+5];
bool he[M+5];
map<int,ll>pma;
map<int,int>mma;
void eular()
{
	mu[1]=phi[1]=1;
	for(i=2;i<=N;i++)
	{
		if(!he[i])
		{
			su[++tot]=i;			
		    phi[i]=i-1;
		    mu[i]=-1;
		}
		for(j=1;j<=tot&&su[j]*i<=N;j++)
		{
		int k=su[j]*i;
		he[k]=1;
		if(i%su[j]==0)
		{phi[k]=su[j]*phi[i];
		mu[k]=0;
		break;
		}
		else
		{
			phi[k]=phi[i]*phi[su[j]];
			mu[k]=-mu[i];
		}	
		}
	}
	for(i=1;i<=N;i++)phe[i]=phe[i-1]+phi[i],mhe[i]=mhe[i-1]+mu[i];	
}
ll dfs1(ll o)
{
	if(o<=N)return phe[o]; 
	if(pma[o])return pma[o];
ll lin=o*(1+o)/2;
ll l;
for(l=2;l<=o;l++)
{
ll oo=o/l;
lin-=((o/oo)-l+1)*dfs1(oo) ;
l=o/oo;	
}	
pma[o]=lin;
return lin;
}
int dfs2(int o)
{	if(o<=N)return mhe[o]; 
	if(mma[o])return mma[o];
int lin=1;
unsigned int l;
for(l=2;l<=o;l++)
{
int oo=o/l;
lin-=((o/oo)-l+1)*dfs2(oo);
l=o/oo;	
}
mma[o]=lin;
return lin;
}
int main()
{scanf("%d",&T);
for(i=1;i<=T;i++)
{
scanf("%lld",&n[i]);
N=max(n[i],N);	
}
N=pow(N,2.000/3) ;
eular();
for(i=1;i<=T;i++)
{
printf("%lld ",dfs1(n[i]));
printf("%d\n",dfs2(n[i]));	
}
}