给定 \(n,m,p\)，其中 \(n,m\) 较大，\(p\) 为质数且不是很大，求 \[\dbinom{n}{m}\bmod p \]\(\rm Lucas\) 定理 \[\dbinom{n}{m}\equiv\dbinom{\left\lfloor\frac{n}{p}\right\rfloor}{\left\lfloor\frac{m}{p}\right\rfloor}\cdot \dbinom{n\bmod p}{m\bmod p}
约数个数和 题目链接：luogu P3327 题目大意 给你 n,m 要你求 ∑i=1~n∑j=1~m d(ij)，d 是约数个数。 多组数据。 思路 首先由一个比较神奇的东西： \(d(ij)=\sum\limits_{x|i}\sum\limits_{y|j}[gcd(x,y)=1]\) 这里证明一下： 先看我们可以只看一个质数的部分。 假设它的因子的出现次数是
小小网格 题目链接：ybtoj高效进阶 21177 题目大意 给你求 ∑i=1~n∑j=1~mφ(gcd(i,j))。 思路 重新写好看点： \(\sum\limits_{i=1}^n\sum\limits_{j=1}^m\varphi(\gcd(i,j))\) \(\sum\limits_{d}\varphi(d)\sum\limits_{i=1}^n\sum\limits_{j=1}^m[\gcd(i,j)=d]\) \(\sum\limits_{d}
T1 数独（WOJ4218）\(\color{green}{100}\) 小模拟/cy，45min左右就调完了。 code: #include<bits/stdc++.h> using namespace std; #define int long long #define in read() inline int read(){ int p=0,f=1; char c=getchar(); while(!isdigit(c)){if(c=='-')f=-1;c
计算一个式子：\(\sum\limits_{i = 1}^n \cfrac{n}{i}\)。 很明显可以直接一个\(for\)循环，\(O(n)\)求出结果，但是我们可以将其优化到\(O(\sqrt n)\)。 例题 AcWing199. 余数之和 给定正整数n和k，计算\((k \mod 1) + (k \mod 1) + ... + (k \mod n)\)的值。 \(1 \leq n, k \leq 10
\(\rm EXGCD\)，即 扩展欧几里得算法，简称 扩欧，是用来求出方程 \[ax+by=\gcd(a,b) \]的整数解的，其中 \(a,b\) 均为整数. 前置芝士：辗转相除法 与 裴蜀定理。 我们考虑辗转相除法的最后一步，当 \(b=0\) 时，要使得 \[ax+0y=\gcd(a,0)=a \]成立，那么只要取 \(x=1\)，\(y\) 取任意整数即可，不妨
G. GCD Festival 推式子 COMPFEST 13 - Finals Online Mirror 题意 给定数组\(a\)，求 \[\sum_{i=1}^n\sum_{j=1}^ngcd(a_i,a_j) \cdot gcd(i,j) \]\[2 \leq n \leq 10^5\\ 1 \leq a_i \leq 10^5 \]分析 \(n = \sum_{d|n} \varphi(d)\) \(gcd(i,j) = \sum_{d|gcd(i,j)
之前在 loj 上写题目的时候碰到的一些问题，在这里记录一下解决方法。 \(\lfloor\dfrac{n}{i^k}\rfloor\) 这种类型的时间复杂度是 \(O(n^{\frac{1}{k+1}})\) 的。 当 \(i< n^{\frac{1}{k+1}}\) 的时候，只有 \(O(n^{\frac{1}{k+1}})\) 种取值。 当 \(i > n^{\frac{1}{k+1}}\)
题目链接 题目链接 题意 求 \(i\perp j\to p_i\perp p_j\) 的排列个数，部分位置钦定。\(n\leq 10^6\) 题解 首先质因子集合相同的数可以连边，其次我们可能会将两个 \(\lfloor \dfrac{n}{p_i}\rfloor\) 的质数 \(p_i\) 进行替换。那么没有钦定任何 \(p_i\) 是，将 \(1\) 看作一个大质数
杜教筛 作用：用来求积性函数前缀和，时间复杂度为\(O(n^\frac{2}{3})\) 积性函数 若对于函数 \(f(x)\) 满足 \(f(x)=f(a)*f(b)\) ，其中 \(x=a*b\) 且 \(gcd(a,b)=1\)，那么 \(f(x)\) 为积性函数。 常见积性函数： \(\mu()\) ——莫比乌斯函数。 \(\phi()\) ——欧拉函数。\(\phi(n)=\su
Problem CF938C Constructing Tests 题目大意： 在一个 \(n \times n\) 的矩阵中填 \(0,1\)，使得每一个 \(m \times m\) 的子矩阵中都包含至少一个 \(0\)，要让 \(1\) 的个数最多。 现在知道最多有 \(x\) 个 \(1\)，问满足条件的 \(n,m\)，输出任意一个。如果不存在则输出 -1。 有多组数据，
数论分块 定义： 数论分块可以在 \(O \sqrt{n}\) 的时间里计算一些有除法下取整的和式。 主要是 将 \(\frac{n}{d}\) 相同的数一起同时计算。 定理： 定理 \(1\): \[a,b,c \in \mathbb{Z}, \lfloor \frac{a}{bc} \rfloor =\lfloor \frac{\lfloor\frac{a}{b}\rfloor}{c}\rfloor \]证明：
题目大意 小 S 在一棵被称为 OJT 的树上刷题。 这棵树上，有 \(n\) 个节点，每个节点上都有一道题目，每个节点上的题目难度可能会不同。 小 S 的能力有限，仅为k个单位能力，在第 \(i\) 个节点上的题目，难度为 \(D_i\)。 由于题目过毒，每做一道题都会杀死小S的脑细胞，使小 S 的能力值下降，做第
洛谷题面传送门 首先推式子： \[\begin{aligned} ans&=\sum\limits_{i=A}^B\sum\limits_{j=1}^i\{\dfrac{i}{j}\} \end{aligned} \]考虑差分，设 \[f(n)=\sum\limits_{i=1}^n\sum\limits_{j=1}^i\{\dfrac{i}{j}\} \]那么 \[ans=f(B)-f(A-1) \]考虑如何计算 \(f(n)\)： \[\be
询问 \(n\mod 1+n\mod 2+n\mod 3+\cdots +n\mod m\) 对 \(10^9+7\) 取模的值. \(1\leq n,m\leq 10^{13}\) 考虑 \(\mathrm O(m)\) 地计算肯定是不可行的 . 是否可以和分解因数一样，只考虑 \(\mathrm O(\sqrt n)\) 之内的内容呢 ？ 是可以的，观察可得对于 \(1\leq x\leq \sqrt n\) ，对
题目链接 题目链接 题意 重排数组，使得相邻两数之差绝对值的最小值尽可能大。\(n\leq 10^5\) 题解 排序，为了方便表述考虑连边。如果某条边的跨度不足 \(\lfloor \dfrac{n}{2}\rfloor\)，我们可以通过调整覆盖它的某条边（变为交叉）来使答案不减。因此我们要连边使得每条边跨度都有 \(\lf
预期 实际 \(\Delta\) A 100 100 0 B 100 100 0 C 100 0 -100 D 10 29 +29 E 50 100 +50 总计 360 329 -31 A Largest Rectangle in a Histogram 单调栈板子 B 没有上司的舞会 树形\(DP\)，\(f_{u,0}\)表示以\(u\)为根的子树在不选\(u\)时最大快乐指数，\(f_{u,1
题目 给定数组\(a\)，对于某个元素\(a_i\)和区间\([l,r]\)满足\(l\le i \le r\)。如果将\(a_l,a_{l+1},...a_r\)排序后，原\(a_i\)在排序后的新位置为\(j\)(值相同是可以任意排序)，那么\(a_i\)的奇异值为\(|j-\lfloor\frac{l+r+1}{2}\rfloor|\)。问每个元素最大的奇异值是多少。 题解
原理 \[a^n=\begin{matrix} \underbrace{ a*a*…*a } \\ n \end{matrix}\\ a^{13}=a^{(1101)_2}=a^8*a^4*a^1 \]应用 矩阵快速幂和多次置换 计算斐波那契数列可以构建\(2*2\)的转移矩阵从\(F_i,F_{i+1}\)到\(F_{i+1},F_{i+2}\)的变换，将转移矩阵用快速幂优化到\(n^3logk\) 把序列置
感觉这几个月做的数学都是把数学当工具的数学题，欧拉函数、莫比乌斯之类的好像上一次做还是在上一次 好几年前，于是在洛谷上找了几个专题来练一练。隔了感觉有两三年了，其实好多东西都忘差不多了。 （待更新） 目录零碎知识欧拉函数和欧拉定理卷积、莫比乌斯函数和莫比乌斯反演一些练习
题面 Binomial 题解 设 ord ( n ) \operatorname{ord}(n) ord(n) 表示
模板题 定义$\lfloor x\rfloor$表示小于等于$x$的最大整数，$\lceil x\rceil$表示大于等于$x$的最小整数 不难发现，若$a\in Z^{+}$且$b,x\in Z$，则$ax\le b\iff x\le \lfloor\frac{b}{a}\rfloor$、$ax\ge b\iff x\ge \lfloor\frac{b-1}{a}\rfloor+1$ 通过拉格朗日插值法，预处理出$S_{k}
Description 问长为 \(n\) ，每个数值域在 \([l,r]\) 的所有 \((r-l+1)^n\) 个序列的每个序列的拆分的方式之和。拆分合法当且仅当，将序列分为若干段，每一段都是一个等比数列。\(1\leq l\leq r\leq 10^7\)，\(n\leq 10^{18}\)。 Solution 容易发现问题可以拆解为求不同长度的等比数列个
数论分块 对于某些数列求和操作，求和的项数过大，直接求和会超时，但是数列中数的种数较少，且相同的数分布连续。此时需要找到一种方式能快速求出每种数的分布边界，然后枚举每一种数快速求和。 基本形式类似于 ∑
会 F 写不来 E 了属于是…… 题意 给定 \(a,b,c,d,n\) ，令 \(f(x)=a+\lfloor\frac{x}{b}\rfloor,f(x)=c+\lfloor\frac{x}{d}\rfloor\) ，求 \(\sum_{i=1}^n[f(x)=g(x)]\) 。 \(T\) 组询问。 \((1\le T\le2\times10^5,1\le a,b,c,d,n\le10^6)\) 题解 数论大讨论水平仍然不行！！！ 首先不考
