avatar
Articles
91
Tags
81
Categories
26

Home
Archive
Tag
Category
Repose
  • Gallery
  • Music
Link
About
Colythme
Search
Home
Archive
Tag
Category
Repose
  • Gallery
  • Music
Link
About

数论

CRT(中国剩余定理)及其扩展
Created2021-07-19|数论
CRT用于求解类似$$\left\{\begin{aligned}&x \equiv a_1 \pmod{m_1} \\&x \equiv a_2 \pmod{m_2} \\&… \\&x \equiv a_n \pmod{m_n}\end{aligned}\right.$$其中,$\forall i, j, ~ (m_i, m_j) = 1$ 设 $M = \prod_{i = 1}^n m_i, ~ M_i = \frac{M}{m_i}$,再设 $M_it_i \equiv 1 \pmod{m_i}$,其中 $1 \le i \le n$ 则可构造一解 $x = \sum\limits_{i = 1}^n a_iM_it_i$,任意解 $x_0 = x + k \times M$ ExCRT用于求解类似$$\left\{\begin{aligned}&x \equiv a_1 \pmod{m_1} \\&x \equiv a_2 \pmod{m_2} \\&… \\&x \equiv a_n \pmod{m_n}\e ...
洛谷3768 - 简单的数学题
Created2020-08-04|数论
根据Crash的数字表格,很容易可以将式子化简为$$\begin{aligned} Ans &= \sum\limits_{i = 1}^n \sum\limits_{j = 1} ij(i, j) \\ &= \sum\limits_{d = 1}^n d^3 \sum\limits_{k = 1}^{\left\lfloor\frac{n}{d}\right\rfloor} \mu(k) k^2 \left( \sum\limits_{i = 1}^{\left\lfloor\frac{n}{kd}\right\rfloor} i \right)^2 \end{aligned}$$感觉 $d, k$ 放在一起式子无法继续化简,主要是有 $kd$ 存在,故令 $T = kd$ ,则有$$Ans = \sum\limits_{T = 1}^n \left( \sum\limits_{i = 1}^{\left\lfloor\frac{n}{T}\right\rfloor} i \right)^2 T^2 \sum\limits_{d | T} d \mu(\frac{T ...
常系数线性递推
Created2020-08-04|数论
前置知识多项式除法|取模该部分见 多项式操作 特征值与特征向量(这部分我也不怎么理解,就大概记个结论) 对 $n$ 阶矩阵 $A$ 的特征向量 $\vec{v}$ 与,有特征值 $\lambda$ 满足$$A\vec v = \lambda\vec v$$移项则有 $(\lambda I - A)\vec v = 0$,其中 $I$ 为单位矩阵 那么有解的充要条件是 $det(\lambda I - A) = 0$ 令 $f(\lambda) = det(\lambda I - A)$,即此时可将 $det (\lambda I - A)$ 看作一个 $n$ 阶多项式,则称 $f(\lambda)$ 为矩阵 $A$ 的特征多项式,对矩阵 $A$,其任意特征值 $\lambda_0$ 都满足 $f(\lambda_0) = 0$ $\text{Hamiton-Cayley}$定理对矩阵 $A$ 及其特征多项式 $f(x)$,满足 $f(A) = O$,其中 $O$ 为零矩阵 证明留坑 常系数齐次线性递推求一个满足 $k$ 阶齐次线性递推数列 $h_i$ 的第 $n$ 项,即$$h_n = ...
多项式操作
Created2020-08-04|数论
多项式求逆该部分见 多项式求逆 多项式除法|取模有 $F(x) = G(x)Q(x) + R(x)$,给定 $F(x)$ 与 $G(x)$,求解 $Q(x), R(x)$,其中 $F(x), G(x)$ 最高次数分别为 $n, m$ 显然地,$Q(x), R(x)$ 的最高次数必定不超过 $n - m, m - 1$,考虑式子变形$$\begin{aligned}&F(x) = G(x)Q(x) + R(x) \\\Rightarrow &F(\frac1x) = G(\frac1x)Q(\frac1x) + R(\frac1x) \\\Rightarrow &x^nF(\frac1x) = x^mG(\frac1x)x^{n - m}Q(\frac1x) + x^{n - m + 1}x^{m - 1}R(\frac1x) \\\end{aligned}$$令 $A_R(x)$ 表示 $A(x)$ 系数翻转后得到的多项式,即 $A(x)$ 系数 $(a_0, a_1, …, a_n)$,$A_R(x)$ 系数 $(a_n, a_{n - 1}, …, a_1) ...
「清华集训2016」你的生命已如风中残烛(卡特兰数的另一种组合意义)
Created2020-08-04|数论
题意将题意转化一下 有 $m = \sum\limits_{i = 1}^n w_i$ 个数,前 $n$ 个中第 $i$ 个为 $w_i − 1$ ,剩下的 $m − n$ 个数都是 $− 1$ 求对这 $m$ 个数 $m!$ 种重排的方案中,有多少种满足重排后的序列任意前缀和不为负 题解卡特兰数的另一种组合意义首先有一个引理 给定 $n$ 个 $1$ 和 $n$ 个 $- 1$,它们组成一个长度为 $2n$ 的序列 $a$,则必定存在一个 $i$ 使得序列 $a_{i + 1…2n} + a_{1…i}$ 满足任意前缀都不小于零 证明也很简单,找到一个使得前缀和最小的位置 $i$,设 $s_{i,j}$ 表示 $a$ 上 $a_i$ 到 $a_j$ 的和,那么显然 $s_{i + 1, 2n} \ge 0$,对 $\forall k \in [1, i]$,$s_{i + 1, 2n} + s_{1, k} \ge 0$ 那么接下来求解 $n$ 个 $1$ 和 $n$ 个 $- 1$ 组成的合法序列(即任意前缀和不小于零)的方案数,也就是卡特兰数 $1$ 的编号为 $1 \sim ...
「国家集训队」Crash的数字表格
Created2020-08-04|数论
题意求 $\sum\limits_{i = 1}^N \sum\limits_{j = 1}^M lcm (i, j)$ Solution易知,原式$$\sum\limits_{i = 1}^N \sum\limits_{j = 1}^M \frac{ij}{\gcd (i, j)}$$枚举 $\gcd (i, j)$ ,且将 $d$ 提出来得$$\sum\limits_{d = 1}^{\min (N, M)} d \sum\limits_{i = 1}^{\left\lfloor\frac{N}{d}\right\rfloor} \sum\limits_{j = 1}^{\left\lfloor\frac{M}{d}\right\rfloor} ij[(i, j) = 1]$$将公式 $\sum\limits_{k | n} \mu(k) = [n = 1]$ 代入,得$$\sum\limits_{d = 1}^{\min (N, M)} d \sum\limits_{i = 1}^{\left\lfloor\frac{N}{d}\right\rfloor} \sum\limits ...
「四校联考」密码
Created2020-08-04|数论
题意给定 $n$ 行 $m$ 列的 $01$ 矩阵,求满足 $\big(\sum_{\forall i \in A, j \in B} value_{i, j}\big) \equiv 0 \pmod{2}$ 的二元组的数量 题解首先考虑 $n = 1$,令 $1$ 的个数为 $t$,可以发现答案$$\sum\limits_{k = 0}^{\lfloor\frac{m}{2}\rfloor} \dbinom{t}{2k} \times 2^{m - t} - 1 = 2^{m - 1} - 1$$至于证明,可以考虑已经选好的数列 $\{a_1, a_2, …, a_n\}$,那么现在考虑加入第 $n +1$ 个,显然加入 $0, 1$ 都是多了两种选择 那么现在考虑将行合并成一列 但是问题是如果需要枚举选的行还是需要 $2^n$ 的复杂度,然后我考场上就想到这就凉了 那么显然每次选好了行,再选列的时候这若干个列会贯穿所有行,那就可以将每行看作一个数,并且将它们合并(异或)起来,这样就变成做 $n = 1$ 的情况了 接下来只要考虑有几种情况会使得其异或和为 $1$ (或不为 $1$) ...
「北京省选集训2019」生成树计数「Matrix-Tree」
Created2020-08-04|数论
前置知识Matrix-Tree定理对于无向图 $G$,定义其度数矩阵 $D$$$\left[\begin{matrix}deg_1 & 0 & 0 & \cdots & 0 \\0 & deg_2 & 0 & \cdots & 0 \\0 & 0 & deg_3 & \cdots & 0 \\\vdots & \vdots & \vdots & \ddots & \vdots \\0 & 0 & 0 & 0 & deg_n\end{matrix}\right]$$定义其邻接矩阵 $C$$$\left[\begin{matrix}0 & a_{1,2} & a_{1,3} & \cdots & a_{1,n} \\a_{2,1} & 0 & a_{2,3} & \cdots & a_{2,n} \\\vdots & \vdots & \vdots &a ...
「NOI2016」循环之美
Created2020-08-04|数论
题意求解 $1 \le x \le n, ~ 1 \le y \le m$ 中满足在 $k$ 进制下 $\frac{x}{y}$ 是有限循环小数的个数 题解很好然后我一开始在计算器上敲了半个小时然后由于眼瞎并且归纳能力差并没有发现规律 设 $\frac{x}{y}$ 的循环节长度为 $l$,则有$$\begin{aligned}\frac{xk^l}{y} - \left\lfloor\frac{xk^l}{y}\right\rfloor &= \frac{x}{y} - \left\lfloor\frac{x}{y}\right\rfloor \\xk^l - \left\lfloor\frac{xk^l}{y}\right\rfloor y &= x - \left\lfloor\frac{x}{y}\right\rfloor y \\xk^l &\equiv x \pmod y \\k &\equiv 1 \pmod y\end{aligned}$$也就是说只要满足 $y \bot k$ 即可满足题意,故可将问题化为$$\begin{aligned ...
「FJOI2016」建筑师
Created2020-08-04|数论
这题一直往 $\text{dp}$ 去想,然后就没了 所以其实并不用注意之前的可看见的建筑的最大值是多少,若把那些能够被看见的建筑之间构成的区间看作一个小组,那么能够被看见的建筑不过是在当前小组中选择一个最大值最为代表排在左(右)边罢了 所以现在问题转化为将 $n - 1$ 个建筑分在 $A + B - 2$ 个盒子,然后盒子内可以随意排列的个数 然后上述问题实际上可以将每个盒子看作一个圆,故问题又转化为 $n - 1$ 个建筑构成的 $A + B - 2$ 个圆排列方案数,即斯特林数 $S (n - 1, A + B - 2)$ 然后要在其中选择一些放左(右)边,故答案即为$$ans = S (n - 1, A + B - 2) \times \dbinom{A + B - 2}{A - 1}$$ 代码12345678910111213141516171819202122232425262728293031323334353637383940414243444546474849505152535455565758596061626364#include <iostream> ...
12
avatar
Colythme
NO GAME NO LIFE
Articles
91
Tags
81
Categories
26
Follow Me
Recent Post
Analysis of Horror Game Design | Creators' Gacha Games - Part One2023-09-10
Analysis of Horror Game Design | Creators' Gacha Games - Part Two2023-09-10
Horror Game闲谈012023-09-05
Operating System2023-08-30
ABOUT ME2023-07-14
Categories
  • 3D3
  • AMP Learning3
  • Essays1
  • Game Design2
  • Machine Learning4
  • Operating System1
  • Social Network1
  • Web31
Tags
2-SAT 矩阵乘法 Game Essays 虚树 数学 网络流 二分图 多项式 点分治 GraphTheory OS Game Design 长链剖分 UE 扫描线 3D 归纳 动态点分治 prufer序列 杜教筛 组合数学 卡特兰数 Blockchain FWT 生成树 RL 博弈论 背包DP LCT 二分答案 ThreeJS 斯坦纳树 iClone 最短路 树上DP 线段树 后缀数组 ML 交互题 AtCoder Tarjan 线性基 树形DP FFT/NTT 斜率优化 概率DP 倍增 NLP 构造 数学期望 CRT/ExCRT 状压DP 容斥原理 集训 可持久化线段树 区间DP 边分治 矩阵树定理 整体二分 启发式合并 分块 杂谈 线性DP Web3D Python DigitalHuman AC自动机 OEIS 最小生成树 C++ 动态规划 Avatar 莫比乌斯反演 思维 最短路DP 贪心 ExGCD 生成函数 随机化 后缀自动机 分层DP
Archives
  • September 20233
  • August 20231
  • July 20239
  • February 20231
  • June 20221
  • May 20221
  • September 20211
  • July 20213
©2020 - 2023 By Colythme
Framework Hexo|Theme Butterfly
Search
Loading the Database