ICode9

精准搜索请尝试: 精确搜索
首页 > 其他分享> 文章详细

Acwing 数论打卡(2)

2021-04-03 16:35:51  阅读:200  来源: 互联网

标签:www eulers 数论 res int 打卡 include 欧拉 Acwing


873. 欧拉函数

https://www.acwing.com/video/298/
欧拉函数定义有f(x)=n 表示1到x中与x互质的数有n个

#include <iostream>
#include<algorithm>
#include<cstdio>
#include<vector>
#include<queue>
#include<stack>
#include<cstring>
using namespace std;
int fun(int x){
    int res=x;
    for(int i=2;i<=x/i;i++){
        if(x%i==0){
            res=res/i*(i-1);//先除后乘 避免数据溢出
            while(x%i==0) x/=i;
        }
    }
    if(x>1) res=res/x*(x-1);
    return res;
}
int main(){
    int n;
    cin>>n;
    while(n--){
        int a;
        scanf("%d",&a);
        cout<<fun(a)<<endl;
    }
  return 0;
}
//  freopen("testdata.in", "r", stdin);

874. 筛法求欧拉函数

https://www.acwing.com/problem/content/876/

#include <iostream>

using namespace std;

typedef long long LL;

const int N = 1000010;


int primes[N], cnt;
int euler[N];
bool st[N];


void get_eulers(int n)
{
    euler[1] = 1;
    for (int i = 2; i <= n; i ++ )
    {
        if (!st[i])
        {
            primes[cnt ++ ] = i;
            //如果p是素数 那么欧拉函数值为p-1
            euler[i] = i - 1;
        }
        for (int j = 0; primes[j] <= n / i; j ++ )
        {
            int t = primes[j] * i;
            st[t] = true;
            if (i % primes[j] == 0)
            {
                euler[t] = euler[i] * primes[j];
                //因为pj是i的一个质因数 所以pj*i的质因数和i的质因数完全一样
                //欧拉函数f(n)=n*(1-1/p1)---(1-1/pn);
                //f(i)=i*(1-1/p1)---(1-1/pn);
                //f(pj*i)=pj*i*(1-1/p1)---(1-1/pn);
                //由于p1-pn完全相同 所以有f(pj*i)=pj*f(i);
                break;
            }
            //i%pj!=0
            //i*pj比i多了一个pj质因数 
            euler[t] = euler[i] * (primes[j] - 1);
        }
    }
}


int main()
{
    int n;
    cin >> n;

    get_eulers(n);

    LL res = 0;
    for (int i = 1; i <= n; i ++ ) res += euler[i];

    cout << res << endl;

    return 0;
}

标签:www,eulers,数论,res,int,打卡,include,欧拉,Acwing
来源: https://www.cnblogs.com/OfflineBoy/p/14614024.html

本站声明: 1. iCode9 技术分享网(下文简称本站)提供的所有内容,仅供技术学习、探讨和分享;
2. 关于本站的所有留言、评论、转载及引用,纯属内容发起人的个人观点,与本站观点和立场无关;
3. 关于本站的所有言论和文字,纯属内容发起人的个人观点,与本站观点和立场无关;
4. 本站文章均是网友提供,不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属;如您发现该文章侵犯了您的权益,可联系我们第一时间进行删除;
5. 本站为非盈利性的个人网站,所有内容不会用来进行牟利,也不会利用任何形式的广告来间接获益,纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

专注分享技术,共同学习,共同进步。侵权联系[81616952@qq.com]

Copyright (C)ICode9.com, All Rights Reserved.

ICode9版权所有