【BZOJ 2705】 [SDOI2012]Longge的问题

hzwer太神啦~~~
k为n的约数 f(k)为 k的个数和
举例: k为2 n为8 那么f(k)=10 分别为 2-1 4-2 6-3 8-4

所以答案 ans=sigma k*f(k); 自己想想就知道是对的!
k可以用sqrt(n) 枚举出来 那么重点就是求f(k)!

因为 k=gcd(m,n) 那么 gcd(m/k,n/k)=1 (*)
所以有多少个m使得 *成立 f(k)就等于多少
又因为 m<=n 所以f(k)=n/k内的质数和 即为欧拉函数值!!

#include <cstdio>
#include <cmath>
#include <iostream>
using namespace std;
#define LL long long
LL f(LL x)
{
	LL 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()
{
	LL n;cin>>n;
	LL ans=0;
	for(int i=1;i*i<=n;i++)
		if(n%i==0)
		{
			ans+=i*f(n/i);
			if(i*i!=n) ans+=n/i*f(i);
		}
	cout<<ans;
	return 0;
}
posted @ 2016-01-08 09:50  sxb_201  阅读(156)  评论(0编辑  收藏  举报