poj 2480 欧拉函数+积性函数+GCD

最新推荐文章于 2020-12-04 21:12:13 发布

Z-Pilgrim

最新推荐文章于 2020-12-04 21:12:13 发布

阅读量1.3k

点赞数 1

分类专栏：数学------------数论

本文链接：https://blog.csdn.net/u011026968/article/details/25485001

版权

数学------------数论专栏收录该内容

20 篇文章 3 订阅

订阅专栏

题目：http://poj.org/problem?id=2480

首先要会欧拉函数：先贴欧拉函数的模板，来源于吉林大学的模板：

//欧拉函数PHI(n)表示的是比n小，并且与n互质的正整数的个数(包括1)。
unsigned euler(unsignedx)
{// 就是公式
    unsigned i, res=x;
    for(i = 2; i < (int)sqrt(x * 1.0) + 1; i++)
        if(x%i==0) {
            res = res / i * (i - 1);
            while(x % i == 0) x /= i; // 保证i一定是素数
    }
    if(x > 1) res = res / x * (x - 1);
    return res;
}

然后欧拉函数在本题中用到的一些性质：
1. 当 n 为素数时， PHI(n) = n - 1 。因为每个比 n 小的正整数都和 n 互素。当 n 为素数 p 的 k 次方时， PHI(n) = p ^ k - p ^ (k - 1) 。因为在 1 到 n 之间的正整数只有 p 的倍数和 n 不互素，这样的数有 (p ^ k / p) 个。
2. 如果 m 和 n 互素，即 GCD(m, n) = 1 ，那么 PHI(m * n) = PHI(m) * PHI(n) 。

然后说本题：
我参考了这个题解，代码也看了o(╯□╰)o，不过人家写的很好啊，膜拜~~

http://blog.csdn.net/terry__j/article/details/7403923

然后贴我的代码：

/****************************************\
欧拉函数+积性函数
    by Pilgrim 2014.5.10
注意：1、//ans=ans*(1+ai*(i-1)/i);--WA,二逼了
         ai/i很可能不是整除，所以这么写
          ans=ans+ans*ai*(i-1)/i;
\****************************************/

#include <iostream>
#include <cstdio>
#include <string>
#include <cstring>
#include <cstdlib>
#include <cmath>
#include <climits>
#include <cmath>

using namespace std;

#define lint long long

int main()
{
    lint n,sqr,tmp,ai,ans,ret;

    while(scanf("%lld",&n)!=EOF)
    {
        tmp=n;
        ans=n;
        sqr = (lint)sqrt(n*1.0)+1;
        for(lint i=2;i<sqr;i++)
        {
            if(tmp%i == 0)
            {
                tmp/=i,ai=1;
                while(tmp%i == 0)tmp/=i,ai++;

                //ans=ans*(1+ai*(i-1)/i);
                ans=ans+ans*ai*(i-1)/i;
            }
        }
        if(tmp!=1)//ans=ans*(1+(tmp-1)/tmp);
            ans=ans+ans*(tmp-1)/tmp;
        printf("%lld\n",ans);
    }

    return 0;
}