自定义博客皮肤VIP专享

*博客头图:

格式为PNG、JPG,宽度*高度大于1920*100像素,不超过2MB,主视觉建议放在右侧,请参照线上博客头图

请上传大于1920*100像素的图片!

博客底图:

图片格式为PNG、JPG,不超过1MB,可上下左右平铺至整个背景

栏目图:

图片格式为PNG、JPG,图片宽度*高度为300*38像素,不超过0.5MB

主标题颜色:

RGB颜色,例如:#AFAFAF

Hover:

RGB颜色,例如:#AFAFAF

副标题颜色:

RGB颜色,例如:#AFAFAF

自定义博客皮肤

-+

原创 实Jordan标准形的推导

众所周知,任意复矩阵有惟一Jordan标准形。对于实矩阵,设其初等因子为p1s1(λ),...,pnsn(λ)p_1^{s_1}(\lambda),...,p_n^{s_n}(\lambda)p1s1​​(λ),...,pnsn​​(λ),影响构造Jordan块的只有形如(x2+px+q)n(x^...

2020-04-25 12:29:58 128 0

原创 拉格朗日乘子法

a brief proof to the Lagrange Multiplier Method

2020-04-14 16:06:26 146 0

原创 《国际最佳数学征解问题分析》P77(单位根反演+生成函数)

问题 给定x1,...,xnx_1,...,x_nx1​,...,xn​,令xi(1)=xi+xi+12x_i^{(1)} = \frac{x_i+x_{i+1}}{2}xi(1)​=2xi​+xi+1​​,i=1,...,ni=1,...,ni=1,...,n,其中xn+1=xix_{n+1}=...

2020-04-05 16:42:17 149 0

原创 Tietze扩张定理

一个延拓函数的方法

2020-04-05 16:36:00 288 0

原创 PCA & SVD & 极分解 简介

A brief intro to PCA and SVD decomposition.

2020-03-28 20:13:21 110 0

原创 Alpha-Beta搜索剪枝最优效率分析

A naive analysis of the best efficiency of Alpha-beta pruning.

2020-03-16 16:17:49 78 0

原创 哈密顿-凯莱(Hamilton-Cayley)定理的两种简要证明

Two brief proofs of Hamilton-Cayley theorem.

2020-03-11 09:48:11 1138 0

原创 《线性代数》 李炯生\查建国\王新茂 中国科学技术大学 第2版 部分习题答案

百度了一圈没有靠谱点的答案,于是便有了这篇博客。。。 文章目录第一章 多项式1.3 整除性与最大公因式习题3习题101.5 实系数与复系数多项式习题1习题2.(4)习题101.6 整系数与有理系数多项式例1例2习题2习题6习题7习题81.8 对称多项式习题3习题6习题8,10习题11. 第一章 多...

2019-10-31 22:37:03 2584 1

原创 Burnside引理的简要证明

A brief proof of Burnside's theorem

2019-09-09 21:58:02 675 1

原创 常系数齐次线性递推优化矩阵快速幂

一般矩阵快速幂的形式 : f(n)=∑i=1k1aif(n−i)+∑i=1k2big(n−i)+bf(n)=\sum_{i=1}^{k_1} a_if(n-i)+\sum_{i=1}^{k_2}b_ig(n-i)+b 可以做到k3lognk^3 \log n的常数递推,不过有更加快速的方法。若...

2019-08-11 13:16:02 2268 0

原创 集训队作业2018: 复读机(生成函数)

题意: 群里有kkk个不同的复读机。为了庆祝平安夜的到来,在接下来的nnn秒内,它们每秒钟都会选出一位优秀的复读机进行复读。非常滑稽的是,一个复读机只有总共复读了ddd的倍数次才会感到快乐。问有多少种不同的安排方式使得所有的复读机都感到快乐(k≤1000,d≤3)(k \le 1000, d \l...

2019-01-18 16:06:02 766 1

原创 集训队作业2018: 喂鸽子(min-max容斥)

题意: 小Z是养鸽子的人。一天,小Z给鸽子们喂玉米吃。一共有n只鸽子,小Z每秒会等概率选择一只鸽子并给他一粒玉米。一只鸽子饱了当且仅当它吃了的玉米粒数量≥k≥k≥k。 小Z想要你告诉他,期望多少秒之后所有的鸽子都饱了。 题解: min-max容斥枚举下集合大小,FFT预处理一下系数算min的期望即...

2019-01-15 16:09:53 545 0

原创 集训队作业2018: 取名字太难了(FFT)

题意: 大概是求∏i=1n(x+i)\prod_{i=1}^n(x+i)∏i=1n​(x+i)系数模ppp意义下的分布。 题解: 分为(∏i=1p−1(x+i))⌊np⌋(\prod_{i=1}^{p-1}(x+i))^{\lfloor\frac{n}{p}\rfloor}(∏i=1p−1​...

2019-01-14 19:23:49 575 0

原创 集训队作业2018: 树(点分治+K短路)

题解: 最近学数分学到意识模糊,做到OI题冷静一下。 联通块强制选根,然后用dfs序转化为一个路径然后就是做K短路了。 用点分治即可在图大小为O(nlog⁡n)O(n \log n)O(nlogn)的图上做K短路,时间复杂度O(nlog⁡2n+klog⁡k)O(n \log ^2 n + k\lo...

2019-01-12 16:50:06 393 2

原创 Codeforces 802C :Heidi and Library (hard)(网络流)

传送门 题解: 比较简单的建图法就是看做小于k条流在一个n∗nn*nn∗n序列上流,其中一些位置是必须流的,然后做个上下界费用流。 不过注意到肯定有一种方案使得这小于kkk条流是不相交的,于是可以直接看做有nnn个点,每个点拆点连−∞-\infty−∞的边,然后规定这个点必须是aia_iai...

2018-11-30 17:33:51 264 0

原创 Codeforces 806F:Test Data Generation(组合数学)

传送门 题解: 相当于是要求: ∑u∑i=1⌊n2u⌋[i为奇数]∑j=1n−1[j为偶数](i−1j)\sum_{u}\sum_{i=1}^{\lfloor\frac{n}{2^u}\rfloor}[i为奇数] \sum_{j=1}^{n-1}[j为偶数] \binom{i-1}{j}u∑​i=...

2018-11-30 15:43:58 213 0

原创 UOJ#433. 【集训队作业2018】串串(循环串/回文串)

传送门 题解: 这道题主要用到的几个性质(具体证明可以看题解): 1.弱双回文串SSS的某个双回文划分ababab,满足aaa是SSS的最长回文前缀,或者bbb是SSS的最长回文后缀。 2.弱双回文串SSS若有两个弱回文划分,则SSS为整周期串。 3.弱双回文串SSS的周期为ttt,则其有∣S∣t...

2018-11-27 20:26:49 735 1

原创 集训队作业2018: 通信(区间DP)

题意: 题解: 想到了大概思路就是不会有相交的非链上的边,然而有一堆细节就鸽了。 然后看题解其实就是暂时忽略两边的值来区间DP就行了,设置的状态还挺妙的,中间记个前缀后缀min优化一下,时间复杂度O(n3)O(n^3)O(n3)。 #include <bits/stdc++.h&a...

2018-11-27 11:36:00 203 0

原创 UOJ#429. 【集训队作业2018】串串划分(循环串)

传送门 题解: 考虑设dpidp_idpi​表示以iii结尾的前缀的划分方案数,因为有2条件的限制,可以得到容斥式子: dpi=∑j(−1)C(Sj+1,i)dpjdp_i = \sum_{j}(-1)^{C(S_{j+1,i})}dp_jdpi​=j∑​(−1)C(Sj+1,i​)dpj​ C(...

2018-11-24 22:41:54 599 0

原创 Atcoder AGC019简要题解

传送门 Reverse and Compare 发现两个串如果翻转后一样且中心不一样,则必定有一个是回文串。然后就只用统计两端字母不同的串的个数了。 #include <bits/stdc++.h> using namespace std; typedef long l...

2018-11-22 09:10:59 202 0

原创 UOJ#428. 【集训队作业2018】普通的计数题(牛顿迭代)

传送门 题解: 把0操作看做是叶子,1操作看做非叶节点,一个操作在另一个操作删除,则另一个操作为这个操作的父亲,于是转化成了满足以下条件的nnn个点的树的计数: 1.父亲标号>儿子。 2.若一个点为非叶节点,记其儿子中叶子节点的数量为TTT,则若其儿子...

2018-11-21 12:00:42 777 0

原创 HDU6172:Array Challenge(Berlekamp-Massey算法)

传送门 题解: 大胆猜想有递推式,然后BM立水之。 #include <bits/stdc++.h> using namespace std; const int RLEN=1<<18|1; inline c...

2018-11-19 18:45:14 151 0

原创 Codechef:Walk on Tree/TREEWALK(Berlekamp-Massey算法)

传送门 题解: O(n3log⁡k)O(n^3 \log k)O(n3logk)的优化应该都会。 然后O(n2)O(n^2)O(n2)用BM求出递推式之后再O(n2log⁡k)O(n^2 \log k)O(n2logk)用特征多项式优化一下就行了。 #include &am...

2018-11-19 18:03:40 237 0

原创 BZOJ1494: [NOI2007]生成树计数(Berlekamp-Massey算法)

传送门 题解: 直接打表+BM算出递推式,BM具体实现可以戳这里 附上一份其丑无比的BM代码: const int L=4e2; namespace bm { int cnt,a[N],fail[N],delta[N]; vector <int&...

2018-11-17 11:39:03 291 0

原创 集训队作业2018:GAME(并查集)

题意: 题解: 把这个DP式子给列出来: fi=si+max⁡j{2∗aj−sj−fj}f_i = s_i + \max_j \{ 2*a_j-s_j- f_j\}fi​=si​+jmax​{2∗aj​−sj​−fj​} 把这个后缀max⁡\maxmax记为mmm的话,每次就是m=max⁡{m,...

2018-11-16 16:12:07 461 0

原创 集训队作业2018: 青春猪头少年不会梦到兔女郎学姐(多限制容斥)

前言: 虽然这道题的名字有点那啥,但是题还是很好的,听说是某道原题的加强版。 题意: 给定nnn种颜色的球,第iii种颜色的球数量为AiA_iAi​个,保证∑i=1nAi≤2∗105\sum_{i=1}^n A_i \le 2*10^5∑i=1n​Ai​≤2∗105,对于这所有(∑Ai)!∏...

2018-11-16 14:19:53 1688 1

原创 Codeforces 809E: Surprise me!(Mobius反演)

传送门 题解: 对于每个iii,处理出:fi=∑a∑b[(vala,valb)==i]φ(vala)φ(valb)dis(a,b)f_i=\sum_a\sum_b[(val_a,val_b)==i]\varphi(val_a)\varphi(val_b)dis(a,b)fi​=a∑​b∑​[(va...

2018-11-13 22:02:47 206 0

原创 Codechef:Binary Tree/COOK82E(Trie)

传送门 题解: 维护每个时刻完整的二叉树信息即可,具体可以用trie树实现,时间复杂度O(n)O(n)O(n)。 #include <bits/stdc++.h> using namespace std; const int RLEN=1<&lt...

2018-11-13 19:34:51 128 0

原创 Codechef:Painting Tree/KILLER(二进制分组)

传送门 题解: 这道题,一开始想到自然是线段树合并了,每个点维护个fi,jf_{i,j}fi,j​表示从jjj开始上到iii,其他内部配对的最小值。不过这样就要支持区间加二次函数求最大值,根据bzoj某道题的经验,显然是不能在线段树上搞的。 然后仔细观察一下,发现这个(h−depx)(h-dep_...

2018-11-13 16:38:36 215 1

原创 Topcoder SRM 710 900pts:Hyperboxes(FMT)

题解: 对于一维是否相交用2(m2)2^{\binom{m}{2}}2(2m​)来表示一下。 然后多维直接FMT做并卷积即可,不过对于一维初始化就需要大力删去重复状态来剪枝了。 #include <bits/stdc++.h> using namespace std; ...

2018-11-13 08:54:47 159 0

原创 NOIP2018:保卫王国(LCT)

传送门 题解: 考场上刚T2导致这道题没写。。。 其实挺傻逼的。 LCT的时候维护一个2*2的矩阵就行了,顺便维护一下虚子树的信息,时间复杂度O(nlog⁡n)O(n \log n)O(nlogn)。 #include <bits/stdc++.h> using nam...

2018-11-12 11:31:21 1617 1

原创 51nod1819 :黑白树 V2(树链剖分)

传送门 题解: NOIP之前当然要做做NOIP题啦! 这道题并不难(想),把一个点的权值定为(u−fau)(u-fa_u)(u−fau​)乘上子树黑点个数,然后就相当于是支持查询链和,链修改,子树修改了,用了链剖分开维护奇偶的情况即可。 #include <bits/...

2018-11-09 18:06:30 212 0

原创 Topcoder SRM 702 1000pts:FindingFriends(分治)

题解: 这道题很妙啊。 显然是要二分之后每个位置找到前后第一个合法位置,分别记为Li,RiL_i,R_iLi​,Ri​,然后要求包含任意一个。 这个时候我们只需要找到1∼n1\sim n1∼n中间第一个非法的位置,然后递归下去做就好了,问题是怎么找到这个非法位置。 显然这是一个区间内的二维偏序,也...

2018-11-09 11:19:18 168 0

原创 Topcoder SRM 700 1000pts:AnyNumber(DP)

题解: 首先放满一行之后我们强制规定他可以继续放,这样总方案数是A(∑∣S∣,n)A(\sum |S|,n)A(∑∣S∣,n),一个方案可能被算很多次。 处理出fi,jf_{i,j}fi,j​表示前iii个第一次放满行为jjj的概率,hi,jh_{i,j}hi,j​表示前iii个种选一些数放满第j...

2018-11-09 09:13:17 130 0

原创 THUPC2017 I :Sum(牛顿恒等式)

题意: 给定数组A1...,AnA_1...,A_nA1​...,An​,对于所有1≤i≤k1 \le i \le k1≤i≤k,求Si=∑jAjiS_i = \sum_{j}A_j^iSi​=∑j​Aji​。 题解: 这道题要用到一个叫牛顿恒等式的玩意儿。 对于nnn次多项式f=∑i=0naix...

2018-11-08 10:05:57 1193 0

原创 计蒜之道2017复赛:商汤智能机器人(组合数学)

传送门 题解: 日常划水。。 枚举一下步数,其实是要求∑t=0b(a+tb)(bt)\sum_{t=0}^b \binom{a+t}{b}\binom{b}{t}∑t=0b​(ba+t​)(tb​)。 不妨记其为S(a,b)S(a,b)S(a,b),发现模数很小,于是可以根据lucas定理...

2018-11-07 22:22:17 119 0

原创 集训队作业2018:Z-function(DP)

传送门 题解: 枚举相同位置的长度是多少,然后可以设计一个DP,fi,jf_{i,j}fi,j​表示第iii位,相同状态为jjj的方案数(注意这里要带个-1的系数方便容斥),然后发现这个状态数很少,就可以过了。 #include <bits/stdc++.h> usi...

2018-11-07 19:10:56 344 0

原创 UOJ#273. 【清华集训2016】你的生命已如风中残烛(组合数学)

传送门 题解: 首先所有位置先-1,然后考虑m!m!m!种排列,如果全部后缀和小等于0(前缀和大等于0)那么是个合法排列,否则不合法。 这个时候有个自然的想法就是把小于0的最小的位置(若有多个则选最靠后的)放到这个排列后面,然后就又对应着一个合法的排列了,不过这样做每个排列对应的方案可能不一样,这...

2018-11-07 12:08:33 288 0

原创 UOJ#420. 【集训队作业2018】矩形(组合数学)

传送门 题解: 这种题再也不想做第二次了,思想很简单,但是细节贼多。 考虑每个aia_iai​(就是题面中的fif_ifi​)的贡献,显然是∑j=in(bhm)j∑k(j+k−ij−i)(ah)k\sum_{j=i}^n (bh^m)^j\sum_k\binom{j+k-i}{j-i} ...

2018-11-06 15:53:39 610 0

原创 UOJ#418. 【集训队作业2018】三角形(线段树合并)

传送门 题解: 加入一个数,相当于是先加上Ai=wiA_i=w_iAi​=wi​,再减去Bi=∑j∈soniwjB_i = \sum_{j \in son_i} w_jBi​=∑j∈soni​​wj​,然后代价就是一个操作序列的前缀最大值。 先考虑一下没有限制的的时候,怎么使得这个前缀最...

2018-11-05 22:45:22 470 0

提示
确定要删除当前文章?
取消 删除