Rose_max的博客

前言在联合省选day2t3中，存在一种使用行列式求导来计算生成树边权和的方法。需要计算出每个位置的代数余子式常见的做法均是套用A×A∗=∣A∣×IA\times A^*=|A|\times IA×A∗=∣A∣×I的等式求解，这种做法在模意义下矩阵rank(A)=n−1rank(A)=n-1rank(A)=n−1时不能正确得出伴随矩阵以下给出正确的求解伴随矩阵的方法代数余子式给nnn阶方阵A=(ai,j)A=(a_{i,j})A=(ai,j​)，定义ai,ja_{i,j}ai,j​的余子式Mi,jM

是拿来给自己抄的哼struct Bignum{ int a[MAXN],len; Bignum(){memset(a,0,sizeof(a));len=0;} friend Bignum operator +(Bignum u1,Bignum u2) { Bignum ret;int len=max(u1.len,u2.len); for(int i=1;i<=len;i...

前言我抄我自己基础的NIM取石子游戏结论：当每堆石子数a1 xor a2 xor a3 xor…xor an>0时，先手必胜，否则先手必败证明：当每堆石子均为0时，必定为必败态不做考虑当石子数异或和不为0时，一定能从某一堆石子中取出相对应的一些石子使得异或和为0。设当前异或和为x，我们从某一堆中取出一些石子使其成为a[i] xor x颗，此时异或和为0当石子数异或和等于0时，无论...

前言抄抄抄%%%干什么的听说是爆踩洲阁筛的东西？？可以用来筛一类积性函数f(x)f(x)f(x)的前缀和必须满足f(p)f(p)f(p)是一个简单多项式并且f(pc)f(p^c)f(pc)能快速计算其中p∈primep\in primep∈prime复杂度空间是n12n^{\frac{1}{2}}n21​时间据说是n34logn\frac{n^\frac{3}{4}}{logn...

前言方便自己看的…做法一般我们可以有一类dp方程譬如f[x]=∑f[y]∗cal[z]f[x]=\sum f[y]*cal[z]f[x]=∑f[y]∗cal[z]，其中cal[z]cal[z]cal[z]代表需要的代价满足y⋃z=xy\bigcup z=xy⋃z=x且y⋂z=∅y\bigcap z=\emptysety⋂z=∅那么朴素做法就是先枚举子集，然后转移，复杂度是3n3^n3n...

前言给我也来一个玩法为啥人家的博客都要放个manachar…感觉这两个没啥关系啊qwq每个状态xxx代表了一个长度为len[x]len[x]len[x]的回文串我们保证他代表的这个回文串是他能代表的里面最大的那个这是一个自动机不是一个树qwq…，他的fail才是一个树son[x][i]son[x][i]son[x][i]表示在状态xxx后加一个字符iii能转移到哪个回文串的状态l...

前言补坑补坑还债还债参考了liu_runda的前置上下界：即一条边存在必须流过的流量与至多流过的流量，分别称为上界和下界流量守恒：即满足一个点的流入流量=流出流量可行流：问一张图是否存在每个点满足流量守恒且满足上下界限制的流将以上条件称为基本条件称一条边(u,v,down,up)(u,v,down,up)(u,v,down,up)为从uuu到vvv，流量下界为downdowndo...

前言博主这个暴力骗分选手get到了人生的本质前置有一类函数，我们要求其的最低点/最高点二分？三分？四五六七八九分？？哦凉凉了…先介绍一个爬山我们随机撒点，然后让这些点去做类似现实的爬山即右边优就去右边，左边优就去左边直到没有比他更优秀的了！也就是相当于到山峰了吧那么这个贪心的算法显然非常容易卡在一个局部最优解出不来了哦怎么办呢…按正常人的思维想，走到了一个最高处，就应该...

前言头都给队友们打烂了啊这玩意还是简单易懂的啊qwq似乎博客已经变成了笔记博客？？用途求一类积性函数的前缀和经典问题有求μ\muμ与ϕ\phiϕ的前缀和，本文将以这两个函数的前缀和为例前置狄利克雷卷积两个函数fff和ggg的狄利克雷卷积f∗gf*gf∗g即为(f∗g)(i)=∑d∣if(i)g(id)(f*g)(i)=\sum_{d|i}f(i)g(\frac{i}{d})(f...

前言我是被逼的不打还真不会打模版题干啥的树上最大独立集？动态改点权？NOIP2018D2T3?…玩法以模版题为例先写出dpdpdp方程，常用的fi,0/1f_{i,0/1}fi,0/1​表示当前点选了/没选的最大方案转移有fi,0=∑max(fson,0,fson,1)f_{i,0}=\sum max(f_{son,0},f_{son,1})fi,0​=∑max(fson...

多项式求逆给你f(x)f(x)f(x)，要求计算g(x)g(x)g(x)满足f(x)g(x)≡1(mod  xn)f(x)g(x)\equiv 1(\mod x^n)f(x)g(x)≡1(modxn)考虑倍增的过程，常数项就直接求逆，否则已经计算完n2\frac{n}{2}2n​的贡献若已知f(x)h(x)≡1(mod&Thin...

搬运自@_rqy不定期更新

符号与约定s[1...x]=s[1]s[2]s[3]...s[x]s[1...x]=s[1]s[2]s[3]...s[x]s[1...x]=s[1]s[2]s[3]...s[x]组成的字符串pre[s,x]=s[1...x],suf[s,x]=s[∣S∣−x+1...∣S∣]pre[s,x]=s[1...x],suf[s,x]=s[|S|-x+1...|S|]pre[s,x]=s[1...x]...

就是个板子。又短又爽我爱了。namespace LCT{ int ch[2][MAXM],fa[MAXM],rev[MAXM]; int sta[MAXM],tp; bool son(int x){return x==ch[1][fa[x]];} bool isroot(int x){return ch[0][fa[x]]!=x&&ch[1][fa[x]]!=x;}...

第一类斯特林数ppp个不同人围着mmm个不同圆桌坐，要求每桌非空，方案数即为S(p,k)S(p,k)S(p,k)递推边界S(p,p)=1(p>=0),S(p,0)=0(p>=1)S(p,p)=1(p>=0),S(p,0)=0(p>=1)S(p,p)=1(p&amp

前言胡乱学了一点点…补在这里可能会慢慢更新吧…定义F(a)=F(a−1)+F(a−2)F(a)=F(a-1)+F(a-2)F(a)=F(a−1)+F(a−2)其中F(0)=0,F(1)=1F(0)=0,F(1)=1F(0)=0,F(1)=1邻项公约数为1gcd(F(a),F(a+1))=1gcd(F(a),F(a+1))=1gcd(F(a),F(a+1))=1数学归纳法一下，首先...

学了一发字典树 发现其实就是一个储存数据结构的东西。。可以省空间 首先有一个根，这个根无实际意义 根往下到叶子节点的每条路径，储存的都是一个串 这样就显而易见了对吧 由于许多个字符串开头有可能相同的，那我们就可以用一个类似树型结构的东西存啦 类似下面这张图 那么abcd就储存在了从左往右的第一条链 abd就存在了第二条链 以此类推…… 查找就更加简单啦。看看对于每个节

update一个原文档的链接方便大家看看 快速傅里叶变换.docxFast Fourier Transformation ——By Rose_max 简单来说，傅里叶变换，在oi里面就一个用途：加速多项式乘法 方法就一个：构造多项式fft点值乘法ifft写在前面 关于学习FFT算法

主要参考的是这篇blog 以及这篇俄文翻译 upd:一些性质参考的是这篇blog，Menci超劲啊！ 因为clj的论文太难懂了其他blog又是特别抽象的对我这种蒟蒻特别不友好 这篇blog是方便自己复习写的所以难免有引用其他blog的地方一些定义定义ch[p][s]表示在状态p上增加一个字母s后能转移到的状态，如果不能转移则为0 定义str(s)表示在原串中走到s这个位置所构...

前言KOI里发现一个我连暴力都不会打的题可以NTT水80分 于是我就下定决心学一手NTT，顺便复习FFT 主要都是推多项式啦一些关于NTT的东西其实NTT也很简单，和FFT差不多 为什么需要NTT呢，因为FFT带了三角函数和浮点数运算，精度误差有时候会比较大 而且对于取模的时候，需要搞个任意膜（膜myy），精度误差更大 所以就出现了NTT关联在FFT中，我们有一个欧拉...

前言之前学了一遍 发现只会无向图生成树 于是就在这里总结一下各种结论吧定义有向树：对于一个有向图，如果无视边的方向是一棵树，那么此有向图就称为有向树 外向树：有向树的特殊情况，下同，所有边的方向都是从根指向叶子 内向树：所有边的方向都是从叶子指向根求法无向图定义一张图的基尔霍夫矩阵为 两个点i,j，如果i=j，那么他在基尔霍夫矩阵中的权为i点在图中的...

N数码判定性问题大意就是给你两个N数码（N为奇数），判定是否能互相转化的问题 首先展开成一维，空格不计 如 1 2 3 4 5 6 7 0 8转化为1 2 3 4 5 6 7 8 然后判断两个序列的逆序对数奇偶性是否相同，相同则有解 证明： 在空格左右移动的时候，展开序列不改变 在空格上下移动的时候，展开序列中有一个数会前移（后移）n-1位，由于n为奇数则n-1为偶数。于是改...

前言补了这个辣鸡算法吧..仙人掌定义：一张图中任意一条边最多只会出现在一个环中DFS树的做法对图进行dfs，记录dfs的编号形成一棵树 显然最后边会分为两类 树边和非树边 环一定是由一些树边和一条非树边构成的 对于仙人掌问题，可以考虑对dfs树进行dp 在遇到环的时候，把环单独拿出来考虑即可…圆方树为了解决大多数静态仙人掌问题，我们就引入了这个算法！...

前言老年颓废选手深感无力…水题一堆堆真是痛苦点点技能树愉悦身心定义有一个多项式A(x)A(x)A(x)，现在要你找一个多项式B(x)B(x)B(x)使得A(x)B(x)≡1(mod  xn)A(x)B(x)\equiv 1(\mod x^n)A(x)B(x)≡1(modxn)mod &Thin...

胡说八道啥都不会啥都想学大概废了拉格朗日插值给你nnn个点(xi,yi)(x_i,y_i)(xi​,yi​)，求经过这些点的n+1n+1n+1次函数及求某个xxx在该函数上的取值由于nnn个点组成的函数是固定的所以你可以考虑构造一种函数，使得当x=xix=x_ix=xi​时取值y=yiy=y_iy=yi​，显然这种函数就是解可以写成∑i=1nΠj!=i(x−xj)yiΠj!=i(...

只是来拖一个板子的…void FWT_or(int *a,int opt){ for(int i=1;i<N;i<<=1) for(int p=i<<1,j=0;j<N;j+=p) for(int k=0;k<i;++k) if(opt==1)a[i+j+k]=

胡思乱想最后的沉沦？？用处求CnmmodpC_n^m mod pCnm​modp的结果，不保证ppp为质数水法将ppp分解质因数，可以知道p=p1k1∗p2k2∗...∗pnknp=p1^{k1}*p2^{k2}*...*pn^{kn}p=p1k1∗p2k2∗...∗pnkn如果我们能知道ans≡a1(mod  p1k1)an...

caioj1178 【题目描述】 给出模板串A和子串B，长度分别为lenA和lenB，要求在线性时间内，对于每个A[i]（1<=i<=lenA)，求出A[i..lenA]与B的最长公共前缀长度 【输入格式】 输入A，B两个串，（lenB<=lenA<=1000000) 【输出格式】 输出lenA个数，表示A[i…lenA]与B的最长公共前缀长度，每个数之前有空格 【样例输入】 aab