ixRic-CSDN博客

原创多项式拉格朗日反演与复合逆

文章目录膜膜膜引入拉格朗日反演扩展拉格朗日反演例题膜膜膜Tiw 永远的神！引入拉格朗日反演及其扩展用于求一类 F(G(x))=H(x)F(G(x)) = H(x)F(G(x))=H(x) 的问题，其中 F(x),H(x)F(x), H(x)F(x),H(x) 为已知多项式函数，G(x)G(x)G(x) 为待求函数。但其只能在 O(nlog⁡n)O(n \log n)O(nlogn) 的时间内求得 G(x)G(x)G(x) 的某一项系数。拉格朗日反演拉格朗日反演用于在 O(nlog⁡n)O(n \

2020-07-11 09:46:30 1237

原创 [集训队作业 2013] 城市规划（指数型生成函数 + 多项式对数） | 错题本

文章目录题目分析错因代码题目[集训队作业 2013] 城市规划分析先求无标号的，最后乘一个 n!n!n! 即可。设 G(x)=∑i=0ngi⋅xiG(x) = \sum_{i = 0}^{n} g_i \cdot x^iG(x)=∑i=0ngi⋅xi 是无标号无向图的个数的生成函数，其中 gi=2Cn2i!g_i = \dfrac{2^{C_n^2}}{i!}gi=i!2Cn2，即 iii 个点构成的无标号无向图个数（考虑每条边选不选，标号去重）。设 F(x)=∑i=0nfi⋅xiF(x

2020-07-09 15:57:46 245

原创 [CodeForces 438E] The Child and Binary Tree（生成函数 + 多项式求逆 + 多项式开根） | 错题本

文章目录题目分析代码题目[CodeForces 438E] The Child and Binary Tree分析首先列出 DP 式，设 dp[i]dp[i]dp[i] 表示用权值和为 iii 的点，形成的合法二叉树数量，于是有 dp[i]={∑k=0i∑j=0i−k(dp[i−k]×dp[i−k−j]×[k∈C])i>01i=0dp[i] = \begin{cases} \sum_{k = 0}^{i} \sum_{j = 0}^{i - k} (dp[i - k] \times dp[i

2020-07-09 15:05:32 180

原创 [CodeForces 757F] Team Rocket Rises Again（最短路径图 + 支配树） | 错题本

文章目录题目分析错因代码题目[CodeForces 757F] Team Rocket Rises Again分析建出最短路径图（是个 DAG）的支配树，找到除起点外支配的点最多的点即可。错因Tid 打成 Dfn；用链式前向星，To 和 Next 数组的大小没有开到 mmm；以为是求要删的点（样例还能过）。代码#include <bits/stdc++.h>typedef long long LL;typedef std::pair<int, int>

2020-07-07 15:36:37 179

原创 [ZJOI 2012] 灾难（支配树） | 错题本

文章目录题目分析错因代码题目[ZJOI 2012] 灾难分析将图中所有边取反，发现要求的就是每个点支配的点数。然而发现一个问题：这个不是一个单源点的有向图。我们只需建一个源点，将其向所有入度为零的点连边，建出支配树即可。当然，这个图是 DAG，事实上拓扑排序就 + LCA 能找到支配点。但我还是用的 LT错因没想到怎么处理多源点的 DAG 的问题；板子打错了两个地方。代码#include <bits/stdc++.h>const int MAXN = 65534;

2020-07-07 09:58:34 165

原创支配树 | 模板

文章目录题目实现思路代码题目【模板】支配树实现思路DFS 一次得到 Dfn 和 Tid 数组；从 Dfn 大的点开始枚举点 uuu：枚举 uuu 的反向邻接点 vvv，则 vvv 所在的并查集的最小 SDom 可以更新 uuu（因为 uuu 的 DFS 树上的父亲还没更新进带权并查集，所以这时候所有的 vvv 都可以在符合定理的条件下更新 uuu），注意不要管不在 DFS 树上的点；在 SDom 树上将 SDom[u] 向 uuu 连边；将 uuu 更新进带权并查集（将自己在并查集上连向

2020-07-06 18:05:09 223

原创 [CodeForces 1109D] Sasha and Interesting Fact from Graph Theory（广义 Cayley 定理 + 组合数学） | 错题本

文章目录题目分析代码题目[CodeForces 1109D] Sasha and Interesting Fact from Graph Theory分析a,ba, ba,b 不影响答案，不妨设 a=1,b=2a = 1, b = 2a=1,b=2，枚举 a,ba, ba,b 之间的边数 i∈[1,n−1]i \in [1, n -1]i∈[1,n−1]，隔板法可得链 a−ba - ba−b 的方案数为 An−2i−1⋅Cm−1i−1A_{n - 2}^{i - 1} \cdot C_{m - 1}

2020-07-06 12:25:13 380

原创 [CodeForces 156D] Clues（Prufer 序列 + 多项式定理） | 错题本

文章目录题目分析错因代码题目[CodeForces 156D] Clues分析记联通块个数为 kkk，枚举每个联通块的度数 did_idi（∑i=1kdi=2(k−1)\sum_{i = 1}^{k} d_i = 2(k - 1)∑i=1kdi=2(k−1)，每条边的贡献为 222，因此度数和为 2(k−1)2(k - 1)2(k−1)），用重排列计算方案数，并且对于一个联通块，取里面一个点的方案数为联通块的大小 sis_isi，可得答案为 ∑∑i=1kdi=2(k−1), di≥

2020-07-06 09:28:34 253

原创 [HNOI2008] 明明的烦恼（Prufer 序列 + 高精度） | 错题本

文章目录题目分析错因代码题目[HNOI2008] 明明的烦恼分析令 ccc 为规定了度数的点的个数， s=∑i=1c(di−1)s = \sum_{i = 1}^{c} (d_i - 1)s=∑i=1c(di−1)，先考虑规定度数的点的方案数，就是一个重排列乘上选位置的方案数：s!∏i=1c(di−1)!⋅Cn−2s\frac{s!}{\prod_{i = 1}^{c} (d_i - 1)!} \cdot C_{n - 2}^{s}∏i=1c(di−1)!s!⋅Cn−2s 再考虑没规定度

2020-07-05 23:10:19 146

原创 [Code+ #2] 白金元首与独舞（建图 + 矩阵树定理推广） | 错题本

文章目录题目分析代码题目[Code+ #2] 白金元首与独舞分析首先不难想到直接用 nmnmnm 个点上下左右连边建图，既然要走出去，就把“外面”建成一个点 rrr，所以问题变为求以 rrr 为根的内向树个数。但这样时间复杂度是 O((nm)3)O((nm)^3)O((nm)3)。观察到 k≤300k \leq 300k≤300，因此考虑只保留 .，把 . 向它上下左右能走到的另一个 . 连边即可。记忆化搜索实现，注意每个点都要搜一遍找环。如果是求有向图以 uuu 为根的个数，我们只需要算 ∣L

2020-07-05 20:34:32 209

原创 [SHOI 2016] 黑暗前的幻想乡（矩阵树定理 + 容斥） | 错题本

文章目录题目分析错因代码题目[SHOI2016] 黑暗前的幻想乡分析把问题转化为“每个公司都要参与”的方案数。假设现在只有 3 个公司，我们先算总方案数，减掉 1 公司不参与的方案数，减掉 2 公司不参与的方案数，减掉 3 公司不参与的方案数，这时候还不是答案，因为我们把 1、2 公司都不参与的方案数，2、3 公司都不参与的方案数和 3、4 公司都不参与的方案数多减了一次（例如：减掉 1 公司不参与的方案数、减掉 2 公司不参与的方案数后，1、2 公司都不参与的方案数被减了两次），所以要把它们加回

2020-07-05 19:20:25 209

原创 [SDOI 2014] 重建（矩阵树定理推广） | 错题本

文章目录题目分析代码题目[SDOI2014] 重建分析问题要求 ∑T(∏e∈Tpe∏e∉T(1−pe))\sum_{T} \left( \prod_{e \in T} p_e \prod_{e \notin T} (1 - p_e)\right)T∑⎝⎛e∈T∏pee∈/T∏(1−pe)⎠⎞ 由于矩阵树只能处理 e∈Te \in Te∈T，所以把式子变成 ∑T(∏e∈Tpe∏e(1−pe)∏e∈T(1−pe))=∏e(1−pe)∑T(∏e∈Tpe1−pe)\begin{aligned

2020-07-04 17:56:40 193

原创 [NOI Online #2 入门组] 魔法（矩阵加速） | 错题本

文章目录题目分析错因代码题目[NOI Online #2 入门组] 魔法分析dp[i][j][k]dp[i][j][k]dp[i][j][k] 表示 iii 到 jjj 用 kkk 次魔法的最小代价，f[i][j]f[i][j]f[i][j] 表示 iii 到 jjj 用一次魔法的最小代价：dp[i][j][k]=min⁡{dp[i][u][k−1]+f[u][j],f[i][u]+dp[u][j][k−1]}dp[i][j][k] = \min \{dp[i][u][k - 1] + f[u][

2020-07-04 16:30:46 156

原创 [NOI Online #2 入门组] 建设城市（组合数学） | 错题本

文章目录题目分析错因代码题目[NOI Online #2 入门组] 建设城市分析正整数 x1,x2,⋯ ,xnx_1, x_2, \cdots, x_nx1,x2,⋯,xn 满足 1≤x1≤x2≤⋯≤xn≤r1 \leq x_1 \leq x_2 \leq \cdots \leq x_n \leq r1≤x1≤x2≤⋯≤xn≤r，则 x1,x2,⋯ ,xnx_1, x_2, \cdots, x_nx1,x2,⋯,xn 的取值方案数为 Cn+r−1r−1C_{n + r - 1}^{

2020-07-04 16:17:34 317

原创 [NOI Online #2 入门组] 荆轲刺秦王（记忆化搜索） | 错题本

文章目录题目分析错因代码题目[NOI Online #2 入门组] 荆轲刺秦王分析记忆化 BFS。记录 T[x][y][c1][c2] (c1≤C1,c2≤C2)T[x][y][c_1][c_2]\ (c_1 \leq C_1, c_2 \leq C_2)T[x][y][c1][c2] (c1≤C1,c2≤C2) 表示能否恰好用 c1c_1c1 次隐身和 c2c_2c2 次瞬移从起始位置走到 xxx 行 yyy 列即可，Num[x][y]Num[x][y]Num[

2020-07-04 16:01:58 279

原创行列式入门与矩阵树定理完整证明

文章目录前置技能行列式定义性质拉普拉斯展开线性性可乘性可加性不重性可倍加性转置不变性可交换性行可交换性列可交换性优化行列式的计算矩阵树定理前置定义一些引理转置引理连通性引理引理 1引理 2引理 3Binet - Cauchy 定理前置技能简单的数学知识高斯消元行列式定义对于一个 n×nn \times nn×n 的矩阵 AAA，记它第 iii 行第 jjj 列的元素为 ai,ja_{i, j}ai,j，以及一个 1∼n1 \sim n1∼n 的排列 ppp，记 λA(p)=(−1)τ(p)

2020-07-03 23:50:38 978 1

原创 C++ 拉格朗日插值法优化 DP

文章目录拉格朗日插值法简介拉格朗日插值法模板DP 优化思路例题一分析代码例题二分析代码A trick拉格朗日插值法简介众所周知，nnn 个点 (xi,yi)(x_i, y_i)(xi,yi)（任意两个点横坐标不相等）可以确定一个 n−1n - 1n−1 次多项式函数 y=f(x)y = f(x)y=f(x)。拉格朗日插值法可以根据这 nnn 个点求出这个多项式 f(x)f(x)f(x)。当然，实际应用中通常求出横坐标为 kkk 的点在该（ nnn 个点确定的）多项式函数上对应的纵坐标的值，代码实现

2020-06-11 22:26:22 1174

原创二进制状态压缩枚举子集

对于二进制状态 SSS，可以用此方法不重不漏地枚举出子状态：for (int sub = S; sub; sub = (sub - 1) & S) { // sub 为 S 的子集}【证明】由 sub=(sub−1)∩Ssub = (sub - 1) \cap Ssub=(sub−1)∩S 可知 subsubsub 每次必然会变小，于是我们只需要证明区间 ((sub−1)∩S,sub)\left((sub - 1) \cap S, sub\right)((sub−1)∩S,sub) 中不

2020-05-30 10:15:40 1784 2

原创洛谷P1117 [NOI2016]优秀的拆分（巧妙的计数方法）

文章目录题目分析代码题目[NOI2016]优秀的拆分分析统计以SiS_iSi开头的形如AA\text{AA}AA的子串的数量，存入L[i]L[i]L[i]；统计以SiS_iSi结尾的形如AA\text{AA}AA的子串，存入R[i]R[i]R[i]。于是把可以把它们拼起来，答案就是∑i=2n(L[i]×R[i−1])\sum \limits_{i = 2}^{n} (L[i] \tim...

2020-05-03 19:57:51 437

原创 NOI.AC170 数数（计数DP）

文章目录题目分析代码题目题目描述求有多少对1∼n1∼n1∼n的排列(a,b)(a, b)(a,b)满足m≤∑i=1nmax⁡{ai,bi}m \leq \sum\limits_{i = 1}^{n} \max\{a_i,b_i\}m≤i=1∑nmax{ai,bi}。两个方案(a,b)(a, b)(a,b)和(a′,b′)(a', b')(a′,b′)不同当且仅当存在iii使得ai≠a...

2020-03-28 15:31:38 233

原创 LOJ530 「LibreOJ β Round #5」最小倍数（二分）

文章目录题目分析代码题目「LibreOJ β Round #5」最小倍数分析令n!=p1a1p2a2⋯pkakn! = {p_1}^{a_1}{p_2}^{a_2} \cdots {p_k}^{a_k}n!=p1a1p2a2⋯pkak，那么ai=∑j=1∞⌊npij⌋a_i = \sum\limits_{j = 1}^{\infty} \left\lfloor \dfrac{n...

2020-03-28 11:32:25 284

原创 LOJ2161 「POI2011 R2 Day1」差值 Difference（细节DP）

文章目录题目分析代码题目「POI2011 R2 Day1」差值 Difference分析考虑枚举两个字符分别作为子序列的出现次数最多和最少的字符。一个性质是，这两个字符到底是不是次数最多或最少的字符并不重要，我们只需要统计最大差值，就能自动避免不符合要求的情况（因为不符合要求的情况一定比大差值小）。简单来说，枚举字符aaa，bbb和区间[l,r][l, r][l,r]，计算[l,r][l...

2020-03-28 11:17:54 341

原创 C++类欧几里得算法

文章目录概述算法版题（本文用绿色\color{green}绿色绿色标注一切函数，用紫色\color{purple}紫色紫色标注部分可以O(1)O(1)O(1)计算的常量式，便于观察）（不知道为什么加了颜色过后括号不能自动变大了，将就看吧= =）概述类欧几里得算法用于求形如f(a,b,c,n)=∑i=0n⌊ai+bc⌋\color{green}{f(a, b, c, n)} \color{o}...

2020-02-20 12:25:48 559

原创 [CTSC2006] 歌唱王国（概率生成函数 + KMP / 哈希）

文章目录题面分析代码题面洛谷P4548 [CTSC2006]歌唱王国分析先了解一下概率生成函数，对于一个随机变量aaa，它的概率生成函数是F(x)=∑i≥0(Pr(a=i)⋅xi)F(x)=\sum\limits_{i\geq0}(\text{Pr}(a=i)\cdot x^i)F(x)=i≥0∑(Pr(a=i)⋅xi)对它求导得F′(x)=∑i≥0(Pr(a=i)⋅i⋅xi−1)F'(...

2020-02-09 17:37:25 497 1

原创 [TopCoder 12984] TorusSailing（高斯消元主元法优化）

文章目录题面分析代码题面Vjudge TorusSailing分析首先得到一个DP方程dp[i][j]=12(dp[(i+1) mod N][j]+dp[i][(j+1) mod M])+1①dp[i][j]=\dfrac{1}{2}\left(dp[(i + 1)\ \text{mod}\ N][j]+dp[i][(j + 1)\ \text{mo...

2020-02-09 14:37:17 545

原创 [SDOI2012]走迷宫（Tarjan + 概率DP + 高斯消元）

文章目录题面分析代码题面洛谷P6030 [SDOI2012]走迷宫分析设dp[u]dp[u]dp[u]表示从uuu走到TTT的期望步数，VuV_uVu表示uuu能到的点的集合，那么有dp[u]=∑v∈Vu1∣Vu∣⋅(dp[v]+1)=1∣Vu∣∑v∈vudp[v]+1\begin{aligned}dp[u]&=\sum\limits_{v\in V_u}\dfrac{1}{|...

2020-02-09 12:58:31 246

原创网络流24题

文章目录1. 餐巾计划问题题面分析代码2.[CTSC1999]家园题面分析代码3.航空路线问题题面1. 餐巾计划问题【洛谷1251】餐巾计划问题题面题目描述一个餐厅在相继的 NNN 天里,每天需用的餐巾数不尽相同。假设第 iii 天需要 rir_iri 块餐巾（i=1,2,...,Ni=1,2,...,Ni=1,2,...,N）。餐厅可以购买新的餐巾，每块餐巾的费用为 ppp 分；或者...

2020-02-07 12:32:32 516

原创斜率优化之凸包优化与李超线段树

文章目录前言例一前言本题的例题可能比较偏，但难度都是正常的（经典例题已经用传统斜率优化AC过了，不想再写）。例一Kalila and Dimna in the Logging Industry不用看题，直接看转移方程即可：dp[i]=min⁡1≤j<i{dp[j]+bj⋅ai}dp[i] = \min\limits_{1 \leq j < i}\{dp[j] + b_j \...

2019-12-28 09:59:41 1275 1

原创 [CodeForces 697F][CodeForces 696D] Legen...（AC自动机+矩阵加速）

文章目录题目分析代码题目time limit per test: 6 secondsmemory limit per test: 256 megabytesinput: standard inputoutput: standard outputDescriptionBarney was hanging out with Nora for a while and now he thin...

2019-11-14 20:07:58 382

原创 C++高斯消元详解

文章目录前言直接来代码前言学了才发现这玩意贼简单虽然据说是复习。直接来高斯消元用于解一个线性方程组，就是：{a1,1⋅x1+a1,2⋅x2+⋯+a1,n⋅xn=a1,n+1a2,1⋅x1+a2,2⋅x2+⋯+a2,n⋅xn=a2,n+1⋯an,1⋅x1+an,2⋅x2+⋯+an,n⋅xn=an,n+1\begin{cases}a_{1,1}\cdot x_1+a_{1,2}\cdot ...

2019-11-13 08:51:45 3955 1

原创 [洛谷P2540]【NOIP2015】斗地主增强版（DP+搜索）

文章目录题目分析代码题目P2540 斗地主增强版分析如果不出顺子，那么怎么出最优是可以DP解决的：dp[i][j][k][l]dp[i][j][k][l]dp[i][j][k][l]表示一副牌有iii个炸弹、jjj个三张、kkk个顺子、lll个单牌的最优出法。在这个状态定义下，原来的4i+3j+2k+l4i+3j+2k+l4i+3j+2k+l张牌变成了i+j+k+li+j+k+li+j+k...

2019-11-11 20:25:12 466

原创 [HNOI2019] 校园旅行（搜索+二分图性质）

文章目录题目分析代码题目注意数据范围是N≤5000N\leq 5000N≤5000！！分析考虑最plain的DP，dp[u][v]=1/0dp[u][v]=1/0dp[u][v]=1/0表示uuu到vvv是否有一条回文路径，如果dp[u][v]=1dp[u][v]=1dp[u][v]=1，那么dp[x][y]=1dp[x][y]=1dp[x][y]=1（xxx为uuu的邻结点，yyy为vv...

2019-11-11 12:21:01 215

原创【牛客CSP-S提高组赛前集训营5】B 十二桥问题（最短路+状压DP）

文章目录题目分析代码题目十二桥问题分析k≤12k\leq12k≤12，果断考虑状压。将所有必须走的边的端点u,vu,vu,v视为关键点，为每个关键点跑Dijkstra，然后状压DP。dp[S][i][0]dp[S][i][0]dp[S][i][0]表示走完集合为SSS的边，最后停在第iii条边的uuu点的最小代价；dp[S][i][1]dp[S][i][1]dp[S][i][1]表示...

2019-11-08 22:30:46 255

原创 C++自动取模的模数类

一直想写，今天终于写出来了。使用方法：ModNumber<MOD> k，可以定义一个模MOD意义下的变量k，它接下来的所有运算都自动模MOD，不需要再手写。具体操作方法请自行探索：#define LL long longtemplate<const int _MOD> struct ModNumber{ int x; ModNumber(){x=0...

2019-11-06 19:02:09 1070

原创浅谈WQS二分、带权二分、凸优化与一类斜率优化DP

幕天题意一个简单的转化二维DP朴素斜率优化带权二分我们先看一道简单题题意分析代码进入正题题目回顾分析加权权与段数的关系单调性两个问题DP朴素DP斜率优化就题论题的分析方法一个普适的结论代码后记除最后一张普适公式外，文中所有图像和公式均为原创，转载请附出处谢谢！！！题意[BZOJ 4518] Sdoi2016 征途求nnn个数a1,a2,⋯ ,ana_1,a_2,\cdots,a_na1,...

2019-11-05 22:53:27 1020 2

原创 [BZOJ 3551][ONTAK2010]Peaks加强版（Kruskal重构树+主席树）

文章目录题目分析代码题目Description在Bytemountains有NNN座山峰，每座山峰有他的高度hih_ihi。有些山峰之间有双向道路相连，共MMM条路径，每条路径有一个困难值，这个值越大表示越难走，现在有QQQ组询问，每组询问询问从点vvv开始只经过困难值小于等于xxx的路径所能到达的山峰中第kkk高的山峰，如果无解输出-1。Input第一行三个数NNN，MMM，QQQ。...

2019-11-02 12:04:58 210

原创 Kruskal重构树入门

Kruskal算法众所周知，Kruskal是一种优秀的基于并查集的最小（最大）生成树算法。它的简要流程为：- 将所有边按边权排序；- 依次枚举所有边的两端点u和v，若它们不在一个联通块内，则将该边加入生成树的边集中，且将u、v合并于一个连通块中。Kruskal重构树定义Kruskal重...

2019-11-02 11:24:59 705 1

原创 [CodeForces 141E]Clearing Up（Kruskal变式）

文章目录题目题目大意分析举个例子原图第一步第二步第三步代码题目Clearing Up题目大意给出一个nnn（1≤n≤1031\leq n\leq 10^31≤n≤103）个点mmm（1≤m≤1051\leq m\leq 10^51≤m≤105）条边的图，每条边有个颜色，要么是S，要么是M，要求一个该图的生成树，使得两种颜色的边各占一半。分析nnn是偶数时无解。举个例子原图原图（红...

2019-10-31 12:11:01 239

原创 [CodeForces 723F]st-Spanning Tree（分类讨论+贪心）

文章目录题目题目大意分析代码题目st-Spanning Tree题目大意给出一个nnn（2≤n≤200 0002\leq n\leq 200\ 0002≤n≤200 000）个点mmm（1≤m≤400 0001\leq m\leq 400\ 0001≤m≤400 000）条边的图（连通，无重边，无自环）和s,t,ds,dts,t,d_s,d_ts,t...

2019-10-31 11:36:56 259

原创 [GYM 101630C]CONNECTIONS（DFS树构造）

文章目录题目题目大意分析代码题目A-L题全题目大意给出一个有向强连通图nnn个点，mmm条边，要求留下2n2n2n(删除m−2nm-2nm−2n)条边使其仍然强连通。保证一定存在解。（2n<m≤1052n<m\leq 10^52n<m≤105）分析边数和nnn有关的那就是树了。所以考虑从某点uuu开始做一个dfs生成树，由于原图强联通，所以一定有这样一个生成树，且...

2019-10-29 17:49:00 213