STcyclone的博客

扩展欧几里得算法详解扩展欧几里得算法详解及模板题（洛谷P5656）解题思路题目链接题目描述数据范围扩展欧几里得算法代码模板时间复杂度算法正确性的证明归纳基础归纳步算法时间复杂度的证明b=0b=0b=0 时b≠0b \neq 0b​=0 且 a<ba < ba<b 时b≠0b \neq 0b​=0 且 a≥ba \ge ba≥b 时引理： ∀a,b∈Z+\forall a,b \in \mathbb{Z}^{+}∀a,b∈Z+，若 a≥ba \ge ba≥b，则有 a mod b≤⌊a2

void phi_table(){ phi[1]=1; fo(i,2,mn){//从2开始计算 if (!phi[i]){//等价于是质数 pri[++pt]=i; phi[i]=i-1; } fo(j,1,pt){ //枚举添加的小质数 te=i*pri

scanf("%d",&n); fo(i,2,n){ if (!isp[i]){ pri[++pt]=i; //printf("%d ",i); } fo(j,1,pt){ t=(long long)i*pri[j]; if (t>n)

在数论函数集上，狄利克雷卷积定义为如下的二元运算 (f*g)(n)=d|nΣf(d)g(n/d)狄利克雷卷积运算满足交换律，结合律，分配律在数论函数集上，存在单位元函数ε(n)满足ε(1)=1 ε(n)=0(n>1) 使得(f*ε)=(ε*f)=f而常值函数1(n)则定义为1(n)=1一个数论函数f(n)与1(n)的狄利克雷卷积并不是f(n)自身，而是一个常见的式子 (f*1)(n)=d|n

乘法逆元对于缩系中的元素，每个数a均有唯一的与之对应的乘法逆元x，使得ax≡1(mod n) 一个数有逆元的充分必要条件是gcd(a,n)=1，此时逆元唯一存在 逆元的含义：模n意义下，1个数a如果有逆元x，那么除以a相当于乘以x。下面给出求逆元的几种方法1 循环找解法给定模m和需要求逆的数x，直接暴力枚举1~m-1 检查是否有x*i=1(mod m)这种算法可以应用与写暴力、对拍、模数较小,

题目大意求∑i=1n∑j=1nijgcd(i,j)\sum_{i=1}^{n}\sum_{j=1}^{n}ijgcd(i,j)对质数p取膜的结果n<=10^10 5*10^8<=p<=1.1*10^9分析看到这种两层求和带gcd的式子 可以考虑直接用莫比乌斯反演求解。 但注意到i j的上界都是n，所以用ϕ\phi推导公式也可以解决这题。(感觉可能会简便一些)∑i=1n∑j=1nijgcd(i,