正题


大意

给出一个nnn,求1∼n" role="presentation" style="position: relative;">1∼n1∼n1\sim n里面和nnn互质的数的和。


解题思路

φ(i)" role="presentation" style="position: relative;">φ(i)φ(i)\varphi(i)表示的就是1∼i1∼i1\sim i里面和iii互质的个数。然后数字和就是φ(i)∗i/2" role="presentation" style="position: relative;">φ(i)∗i/2φ(i)∗i/2\varphi(i)*i/2。


代码

#include<cstdio>
#include<cmath>
#include<algorithm>
#define ll long long
using namespace std;
ll a,ans,n;
int main()
{scanf("%lld",&a);n=a;ans=a;for (ll i=2;i<=sqrt(a);i++)if (a%i==0){ans=ans/i*(i-1);//计算while(a%i==0) a/=i;//取出质数}if (a>1) ans=ans/a*(a-1);//处理a=2的特判printf("%lld",n*ans/2);
}

jzoj1164-求和【欧拉函数,数论】(筛欧拉函数模板)相关推荐

  1. matlab狄利克雷函数,数论入门1——莫比乌斯函数,欧拉函数,狄利克雷卷积,线性筛,莫比乌斯反演,杜教筛...

    数论入门1 一个菜鸡对数论的一点点理解... 莫比乌斯函数 定义函数$\mu(n)$为: 当n有平方因子时,$\mu(n)=0$. 当n没有平方因子时,$\mu(n)=(-1)^{\omega(n)} ...

  2. 中高级数论 [欧拉函数线性筛,二次剩余]

    欧拉函数线性筛 对于素数ppp, φ(p∗i)={p−1i=1p∗φ(i)p∣i(p−1)∗φ(i)p∤i\varphi (p*i)= \begin{cases} p-1& i=1\\ p*\ ...

  3. 数学/数论专题:莫比乌斯函数与欧拉函数

    数学/数论专题:莫比乌斯函数与欧拉函数(进阶) 0. 前言 1. 前置知识 2. 正文 3. 总结 4. 参考资料 0. 前言 本篇文章会从狄利克雷卷积的角度,讨论莫比乌斯函数与欧拉函数的相关性质. ...

  4. 数论一之定理证明——裴蜀/威尔逊/费马/扩展欧几里得/[扩展]欧拉/[扩展]中国剩余定理,欧拉函数,逆元,剩余系,筛法

    打死没想到会在H老师处学懂数论 同余,整除 模运算 埃式筛法 欧拉筛法 最大公约数和最小公倍数 辗转相除法 更相减损术 裴蜀定理 威尔逊定理 费马定理 同余等价类.剩余系.缩系 欧拉函数 欧拉定理 扩 ...

  5. 筛表合集(素数筛 欧拉函数筛 莫比乌斯函数筛)

    [目录] 一.素数筛 1.素数判断 2.素数普通筛 3.素数线性筛 4.素数区间筛 二.欧拉函数筛 三.莫比乌斯函数筛 [素数筛] 1.直接判定质数 bool judgePrime( int num ...

  6. 【数学知识】三种方法求 [1,n] 中所有数欧拉函数(线性筛欧拉函数优化至 O(n) )

    整理的算法模板合集: ACM模板 ①直接求小于或等于n,且与n互质的数个数(求[1,n]中所有数的欧拉函数时间复杂度:O(nn)O(n\sqrt{n})O(nn​)) ②求[1,n]之间每个数的质因数 ...

  7. 数学--数论-数论函数-欧拉函数

    **欧拉函数定义 对正整数n,欧拉函数是少于或等于n的数中与n互质的数的数目.例如euler(8)=4,因为1,3,5,7均和8互质. Euler函数表达通式: 其中p1,p2--pn为x的所有素因数 ...

  8. [模板] 线筛欧拉函数

    今天复习了一些数学方面的东西,线性筛欧拉是个很重要的东西,先贴代码,免得以后又忘了. #include <iostream> #include <cstdlib> typede ...

  9. 埃氏筛 线性筛(欧拉筛) 算法解析

    埃氏晒 埃拉托斯特尼筛法,简称埃氏晒,是一种用来求自然数n以内的全部素数. 他的基本原理是,如果我们要获得小于n的所有素数,那就把不大于根号n的所有素数的倍数剔除. 埃氏晒的原理很容易理解,一个合数, ...

  10. 素数的线性筛法java,埃氏筛 线性筛(欧拉筛) 算法解析

    埃氏晒 埃拉托斯特尼筛法,简称埃氏晒,是一种用来求自然数n以内的全部素数. 他的基本原理是,如果我们要获得小于n的所有素数,那就把不大于根号n的所有素数的倍数剔除. 埃氏晒的原理很容易理解,一个合数, ...

最新文章

  1. 快速原型工具 原型可视化
  2. linux恢复出厂设置_怎么恢复tp-link路由器出厂设置 恢复tp-link出厂设置方法【详解】...
  3. 26Java中的多态
  4. Vue 全家桶 + Express 实现的博客(后端API全部自己手写)
  5. php js登录,php+js实现单点登录
  6. TCP、UDP相关协议使用的端口号
  7. parallels desktop 缺少组件_德国警方扣留光伏组件盗贼
  8. Careercup - Microsoft面试题 - 5752271719628800
  9. 判断release模式_AbstractQueuedSynchronizer共享模式与基于Condition的等待/通知
  10. 关于VS2017编译UE4项目,输出窗口有乱码的解决方案
  11. prometheus命令_Prometheus
  12. OpenCV图像处理(18)——文件夹下所有图像转灰度(14-15综合)
  13. VS Code 的粘性滚动预览 - 类似于 Excel 的冻结首行
  14. ERP-非财务人员的财务培训教(一.二)------财务基础知识
  15. HYSBZ 3991 寻宝游戏
  16. 三个灭点来衡量一个立方体
  17. 给 perl hacker 们看的(ReShip)
  18. java url解码解不了_java – 为什么URL没有完全解码?
  19. 【手把手带你刷Leetcode力扣】10.数据结构 -图
  20. blender使用小技巧

热门文章

  1. java .this的用法_JAVA中this用法小结
  2. postgis数据库优化_国内首个 Serverless 数据库来了,技术架构全揭秘!
  3. eval() php,js-eval编码,js-eval解码
  4. 华为的型号命名规则_华为最实惠5G手机来了!畅享Z 5G宣布:5月24日发
  5. sql datetime 排序_超全的数据库建表/SQL/索引规范,建议贴在工位上!
  6. java实用教程——组件及事件处理——设置组件的位置(相对于窗口具体位置和布局)
  7. 7-27 家谱处理 (30 分)(详解+map做法)map真香啊
  8. C++ 学习之旅(7)——指针pointer
  9. 「3.4w字」超保姆级教程带你实现Promise的核心功能
  10. 一份【软件工程】的学习指南已到达,请注意查收!!