author: cust-- ZKe --------------------- 这里以连乘积加括号问题为背景: 由于矩阵的乘积满足结合律,且矩阵乘积必须满足左边矩阵的列数的等于右边矩阵的行数,不同的计算顺序,需要的乘法运算次数不一样.加括号可以改变计算顺序,合理安排计算顺序可以大大降低计算次数. 给乘积算式加括号的方法数是一个计数问题.它的模型是卡特兰数. 比如有矩阵A,B,C,D,有五种加括号方式 ((A*B)*C)*D (A*(B*C))*D (A*B)*(C*D) A*(B*(C*D))

前记 最近,公司一位挺优秀的总务离职,欢送宴上,她对我说“你是一位挺优秀的程序员”,刚说完,立马道歉说“对不起,我说你是程序员是不是侮辱你了?”我挺诧异,程序员现在是很低端,很被人瞧不起的工作吗?或许现在连卖盗版光盘的,修电脑的都称自己为搞IT的,普通人可能已经分不清搞IT的到底是做什么的了.其实我想说,程序员也分很多种的,有些只能写if-then-else,有些只能依葫芦画瓢,但真正的程序员我想肯定是某个领域的专家,或许他是一位数学家,或许他是一位物理学家,再或许他是计算机某个细分领域的专家,

思路来源:http://blog.csdn.net/u013654696/article/details/23037407#comments [做浙大校赛的时候没有看这道题,事后做的.思路不是自己的,但代码是自己敲的,由于伦家不懂如何用TeX敲出如此优美的公式,所以具体请看上面的博客链接(づ￣3￣)づ╭.虽然说思路对应下的代码很好敲,但如果在比赛中我肯定不一定想得到这么做.在这道题中,线段树的节点与区间存的不是具体的值,是对应的一个矩阵.做了这道题也让我明白了线段树原来还可以那么好用吖(づ￣3￣

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6470 题意:求,直接矩阵快速幂得f(n)即可 构造矩阵如下: n^3是肯定得变换的,用二项式展开来一点点配吧 我们会发现中间6*6的矩阵是个常数矩阵,则可以化为A=B^(N-2)*C(n-2次幂是因为我们求解是从N=3开始的),根据矩阵快速幂算出B^(N-2)次幂即可以了 矩阵快速幂的时间复杂度是logn 注意:我们中间构造的矩阵必须是一个方阵,矩阵快速幂的难点就在于构建中间的方阵 比如矩阵A,B,

题目描述 一张N个点M条边的无向图,从A走到B,要求:每一次不能立刻沿着上一次的边的反方向返回.求方案数. 输入 第一行:五个整数N,M,t,A,B. N表示学校里的路口的个数 M表示学校里的路的条数 t表示HH想要散步的距离 A表示散步的出发点 B则表示散步的终点. 接下来M行 每行一组Ai,Bi,表示从路口Ai到路口Bi有一条路. 数据保证Ai != Bi,但不保证任意两个路口之间至多只有一条路相连接.  路口编号从0到N -1.  同一行内所有数据均由一个空格隔开,行首行尾没有多余空格.没

题意: 冈伦从二维平面上(0,0)走到(k,0),(k<=1e18),每次有三个行动方向:右上一格.右方一格.右下一格,问一共有多少种走的方案 限制:每段x都有一个天花板,一共有n段天花板(n<=100),每段天花板都有一个高度(<=15),不能走得比该段天花板高 分析: 将其分成n段 发现就是执行n段矩阵快速幂,这n段矩阵快速幂的转移矩阵不同,转移次数也不同 首先预处理出每段天花板的转移次数(长度) 为了方便实现,建个结构体mat将矩阵的初始化.转移矩阵的构造以及矩阵乘法全部写入结构体

