一位蒟蒻的小博客

一道字符串结论+数论的题

一个经典计算几何问题辛普森积分，圆的面积并

Min25筛基本应用。

有关斐波那契拆分的一道题目

题意给定一个长度为n的序列和一个常数p%，有三种操作：区间加区间赋值输出所有区间内占比>=p%的所有元素。但是可以包含<=p%的元素。n≤105,p>=20n\leq10^5,p>=20n≤105,p>=20分析找一个最小的k使得p%>1/kp\%>1/kp%>1/k。将限制放宽到占比>1/k，即出现次数大于n/k下取整n/...

题意给你一个仙人掌，随机选一个点排列。按照这个排列的顺序将所有点删掉，求每个点被删掉时连通块大小之和的期望。n≤400n\leq 400n≤400分析树上一条路径可行的概率是1/len.环上做一个容斥：路径A可行+路径B可行-路径AB均可行。仙人掌上做dp容斥即可。仙人掌上有几个处理的技巧：仙人掌不需要tarjan. 暴力建出所有点向环顶连边的树。做dp的时候，所有边可以被...

ETT维护边连通块括号序解决同色连通块问题。

分析猜想一下这个函数是分奇偶凸的。事实上也的确是这样。考虑分治之后合并的过程。已经处理出了左右两边选多少个的最大值和最小值。然后按照奇偶分类地合并即可得到新的值。可以归纳地证明他是凸的。考虑更新的最优解如何求得即可。每次必定是左边多选两个或者右边多选两个，取最优的即可。这个东西可以推广到求闵可夫斯基和的凸包上，具体证明大致是差不多的。#include <bits/stdc+...

本文介绍了多项式求逆，开方，exp，ln，次幂，复合逆的基本概念与思想。并附有简略求解方法与实现。

树分块区间众数

简略题意给一个序列，支持三种操作：查询某个位置的值区间赋值不断给某个区间加一个数，直到其中不存在42的幂为止。n≤105，所有输入的数在109范围内n\leq 10^5，所有输入的数在10^9范围内n≤105，所有输入的数在109范围内思路看到这种操作基本上就是个势能分析线段树题构不出太大的数，所以每个数能加的次数就是log次。将每个位置设为他到下一个幂的距离（真实值为ne...

2-sat问题的求解对称性2 satO(n+m)

