欧拉函数和遗忘自动机 SX 的故逝
欧拉函数 \(\varphi(n)\) 定义为小于 \(n\) 与 \(n\) 互质的数字,炒个例子,\(\varphi(10) = 4\),因为 \(1,3,7,9\) 与 \(10\) 互质。
怎么求?\(n\) 唯一分解成 \(p_1^{c_1}p_2^{c_2}\dots p_k^{c_k}\),则 \(\varphi(n) = n\prod\limits_{i = 1}^k \dfrac{p_i - 1}{p_i}\),证明可以用容斥做(懒得写了┓( ´∀` )┏)
它显然可以用分解质因数求单个的,也可以用埃筛求多个的,长这样
void getphi() {
for(int p = 1; p <= MAXN; p++) phi[p] = p;
for(int p = 2; p <= MAXN; p++)
if(phi[p] == p)
for(int i = p; i <= MAXN; i += p)
phi[i] = phi[i] / p * (p - 1);
}
原因显而易见。
最重要的是几条性质(也就是说我水了一半是吧 QAQ)
与 \(n\) 互质的数之和为 \(\dfrac{\varphi(n)n}{2}\)。
这是因为 \(\gcd(n - x, n) = \gcd(n, x)\),则如果 \(x\) 与 \(n\) 互质,那么 \(n - x\) 亦与 \(n\) 互质,成对出现,和为 \(n\),有 \(\dfrac{\varphi(n)}{2}\),总和为 \(\dfrac{\varphi(n)n}{2}\)。欧拉函数是积性函数,即 \(\varphi(ab) = \varphi(a)\varphi(b),\gcd(a, b) = 1\)
显而易见,分解质因数即可证明。若 \(p\) 为质数,且 \(p|n, p^2|n\),则有 \(\varphi(n) = \varphi(n/p) \times p\)。
由于 \(p|n, p^2|n\),因此 \(n/p\) 与 \(n\) 有同样的质因子,\(\varphi(n)\) 与 \(\varphi(n/p)\) 用公式写一下显然成立(若 \(p\) 为质数,且 \(p|n, p^2\not | n\),则有 \(\varphi(n) = \varphi(n/p) \times (p - 1)\)。
由于 \(p|n, p^2\not|n\),则 \(n/p\) 没有质因子 \(p\),即 \(\gcd(n/p, p) = 1\),由于这是积性函数则 \(\varphi(n) = \varphi(n / p)\times \varphi(p) = \varphi(n/ p) \times (p - 1)\)。
运用 3 4 可以用欧拉筛筛这个 qwq
- \(\sum\limits_{d|n}\varphi(d) = n\)
欧拉。。欧拉反演。。。?
不会证明,懒得去学,跑路,有时间填坑 qAq
欧拉筛怎么筛这个(欧拉超级加倍 2333333)
欧拉筛的核心在于将一个合数拆成 \(p\times k\),其中 \(p\) 是这个数字的最小质因数。如果 \(p\) 与 \(k\) 互质,则有 \(\varphi(pk) = \varphi(p)\times (k - 1)\),否则 \(varphi(pk) = \varphi(p/k)\times p\)。原因是因为上面的 3 4 两条性质 QWQ
代码不放了,都讲到这个份上了就都会了~(主要是暂时找不到懒得写了┓( ´∀` )┏(懒狗。。。))
式子的话,一般是求和,然后我们每次是枚举 \(\gcd\),最后一波乱搞搞成欧拉函数的形式 qwq
因题而异,不写了 qwq
欧拉函数和遗忘自动机 SX 的故逝的更多相关文章
- 【bzoj2401】陶陶的难题I “高精度”+欧拉函数+线性筛
题目描述 求 输入 第一行包含一个正整数T,表示有T组测试数据.接下来T<=10^5行,每行给出一个正整数N,N<=10^6. 输出 包含T行,依次给出对应的答案. 样例输入 7 1 10 ...
- hdu2588 GCD (欧拉函数)
GCD 题意:输入N,M(2<=N<=1000000000, 1<=M<=N), 设1<=X<=N,求使gcd(X,N)>=M的X的个数. (文末有题) 知 ...
- BZOJ 2705: [SDOI2012]Longge的问题 [欧拉函数]
2705: [SDOI2012]Longge的问题 Time Limit: 3 Sec Memory Limit: 128 MBSubmit: 2553 Solved: 1565[Submit][ ...
- BZOJ 2818: Gcd [欧拉函数 质数 线性筛]【学习笔记】
2818: Gcd Time Limit: 10 Sec Memory Limit: 256 MBSubmit: 4436 Solved: 1957[Submit][Status][Discuss ...
- COGS2531. [HZOI 2016]函数的美 打表+欧拉函数
题目:http://cogs.pw/cogs/problem/problem.php?pid=2533 这道题考察打表观察规律. 发现对f的定义实际是递归式的 f(n,k) = f(0,f(n-1,k ...
- poj2478 Farey Sequence (欧拉函数)
Farey Sequence 题意:给定一个数n,求在[1,n]这个范围内两两互质的数的个数.(转化为给定一个数n,比n小且与n互质的数的个数) 知识点: 欧拉函数: 普通求法: int Euler( ...
- 51Nod-1136 欧拉函数
51Nod: http://www.51nod.com/onlineJudge/questionCode.html#!problemId=1136 1136 欧拉函数 基准时间限制:1 秒 空间限制: ...
- 欧拉函数 - HDU1286
欧拉函数的作用: 有[1,2.....n]这样一个集合,f(n)=这个集合中与n互质的元素的个数.欧拉函数描述了一些列与这个f(n)有关的一些性质,如下: 1.令p为一个素数,n = p ^ k,则 ...
- FZU 1759 欧拉函数 降幂公式
Description Given A,B,C, You should quickly calculate the result of A^B mod C. (1<=A,C<=1000 ...
- hdu 3307 Description has only two Sentences (欧拉函数+快速幂)
Description has only two SentencesTime Limit: 3000/1000 MS (Java/Others) Memory Limit: 65536/32768 K ...
随机推荐
- ArcEngine要素编辑遇到的一些问题
1.如何开启编辑 IMap myMap = this._Aplication.ActiveView.FocusMap; IWorkspace myWorkspace = (myMap25Sheet.P ...
- 【Java SE】Day05数组
一.数组的定义和访问 1.初始化 动态new int[10];--默认值 静态new int[]{1,2,3};,省略为{1,2,3}; 2.访问 长度arr.length属性(数组的属性) 打印数组 ...
- 【Day03】Spring cloud:源码讲解与容器化初探
今日内容 原理和源码 容器化过度 一.Naocs 1.介绍 server端 启动入口类(Spring Boot项目,提供8848端口的监听访问) 源码包含InstanceController类(ser ...
- 单节锂电池充电管理芯片,IC电路图
PW4054 是一款性能优异的单节锂离子电池恒流/恒压线性充电器.PW4054 适合给 USB 电源以及适配器电源供电.基于特殊的内部 MOSFET 架构以及防倒充电路, PW4054 不需要外接检测 ...
- JavaScript冒泡排序+Vue可视化冒泡动画
冒泡排序(Bubble Sort)算是前端最简单的算法,也是最经典的排序算法了.网上JavaScript版本的冒泡排序很多,今天用Vue实现一个动态的可视化冒泡排序. 01.JavaScript冒泡排 ...
- element-ui中table表格表头和表格内容都水平居中,以及斑马纹背景颜色修改
<el-table :data="detalData" stripe //斑马纹 border :header-cell-style="{textAlign: 'c ...
- 一个开放源代码,实现动态IL注入(Hook或补丁工具)框架:Lib.Harmony(Patch,PatchAll,Prefix,Postfix,Transpiler)
详情请参考原文:一个开放源代码,实现动态IL注入(Hook或补丁工具)框架:Lib.Harmony
- [深度学习] tf.keras入门1-基本函数介绍
目录 构建一个简单的模型 序贯(Sequential)模型 网络层的构造 模型训练和参数评价 模型训练 模型的训练 tf.data的数据集 模型评估和预测 基本模型的建立 网络层模型 模型子类函数构建 ...
- 经典问题 1 —— DAG 上区间限制拓扑序
问题描述 给定一个 DAG,求一个拓扑序,使得节点 \(i\) 的拓扑序 \(\in [l_i,r_i]\). 题解 首先进行一个预处理:对于所有 \(u\),令 \(\forall (v,u)\in ...
- 《深入理解Java虚拟机》第三章读书笔记(一)——垃圾回收算法
参考书籍<深入理解java虚拟机>周志明著 系列文章目录和关于我 本文主要介绍垃圾回收理论知识 1.jvm哪些区域需要进行垃圾回收 虚拟机栈,本地方法栈,程序计数器都是线程私有的,随线程而 ...