手动博客搬家: 本文发表于20181212 09:37:21, 原地址https://blog.csdn.net/suncongbo/article/details/84962727 呜啊怎么又是数学了啊...数学比例\(\frac{16}{33}=0.4848\) orz yhx-12243神仙 题目链接: https://codeforces.com/contest/947/problem/E 题意: 有一个\([0,n]\)的随机数\(x\)初始为\(i\)的概率为\(p_i\). \(m

题目链接:http://acm.hdu.edu.cn/showproblem.php?pid=6395 因为题目数据范围太大,又存在递推关系,用矩阵快速幂来加快递推. 每一项递推时  加的下取整的数随着n变化,但因为下取整有连续性(n一段区间下取整的数是相同的),可以分块,相同的用矩阵快速幂加速 想了好久..如果最小的开始的值是[p/i]的数为i,那连续的一段长度是[p/(p/i)]-i+1,但为什么分段数是根号n级别啊?... 套矩阵快速幂,时间复杂度O(sqrt(n) * log(n)) ⎧

公交线路 Description 小Z所在的城市有N个公交车站,排列在一条长(N-1)km的直线上,从左到右依次编号为1到N,相邻公交车站间的距离均为1km. 作为公交车线路的规划者,小Z调查了市民的需求,决定按下述规则设计线路: 1.设共K辆公交车,则1到K号站作为始发站,N-K+1到N号台作为终点站. 2.每个车站必须被一辆且仅一辆公交车经过(始发站和终点站也算被经过). 3.公交车只能从编号较小的站台驶往编号较大的站台. 4.一辆公交车经过的相邻两个 站台间距离不得超过Pkm. 在最终设计

明明是水题结果没切掉……降智了…… 首先令 $c$ 为序列中 $0$ 的个数,那么排序后序列肯定是前面 $c$ 个 $0$,后面 $n-c$ 个 $1$. 那么就能上 DP 了.(居然卡在这里……) $f[i][j]$ 表示经过 $i$ 次操作后,前 $c$ 个数中有 $j$ 个 $0$ 的方案数.答案就是 $\dfrac{f[k][c]}{\sum f[k][i]}$. 这个状态的好处就是可以直接求出以下这些值: 前 $c$ 个数中 $1$ 的个数为 $c-j$ 后 $c$ 个数中 $0$ 的

原文链接www.cnblogs.com/zhouzhendong/p/Cayley-Hamilton.html Cayley-Hamilton定理与矩阵快速幂优化.常系数线性递推优化 引入 在开始本文之前,我们先用一个例题作为引入. 给定一个 \(n \times n\) 的矩阵 \(M\) , 求 \(M ^ k\) . \(n\leq 50, k\leq 10 ^ {50000}\) . 注意到 \(n\) 十分小,但是 $ \log k$ 非常大.如果使用传统的矩阵快速幂,时间复杂度为 \

