关闭

UVA 11426 (欧拉函数)

169人阅读 评论(0) 收藏 举报
分类:

白书125页例题。

gcd(x,n)所有可能的值是n的约数,要满足gcd(x,n) = i(i是n的约数)必然有gcd (x/i, n/i) = 1,也就是ph(n/i)的值。也就是对于每一个数i,如果i是n的约数,那么i在以n结尾的gcd中出现phi(i)次,这就是i对答案的贡献。

 

#include <bits/stdc++.h>
using namespace std;
#define maxn 4000000

long long phi[maxn+5];
long long f[maxn+5], sum[maxn+5];

void init () { //打表
    for (int i = 2; i <= maxn; i++)
        phi[i] = 0;
    phi[1] = 1;
    for (int i = 2; i <= maxn; i++) if (!phi[i]) {
        for (int j = i; j <= maxn; j += i) {
            if (!phi[j])
                phi[j] = j;
            phi[j] = phi[j] / i * (i-1);
        }
    }
    memset (f, 0, sizeof f);
    for (int i = 1; i <= maxn; i++) {
        for (int j = i+i; j <= maxn; j += i)
            f[j] += i*phi[j/i];
    } 
    sum[2] = 1;
    for (int i = 3; i <= maxn; i++)
        sum[i] = sum[i-1]+f[i];
    return ;
}

int n;

int main () {
    //freopen ("in", "r", stdin);
    init ();
    while (cin >> n && n) {
        cout << sum[n] << endl;
    }
    return 0;
}

0
0

猜你在找
【直播】机器学习&数据挖掘7周实训--韦玮
【套餐】系统集成项目管理工程师顺利通关--徐朋
【直播】3小时掌握Docker最佳实战-徐西宁
【套餐】机器学习系列套餐(算法+实战)--唐宇迪
【直播】计算机视觉原理及实战--屈教授
【套餐】微信订阅号+服务号Java版 v2.0--翟东平
【直播】机器学习之矩阵--黄博士
【套餐】微信订阅号+服务号Java版 v2.0--翟东平
【直播】机器学习之凸优化--马博士
【套餐】Javascript 设计模式实战--曾亮
查看评论
* 以上用户言论只代表其个人观点,不代表CSDN网站的观点或立场
    个人资料
    • 访问:91656次
    • 积分:5349
    • 等级:
    • 排名:第5013名
    • 原创:463篇
    • 转载:0篇
    • 译文:0篇
    • 评论:11条
    最新评论