扩展の灰(Extended Ash/Cooevjnz)

原创 一些基础算法的模板(持续更新)

更新中//Templates From Extended_Ash/Cooevjnz/JacaJava/Tubbcrafft//To be continued...//Suffix Automationchar str[N]; int s[N][26],mx[N],f[N],sz[N]; int last=1,cnt=1,n,v[N],r[N],ans=0; inlin

原创 OI一些奇怪的优化手段和常用技巧

经典的O3优化（一般写在开头）#pragma GCC optimize("O3")#pragma G++ optimize("O3")G++手动扩大栈（写在main的开始）int size = 256 << 20; // 256MB char *p = (char*)malloc(size) + size; __asm__("movl %0, %%esp\n" :: "r"(

原创 容斥原理基本思想

基本原理：对于f(n),g(n)f(n),g(n)f(n),g(n)如果有以下关系g(i)=∑j=inf(j)∗C(j,i)g(i)=\sum_{j=i}^nf(j)*C(j,i)g(i)=j=i∑n​f(j)∗C(j,i)而且存在系数F(n),G(n)F(n),G(n)F(n),G(n)满足F(n)=∑i=1nC(n,i)∗G(i)F(n)=\sum_{i=1}^nC(n,i)*G(i)F(n)=i=1∑n​C(n,i)∗G(i)那么就可以得到结果∑i=1nf(i)∗F(i)=∑i=1ng(i)∗G(

原创 一些反演求和过程？

.问题∑i,j&amp;lt;=niaja(i,j)b\sum_{i,j&amp;lt;=n}i^aj^a(i,j)^bi,j&lt;=n∑​iaja(i,j)b记f(d)=∑(i,j)=diajaF(d)=∑d∣i,d∣jiajaf(d)=\sum_{(i,j)=d}i^aj^a \\ F(d)=\sum_{d|i,d|j}i^aj^af(d)=(i,j)=d∑​iajaF(d)=d∣i,d∣...

原创 UOJ62 怎样跑的更快

bi=∑j=1n(i,j)c−d∗idjdzjb_i=\sum_{j=1}^n(i,j)^{c-d}*i^dj^dz_jbi​=j=1∑n​(i,j)c−d∗idjdzj​biid=∑j=1nzjjd(i,j)c−d\frac{b_i}{i^d}=\sum_{j=1}^nz_jj^d(i,j)^{c-d}idbi​​=j=1∑n​zj​jd(i,j)c−d记(i,j)=x,Bi=biid,Zi=...

原创 LOJ 6495~6497「雅礼集训 2018 Day1」

.真是看懂std就秒的三道题，但是看懂就很难啊（感谢superguymj和memset0两个大佬的代码）稍微说一下题解-------------------------------------------------------- 「雅礼集训 2018 Day1」树这道题的正解被superguymj大佬吊打了，原来O(n⋅2n)O(n·2^n)O(n⋅2n)变成了O(n4)O(n^4)...

原创 一般做题流程(雾)

.如果一道题没有题解怎么办呢介绍一个非常有效的办法：首先打开LOJ点标签雅礼集训点开没做的第1道题，看题，并暗示自己这题很简单，稍微看一下std就会了点“统计”点开superguymj的代码看（如果没有就随便找一个）看完并理解后开始写代码多次尝试提交然后AC很快你就会成为神仙辣我太菜了永远不能成为神仙距离THUWC还有一个月吧，共勉...

原创 PKUSC2018部分题解

.LOJ的std+数据+当时讲题的记忆硬是刚了五道题（D1t3咕咕咕）稍微写一下题解LOJ6432 「PKUSC2018」真实排名一道比较简单的题，不过需要注意很多细节考虑排名不变的两种情况1.自己分数xxx不变，那么所有分数在[[x2],x)[[\frac{x}{2}],x)[[2x​],x)中的其他人的分数也不能变，否则必然有在自己后面的人到自己前面2.自己分数xxx翻倍，那么...

原创 Bzoj1185 [HNOI2007]最小矩形覆盖 + Bzoj1069 [SCOI2007]最大土地面积

.题目链接1题目链接22道旋转卡壳算法+凸包模板题应该说是一个比较标准的算法了，凸包就是排序+单调栈搞一下就好旋转卡壳就是一个运用了凸包单调性的思想，每次逆时针走一条边，就把对点移动一个第一题稍微麻烦，需要枚举矩形的一条边，让后找水平方向和垂直方向最远距离，分别用点积和叉积就可以判第二题直接枚举对角线，让后找左右点到直线最长距离稍微注意基本运算和函数的实现就可以了，剩下的都是可以现...

原创 计算几何模板

超级大毒瘤目前有仅有基本运算和求凸包以及面积估计很快会加上半平面交&amp;旋转卡(qia)壳(ke)动态凸包……这辈子都不会有的（真香）Debugging…#include&lt;math.h&gt;#include&lt;stdio.h&gt;#include&lt;string.h&gt;#include&lt;algorithm&gt;#define eps 1e-7...

原创 KD-Tree模板

.题：Bzoj2850参考资料好像其实挺好写的（就是慢了很多很多），复杂度是O(kn1−1k)O(kn^{1-\frac{1}{k}})O(kn1−k1​)只能算是一个“分块级别”的暴力，离线段树这种还差很多#include&lt;stdio.h&gt;#include&lt;string.h&gt;#include&lt;algorithm&gt;#define N 50010...

原创 Bzoj2434 阿狸的打字机

纪念一次AC，水一篇blog又是喜闻乐见的广义SAM题辣题意：给一个trie，每次问询一个串在另一个串出现的次数思路：建立广义SAM，离线所有询问，每次讲一个主串在fail树上的节点+1，处理询问就用fenwick维护子树和这里注意，1.广义SAM要用bfs建2.每次从一个主串到下一个时，不要整个都减掉而是要在LCA处会和，这样复杂度就是对的代码有点长但是很好打#include&...

原创 LCT模板

此题是bzoj2631Lct真的是看一眼就懂的数据结构，Tarjan太强了考虑轻重链剖分，我们可以将一条路径分成log个区间而在Lct上，用splay来完成这一过程，每次将节点到根的路径强制变为重路径(Access操作)这样，就可以很方便的来完成链上的查询和修改了而换根，需要维护一个Reverse标记，每次access(x)access(x)access(x)之后将整颗splay翻转反...

原创 ETT模板

其实这个是弹飞绵羊辣ETT比lct好写很多简单来说就是维护树的dfn序列，因为其长度为2n所以效率会不如lct这里只有换父亲的操作（用一个超级根作为所有树的父亲，就可以支持link和cut），维护的信息是节点的深度，换根到时候再学说一下几个实现的技巧可以用splay来维护，由于需要用到抽出一段的操作，建议开一个数组存每个节点的前驱和后继，这样可以省去找前驱后继的常数，但是依然没有LCT的...

原创 Latex代码测试

#include&amp;amp;lt;queue&amp;amp;gt;\#include&amp;amp;lt;queue&amp;amp;gt;#include&amp;lt;queue&amp;gt;#include&amp;amp;lt;stdio.h&amp;amp;gt;\#include&amp;amp;lt;stdio.h&amp;amp;gt;#include&amp;lt;stdio.h&

原创 字符串专题1

.都是些bzoj原题辣，这几天刚做的 Bzoj4032有趣的dp题，也要用到各种自动机注意到题目的两个关键词“子串”和“子序列”考虑对A和B串建立后缀自动机和序列自动机序列自动机：可以识别一个序列所有子序列的自动机 想必学过自动机的各位都知道这个玩意怎么建，这里不再阐述让后我们考虑这些询问询问1：直接用SAM做类似LCS的做法求出A每个前缀和B匹配的最长距离，取最小+1询问2：...

原创 后缀数组SA-IS模板

.SAIS真的好快啊，在UOJ能卡进第二页推荐几个讲的比较好的博客链接1链接2链接3和一篇讲了各种后缀数组算法的论文Code:#pragma GCC optimize("O3")#pragma G++ optimize("O3")#include&lt;stdio.h&gt;#include&lt;string.h&gt;#include&lt;algorithm&gt;...

原创 自适应simpson积分模板

.好像还挺好写的具体思路就是一个区间，判断一下两半积分和整体差值是否过大，如果很大那么分成两边继续递归，否则直接返回（复杂度未知，常数很大）#include&lt;math.h&gt;#include&lt;stdio.h&gt;#include&lt;string.h&gt;#define db doubleusing namespace std;db A,B,C,D,l,r,...

原创 超级读入&输出优化

.mmapmmapmmap真的是黑科技，用了之后时间减半但是输出还是只能靠writewritewrite，不知道为什么用mmapmmapmmap会RERERE应该算是一个比较好的快速读写模板了吧，到时候再把输出也换成mmapmmapmmap（以下代码均在LinuxLinuxLinux环境下编译，windowswindowswindows的g++g++g++不行）#include&lt;s...

原创 51Nod1680 区间求和

.题目点这里大概是这么一个问题∑k=1n∑i=1n−k+1∑j=i+k−1n区间前k大元素和\sum_{k=1}^n\sum_{i=1}^{n-k+1}\sum_{j=i+k-1}^n区间前k大元素和k=1∑n​i=1∑n−k+1​j=i+k−1∑n​区间前k大元素和观察式子并加以变形=∑i=1n∑j=in∑k=1j−i+1区间前k大元素和=\sum_{i=1}^n\sum_{j=i}^n...

原创 序列求和V2

.题目：计算∑i=1nri∗ik\sum_{i=1}^nr^i*i^ki=1∑n​ri∗ik其中n≤109,k≤103n\le10^9,k\le10^3n≤109,k≤103一个非常巧妙的转化，设S(x)=∑i=1nri∗ixS(x)=\sum_{i=1}^nr^i*i^xS(x)=i=1∑n​ri∗ix我们要求的就是S(k)S(k)S(k)考虑将式子变形S(k)=∑i=1nri∗ik...

原创 配对堆模板

配对堆是一种可并堆题意：两种操作，合并两个堆或者查询一个堆的最小值，n≤106n\leq 10^6n≤106Pairing−HeapPairing-HeapPairing−Heap还挺好写的，不过并没有传说中那么快这里没有Decrease−KeyDecrease-KeyDecrease−Key的代码，说一下大概怎么实现额外记录一个FatherFatherFather域，执行Decrease...

原创 各种杂题合集

.好久没有往blog搬运题解了，稍微总结一下1.基环树一个非常套路的题目类型，都是可以直接复制代码的那种基本操作就是，dfs找环，让后拉出来倍长，最后在上面做dp即可CodeForces 835F Roads in the Kingdom大概就是要找删掉一条边后使得树最小的直径先dp环旁边的数，找到子树的到根最长链，设为fff倍长后令Fi=Si+fi,Gi=fi−SiF_i=S...

原创 伯努利数

.生成函数B(z)=zez−1=∑Bnznn!B(z)=\frac{z}{e^z-1}=\sum B_n\frac{z^n}{n!}B(z)=ez−1z​=∑Bn​n!zn​递推式(n2n^2n2)B0=1B_0=1B0​=1Bn=−∑k=0nBkn−k+1B_n=-\sum_{k=0}^n\frac{B_k}{n-k+1}Bn​=−k=0∑n​n−k+1Bk​​也可以写成这样∑k=0nC...

原创 多项式相关模板#2

.最近出了几个这样的题于是就把模板库完善了一下，依然是全程NTT实现，模数为998244353，依然没有多点插值和多点求值不过这次画风大变，因为结构体化的多项式模板其实无论用什么实现方式（数组，指针，vector）都不能同时兼顾效率和灵活性两方面，所以整个代码都是用数组实现稍微写一下一些需要注意的点：多项式求逆大概是基于一个思想:倍增首先当deg(A)=1deg(A)=1deg(A)=...

原创 FWT模板

.稍微学了一下FWT发现挺好写这个东西和FFT其实区别比较大，所以比较适合理性理解比如对于or运算的FWT，我们考虑这样一个式子FWTor(A)[i]=∑j∣i=iA[j]FWT_{or}(A)[i]=\sum_{j|i=i}A[j]FWTor​(A)[i]=j∣i=i∑​A[j]以及我们要求的东西(A∣B)[k]=∑i∣j=kA[i]∗B[j](A|B)[k]=\sum_{i|j=k}...

转载 有没有比图灵机能力更强的计算模型？

有，而且还不少。他们被称为超计算（Hyper computation）模型。超计算，是一个研究比图灵机计算能力更强的计算能力的计算机器的理论计算机科学分支。主要有以下部分模型： A.谕示机. （Oracle Machine） 带“黑箱”的图灵机。由图灵本人亲自提出，“黑箱”就是一个谕示，经过一个谕示就可以得到一个问题的判定结果。所有 Hypercomputation 的“原型机...

原创 音乐测试

原创 几种排序算法的对比

!最近干了一些有意思的事情:拿出来写一写关于各种sort!比如，c++自带的快排，归并排序，还有stdlib库里面的qsort还有我自己实现的二进制基排，65536进制基数排序，普通归并排序，FFT式归并排序以及基于队列实现的queue-mergesort（其实本质还是mergesort啦）目前测试结果来看，表现最好的，应该是普通归并排序（65536基排不算），排序 10710^710...

原创 CSDN-markdown编辑器使用帮助

这里写自定义目录标题欢迎使用Markdown编辑器新的改变功能快捷键合理的创建标题，有助于目录的生成如何改变文本的样式插入链接与图片如何插入一段漂亮的代码片生成一个适合你的列表创建一个表格设定内容居中、居左、居右SmartyPants创建一个自定义列表如何创建一个注脚注释也是必不可少的KaTeX数学公式新的甘特图功能，丰富你的文章UML 图表FLowchart流程图导出与导入导出导入欢迎使用Ma...

原创 SPFA的两个优化

最近看到知乎上面有人讨论怎么卡spfa并给出了这道题作为实战地点 于是我就搞了两个优化把出台人精心构造的数据过了 优化1： 这个优化不是我想出来的，如果没错应该是Menci的做法 本来这个优化已经可以过掉原来的五个数据的，让后fstqwq（出题人）搞了一个新数据把这个优化卡T了 其实就是每次入队，如果队头比队尾大那么就交换（感觉很像SLF但是比SLF强） 让后为了Ac我就加了第二...

原创 状压DP十六连测

DP应该先学插头 状压DP十六连 A[寿司晚宴] 开始想到了反演，设f(d)f(d)f(d)表示两边选出来的公因数为d的方案数，设F(d)F(d)F(d)… 好的介绍一个套路，按照最大的质因子分类（没错这不是洲阁筛） 我们发现小于500−−−√500\sqrt{500}的质数只有8个，可以搞一个2162162^{16}的状态表示一种方案 让后，我们设T(x)T(x)T(x)表示x...

原创 高斯消元模板

高斯消元还是很好理解的，于是就找了个模板题做了下:Bzoj1013 不过好像用在那些dp题上面还是很不方便，感觉还是自己太菜了 这个题就很简单了，设出球心坐标O(O1,O2...On)O(O1,O2...On)O(O_1,O_2...O_n) 写出距离方程式 ∑ni=1(X0,i−Oi)2=∑ni=1(X1,i−Oi)2=...=∑ni=1(Xn,i−Oi)2∑i=1n(X0,i−Oi...

原创 CDQ分治&整体二分九连测

整体二分好难a CDQ分治&amp;amp;amp;amp;整体二分九连： A[适者] 先来一个不是cdq分治的题（虽然也可以用分治做） 首先先来考虑按照什么顺序来搞掉这些机甲 对于两个相邻的机甲i,j，先i后j的代价是Ai∗Di+(Di+Dj)∗AjAi∗Di+(Di+Dj)∗AjA_i*D_i+(D_i+D_j)*A_j 那么如果i应该比j先消灭就有 Ai∗Di+(Di+Dj)∗Aj&amp;amp;amp;lt;A...

原创 E文本加密器0.2 Aplha版

自己做的一个有趣的东西，可以将一段文字加密 目前暂不支持中文和非文本文件（目前主要就是加密code） 测试效果： 原文After 4 years' waiting, the game "Chinese Paladin 5" finally comes out. Tomato is a crazy fan, and luckily he got the first release....

原创 斜率DP十连测

最近学校里很多题目都没时间做，顺便来写一下博客 斜率DP十连： A[征途] 注意到，原题的式子，等价于m∑a2i−(∑ai)2m∑ai2−(∑ai)2m\sum a_i^2-(\sum a_i)^2后面部分是常数 那么我们就可以写出dp方程f[i,j]=min{f[i−1,k]+(sj−sk)2},k&amp;amp;amp;amp;amp;lt;jf[i,j]=min{f[i−1,k]+(sj−sk)2},k&amp;amp;amp;amp;amp;lt;...

转载 数学公式测试

Testing 癌症，而且是恶性程度很高的哪一种。 癌症，而且是恶性程度很高的哪一种。 一开始，全家人都不相信，妈妈带着我把本地所有三甲医院都转遍了。一开始，全家人都不相信，妈妈带着我把本地所有三甲医院都转遍了。 确诊，我爸说，瞬间天昏地暗，站不稳。确诊，我爸说，瞬间天昏地暗，站不稳。 手术前一天，妈妈做梦去世的奶奶要带我走，梦里哭着喊着把我拉了回来。手术前一天，妈妈做梦去世的奶奶要...

原创 51Nod1038 X^A Mod P

题目看这里 经典的n次剩余问题，用到很多数论知识点 1.扩展gcd 2.原根 3.离散对数 4.n次剩余 说一下这个算法的流程 首先，我们的方程为xn=a(mod&amp;amp;nbsp;m)xn=a(mod&amp;amp;nbsp;m)x^n=a(mod\ m) m为质数 那么，我们首先要找m的原根g，这里g要满足的性 质就是对于任意&amp;amp;lt;fontsize=8&amp;amp;gt;&amp;amp;lt

原创 51Nod2026 Gcd and Lcm

题目看这里 一个非常好的题！ 好的，看到题目就很懵逼 首先这个f不就是ϕϕ\phi吗，认真一看才发现不对 让后问题？f(lcm)*f(gcd)? 肯定有问题，推了一会没有结论，去看看题解： 有这么一个神奇结论f(gcd(x,y))∗f(lcm(x,y))=f(x)∗f(y)f(gcd(x,y))∗f(lcm(x,y))=f(x)∗f(y)f(gcd(x,y))*f(l...

原创 51Nod1192 Gcd表中的质数

题目看这里 又到了推式子的时候了，莫比乌斯反演入门题 ∑i=1n∑j=1m[prime(i,j)]∑i=1n∑j=1m[prime(i,j)]\sum_{i=1}^n\sum_{j=1}^m[prime(i,j)] 首先搞一个f(d)f(d)f(d)表示有多少对(i,j)的gcd就是d 那么原式=∑nd=1[prime(d)]∗f(d)=∑d=1n[prime(d)]∗f(d)=...