jzoj1164-求和【欧拉函数,数论】(筛欧拉函数模板)
正题
大意
给出一个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-求和【欧拉函数,数论】(筛欧拉函数模板)相关推荐
- matlab狄利克雷函数,数论入门1——莫比乌斯函数,欧拉函数,狄利克雷卷积,线性筛,莫比乌斯反演,杜教筛...
数论入门1 一个菜鸡对数论的一点点理解... 莫比乌斯函数 定义函数$\mu(n)$为: 当n有平方因子时,$\mu(n)=0$. 当n没有平方因子时,$\mu(n)=(-1)^{\omega(n)} ...
- 中高级数论 [欧拉函数线性筛,二次剩余]
欧拉函数线性筛 对于素数ppp, φ(p∗i)={p−1i=1p∗φ(i)p∣i(p−1)∗φ(i)p∤i\varphi (p*i)= \begin{cases} p-1& i=1\\ p*\ ...
- 数学/数论专题:莫比乌斯函数与欧拉函数
数学/数论专题:莫比乌斯函数与欧拉函数(进阶) 0. 前言 1. 前置知识 2. 正文 3. 总结 4. 参考资料 0. 前言 本篇文章会从狄利克雷卷积的角度,讨论莫比乌斯函数与欧拉函数的相关性质. ...
- 数论一之定理证明——裴蜀/威尔逊/费马/扩展欧几里得/[扩展]欧拉/[扩展]中国剩余定理,欧拉函数,逆元,剩余系,筛法
打死没想到会在H老师处学懂数论 同余,整除 模运算 埃式筛法 欧拉筛法 最大公约数和最小公倍数 辗转相除法 更相减损术 裴蜀定理 威尔逊定理 费马定理 同余等价类.剩余系.缩系 欧拉函数 欧拉定理 扩 ...
- 筛表合集(素数筛 欧拉函数筛 莫比乌斯函数筛)
[目录] 一.素数筛 1.素数判断 2.素数普通筛 3.素数线性筛 4.素数区间筛 二.欧拉函数筛 三.莫比乌斯函数筛 [素数筛] 1.直接判定质数 bool judgePrime( int num ...
- 【数学知识】三种方法求 [1,n] 中所有数欧拉函数(线性筛欧拉函数优化至 O(n) )
整理的算法模板合集: ACM模板 ①直接求小于或等于n,且与n互质的数个数(求[1,n]中所有数的欧拉函数时间复杂度:O(nn)O(n\sqrt{n})O(nn)) ②求[1,n]之间每个数的质因数 ...
- 数学--数论-数论函数-欧拉函数
**欧拉函数定义 对正整数n,欧拉函数是少于或等于n的数中与n互质的数的数目.例如euler(8)=4,因为1,3,5,7均和8互质. Euler函数表达通式: 其中p1,p2--pn为x的所有素因数 ...
- [模板] 线筛欧拉函数
今天复习了一些数学方面的东西,线性筛欧拉是个很重要的东西,先贴代码,免得以后又忘了. #include <iostream> #include <cstdlib> typede ...
- 埃氏筛 线性筛(欧拉筛) 算法解析
埃氏晒 埃拉托斯特尼筛法,简称埃氏晒,是一种用来求自然数n以内的全部素数. 他的基本原理是,如果我们要获得小于n的所有素数,那就把不大于根号n的所有素数的倍数剔除. 埃氏晒的原理很容易理解,一个合数, ...
- 素数的线性筛法java,埃氏筛 线性筛(欧拉筛) 算法解析
埃氏晒 埃拉托斯特尼筛法,简称埃氏晒,是一种用来求自然数n以内的全部素数. 他的基本原理是,如果我们要获得小于n的所有素数,那就把不大于根号n的所有素数的倍数剔除. 埃氏晒的原理很容易理解,一个合数, ...
最新文章
- 快速原型工具 原型可视化
- linux恢复出厂设置_怎么恢复tp-link路由器出厂设置 恢复tp-link出厂设置方法【详解】...
- 26Java中的多态
- Vue 全家桶 + Express 实现的博客(后端API全部自己手写)
- php js登录,php+js实现单点登录
- TCP、UDP相关协议使用的端口号
- parallels desktop 缺少组件_德国警方扣留光伏组件盗贼
- Careercup - Microsoft面试题 - 5752271719628800
- 判断release模式_AbstractQueuedSynchronizer共享模式与基于Condition的等待/通知
- 关于VS2017编译UE4项目,输出窗口有乱码的解决方案
- prometheus命令_Prometheus
- OpenCV图像处理(18)——文件夹下所有图像转灰度(14-15综合)
- VS Code 的粘性滚动预览 - 类似于 Excel 的冻结首行
- ERP-非财务人员的财务培训教(一.二)------财务基础知识
- HYSBZ 3991 寻宝游戏
- 三个灭点来衡量一个立方体
- 给 perl hacker 们看的(ReShip)
- java url解码解不了_java – 为什么URL没有完全解码?
- 【手把手带你刷Leetcode力扣】10.数据结构 -图
- blender使用小技巧
热门文章
- java .this的用法_JAVA中this用法小结
- postgis数据库优化_国内首个 Serverless 数据库来了,技术架构全揭秘!
- eval() php,js-eval编码,js-eval解码
- 华为的型号命名规则_华为最实惠5G手机来了!畅享Z 5G宣布:5月24日发
- sql datetime 排序_超全的数据库建表/SQL/索引规范,建议贴在工位上!
- java实用教程——组件及事件处理——设置组件的位置(相对于窗口具体位置和布局)
- 7-27 家谱处理 (30 分)(详解+map做法)map真香啊
- C++ 学习之旅(7)——指针pointer
- 「3.4w字」超保姆级教程带你实现Promise的核心功能
- 一份【软件工程】的学习指南已到达,请注意查收!!