51NOD 1040 1040 最大公约数之和 数论 欧拉函数

原创 2016年12月13日 14:39:28
1040 最大公约数之和
题目来源: rihkddd
基准时间限制:1 秒 空间限制:131072 KB 分值: 80 难度:5级算法题 收藏  关注
给出一个n,求1-n这n个数,同n的最大公约数的和。比如:n = 6
1,2,3,4,5,66的最大公约数分别为1,2,3,2,1,6,加在一起 = 15
Input
1个数N(N <= 10^9)
Output
公约数之和
Input示例
6
Output示例
15

设整数d,满足n%d=0;
gcd(n,k)=d —> nx + ky = d —-> n/d * x + k/d *y = 1
显然 对任意k<=n,满足gcd(n/d,k/d)=1的k的个数 就是gcd(n,k)=d的k个数
而这个k的个数就等于phi(n/d) (n的欧拉函数)
所以枚举d,将d*phi(n/d)累加就是最终答案


#include<iostream>
#include<stdlib.h>
#include<stdio.h>
#include<string>
#include<vector>
#include<deque>
#include<queue>
#include<algorithm>
#include<set>
#include<map>
#include<stack>
#include<time.h>
#include<math.h>
#include<list>
#include<cstring>
#include<fstream>
//#include<memory.h>
using namespace std;
#define ll long long
#define ull unsigned long long
#define pii pair<int,int>
#define INF 1000000007
#define pll pair<ll,ll>
#define pid pair<int,double>

void caulate_factor(int n,vector<int>&factor){//分解n的质因子
    factor.clear();
    for(int i=2,end=sqrt(n+1);i<=end;++i){
        if(n%i==0){
            factor.push_back(i);
            while(n%i==0){
                n/=i;
            }
        }
        if(i>=n)
            break;
    }
    if(n!=1)
        factor.push_back(n);
}

ll phi(ll n,vector<int>&factor){
    ll ans=n;
    for(int i=0;i<factor.size();++i){
        ans=ans*(factor[i]-1)/factor[i];
    }
    return ans;
}

ll slove(int n){
    ll ans=0;
    vector<int>factor_fac;
    for(int i=1,end=sqrt(n);i<=end;++i){//枚举d 
        if(n%i==0){
            int factor=i;
            caulate_factor(n/factor,factor_fac);
            ans+=factor*phi(n/factor,factor_fac);
            if(n/i>end){
                factor=n/i;
                caulate_factor(n/factor,factor_fac);
                ans+=factor*phi(n/factor,factor_fac);
            }
        }
    }
    return ans;
}

int main()
{
    //freopen("/home/lu/文档/r.txt","r",stdin);
    //freopen("/home/lu/文档/w.txt","w",stdout);
    int n;
    while(~scanf("%d",&n)){
        printf("%lld\n",slove(n));
    }
    return 0;
}
版权声明:本文为博主原创文章,未经博主允许不得转载。

51nod-1040-最大公约数之和(欧拉函数)

给出一个n,求1-n这n个数,同n的最大公约数的和。比如:n = 6 1,2,3,4,5,6 同6的最大公约数分别为1,2,3,2,1,6,加在一起 = 15 Input ...
  • qq978874169
  • qq978874169
  • 2016年05月06日 14:53
  • 477

51nod oj 1040 最大公约数之和 1060 最复杂的数【打表+筛重】

1040 最大公约数之和 题目来源: rihkddd 基准时间限制:1 秒 空间限制:131072 KB 分值: 80 难度:5级算法题  收藏  关注 给出一个n,求1-n这...
  • leibniz_zhang
  • leibniz_zhang
  • 2016年09月15日 22:18
  • 821

N个整数的最大公约数

经过一年的学习,发觉C++实在太难了,所以暂且放一边,学习一门相对容易入手的语言Python 经过一段时间学习发现Python真的很简洁,和适合我。哈哈 #——————————————————————...
  • IOS9527o
  • IOS9527o
  • 2016年09月06日 22:43
  • 429

最大公约数之和

描述 题目很简单,求出: 输入每行一个数n(n输出每个结果占一行。样例输入 12 样例输出 40 这道题采用了欧拉函数的性质, 1-12与12的最大公约数和就为...
  • wsnbb123456789
  • wsnbb123456789
  • 2015年11月28日 16:15
  • 329

51NOD 1040 最大公约数之和(分析 + 欧拉函数)

传送门 1040 最大公约数之和 题目来源: rihkddd 基准时间限制:1 秒 空间限制:131072 KB 分值: 80 难度:5级算法题 给出一个n,求1-n这n个数,同n的最大公约数...
  • qingshui23
  • qingshui23
  • 2016年05月30日 10:53
  • 796

51nod1040 最大公约数之和 (欧拉函数 )

51nod1040 最大公约数之和 (欧拉函数 )
  • u013790563
  • u013790563
  • 2015年04月20日 17:10
  • 863

51nod 1040 最大公约数之和

题解这题是欧拉函数的经典题。欧拉函数φ(p)表示的是与p互质的数的个数(包括1)。这和本题好像没有什么关系,这是就需要我们来转化了。gcd(n, i) = x(x为这个因子)的个数,可以表示成gcd(...
  • DYT_B
  • DYT_B
  • 2018年02月09日 16:33
  • 17

三种方法求n个数的最大公约数

/*****************************************************  File name:穷举法求最大(小)公约(倍)数  // 文件名   Auth...
  • x_i_xw
  • x_i_xw
  • 2017年03月22日 18:31
  • 1132

【51nod】1040 最大公约数之和 欧拉函数

题目传送门我连这种水题都不会做啦……感觉要炸啊……nn的最大公约数必定是nn的因数xx,如果有gg个数和nn的最大公约数为xx,那么这gg个数对答案的贡献显然为x∗gx*g。那么问题就转化成了如何快速...
  • lyfsb
  • lyfsb
  • 2017年10月16日 19:03
  • 40

51nod 1040 最大公约数之和 欧拉函数

题目连接: https://www.51nod.com/onlineJudge/questionCode.html#!problemId=1040 Description 给出一个n,求1-n这...
  • hjf1201
  • hjf1201
  • 2017年11月03日 11:51
  • 63
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:51NOD 1040 1040 最大公约数之和 数论 欧拉函数
举报原因:
原因补充:

(最多只允许输入30个字)