题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=5118 题解 这个题一看就是不可做的样子. 求斐波那契数列的第 \(n\) 项,\(n \leq 2^{10^{15}}\)??? 这样人怎么矩阵快速幂啊. 等等这个模数很神奇啊. \(1125899839733759\) 好像是一个质数,还以 \(9\) 结尾. 那么 \(5\) 对于 \(1125899839733759\) 一定有二次剩余咯. 那么根据 Fib 的通项公式 \[ f(n)

传送门 解题思路 以前bpw讲过的一道题,顺便复习一下矩阵乘法.做法就是拆点,把每个点拆成\(9\)个点,然后挨个连边.之后若\(i\)与\(j\)之间的边长度为\(x\),就让\(i\)的第\(x\)个点和\(j\)的第\(1\)个点连边.然后就是一个矩阵快速幂,时间复杂度\(O((n*9)^3log(n*9))\). 代码 #include<iostream> #include<cstdio> #include<cstring> #include<algori

题目描述 分析 定义出\(dp[i][j]\)为第\(i\)列涂\(j\)种颜色的方案数 然后我们要解决几个问题 首先是求出某一列涂恰好\(i\)种颜色的方案数\(d[i]\) 如果没有限制必须涂\(i\)种,而是有的颜色可以不涂,那么方案数为\(i^n\) 为了避免少涂的情况,我们减去只涂\(1 \sim i-1\)种颜色的方案数 即\(d[i]=i^n-\sum_{j=1}^{i-1}C_i^j \times d[j]\) 初始化为\(d[1]=1\) 接下来考虑转移 \(f[i][j]=f

[ 矩 阵 乘 法 ] 裴 波 拉 契 数 列 I I [矩阵乘法]裴波拉契数列II [矩阵乘法]裴波拉契数列II Description 形如 1 1 2 3 5 8 13 21 34 55 89 144-的数列,求裴波拉契数列的第n项. Input n (1< n <2^31) Output 一个数为裴波拉契数列的第n项mod 10000; Sample Input 123456789 Sample Output 4514 题目解析 首先看题面,是斐波那契数列.首先想到递归,但考虑到N的值

Atcoder 题面传送门 & 洛谷题面传送门 这是一道难度 Cu 的 AGC E,碰到这种思维题我只能说:not for me,thx 然鹅似乎 ycx 把题看错了? 首先这个平方与乘法比较喜闻乐见,很容易与组合联系在一起,于是我们不妨把题目条件翻译成组合的语言: 有一排 \(n\) 个格子,你要在其中插入若干个隔板将其隔成若干段 有 \(m\) 个特殊格子 \(a_1,a_2,\dots,a_m\),\(\forall i\in [1,m]\) 你禁止在 \(a_i\) 与 \(a_{i}+

A,B,C傻逼题,就不说了. E题: #include <iostream> #include <cstdio> #include <cstring> #include <cmath> #include <algorithm> using namespace std; ,mod=; typedef long long int64; ]; int64 a[maxn],s1[maxn<<],s2[maxn<<]; struct

fibonacci 数列及其延展 fibonacci计算 fibonacci数列是指 0,1,1,2,3,5,8,13,21……这样自然数序列,即从第3项开始满足f(n)=f(n-1)+f(n-2): 递归实现非常简单: long long fibonacci(unsigned int n) { ] = {, }; ) return result[n]; ) + fibonacci(n-); } 以计算f(10)为例,必须先求得f(9)和f(8),要计算f(9),又必须先求得f(8)和f(7),

首先这道题是GT考试的加强版本QAQ 当n<k的时候,答案显然是10^n 当n=k的时候,答案显然是10^n-1 这样就有20分辣 之后我们考虑k<=20的做法 显然设f(i,j)表示前i位匹配到了第j位 之后用KMP来转移就可以了 如果不用矩阵乘法优化的话时间复杂度O(n*k*10) 如果用矩阵乘法的话时间复杂度O(k^3logn) 其实就是GT考试的做法辣 这样就有50分了 之后我们考虑用dp来容斥 设fi表示前i位且第一次匹配的位置是[i-k+1,i]的方案 不难发现fi=m^(i-k)

啦啦啦,今天的考试题 不过原来考试题的n<=10w 由于我有更好的做法,所以我就改成20亿辣 本来先说一说考试题的正解做法的 但是复杂度是O(nlogm),实在是太渣了 所以还是说一说我的做法吧 首先假定都会写裸的DP 我们考虑A,B,如果B不能转移到A,当且仅当A不等于B且A%B=0 很容易发现当A是素数时,A只不能从1那里转移过来,也就是说素数的转移都是一样的 换句话说,在状态里当长度一定时以每个素数结尾的方案一定是一样的 这就给了我们一些启发,很容易得到结论若两个数的唯一分解式的指数排序后

题目链接:http://poj.org/problem?id=3318 http://meizhe143.blog.163.com/blog/static/38938362007102995121360/     介绍随机数产生方法. http://blog.csdn.net/hanxiang_acm/article/details/3331113 上面这个单纯的随机产生数,比如产生50000,我试了多次结果都是WA. 正解应该是左乘或右乘一个一行的矩阵,原来的时间复杂度是O(n3) ,现在直接