Comet OJ - Contest #8 E.神奇函数（欧拉函数性质）

题目描述

BB

其实是一道sb题
10¹³+10组数据足以把
杜教筛/min25/洲阁筛/反演+筛μ/分块+筛质数
给送上天了
所以正解肯定是T√n的做法

性质

欧拉函数有一个著名的性质：
$n={\sum }_{d|n}\phi (d)$
证明：
设$F(n)={\sum }_{d|n}\phi (d)$，则
$F(n)\ast F(m)={\sum }_{i|n}\phi (i)\ast {\sum }_{j|m}\phi (j)$（nm互质）
$={\sum }_{i|n}{\sum }_{j|m}\phi (i\ast j)$
$=F(n\ast m)$
所以证得F(n)是积性函数
求$F(p^k)$（p为质数）
$F(p^k)={\sum }_{i=0}^k\phi (p^i)$
$=({\sum }_{i=1}^k{p}^i\ast (1-\frac{1}{p}))+1$
$=({\sum }_{i=1}^k{p}^{i-1}\ast (p-1))+1$
$=({\sum }_{i=1}^k{p}^i-p^{i-1})+1$
$=p^k-p^0+1$
$=p^k$
由于F(n)是积性函数，且F(p^k)=p^k，所以可以推得F(n)=n（对于任意n）
所以
$F(n)={\sum }_{d|n}\phi (d)$
$n={\sum }_{d|n}\phi (d)$
参考：https://blog.csdn.net/liuzibujian/article/details/81086324

题解

题目的$f(n)=\prod {a_i}^{\lfloor \frac{p_i}{2}\rfloor }$
可以发现，每个质因子的质数除了2
那么
$f(n)={\sum }_{d|f(n)}\phi (d)$
本质上，枚举d其实是枚举n中指数为1的倍数的约数
由于先前除了2，考虑把d和f(n)的指数都乘以2，这样得到的实际上是一样的
所以
$f(n)={\sum }_{d^2|n}\phi (d)$
$ans={\sum }_{i=1}^{n}f(i)={\sum }_{i=1}^n{\sum }_{d^2|i}\phi (d)$
$={\sum }_{d=1}^{\lfloor \sqrt{n}\rfloor }{\sum }_{d^2|iandi⩽n}\phi (d)$
$={\sum }_{d=1}^{\lfloor \sqrt{n}\rfloor }\phi (d)\ast \lfloor \frac{n}{d^2}\rfloor$
没了

code

#include <algorithm>
#include <iostream>
#include <cstdlib>
#include <cstring>
#include <cstdio>
#include <cmath>
#define fo(a,b,c) for (a=b; a<=c; a++)
#define fd(a,b,c) for (a=b; a>=c; a--)
#define Len 3162277
using namespace std;

bool f[Len+1];
int p[Len+1];
int phi[Len+1];
int T,N,i,j,k,l,len;
long long n,ans;

void init()
{
    
    
	int i,j;
	
	memset(f,0,sizeof(f));
	len=0;
	
	phi[1]=1;
	fo(i,2,Len)
	{
    
    
		if (!f[i])
		{
    
    
			p[++len]=i;
			phi[i]=i-1;
		}
		
		fo(j,1,len)
		if ((long long)i*p[j]<=Len)
		{
    
    
			f[i*p[j]]=1;
			phi[i*p[j]]=phi[i]*p[j];
			
			if (!(i%p[j]))
			break;
			
			phi[i*p[j]]=phi[i*p[j]]/p[j]*(p[j]-1);
		}
		else
		break;
	}
}

int main()
{
    
    
//	freopen("e.in","r",stdin);
	
	init();
	
	scanf("%d",&T);
	for (;T;--T)
	{
    
    
		scanf("%lld",&n);
		N=floor(sqrt(n));
		
		ans=0;
		fo(i,1,N)
		ans+=phi[i]*(n/i/i);
		
		printf("%lld\n",ans);
	}
}