【BZOJ 2705】 [SDOI2012]Longge的问题

来源:互联网 发布:淘宝如何设置冰点促销 编辑:程序博客网 时间:2024/06/08 09:43

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 longLL 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;}
0 0
原创粉丝点击