矩阵乘法
文章平均质量分 94
Coco_T_
东南大学本科毕业生,浙江大学在读研究生
展开
-
矩阵快速幂模板
题目描述给定n*n的矩阵A,求A^k输入输出格式输入格式: 第一行,n,k 第2至n+1行,每行n个数,第i+1行第j个数表示矩阵第i行第j列的元素输出格式: 输出A^k 共n行,每行n个数,第i行第j个数表示矩阵第i行第j列的元素,每个元素模10^9+7输入输出样例输入样例#1: 2 1 1 1 1 1输出样例#1: 1 1 1 1这里写代码片#include<cstdio>原创 2017-05-31 07:55:38 · 281 阅读 · 0 评论 -
可做题2(Code+12月网络赛)(扩欧+矩阵快速幂)
题目背景“codeplus比赛的时候在做什么?有没有空?能来解决丢番图方程问题吗?”sublinekelzrip这样问qmqmqm。当然,qmqmqm并不会丢番图方程问题,所以sublinekelzrip改为提出了另一个题目,现在请你帮助qmqmqm解决这个题目。题目描述这个问题是这样的:若一个数列aa满足条件an=an−1+an−2,n>3,而a1,a2a_n=a_{n-1}+a_{n-2},n>原创 2017-12-25 19:14:57 · 790 阅读 · 0 评论 -
18.1.7 T3(AC+dp+矩阵)
问题描述 Magic Land上的人们总是提起那个传说:他们的祖先 John 在那个东方岛屿帮助 Koishi 与其姐姐 Satori 最终战平。而后,Koishi 恢复了读心的能力 如今,在 John 已经成为传说的时代,再次造访那座岛屿的人们却发现 Koishi遇到了新麻烦。这次她遇到了 Flandre Scarlet——她拥有可以使用禁忌魔法而不会受到伤害的能力。为了说明什么是禁忌魔原创 2018-01-07 19:16:11 · 262 阅读 · 0 评论 -
bzoj4861 [Beijing2017]魔法咒语(AC+矩阵优化dp)
DescriptionChandra 是一个魔法天才。从一岁时接受火之教会洗礼之后,Chandra 就显示出对火元素无与伦比的亲和力,轻而易举地学会种种晦涩难解的法术。这也多亏 Chandra 有着常人难以企及的语言天赋,让她能轻松流利地说出咒语中那些极其拗口的魔法词汇。直到十四岁,开始学习威力强大的禁咒法术时,Chandra 才遇到了障碍。根据火之魔法规则,禁咒的构成单位是 N 个基本词汇。施法时原创 2018-01-22 17:31:11 · 491 阅读 · 0 评论 -
bzoj4180 字符串计数(SAM+二分+矩阵加速dp)
题目链接分析: 问题建立在字符串的子串上 而SAMSAMSAM包含了字符串的所有子串,所以可以往这个方向考虑一开始我设计的状态:f[i][j]f[i][j]f[i][j]表示第iii位是字符jjj,用的最大步数 然后我就很直接的想到,在转移的时候可以找到以jjj为起点的子串(j,...,k):len(j,...,k):len(j,...,k):len f[i+len−1][k]=m...原创 2018-01-30 20:44:46 · 393 阅读 · 0 评论 -
hdu2865 Birthday Toy(置换+矩阵优化dp)
题目链接分析: 首先中间的大珠子需要一种颜色,而这种颜色是不能为周围珠子使用的:m=k-1如果没有“相邻颜色不同”这个限制,这道题就和poj2154一样了 poj2154 旋转i个珠子,轮换个数就是gcd(n,i)" role="presentation" style="position: relative;">gcd(n,i)gcd(n,i)gcd(n,i) 不原创 2018-02-05 20:56:56 · 382 阅读 · 0 评论 -
poj2888 Magic Bracelet(置换+矩阵)
题目链接分析: 做过hdu2865之后就不虚了n个珠子的手镯,n种置换,每个置换的轮换个数gcd(i,n)" role="presentation" style="position: relative;">gcd(i,n)gcd(i,n)gcd(i,n) (比较基础的东西)对于一个置换(转动i" role="presentation" style="position: relati原创 2018-02-06 08:27:39 · 307 阅读 · 0 评论 -
全场六零赛(1.29)
T1【题目描述】 吾乃闪耀知识的灯塔。 B 君有个n个点,m个边的仙人掌。所谓仙人掌,就是任何一个点至多属于一个环。 每条边有 1/2 的概率被删掉。问期望剩下多少个边联通块。所谓边联通块,就是问剩下的边,构成多少个联通块,单独一个点不算做联通块。 B君不喜欢实数,B君想知道答案乘以2m" role="presentation" style="position: relative;">原创 2018-01-29 16:45:07 · 328 阅读 · 0 评论 -
poj3744 Scout YYF I(概率dp+矩阵加速)
题目链接分析: 可以看出一个很naive的转移方程: f[i]=p∗f[i−1]+(1−p)∗f[i−2]f[i]=p∗f[i−1]+(1−p)∗f[i−2]f[i]=p*f[i-1]+(1-p)*f[i-2](i−1,i−2i−1,i−2i-1,i-2都不是地雷)在两个相邻的地雷之中,有一段安全区域 我突然发现,通过这个安全区域的方案数恰好是斐波那契数列(f[i]=f[i−1]+f...原创 2018-03-06 18:44:55 · 244 阅读 · 0 评论 -
打瞌睡hu测(2.23)
T1T2分析: 考场上我想了一个奇技淫巧,竟然过了40% 考虑一棵树,以每个点为根记录距离ta最远的叶节点的距离(dis) 找到这棵树中dis值最大的结点,暂且称这个结点为此树的重心根据贪心的想法,我们直接把重心设置为黑色 但是重心有可能存在某一子结点XXX,使得dis[X]+w(重心,X)>mdis[X]+w(重心,X)>mdis[X]+w(重心,X)>m...原创 2018-02-23 17:39:20 · 195 阅读 · 0 评论 -
loli终于又来hu测了 T3.回文串(dp+矩阵加速【difficult)
分析: 神题 心里默念着不能放弃,于是10分到手30%题意相当于问有多少回文串包含sss这个子序列 考虑如果判断一个回文串中是否包含子序列sss 我们可以从小到大枚举iii,分别用aiaia_i和an−i+1an−i+1a_{n-i+1}去匹配当前sss两端的字母,并将sss中被成功匹配的字符删去, 我们可以据此进行dp,设f[i][l][r]f[i][l][r]...原创 2018-03-26 16:33:44 · 230 阅读 · 0 评论 -
爆零赛(12.24)
平安夜的hu测 mdzz最近就是有点颓,明知道这次成绩要上报,还是妥妥的爆了零。。。T1 分析:%60:做这道题之前,要了解一些复数运算 手玩了r的1~8次方之后,发现i前面一定会有一个sqrt(7),而常数项中一定没有sqrt(7) 实际上这就是复数的计算法则: 而我们需要求的答案也是符合这个形式的(只有i前面有sqrt(7)) 那我们就可以把i前面的sqrt(7)抹掉 (但是如果原创 2017-12-24 18:25:09 · 329 阅读 · 0 评论 -
爆零赛T2(矩阵加速递推式)
前言: 曲神:今天中午做了个噩梦,梦见给DP讲题,讲了半个小时。。。分析: 终于在曲神,ZYXZYXZYX和舒老师的共同努力下,zz博主终于学会了70%做法动态规划 我们首先就要设计状态:f[i][j]:表示到第i位为止,[i-j+1,i]这j位每一位上的颜色都不同(一共有j种颜色),这种情况下的方案数我们需要从f[i-1]转移到f[i] 有两种情况: (色差盲福利) 得到的状态转移:原创 2017-12-25 15:31:41 · 377 阅读 · 0 评论 -
luogu1962斐波那契数列
题目背景大家都知道,斐波那契数列是满足如下性质的一个数列:• f(1) = 1 • f(2) = 1 • f(n) = f(n-1) + f(n-2) (n ≥ 2 且 n 为整数)题目描述请你求出 f(n) mod 1000000007 的值。输入输出格式输入格式: ·第 1 行:一个整数 n输出格式: 第 1 行: f(n) mod 1000000007 的值输入输出样例输入样例#1:原创 2017-05-31 18:44:04 · 561 阅读 · 0 评论 -
nyoj301递推求值
描述 给你一个递推公式:f(x)=a*f(x-2)+b*f(x-1)+c并给你f(1),f(2)的值,请求出f(n)的值,由于f(n)的值可能过大,求出f(n)对1000007取模后的值。注意:-1对3取模后等于2输入 第一行是一个整数T,表示测试数据的组数(T<=10000) 随后每行有六个整数,分别表示f(1),f(2),a,b,c,n的值。 其中0<=f(1),f(2)<100,-10原创 2017-05-31 18:48:41 · 293 阅读 · 0 评论 -
luogu2044 随机数生成器
题目描述栋栋最近迷上了随机算法,而随机数是生成随机算法的基础。栋栋准备使用线性同余法(Linear Congruential Method)来生成一个随机数列,这种方法需要设置四个非负整数参数m,a,c,X[0],按照下面的公式生成出一系列随机数{Xn}: X[n+1]=(aX[n]+c) mod m其中mod m表示前面的数除以m的余数。从这个式子可以看出,这原创 2017-05-31 18:58:50 · 244 阅读 · 0 评论 -
luogu2461递归数列
题目描述一个由自然数组成的数列按下式定义:对于i <= k:ai = bi 对于i > k: ai = c1ai-1 + c2ai-2 + … + ckai-k 其中bj 和 cj (1<=j<=k)是给定的自然数。写一个程序,给定自然数m <= n, 计算am + am+1 + am+2 + … + an, 并输出它除以给定自然数p的余数的值。输入输出格式输入格式: 输入文件spp.in由四原创 2017-05-31 18:54:37 · 321 阅读 · 0 评论 -
bzoj1875 [SDOI2009]HH去散步(矩阵)
Description HH有个一成不变的习惯,喜欢饭后百步走。所谓百步走,就是散步,就是在一定的时间 内,走过一定的距离。 但是同时HH又是个喜欢变化的人,所以他不会立刻沿着刚刚走来的路走回。 又因为HH是个喜欢变化的人,所以他每天走过的路径都不完全一样,他想知道他究竟有多 少种散步的方法。 现在给你学校的地图(假设每条路的长度都是一样的都是1),问长度为t,从给定地 点A走到给定地点B共有多少原创 2017-08-09 09:31:10 · 427 阅读 · 0 评论 -
bzoj1297 [SCOI2009]迷路(矩阵优化)
Description windy在有向图中迷路了。 该有向图有 N 个节点,windy从节点 0 出发,他必须恰好在 T 时刻到达节点 N-1。 现在给出该有向图,你能告诉windy总共有多少种不同的路径吗? 注意:windy不能在某个节点逗留,且通过某有向边的时间严格为给定的时间。Input 第一行包含两个整数,N T。 接下来有 N 行,每行一个长度为 N 的字符串。 第i行第j列为’0’原创 2017-09-28 20:26:16 · 316 阅读 · 0 评论 -
bzoj1009 [HNOI2008]GT考试(AC+矩乘优化dp)
Description 阿申准备报名参加GT考试,准考证号为N位数X1X2….Xn(0<=Xi<=9),他不希望准考证号上出现不吉利的数字。 他的不吉利数学A1A2…Am(0<=Ai<=9)有M位,不出现是指X1X2…Xn中没有恰好一段等于A1A2…Am. A1和X1可以为 0Input 第一行输入N,M,K.接下来一行输入M位的数。 N<=10^9,M<=20,K<=1000Out原创 2017-09-28 10:41:37 · 316 阅读 · 0 评论 -
随 (rand)(校内hu测10.6T1)(dp+矩阵+数论)
【题目描述】 给出n个正整数a1,a2…an和一个质数mod.一个变量x初始为1.进行m次操作.每次在n个数中随机选一个ai,然后x=x*ai%mod.问m次操作之后x的取值的期望. 答案一定可以表示成a/b的精确分数形式.a和b可能很大,所以只需要输出a*(b^(10^9+5))模10^9+7的结果.【输入格式】 第一行三个整数n,m,mod. 接下来一行n个空格隔开的正整数a1,a2…原创 2017-10-06 15:01:07 · 512 阅读 · 0 评论 -
bzoj4818 [Sdoi2017]序列计数(矩阵)
DescriptionAlice想要得到一个长度为n的序列,序列中的数都是不超过m的正整数,而且这n个数的和是p的倍数。Alice还希望,这n个数中,至少有一个数是质数。Alice想知道,有多少个序列满足她的要求。 Input 一行三个数,n,m,p。 1<=n<=10^9,1<=m<=2×10^7,1<=p<=100 Output 一行一个数,满足Alice的要求的序列数量,答案对201原创 2017-08-07 20:34:13 · 459 阅读 · 0 评论 -
UVa10870 - Recurrences(矩阵)
题目链接分析: 矩阵快速幂求递推式 练模板//这里写代码片#include<cstdio>#include<cstring>#include<cstdio>#define ll long longusing namespace std;ll n,p,f[20],a[20];int d;struct node{ ll H[20][20]; node operator *原创 2017-10-08 19:40:08 · 235 阅读 · 0 评论 -
矩阵初步(线性代数在信息学中的体现)
什么是矩阵:矩阵(Matrix)在数学上,矩阵是指纵横排列的二维数据表格,最早来自于方程组的系数及常数所构成的方阵。这一概念由19世纪英国数学家凯利首先提出。矩阵是高等代数学中的常见工具,也常见于统计分析等应用数学学科中。在物理学中,矩阵于电路学、力学、光学和量子物理中都有应用;计算机科学中,三维动画制作也需要用到矩阵。 矩阵的运算是数值分析领域的重要问题。将矩阵分解为简单矩阵的组合可以在理...原创 2019-08-06 17:20:35 · 671 阅读 · 0 评论