【BZOJ 2705】[SDOI2012]Longge的问题 欧拉函数

水水的数学题

枚举gcd所以

ans=∑d*∑gcd(i,n)==d=∑d*∑gcd(i/d,n/d)==1

第二个求个可以显然发现就是欧拉函数，问题解决

可是。。。由于最近写数学题全部都是O(n)用积性函数性质来筛欧拉函数。。。结果想了半天没法筛啊，数据范围那么大，结果。。。。只用每一次根号n的复杂度暴力求就好了脑残。

#include<cstdio>
#include<cstring>
#include<iostream>
#include<cmath>
#define maxn 100021
#define LL long long
using namespace std;
int n,m;
int phi(int x){
	int ans=x;
	for(int i=2;i*i<=x;i++){
		if(x%i==0){
			ans=ans/i*(i-1);
			while(x%i==0)x/=i;
		}
	}
	if(x>1)ans=ans/x*(x-1);
	return ans;
}

int main(){
	scanf("%d",&n);m=sqrt(n);
	LL ans=0;
	for(int i=1;i<=m;i++){
		if(n%i)continue;
		ans+=(LL)i*phi(n/i);
		if(i*i!=n)ans+=(LL)(n/i)*phi(i);
	}
	printf("%lld",ans);
	return 0;
}