数学
文章平均质量分 95
Wangduoyu
这个作者很懒,什么都没留下…
展开
-
容斥原理与Mobius函数
容斥原理 百度百科:容斥原理 设S1,S2,⋯ ,SnS_1, S_2, \cdots ,S_nS1,S2,⋯,Sn为有限集合,∣S∣|S|∣S∣表示集合SSS的大小,则: ∣⋃i=1nSi∣=∑i=1n∣Si∣−∑1≤i<j≤n∣Si∩Sj∣+∑i≤i<j<k≤k∣Si∩Sj∩Sk∣+⋯+(−1)n+1∣S1∩⋯∩Sn∣\left | \bigcup_{i=1}^{n} S_i \right | = \sum_{i=1}^{n}\left | S_i \right | - \sum原创 2021-11-01 18:05:40 · 348 阅读 · 0 评论 -
古代猪文--(全是数论QwQ)
题目链接 题意: Sol: 首先当qqq为999911659999911659999911659时,答案为000. 否则的话,因为999911659999911659999911659是质数,根据欧拉定理的推论有: q∑d∣nCnd≡q∑d∣nCnd mod 999911658( mod 999911659)q^{\sum_{d|n} C_n^d} \equiv q^{\sum_{d|n} C_n^d \bmod 999911658} ( \bmod \ 999911659)q∑d∣nC原创 2021-10-31 19:27:58 · 168 阅读 · 0 评论 -
高斯消元解线性方程组
高斯消元:O(n3)O(n^3)O(n3) 写出增广矩阵,将增广矩阵通过初等行变换变成阶梯型,然后往回带 适用于求解包含nnn个方程,nnn个未知数的多元线性方程组 例如下式多元方程组: {a1,1∗x1+a1,2∗xi+⋯+a1,n∗xn=b1a2,1∗x1+a2,2∗xi+⋯+a2,n∗xn=b2⋮an,1∗x1+an,2∗xi+⋯+an,n∗xn=bn\begin{cases} a_{1,1}*x_1 + a_{1,2}*x_i +\quad \cdots \quad+ a_{1, n}*x_n =原创 2021-10-30 16:11:44 · 139 阅读 · 0 评论 -
矩阵乘法加速递推
矩阵快速幂加速递推式原创 2021-10-30 13:41:24 · 468 阅读 · 0 评论 -
中国剩余定理以及扩展中国剩余定理
中国剩余定理 设m1,m2,⋯ ,mnm1, m2,\cdots, m_nm1,m2,⋯,mn是两两互质的整数,对于下式的一元线性同余方程组有整数解 {x≡a1( mod m1)x≡a2( mod m2)⋮x≡an( mod mn)\begin{cases} x \equiv a_1(\bmod m_1) \\ x \equiv a_2(\bmod m_2) \\ \qquad \vdots \\ x \equiv a_n (\bmod m_n) \end{cases}⎩⎪⎪⎪⎪⎨⎪⎪⎪⎪⎧x≡a1原创 2021-10-29 20:38:21 · 258 阅读 · 0 评论 -
扩展欧几里得求逆元以及逆元的应用
ax≡1( mod m)ax \equiv 1( \bmod m)ax≡1(modm) 若aaa与mmm互质,那么xxx则为aaa的逆元 可推得式子: ax+my=1ax + my = 1ax+my=1 利用扩展欧几里得求: #include <bits/stdc++.h> using namespace std; int exgcd(int a, int b, int &x, int &y) { if(!b) { x = 1; y =原创 2021-10-29 14:03:15 · 609 阅读 · 0 评论 -
数学知识——数论分块
对于一个式子: ∑i=1n⌊ni⌋\sum_{i=1}^{n} \left \lfloor \frac{n}{i} \right \rfloori=1∑n⌊in⌋ 我们容易发现在某段区间[l,r][l, r][l,r], ⌊ni⌋\left \lfloor \frac{n}{i} \right \rfloor⌊in⌋的值是一样的。那么我们如何获取这个区间呢。 我们假设区间左端点为lll,同时令k=⌊nl⌋k = \left \lfloor \frac{n}{l} \right \rfloork=⌊l原创 2021-10-28 01:10:21 · 529 阅读 · 1 评论