dashu497731727-CSDN博客

转载 Luogu 4602 [CTSC2018]混合果汁

BZOJ 5343福利题。对于每一个询问可以二分$d$，然后把满足条件的果汁按照$p$从小到大排序贪心地取$L$升看看满不满足价格的条件。那么按照$p$建立权值主席树，$chk$的时候在主席树上走一走算出价格即可。当然也可以整体二分。时间复杂度都是$O(nlog^2n)$。Code：#include <cstdio>#include &...

2019-02-24 09:20:00 121

转载 Luogu 4433 [COCI2009-2010#1] ALADIN

LOJ 2958线段树 +类欧主要矛盾就是$\sum\limits_{i = l}^{r}Ai \mod B$怎么做。拆成前缀和相减的形式，有$$\sum_{i = 0}^{r}(Ai \mod B) - \sum_{i = 0}^{l - 1}(Ai \mod B)$$我们知道$$a \mod b = a - \left \lfloor \frac{a}{b}...

2019-02-16 16:56:00 146

转载 Luogu 5170 【模板】类欧几里得算法

原理不难但是写起来非常复杂的东西。我觉得讲得非常好懂的博客。传送门我们设$$f(a, b, c, n) = \sum_{i = 0}^{n}\left \lfloor \frac{ai + b}{c} \right \rfloor$$$$g(a, b, c, n) = \sum_{i = 0}^{n}i\left \lfloor \frac{ai + b}{c} \...

2019-02-12 14:38:00 151

转载 COGS 2189 帕秋莉的超级多项式

放模板啦！以后打比赛的时候直接复制过来。说句实话vector的效率真的不怎么样，但是似乎也还行，最主要是……写得比较爽。#include <cstdio>#include <cstring>#include <algorithm>#include <vector>#include <cmath&gt...

2019-02-09 22:13:00 111

转载 Luogu 4512 【模板】多项式除法

高级操作，感觉非常神仙。题目中的字母太难懂了，重新定义一下。$$A(x) = B(x) * C(x) + D(x)$$其中，$A(x)$的次数是$n$，$B(x)$的次数是$m$，$A, B$都已知，要求$C$的次数是$n - m$，$D$的次数小于$m$。定义一种操作:如果$A$的次数为$n$，那么$$A_R(x) = x^nA(\frac{1}{x})$$...

2019-01-26 11:18:00 107

转载 PKUWC2019 真·游记

由于我这个大傻逼考试结果实在是没法看，所以这篇游记将尽可能略去和考试有关的内容，变成一篇真正的游记……接下来的内容中将会出现各种旅游攻略，寻求干货的同学可以提前左转了。Day -7学考终于结束了，可以停课了。感觉停课和闹着玩一样，短短一个星期，又能做什么事情呢？NOIP之后就一直颓废，原来觉得再也不可能有这种级别的机会了，后来发现竟然在这边卡线了，作为一个什么东西都不...

2019-01-22 20:35:00 118

转载 Luogu 4491 [HAOI2018]染色

BZOJ 5306考虑计算恰好出现$s$次的颜色有$k$种的方案数。首先可以设$lim = min(m, \left \lfloor \frac{n}{s} \right \rfloor)$，我们在计算的时候只要算到这个$lim$就可以了。设$f(k)$表示出现$s$次的颜色至少有$k$种的方案数，则$$f(k) = \binom{m}{k}\binom{n}{ks}\f...

2019-01-19 09:53:00 94

转载 Luogu 3321 [SDOI2015]序列统计

BZOJ 3992点开这道题之后才发现我对原根的理解大概只停留在$998244353$的原根是$3$……关于原根：点我首先写出$dp$方程，设$f_{i, j}$表示序列长度为$i$当前所有数乘积模$m$为$j$的方案数，有转移$$f_{i, x * y \mod m} = \sum_{y \in s} f_{i - 1, x}$$把$x$和$y$取个对数就可以变成...

2019-01-18 20:02:00 69

转载 Luogu 3702 [SDOI2017]序列计数

BZOJ 4818感觉不难。首先转化一下题目，“至少有一个质数”$=$“全部方案”$ - $“一个质数也没有”。注意到$m \leq 2e7$，$[1, m]$内的质数可以直接筛出来。设$f_{i, j}$表示当前长度序列为$i$，当前和模$p$的值是$j$的方案数，直接无脑枚举$m$转移复杂度是$O(nmp)$的，但是发现每一次转移形式都是相同的。$$f_{i, x...

2019-01-18 08:26:00 95

转载 CF 990 Educational Codeforces Round 45

既然补了就简单记录一下。感觉还算有一点营养。官方题解传送门：点我ACommentary Boxes对拆掉$n \mod m$个和新建$m - (n \mod m)$求个最小。#include <cstdio>#include <cstring>#include <algorithm>using namespac...

2019-01-17 21:26:00 95

转载 Luogu 4705 玩游戏

看见这个题依稀想起了$5$月月赛时候的事情，到现在仍然它感觉非常神仙。游戏$k$次价值的期望答案$$ans_k = \frac{1}{nm}\sum_{i = 1}^{n}\sum_{j = 1}^{n}(a_i + b_j)^k$$二项式定理展开$$ans_k=\frac{1}{nm}\sum_{i = 1}^{n}\sum_{j = 1}^{m}\sum_{t = 0...

2019-01-17 14:11:00 71

转载 CF 438E The Child and Binary Tree

BZOJ 3625吐槽BZOJ上至今没有卡过去，太慢了卡得我不敢交了……一件很奇怪的事情就是不管是本地还是自己上传数据到OJ测试都远远没有到达时限。本题做法设$f_i$表示权值为$i$的二叉树的个数，因为一棵二叉树可以通过左右儿子构建起来转移，我们可以得到转移：$$f_w = \sum_{x, y, w - (x + y) \in c} f_x * f_y$$...

2019-01-17 08:34:00 94

转载 Luogu 4726 【模板】多项式指数函数

补补补……这个题的解法让我认识到了泰勒展开的美妙之处。泰勒展开泰勒展开就是用一个多项式型的函数去逼近一个难以准确描述的函数。有公式$$f(x)\approx g(x) = g(x_0) + \frac{g'(x_0)}{1!}(x - x_0) + \frac{g^{(2)}(x_0)}{2!}(x - x_0)^2 + \cdots + \frac{g^{(n)}(...

2019-01-16 19:13:00 115

转载 Luogu 4725 【模板】多项式对数函数

继续补全模板。要求$$g(x) = ln f(x)$$两边求导，$$g'(x) = \frac{f'(x)}{f(x)}$$然后左转去把多项式求导和多项式求逆的模板复制过来，就可以计算出$g'(x)$，接下来再对$g'(x)$求不定积分即可。虽然我也不是很会不定积分，但是这就是求导的逆过程，相当于把求完导之后的函数搞回去。因为$(a_ix^i)' = ia_i...

2019-01-16 15:10:00 90

转载 Luogu 4245 【模板】任意模数NTT

这个题还有一些其他的做法，以后再补，先记一下三模数$NTT$的方法。发现这个题不取模最大的答案不会超过$10^5 \times 10^9 \times 10^9 = 10^{23}$，也就是说我们可以取三个满足$NTT$性质的模数先算然后再合并起来。比如三个模数可以分别取$998244353, 1004535809, 469762049$。那么我们现在要做的就是合并三个同余方...

2019-01-16 14:19:00 128

转载 Luogu 4238 【模板】多项式求逆

疯狂补板中。考虑倍增实现。假设多项式只有一个常数项，直接对它逆元就可以了。现在假如要求$G(x)$$$F(x)G(x) \equiv 1 (\mod x^n)$$而我们已经求出了$H(x)$$$F(x)H(x) \equiv 1(\mod x^{\left \lceil \frac{n}{2} \right \rceil})$$两式相减，$$F(x)(G(...

2019-01-16 13:52:00 126

转载 Luogu 2000 拯救世界

从胡小兔的博客那里过来的，简单记一下生成函数。生成函数数列$\{1, 1, 1, 1, \cdots\}$的生成函数是$f(x) = 1 + x + x^2 + x^3 + \cdots$，根据等比数列求和公式，可以得到$f(x) = \frac{1}{1 - x}$。把两边分别平方，得到$$\frac{1}{(1 - x)^2} = (1 + x + x^2 + x^3...

2019-01-16 08:21:00 71

转载 LOJ 2537 「PKUWC2018」Minimax

BZOJ 5461。线段树合并优化$dp$。假设所有离散之后的权值$\in [1, m]$，对于一个点$x$它的权值是$i$的概率是$f(x, i)$，那么1、假如这个点只有一个儿子$y$，那么$f(x, i) = f(y, i)$。2、假如这个点有两个儿子$y, z$，那么$$f(x, i) = f(y, i)\sum_{j = 1}^{m}f(z, j)(p_x[...

2019-01-15 10:51:00 100

转载 CF 1097D Makoto and a Blackboard

算是记一下昨天晚上都想了些什么官方题解点我简单题意给定两个正整数$n$和$k$，定义一步操作为把当前的数字$n$等概率地变成$n$的任何一个约数，求$k$步操作后的期望数字，模$1e9 + 7$。$$n \leq 10^{15}, k \leq 10^4$$我的思路设$f(n, k)$表示$n$在$k$步操作之后的期望数字，假设$n$的约数有$m$个，分别为...

2019-01-05 10:07:00 79

转载 CF 662C Binary Table

用FWT优化计算。首先发现行数很小，想到一个暴力的方法，就是以一个二进制位$0$表示这一行不翻转而二进制位$1$表示这一行翻转，然后$2^n$枚举出所有行的翻转情况，再$O(m)$计算所有的结果。用$a_i$表示第$i$列的原来的情况，有计算式：$$ans_s = \sum_{i = 1}^{m}(a_i \oplus s) * min(bit_{a_i \oplus s}...

2019-01-05 09:12:00 89

转载 Luogu 3768 简单的数学题

莫比乌斯反演 + 杜教筛$$\sum_{i = 1}^{n}\sum_{j = 1}^{n}ijgcd(i, j)$$$$= \sum_{d = 1}^{n}d\sum_{i = 1}^{n}\sum_{j = 1}^{n}ij[gcd(i, j) == d]$$$$= \sum_{d = 1}^{n}d^3\sum_{i = 1}^{\left \lfloor \frac{...

2019-01-02 22:31:00 85

转载 Luogu 4213 【模板】杜教筛（Sum）

当作杜教筛的笔记吧。杜教筛要求一个积性函数$f(i)$的前缀和，现在这个东西并不是很好算，那么我们考虑让它卷上另外一个积性函数$g(i)$，使$(f * g)$的前缀和变得方便计算，然后再反推出这个$f$函数的前缀和。$$\sum_{i = 1}^{n}(f * g)(i) = \sum_{i = 1}^{n}\sum_{d | i}g(d)f(\frac{i}{d}) = ...

2019-01-02 21:12:00 55

转载 Luogu 4781 【模板】拉格朗日插值

模板题。拉格朗日插值的精髓在于这个公式$$f(x) = \sum_{i = 1}^{n}y_i\prod _{j \neq i}\frac{x - x_i}{x_j - x_i}$$其中$(x_i, y_i)$是给定的$n$个点值。代入任何一个给定的点值坐标$x_k$，都会发现这个式子等于$y_k$成立，因为对于任何$i \neq k$，后面的系数都至少有一项为$0$，而...

2019-01-01 20:53:00 97

转载 CF 932E Team Work

推式子题。首先有公式：$$n^k = \sum_{i = 0}^{k}\binom{n}{i}S(k, i)*i!$$其中$S$表示第二类斯特林数。左边表示$k$个不同的小球放$n$个不同的盒子允许有空盒的方案数，而右边先枚举非空盒的数量，选择非空的盒子，然后再把$k$个小球全部放入，因为盒子是不同的，所以还要乘上$i!$。把公式代进去：$$\sum_{i = 1...

2019-01-01 14:45:00 81

转载 CF 961G Partitions

推不动式子我们考虑每一个$w_i$对答案的贡献，因为题目中定义集合的价值为$W(S) = \left | S \right |\sum_{x \in S}w_x$，这个系数$\left | S \right |$可以看作集合中所有的元素（包括$i$自己）对$i$产生了一次贡献，那么我们考虑一个元素$j$对$i$的贡献：1、$j == i$的时候，相当于求把$n$个小球放到$k$个...

2019-01-01 12:25:00 92

转载 Luogu 4844 LJJ爱数数

LOJ 6482设$d = gcd(a, b)$，$xd = a$，$yd = b$，因为$\frac{1}{a} + \frac{1}{b} = \frac{a + b}{ab} = \frac{1}{c}$，所以$c(x + y)= xyd$。因为$d$不整除于$c$，那么$d | (x + y)$，把$d$除过去，$$\frac{x + y}{d} = \frac{xy...

2019-01-01 12:09:00 74

转载 CF 1091E New Year and the Factorisation Collaboration

昨晚Good Bye 2018D题没做出来，车翻大了……官方题解传送门初赛知识：一个无向图所有顶点度数之和为偶数。然而这东西还有一个高端的名字：Handshaking lemma但是这并不是本题的重点，另外一个看上去很高端的东西才是本题的重点：Erdős–Gallai theorem对于一个无向图的度数序列$d$，先从大到小排序，即满足$d_1 \geq d_...

2018-12-31 12:17:00 127

转载 Luogu 4841 城市规划

BZOJ 3456 权限题太菜了推不出式子我们设$f(n)$表示$n$个点的无向连通图的数量，那么有$$f(n) = 2^{\binom{n}{2}} - \sum_{i = 1}^{n - 1}\binom{n - 1}{i - 1}f(i)2^{\binom{n - i}{2}}$$思路就是全部减去不合法的，枚举$1$号点所在的联通块的大小，剩下随便生成一张无向图。...

2018-12-30 19:38:00 56

转载 Luogu 4721 【模板】分治 FFT

还不会这题的多项式求逆的算法。发现每一项都是一个卷积的形式，那么我们可以使用$NTT$来加速，直接做是$O(n^2logn)$的，我们考虑如何加速转移。可以采用$cdq$分治的思想，对于区间$[l, r]$中的数，先计算出$[l, mid]$中的数对$[mid + 1, r]$中的数的贡献，然后直接累加到右边去。容易发现，这样子每一次需要用向量$[l,l + 1, l + ...

2018-12-29 19:46:00 69

转载 Luogu 4091 [HEOI2016/TJOI2016]求和

BZOJ 4555一道模板题。第二类斯特林数有公式：$$S(n, m) = \frac{1}{m!}\sum_{i = 0}^{m}(-1)^i\binom{m}{i}(m - i)^n$$考虑它的组合意义：$S(n, m)$表示$n$个不相同的小球放到$m$个相同的盒子里而且不能有空盒的方案数。我们枚举空盒有$i$个，然后进行容斥。因为盒子没有区别，所以最后得到的值...

2018-12-27 10:53:00 75

转载 Luogu 3723 [AH2017/HNOI2017]礼物

BZOJ 4827$$\sum_{i = 1}^{n}(x_i - y_i + c)^2 = \sum_{i = 1}^{n}(x_i^2 + y_i^2 + c^2 - 2 * x_iy_i + 2c * x_i - 2c * y_i) = \sum_{i = 1}^{n}x_i^2 + \sum_{i = 1}^{n}y_i^2 + nc^2 + (2\sum_{i = 1}^{...

2018-12-25 21:13:00 98

转载 FFT笔记

打开了政治书又合上了。学考爆炸预定。趁着自己还能记得一点东西稍微做一下笔记，要不然又会跟联赛前一样忘得一干二净。attack大佬的blog复数定义$i^2 = -1$，那么所有复数可以表示成$a + bi$的形式（$a, b$为实数）。在复平面中，$x$轴代表实数，$y$轴（除了原点外）代表虚数。那么复数$(a + bi)$可以写成形如$(a, b)$的向量，...

2018-12-22 20:22:00 87

转载 Luogu 4900 食堂

一道把很多东西放在一起的练手题。$$\sum_{i = A}^{B}\sum_{j = 1}^{i}\left \{ \frac{i}{j} \right \} = \sum_{i = A}^{B}\sum_{j = 1}^{i}\frac{n - \left [ \frac{n}{j}\right] * j}{i} = \sum_{i = A}^{B}(i * \sum_{j = ...

2018-12-12 20:28:00 48

转载 Luogu 4155 [SCOI2015]国旗计划

BZOJ 4444倍增 + 贪心。发现是一个环，先按照套路把环断开复制一倍，这样子的话覆盖完整个环就相当于覆盖一条长度不小于$m$的链，注意这样子有一些区间在新的这条链上会出现两次。我们为了找到最小的满足要求的答案，在选择完一个区间$[l, r]$之后会选择左端点不超过$r$但是右端点尽量大的区间，因为题目保证了所有的区间不相互包含，这样子的话我们只要找到左端点最靠右的区间就...

2018-12-12 14:54:00 107

转载 Luogu 4069 [SDOI2016]游戏

BZOJ 4515树链剖分 + 李超线段树要求支持区间插入一条线段，然后查询一个区间内的最小值。可以使用李超线段树解决，因为要维护一个区间内的最小值，所以每一个结点再维护一个$res$表示这个区间内的最小值。因为本题把问题搬到了树上，剖一下就可以了。我们可以把一个点$x$到根的距离$dis_x$记为一个点在二维平面上的横坐标，这样子可以保证一条重链上的点的距离递增，并且不改变...

2018-12-10 16:58:00 62

转载 Luogu 4254 [JSOI2008]Blue Mary开公司

BZOJ 1568学习了一波李超线段树。大佬blog这个东西专门用来维护插入一条线段和区间/单点的最大/最小值。插入的时候讨论：1、如果当前结点上没有线段，那么直接插入。2、如果当前结点上的线段一定比要插入的线段优/劣，那么直接覆盖或者返回。3、如果当前结点上的线段和要插入的线段有交点，那么把优的部分比劣的部分多的线段放在当前结点上，然后把另一条线段下放...

2018-12-05 19:58:00 67

转载 Luogu 4251 [SCOI2015]小凸玩矩阵

BZOJ 4443二分答案 + 二分图匹配外层二分一个最小值，然后检验是否能选出$n - k + 1$个不小于当前二分出的$mid$的数。对于每一个$a_{i, j} \geq mid$，从$i$向$j + n$连一条边，然后跑二分图最大匹配即可。菜的很，二分图匹配都写不对……注意数组要开到双倍$n$。时间复杂度$O(nmlogn)$。Code：#in...

2018-12-02 08:54:00 64

转载 Luogu 5043 【模板】树同构（[BJOI2015]树的同构）

BZOJ 4337简单记录一种树哈希的方法：以$x$为根的子树的哈希值为$\sum_{y \in son(x)}f_y*base_i$，$f_y$表示以$y$为根的树的哈希值，其中$i$表示$f_y$在若干个儿子中的排名，$base$是$rand$出的对一个质数取模之后的很大的数。对于本题这样的情况，可以每一个结点都拿出来作为根计算一下，然后再把所有的结果排个序，如果两棵树同构那...

2018-12-01 17:34:00 80

转载 NOIP2018 解题笔记

D1T1 铺设道路　　在场上并没有想到积木大赛这道原题。　　差分之后可以把在$[l, r]$这段区间$ - 1$变成在$l$处$ - 1$，在$r + 1$处$ + 1$，然后最终目标是使$\forall i \in [1, n] \ \Delta d_i == 0$成立。就想着把正负数配一配对，然后输出了正数绝对值和负数绝对值的$max$，这导致了我的代码非常鬼畜。　　出来...

2018-11-24 18:03:00 100

转载 CF916E Jamie and Tree

这大概是NOIP前最后一道题了，顺便复习一些模板。考虑到不是真正要进行“换根”这个操作，我们发现在查询子树和修改子树的时候只要按照遥远的国度这题的方法分类讨论一下就好了，我们考虑一下如何换根意义下的$lca$。先分类讨论，假设当前结点是$(u, v)$根是$rt$。1、$u, v$都在$rt$的子树内，那么$lca$不变。2、$u, v$有一个在$rt$的子树内（假设为$...

2018-11-08 14:35:00 85

空空如也

空空如也