繁凡さん的博客

原创 【学习笔记】超简单的多项式求逆（含全套证明）

整理的算法模板合集： ACM模板目录多项式求逆一、分治FFT二、倍增法及其证明三、多项式求逆例题[P4238 【模板】多项式乘法逆](https://www.luogu.com.cn/problem/P4238)点我看多项式全家桶(●^◡_◡◡​^●)多项式求逆一、分治FFT二、倍增法及其证明三、多项式求逆例题P4238 【模板】多项式乘法逆直接按照上述的思路实现一下模板即可。注意要中间 h(x)h(x)h(x) 要取的是模 x⌈n⌉2x^{\frac{\lceil n \rceil}

原创 【学习笔记】多项式全家桶（包含全套证明）

目录0x00.多项式基本概述0x01.多项式的基础概念0x02. 多项式的度0x03. 多项式的乘法0x04. 多项式的逆元0x10.拉格朗日插值0x20.多项式乘法0x21.快速傅里叶变换（FFT）0x22.快速数论变换（NTT）0x23.快速沃尔什变换（FWT）0x24.快速莫比乌斯变换（FMT）0x30.多项式求逆0x40.多项式开方0x50.多项式除法 / 取模0x60.多项式牛顿迭代0x70.多项式多点求值 / 快速插值0x80.多项式三角函数0x90.多项式反...

原创 【学习笔记】超简单的快速数论变换（NTT）（FFT的优化）（含全套证明）

整理的算法模板合集： ACM模板目录一、前置知识二、快速数论变换（NTT）三、NTT证明（和FFT的关系）四、NTT模板点我看多项式全家桶(●’◡’●)（全家桶还没写完，待更链接hhh ）一、前置知识快速傅里叶变换FTT（学完FFT再看NTT哈，）原根（链接待更，下面有原根的基础说明，足够今天的NTT使用了）二、快速数论变换（NTT）三、NTT证明（和FFT的关系）我们发现FFT里的大多操作，都跟单位根没有什么关系，（例如选点插值，奇偶分治），我们随便选择一个n个点v,v2

原创 【几何专题】计算几何基础

整理的算法模板合集： ACM模板目录AcWing 2984. 线段 / POJ3304（计算几何、常用技巧转换）AcWing 2983. 玩具 / POJ 2318.toys（计算几何基础、二分、判断点和直线的位置关系）AcWing 2984. 线段 / POJ3304（计算几何、常用技巧转换）常用技巧转换，寻找一个直线旋转我们发现如果可以找到一个直线与所有的线段都相交，那么我们做一个垂直这个直线的垂线，所有的线段的投影一定都交于这个直线与垂线的垂足处，我们可以很轻松地画图验证这一定理。所

原创 【学习笔记】拉格朗日插值

整理的算法模板合集： ACM模板目录P4781 【模板】拉格朗日插值点我看多项式全家桶(●’◡’●)（全家桶还没写完，待更链接hhh ）P4781 【模板】拉格朗日插值我们按照上面的思路，带入n个点求一下f(k)f(k)f(k)即可。时间复杂度O(n2)O(n^2)O(n2)。注意本题还要求逆元，为了防止求逆元的时间复杂度影响整体的时间复杂度，所以我们分别计算出分子和分母，再将分子乘进分母的逆元，累加进最后的答案，时间复杂度的瓶颈就不会在求逆元上，总体的时间复杂度为O(n2)O(n^2)O

原创 【学习笔记】超简单的快速傅里叶变换（FFT）（含全套证明）

整理的算法模板合集： ACM模板目录多项式复数复数的单位根 / 单位向量离散傅里叶变换（DTF）离散傅里叶逆变换（IDTF）FFT算法整体流程代码实现递归版迭代版“三步变两步”优化P1919 【模板】FFT快速傅里叶变换学习笔记（用ipad记录的笔记）：超级简单的快速傅里叶变换！只要基础够扎实，顺着推一遍没有什么难以理解的，我学的整个过程没有一点卡壳，真的很爽，整整写了三个小时，写了满满6页的笔记。主要是内容太多了。首先是一些基础概念：多项式复数复数的单位根 / 单位向量离散傅里叶变

原创 构造莫比乌斯函数

原创 【学习笔记】杜教筛

学习笔记：杜教筛其实特别简单，就是一个构造+和式转换，利用Dirichelt卷积构造两个积性函数卷起来，将要求的前缀和s(n)s(n)s(n)构造成s(⌊ni⌋)s(⌊ \frac{n}{i}⌋)s(⌊in​⌋)的形式，这样我们就可以用整除分块来优化复杂度，可以快速解决一类积性函数的前缀和，n可以达到 1e9 ~ 1e10 ，积性函数比如莫比乌斯函数μ(x)μ(x)μ(x)，欧拉函数ϕ(x)\phi(x)ϕ(x)，约数和函数σk(x)σk(x)σk(x)，约数个数函数σ(x)σ(x)σ(x)，黎曼函数ζ

原创 P1444 [USACO1.3]虫洞 wormhole，新手向超详细讲解（搜索、判环、拆点、好题）

题目传送门首先我们根据题意发现奶牛一共只会出现两种动作（两种情况），一是向右走，二是传送。我们并不知道奶牛的位置，而题目所求的是有多少种分配虫洞传送的分配方案，使得奶牛可以陷入死循环，所以奶牛的位置是没有用的，或者说奶牛的位置可以在任何位置，只要分配的虫洞方案存在死循环的可能。所以我们只需要考虑所有的分配方案即可。我们发现奶牛除了传送以外，只能往右走。也就是说我们所有的y相同，而x不相同的点之间是可以走到的。所以我们先考虑如何得到所有的可以通过往右走连接起来的点。我们只需要将y作为第一关键字，x作

原创 AcWing 2984. 线段 / POJ3304（计算几何、常用技巧转换）

常用技巧转换，寻找一个直线旋转我们发现如果可以找到一个直线与所有的线段都相交，那么我们做一个垂直这个直线的垂线，所有的线段的投影一定都交于这个直线与垂线的垂足处，我们可以很轻松地画图验证这一定理。所以现在问题就变成了找一个可以穿过所有线段的直线。我们肯定不能暴力找所有的坐标判断是否相交。这里我们使用一个计算几何常用的技巧：我们发现所有可能的答案，至少从与一个线段有交点开始。所以我们可以先找一个过一个线段的直线，然后我们旋转他（传统艺能）。我们把这个直线的定点定为这个线段的端点，然后旋转。因为我们.

原创 AcWing 2983. 玩具 / POJ 2318.toys（计算几何基础、二分、判断点和直线的位置关系）

计算几何基础题题目大意就是一个盒子被分成了若干个区域，有m个小球，问每个区域里分别有多少个小球。我们首先考虑暴力。我们如何判断小球是否在一个区域内呢，我们发现一个小球在区域x，也就是说所有小于x的所有隔板都在小球坐标的左边，所有大于x的隔板都在小球坐标的右边，所以我们可以暴力n2n^2n2枚举每个小球的坐标，判断有多少个隔板（线段）在小球的左边，个数就是答案（区域从0到n）。数据5000，10组数据，n2n^2n2暴力数据达到了2∗1082*10^{8}2∗108，2亿，不可过，所以我们考虑优化。.

原创 组合恒等式

组合数恒等式常用组合数计算公式（转载）组合数恒等式

原创 Dirichlet前缀和及其拓展

整理的算法模板合集： ACM模板目录DirichletDirichletDirichlet 前缀和DirichletDirichletDirichlet 后缀和倒推 DirichletDirichletDirichlet 前缀和倒推 DirichletDirichletDirichlet 后缀和DirichletDirichletDirichlet 前缀和题目传送门题目大意：给定一个长度为 nnn 的数列 a1,a2,a3,…,ana_1,a_2,a_3,\dots,a_na1​,a2​,a3​

原创 SP5971 LCMSUM - LCM Sum（莫比乌斯反演 ，推柿子，经典）

整理的算法模板合集： ACM模板题目传动门#include <iostream>#include <cstdio>#include <algorithm>#include <queue>#include <cmath>#include <cstring>#include <vector>using namespace std;typedef long long ll;typedef pair&l

原创 模板 - 快速求n所有因数（快于暴力sqrt(n) ）

整理的算法模板合集： ACM模板P1072 Hankson 的趣味题题目大意：已知正整数a0,a1,b0,b1a_0,a_1,b_0,b_1a0​,a1​,b0​,b1​ ，设某未知正整数 xxx 满足：1． xxx 和 a0a_0a0​的最大公约数是 a1a_1a1​ ；2． xxx 和 b0b_0b0​ 的最小公倍数是 b1b_1b1​求满足条件的xxx的个数。优化做法可以在O(nlog n)O(\frac{\sqrt{n}}{log\ n})O(log nn​​)的

原创 【题单 - 数学专题】最大公约数

整理的算法模板合集： ACM模板目录Part 6.2.2 最大公约数P5435 基于值域预处理的快速 GCD算法P5436 【XR-2】缘分P1029 最大公约数和最小公倍数问题P1414 又是毕业季IIP2152 [SDOI2009]SuperGCDP1072 Hankson 的趣味题Part 6.2.2 最大公约数如果两个数有一个共同的约数，那么这个约数就被称为公约数。最大公约数就是指这两个数的所有公约数中，最大的一个。求解两个数的最大公约数，可以采用欧几里得算法解决。P5435 【模

原创 【数学专题】整除相关 - 素数

Part 6.2 整除相关与整除相关的概念有很多，比较常用的有素数，最大公约数和欧拉函数。Part 6.2.1 素数素数，指的是除 1 和它本身之外没有其他约数的数。P4718 【模板】Pollard-Rho算法P1075 质因数分解P2441 角色属性树P5535 【XR-3】小道消息.P4718 【模板】Pollard-Rho算法直接上模板#include <bits/stdc++.h>using namespace std;#define rg reg

原创 模板 - 莫比乌斯反演（常用技巧）

整理的算法模板合集： ACM模板目录莫比乌斯反演常用技巧经典模板例题莫比乌斯反演莫比乌斯函数：μ(n)={1n=1(−1)kn=p1p2p3…pk,k为n的本质不同的质因子的个数0其他情况n含有平方因子\mu(n) = \left\{ \begin{matrix} 1 & n=1\\ (-1)^k & n = p_1p_2p_3\dots p_k,k为n的本质不同的质因子的个数\\ 0 & \text{其他情况n含有平方因子} \end{m

原创 luoguP2522 [HAOI2011]Problem b（莫比乌斯反演 + 数论分块 + 容斥原理）

#include<iostream>#include<cstdio>#include<algorithm>#include<queue>#include<cmath>#include<cstring>using namespace std;const int N = 500007, M = 500007,INF = 0x3f3f3f3f;typedef long long ll;int read(){ ..

原创 第十二届西南石油大学程序设计新生赛官方题解

比赛地址：https://ac.nowcoder.com/acm/contest/9681#question邀请码：swpu2020大家可以点开上面的比赛地址进入比赛页面提交代码补题。A、杰杰国王的平等节（easy）签到题，只需要找到序列最大值计算即可。#include <iostream>using namespace std;int a[1000005];int main() { int n, m = 0, sum = 0; scanf("%d", &amp

原创 积性函数、狄利克雷卷积、莫比乌斯反演

积性函数、狄利克雷卷积、莫比乌斯反演

原创 【数学专题】组合数学与计数

整理的算法模板合集： ACM模板目录#include <cstdio>#include <algorithm>#include <iostream>#include <cstring>using namespace std;const int N = 100007, mod = 5000011;int n, k;int f[N];int s[N];int main(){ scanf("%d%d", &n,