摘要: 求并是容易的,求交相对复杂一点。 点击查看代码 void solve() { SP x; INT(n); FOR(n) { U32(v); x.add_element(v | static_cast<u64>(v) << 32); } INT(m); FOR(m) { U32(v); x.add_e 阅读全文
posted @ 2024-05-08 18:06 皮皮的橙子树 阅读(19) 评论(0) 推荐(0)
摘要: 给 u32 求逆的,可能是牛迭 点击查看代码 constexpr u32 inv(u32 x){ u32 iv=2-x; for(int i=0;i<4;++i){ iv*=2-x*iv; } return iv; } 阅读全文
posted @ 2024-05-06 15:04 皮皮的橙子树 阅读(14) 评论(0) 推荐(0)
摘要: 引言 朴素的 Ln Exp 在模数非常小(小于长度)的情况下会比较拟人,于是我们需要扩域,于是经过我一晚上+一早上的研究现在至少得到了 \(O(n\log nF(\log^2n))\) 的保证正确的 Exp,\(F(x)\) 表示我们要保持 \(x\) 位精度所要付出的代价,这个是通过牛迭得到的精度 阅读全文
posted @ 2024-05-06 00:21 皮皮的橙子树 阅读(143) 评论(0) 推荐(0)
摘要: 25/4/2024 我刚刚发现一个神奇的事情,矩阵树定理的矩阵消元的时候,矩阵的两行是不需要交换的,只要交换就意味着答案=0。 目前不确定正确性,我尝试在lg矩阵树定理模版、lg BEST定理模板、loj3304作业题提交这种逻辑的做法,都通过了。我本地也对这个问题做了大量对拍,对于有向和无向的情况 阅读全文
posted @ 2024-04-25 20:47 皮皮的橙子树 阅读(61) 评论(0) 推荐(0)
摘要: 矩阵树构建(记得删去一行一列,有向树删根): 无向:a[x][y]--,a[y][x]--,a[x][x]++,a[y][y]++; 内向:a[x][y]--,a[x][x]++;(或者a[y][x]--,a[x][x]++;) 外向:a[x][y]--,a[y][y]++;(或者a[y][x]-- 阅读全文
posted @ 2024-04-25 16:38 皮皮的橙子树 阅读(18) 评论(0) 推荐(0)
摘要: 这两个问题都是卡满上界的,非常深刻。 https://codeforces.com/contest/1365/problem/G https://codeforces.com/contest/1705/problem/F 阅读全文
posted @ 2024-04-23 21:32 皮皮的橙子树 阅读(12) 评论(0) 推荐(0)
摘要: 从 https://codeforces.com/gym/101380/problem/C 谈起 关于本题: 1、题意:构造一个 \(n\times n\) 的 01 矩阵,让他的每一个顶着上顶边的子方阵的行列式模二意义下不为 0。 2、做法: 我们手搓第一行和第二行是唯一的,然后我们可以想到构造主 阅读全文
posted @ 2024-04-23 17:35 皮皮的橙子树 阅读(22) 评论(0) 推荐(0)
摘要: 常见的代数结构有:poly(形式幂级数),soly(集合幂级数,多元多项式的特例,不过这个名字是我自己起的),矩阵。 乘法:将其余维抽象,视作参数即可,多项式可以使用bluestein做到常数维度无关变换。有趣的是,矩阵同样可以作为参数。 poly 白兔之舞 matrix 整式递推 求逆: 当出现带 阅读全文
posted @ 2024-04-23 15:19 皮皮的橙子树 阅读(24) 评论(0) 推荐(0)
摘要: 1)最简单的,就是这个式子 \(\sum_\limits{b=0}^{n} \binom{a+b}{a} = \binom{n+a+1}{a+1}\) GF 证明: \[\begin{aligned} &Ans\\ =&\sum_{b=0}^{n}[x^b]\frac{1}{(1-x)^{a+1}} 阅读全文
posted @ 2024-04-15 00:21 皮皮的橙子树 阅读(28) 评论(0) 推荐(0)
摘要: 0、随机实变量的性质 连续性,对称性。 利用对称性,我们可以化简很多形式。 对于 \(x\in [0,a],x\ge y\) 可以化简成 \(x+y\le a\) 利用对称性,我们可以构造很多巧妙的双射。 1、连续问题离散化技巧 [0,n) 的随机实变量可以用一个 [0,n) 的整数+一个 这种做法 阅读全文
posted @ 2024-04-04 23:12 皮皮的橙子树 阅读(23) 评论(0) 推荐(0)