欧拉函数应用

原创 2017年07月27日 21:50:19

题目:https://vjudge.net/contest/173827#problem/E(UVA:11426)

这题真是值得反省自己,之前做过一遍,比赛的时候硬是没想起来怎么做的,完全没有思路,这就很难受,很绝望,所以这回认认真真的把这个题又补了一遍。。。

题意:

求sum(gcd(i,j),1<= i < j < =n) 1 < n< 4000001

问题转化成怎么求f(n),对于一个n来说,枚举因子乘上个数即可。

我们假设b[n]表示1到n-1与n的gcd的和,那么G[n]=G[n-1]+b[n];

a[i]表示与gcd(n, x)= i 的x的个数;b[n]=sum( a[i] * i ) , 所以我们只需求a[i]即可;根据gcd(n, x)=i —–>gcd(n/i, x/i) = 1,

因此仅仅要求出欧拉函数phi(n / i),就能够得到与n / i互质的个数,从而求出gcd(x , n) = i的个数,这样总体就能够求解了

详见博客:(1)http://www.cnblogs.com/zhengguiping–9876/p/4998848.html

(2)http://blog.csdn.net/hyogahyoga/article/details/8520895

补充:对于欧拉函数,根据分解定理可以求得phi【n】,但是要求每个人数的欧拉函数值的时候,可以用欧拉打表法,就不需要一个一个的求了,打表法如下:

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);
            }
        }
    }

关于本题的博客如下:

#include <iostream>
#include <cstring>
#include <cstdio>

using namespace std;
typedef long long ll;
const ll maxn=4000001;
ll phi[maxn],s[maxn],b[maxn];
int main()
{
    ll n;
    memset(phi,0,sizeof(phi));
    memset(s,0,sizeof(s));
    memset(b,0,sizeof(b));
    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);
            }
        }
    }
    for(int i=1;i<maxn;i++)
        for(int j=i+i;j<maxn;j+=i)
          b[j]+=i*phi[j/i];
    for(int i=2;i<maxn;i++)
        s[i]=s[i-1]+b[i];
    while(scanf("%lld",&n)&&n)
    {
        printf("%lld\n",s[n]);
    }
    return 0;
}

(今天被教主说除了数学啥都不会,,,就很难受,,于是我决定最近去学学图论和数据结构,先去学一个周的线段树和后缀数组,树状数组,,,,)

版权声明:本文为博主原创文章,未经博主允许不得转载。 举报

相关文章推荐

codeforces round 17 D(扩展欧拉函数的应用)

链接:http://codeforces.com/contest/17/problem/D 题意,给你三个数,b,n,c 求b*(b-1)^(n-1)%c,结果等于0输出c,不等于0输出结果; ...

扩展欧几里得、求乘法逆元及其应用、中国剩余定理(互质版和非互质版)、欧拉函数、快速判素数模板

互质版: #include #include #include using namespace std; typedef __int64 int64; int64 a[15...

我是如何成为一名python大咖的?

人生苦短,都说必须python,那么我分享下我是如何从小白成为Python资深开发者的吧。2014年我大学刚毕业..

poj 2480 (欧拉函数应用)

点击打开链接 //求SUM(gcd(i,n), 1<=i<=n) /* g(n)=gcd(i,n),根据积性定义g(mn)=g(m)*g(n)(gcd(m,n)==1) 所以gcd(i,...

poj2773(欧拉函数的应用 或者 用容斥原理+dfs+二分)

题意:给出m和k,求与m互质的第k个数 要用到的zh

HDU 1395 欧拉函数的基本应用

题目给的很明显就是欧拉函数的表达式,由于底是2,所以被取模的数只要是奇数,就能保证与2互质,而定义在一数集上的取模,是一个循环群, 若2^n%p==1, phi(p)是一个循环节,单不能保证它是最...

HDU 2588 NOJ 1126 GCD 欧拉函数应用问题

HDU 2588 NOJ 1126 GCD 欧拉函数应用问题 欧拉函数变形

HDU 2588(欧拉函数的应用)

GCD Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Total Submis...

poj 2773 Happy 2006(欧拉函数应用)

http://poj.org/problem?id=2773 题意:输入n,k,求与n不互素的第k个数,k可能大于n。 思路:以n=6为例,与6互素的数有一定规律。{1,5},{7...

欧拉函数应用

在n*n平面上有n*n-m-1(除原点与另m个点),可以从原点引k条线,询问最多能穿过多少点。 设法最简表示每种方案(i,j)->i,j互质,容易想到最多有(phi[2]+phi[3]+...+ph...

2773 Happy 2006 //欧拉函数的应用

Happy 2006Time Limit: 3000MS Memory Limit: 65536KTotal Submissions: 5258 Accepted: 1576DescriptionTw...
返回顶部
收藏助手
不良信息举报
您举报文章:深度学习:神经网络中的前向传播和反向传播算法推导
举报原因:
原因补充:

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