2019年11月_「已注销」

原创数学杂记

记录一些琐碎的数学知识。以下内容重点放在记结论，以便作者复习。慢慢补充吧。最大公约数由公约数的基本性质可得： inline int gcd(int a,int b){return b?gcd(b,a%b):a;} 扩展欧几里得求解关于x,yx,yx,y的二元一次方程：ax+by=gcd(a,b)ax+by=gcd(a,b)ax+by=gcd(a,b) 当等式右端不是gcd(a,b)gcd(...

2019-11-11 21:09:56 260

原创基础向模板收集(字符串篇)

KMP 模板在循环中，iii表示当前要去匹配第iii位，jjj表示TTT的前jjj位和SSS的前i−1i-1i−1位可以匹配上。现在就看T[j+1]T[j+1]T[j+1]和S[i]S[i]S[i]是否匹配。读入需要判断EOF。 #include<bits/stdc++.h> #define VI std::vector<int> #define ll long lon...

2019-11-15 09:40:23 163

原创基础向模板收集(图论篇)

最短路计数要记录visvisvis数组避免算重。 if(dis[v]>dis[u]+W[i]) dis[v]=dis[u]+W[i],num[v]=num[u],Q.push(mp(-dis[v],v)); else if(dis[v]==dis[u]+W[i]) num[v]+=num[u]; 这里必须加elseelseelse，因为上面那个算完之后下面那个就成立了。模板 #incl...

2019-11-11 22:02:03 192

原创 HDU6172 Array Challenge

传送门佛辣。推出前几项后可以用BMBMBM算法得出⌊an⌋\lfloor \sqrt{a_n} \rfloor⌊an⌋的最短线性递推式。注意是直接算出来不取模，然后开根，并不是CipollaCipollaCipolla。嗯就是这个东西。。很神奇。它居然可以线性递推。所以要勇敢大胆打表蒙结论啊。 #include<bits/stdc++.h> #define ll long...

2019-11-08 20:24:31 168

原创 bzoj4162 shlw loves matrix II

传送门对于一个k∗kk*kk∗k的矩阵AAA，它的特征多项式是一个kkk次多项式。下面我们把多项式的自变量xxx看做一个矩阵。首先我们假设已经求出了特征多项式f(x)f(x)f(x) 特征多项式满足：f(A)=0f(A)=0f(A)=0。(Cayley-Hamilton定理) 我们令G(x)=xnG(x)=x^nG(x)=xn，那么AnA^nAn即为G(A)G(A)G(A)，也就是要求的东西...

2019-11-07 22:06:26 202

原创 bzoj4161 Shlw loves matrixI

传送门常系数齐次线性递推学习学习了任意模数多项式乘除法O(n2)O(n^2)O(n2)。除法就是手动模拟大除法。需要的东西就是一个转移矩阵AAA的特征多项式。要把AnA^{n}An看做是一个只有一项的nnn次多项式。然后知道初始的那几个就可以推出来了。 #include<bits/stdc++.h> #define ll long long #define re regist...

2019-11-07 11:35:36 151

原创 SDOI2017 文本校正

传送门分为6种情况。 ABC:ABC:ABC:哈希O(1)O(1)O(1)判断 CAB:CAB:CAB:枚举CCC的位置，哈希判断ABABAB，O(n)O(n)O(n) BCA:BCA:BCA:枚举AAA的位置，哈希判断BCBCBC，O(n)O(n)O(n) CBA:CBA:CBA:金策字符串算法选讲如果s=abs=abs=ab，a,ba,ba,b都是回文串，称sss为一个双回文串。引理：如...

2019-11-01 19:19:52 665

g21wcr的博客