题意有n个点，n-1条边。每条边是(i,i+1)(i,i+1)(i,i+1)，长度为w[i]w[i]w[i]。每走1长度要消耗1的油，每走到这个点上可以获得g[i]g[i]g[i]的油。你有K次机会给某个点的g[i]+1g[i]+1g[i]+1。求最大的满足下述条件的区间长度：满足存在一个给油的方案，使可以从L不回头的走到R，也可以从R不回头的走到L。思路 O(nlog⁡2n)O(n\l...

好像是一种比较新的筛法，网上资料都是18年的 赶上时代潮流了？？用途筛一些比较神奇的函数前缀和。（或者询问区间内什么的） 筛的东西需要是积性函数。（题目所求可以转换为积性函数啥的）复杂度空间n^0.5，时间据说为n0.75/lognn0.75/lognn^{0.75} / logn 2s大概可以做两次1e10吧。思想递推模拟线筛。具体步骤设g(i,j)...

题解考虑多组询问，我们需要找到一种对于任意一个起点都可以快速计算的方法。 假如设FiFiFi表示从i走到终点的答案，那么对于不同的时间，难以将参数T分离。 于是，假如我们可以计算GiGiGi表示点i的期望经过次数，一条边(x,y)的贡献就是Gx⋅P(x,y)⋅LenGx⋅P(x,y)⋅LenGx \cdot P(x,y) \cdot Len，那么就可以根据边的p分成两部分，统计...

描述一棵树，每次询问某个点与他距离第k小的点。询问数，点数1e5.要点先建点分树。在点分树上统计。每个分治中心开数组记录其范围内所有到他的距离，排序之后就可以快速统计了。会算重，每次要减去当前分治中心范围在上一级分治中心中的贡献（不能来自同一颗子树）。多维护一个距离数组表示其到上一级的距离就可以计算了。#include &amp;lt;cstdio&amp;gt;#include &amp;lt;i...

描述n&amp;lt;=300，给定有权边，求生成树大小和所有生成树边权乘积和。要点基尔霍夫矩阵：c[i][i]c[i][i]c[i][i]为点i的度数，c[i][j]=−(i,j之间边数)c[i][j]=-(i,j之间边数)c[i][j]=−(i,j之间边数)行列式：枚举每一个1…n的排列，将每行对应的列乘起来， 再乘上(−1)逆序对个数(-1)^{逆序对个数}(−1)逆序对个数之和。PTY...

Description一个无向图，求从i到j简单路径（无重点）上的最小点权。n,m&lt;=1e5分析考虑走到某个点双上时，假如不是同点进同点出，就可以获得这个点双上的所有贡献。于是建出圆方树，在方点开个multiset维护点权最小值。每次询问就是链上最值。但是考虑修改，因为一个点可能在多个点双内，要摊一下修改和查询的复杂度。我们使得每个方点只管儿子，不管父亲。 这样假如lca是方...

Description懒得写背景了，给你一个字符串init，要求你支持两个操作(1):在当前字符串的后面插入一个字符串(2):询问字符串s在当前字符串中出现了几次？(作为连续子串)你必须在线支持这些操作。毒瘤显然就是sam求endpos大小，强套LCT。时隔多年又打了一次维护子树大小。大致就是维护：fsum[x]表示虚父亲为x的tsum和。tsum[x]表示x的splay和，=...

题面 n&lt;=50巧妙构图想了好久流还是不会 注意题目给的奇偶性质 ： 一个L中间必定X+Y=奇数。 显然L除了中间之外的旁边的两个点一定是X+Y=偶数的。 仔细观察bing就bu能发现，这两个点一个是偶+偶，一个是奇+奇。 根据这个构图，就变成了简单的最大费用可行流。 别忘了最大流量限制为m. 关于可行流，有一个结论是，当最长路小于0时，不再进行增广。 感性理...

水平有限，此处只记录一些OI中线性基的思路，基本与线性代数无关了。定义线性基一般用于求解异或问题。 对n个数求其线性基，得到最多 二进制下位数m 个数。其满足三个性质： 1. 线性基能相互异或得到原集合的所有相互异或得到的值。 2. 线性基是满足性质1的最小的集合 3. 线性基没有异或和为0的子集。注意一点，线性基内不一定是原数。 但我们称被成功插入线性基的数为在线性基中...

例题jzoj3658文本编辑器无旋Treap有旋Treap由于splay的存在而无用武之地了。 优点：代码较短，原生支持区间分裂合并，并支持可持久化。 缺点：较splay大概有2的常数。核心操作定义请戳split(root,x) 将树root分解为两颗树，第一颗为其中[1,x]，第二颗为[x+1,size] 返回是保存上述两个根的pair。merge(a,...

题意四边形不等式对于二元函数w(i,j)，若任意a&amp;lt;=b&amp;lt;=c&amp;lt;=da&amp;lt;=b&amp;lt;=c&amp;lt;=daw(a,c)+w(b,d)&amp;lt;=w(a,d)+w(b,c)w(a,c)+w(b,d)&amp;lt;=w(a,d)+w(b,c)w(a,c) + w(b,d) w(i,j)&amp;gt;=w(i′,j′)，[i,j]包含

mk&lt;=5e6，m&lt;=5e4mk&lt;=5e6，m&lt;=5e4mk[n=1]=∑m的集合划分A&nbsp;&nbsp;&nbsp;&nbsp;∏(ai−1)!⋅(−1)ai−1[n=1]=∑m的集合划分A&nbsp;&nbsp;&nbsp;&nbsp;∏(ai−1)!⋅(−1)ai−1[n=1] = \sum_{m的集合划分A} ~~~~\prod(ai-1)! \cdot ...

题面梦游中的你来到了一棵 N 个节点的树上. 你一共做了 Q 个梦, 每个梦需要你从点 u 走到 点 v 之后才能苏醒, 由于你正在梦游, 所以每到一个节点后,你会在它连出去的边中等概率地 选择一条走过去, 为了确保第二天能够准时到校, 你要求出每个梦期望经过多少条边才能苏 醒. 为了避免精度误差, 你要输出答案模10^9 + 7的结果. 对于 100%的数据, N &lt;= 100000...

前置知识1.群 2.置换 3.置换群问题背景求解等价染色问题。（万恶的本质相同） 比如用k种颜色给一个2*2的方格染色，并认为旋转后相同的方案是本质相同的。 求有多少本质不同的方案。Burnside定理上面的问题也就是在求等价类个数了。 这个问题中，置换群就是恒等置换，转90，转180，转270。 先给出式子，将在后面证明。 等价类个数=∑C(f)|G|等价类个...

分析死因：思路错了一开始在考虑欧拉序和原序列单调栈的问题，这样想其实可以分治（超麻烦）。注意到一个结论，假如你要求n个点的lca，那么你可以以任意顺序排序，然后对相邻求lca，再求深度最小的即可。证明很显然，考虑欧拉序，答案肯定会至少被一组相邻的点盖到，这样问题就很简单了，做两遍单调栈再枚举答案即可。其实还有两种做法：将dep转化为个数，这样其实就是查询子树内在给出序列内所有连续段...

欧几里得算法又名辗转相除法形式gcd(a,b) = gcd(b,a mod b)证明1.设a=kb+r,r=amodb=a−kba=kb+r,r=a mod b=a-kb2.设dd为a,ba,b的公约数，则amoda mod d=0，bmodd=0，b mod d=0d=0 记作d|a、d|bd|a、d|b3.r=a−kbr=a-kb由2得，rd\dfrac{r}{d}是一个整数，所以d|rd|

并不是什么新方法，老早就有了。背景二进制的异或FWT实际上在做这样一件事：构造一个行列式不为0的Len阶方阵T，使得TA⋅TB=TCTA\cdot TB=TCTA⋅TB=TC将矩阵乘法拆开看，这个矩阵满足：∀x,i,j    T(x,i)T(x,j)=T(x,i⊕j)\forall {x,i,j}~~~~T(x,i)T(x,j)=T(x,i\oplu...

题意A,B,CA, B, CA,B,C都是n阶01方阵，给定CCC，现求有多少对有序(A,B)(A,B)(A,B)满足A∗B=CA * B = CA∗B=C.乘法是模2意义下的矩阵乘法。n≤2000n \leq2000n≤2000？迷得很这个题模2意义就相当于加法是xor，数量乘法是and，仍然能组成一个线性空间。因此下面的讨论是在数域0,1数域{0,1}数域0,1下的进行的。首先...

小声bb去年做了保卫王国之后感觉自己会动态dp了，今天看下全都还给出题人了…省选前临急抱佛脚吧。动态dp是什么就是先给一个dp问题，然后修改里面的一些值 / 给出一些限制，要你快速求出新的dp答案。方法大概分三种：倍增（适用于无权值修改，只有限制的情况）树链剖分LCT其中树剖是比较不推荐的，太长了。但是树剖比较好理解，可以用这个来学。注意不要拘泥于原有的完整树形结构。为了...

对于长度为x的平方串，只需要每隔x做一个关键点，然后对相邻关键点做lcp和lcs就可以找出每一个平方串。sa或者hash+二分都是可以的找出平方串的区间，下面要实现的操作就是实现区间中x和x+le合并。用倍增维护一下并查集，开log个并查集，若x,y在第i个并查集里被并起来，意味着x..x+2i−1x..x+2^i-1x..x+2i−1与y..y+2i−1y..y+2^i-1y..y+2i...

库默尔定理：C(n,m)C(n,m)C(n,m)中含有的质数p幂次= n−mn - mn−m在p进制下的借位次数= m+(n−m)m + (n - m)m+(n−m)在p进制下的进位次数。原理：Ans=∑k&gt;0npk−(n−mpk+mpk)Ans = \sum_{k&gt;0}\frac n {p^k}-(\frac{n-m}{p^k} + \frac{m}{p^...

例题给定G(次数界3w左右)，求F(x)=11−G(x)F(x)=\frac 1 {1-G(x)}F(x)=1−G(x)1​的第n项，保证G(x)[x0]=0G(x)[x^0]=0G(x)[x0]=0。n在1e9以上。有毛病毒瘤毒瘤毒瘤式子就是F(x)=1+F(x)G(x)F(x)=1+F(x)G(x)F(x)=1+F(x)G(x)这是什么意思呢其实就相当于一个齐次线性递推。f...

简单模板题：6121. 【GDOI2019模拟2019.4.13】数据结构首先建树，一层按x二分，一层按y二分交替。用algorithm中的nth_element(开始，中间，结束（开），比较函数)就可以把点分成两边。但是两边是乱序的。这是线性的.然后把每个点控制的矩形区域更新上来，接下来就类似二维线段树了。据说询问一个矩阵是O(n0.5)O(n^{0.5})O(n0.5)的。证明不显然...

存个板子首先加限制的四条边，然后先去平行，然后排极角序。每次加入一条边的时候，若队头两条线交点不在新半平面内，就出掉队尾。然后队头类似最后记得去掉尾部多加的半平面。jzoj6093#include <cstdio>#include <iostream>#include <cstring>#include <algorithm>#in...

迭代法求零点已知fff是一个多项式对多项式的函数，现要逼近f的零点，采用倍增法：假如要求A模xnx^nxn意义下的值，则预先求出A0为A模xn/2x^{n/2}xn/2意义下的值由泰勒展开，f(A)=f(A0)+f′(A0)(A−A0)f(A)=f(A0)+f&amp;#x27;(A0)(A-A0)f(A)=f(A0)+f′(A0)(A−A0)在模xnx^nxn意义下成立。令f(A)=0f(...

很常见的推式子套路，有个组合数取模的黑科技：现在要求C(n,1),C(n,2)…C(n,m) (n&lt;=1e9,m&lt;=1e6)对一个大数p取模的结果，可以先将p分解质因数，然后将出现在p中的质因子单独计算指数以避免除法，其余的部分是与p互质的，直接使用逆元计算。（分解质因数时顺便求phi即可）当只有乘除的时候就可以使用这种规避逆元的方法，计算一次乘除的复杂度是质因子个数的。最后...

看到gcd，直接求也不好求，就可以先考虑一下反演。然后发现答案就是∑f(u)⋅ϕ\sum f(u)\cdot \phi∑f(u)⋅ϕ，f(u)是gcd是u的倍数的方案数。有个结论是∑d∣xμ(d)⋅(x/d)=ϕ(x)\sum_{d|x}\mu(d)\cdot(x/d)=\phi(x)∑d∣x​μ(d)⋅(x/d)=ϕ(x)，就是一个容斥，最终只有与x无公共质因子的数会被算到。首先将所有数...

失智又被题目吓到了考虑树的情况，连通块数目就是总点数-存在的边。考虑仙人掌的情况，连通块数目就是总点数-存在的边+存在的环因为是0/1分开考虑，所以所谓存在的X就是要求相关联的点颜色一样。使用lct维护圆方树就可以了。算一个环都变黑的概率可以容斥，因为总环中点数是不超过n+m的，暴力就可以了。#include &lt;cstdio&gt;#include &lt;iostrea...