![](https://img-blog.csdnimg.cn/20201014180756928.png?x-oss-process=image/resize,m_fixed,h_64,w_64)
矩阵快速幂
文章平均质量分 83
wjw1340
这个作者很懒,什么都没留下…
展开
-
HDU4565 So Easy!【矩阵快速幂】
题意:给出a,b,n,m。求思路:(a + sqrt(b))^n = Xn + Yn * sqrt(b)(a+sqrt(b))^(n+1) = (Xn + Yn * sqrt(b)) *(a + sqrt(b)) = aXn+bYn + (Xn+aYn) * sqrt(b)Xn+1 = aXn+bYn; Yn+1 =(Xn+aYn);X0 = 1;Y0 = 0原创 2017-07-29 18:59:12 · 339 阅读 · 0 评论 -
HDU4549 M斐波那契数列【矩阵快速幂】
题意:F0 = a ,F1 = b ,Fn = Fn-1 * Fn-2 (n > 1),求Fn思路:F2 = a*b F3 = a*b^2 F4 = a^2*b^3 F5 = a^3*b^5 F6 = a^5*b^8 观察易得: Fn原创 2017-07-31 16:54:23 · 223 阅读 · 0 评论 -
HDU4990 Reading comprehension【矩阵快速幂】
题意:偶数项f[i] = 2f[i-1],奇数项f[i] = 2f[i-1]+1思路:f[i] = f[i-1] + 2f[i-2] + 1 直接一个公式,矩阵快速幂开始我没有发现这个通项公式,我看到的是偶数项f[i] = 4f[i-2] + 2,快速幂偶数项,奇数项f[i] = f[i+1] / 2,好像有点道理,但是仔细想想这里可是取模的。没事瞎搞什么除法,奇数项f原创 2017-07-27 09:29:58 · 320 阅读 · 0 评论 -
HDU4965 Fast Matrix Calculation【矩阵快速幂】
题意:C = A * B, M = C ^ (n*n),求M中元素MOD 6以后的SUM思路:A*B最大1000*1000,快速幂要炸,我们转换一下比如(A*B)^4 = ABABABAB = A(BA)(BA)(BA)B = A*(BA)^3*BB*A最大6*6 这样就快了为什么我结构体里开数组s[1005][1005],程序一运行就停,照理说这也不大啊,不清楚怎原创 2017-07-28 17:45:34 · 227 阅读 · 0 评论 -
HDU3411 Snail Alice【矩阵快速幂】
题意:,求f(n), q=x1^y1+z1, n=2^y2+z2 ,0思路:推公式就算了,打个表,找规律,发现f(n) = (q-1)*f(n-1) + q * f(n-2),矩阵快速幂,我们要求一个矩阵的n-1次,n =2^y2+z2,y2可以达到几万,2^y2存不下,但是底数是2,2^y2相当于1个1后面跟了y2个0,把那个1减过来,就是连续的y2个1,模拟下y2个1的快速幂,原创 2017-08-30 08:56:55 · 259 阅读 · 0 评论 -
Light1052 String Growth【矩阵快速幂】
题意:一个字符串由a,b组成,字符串转化为下一个的规则为:a->b,b->ab,告诉你,第n个、第m个字符串的长度,算第k个的长度思路:这个变化规则,我们可以发现长度的变化是斐波拉契数列,我们设第1项有q个a,p个b,往下推,q、p前面的系数也是斐波拉契数列,我们快速幂算第n、m项q、p的前的系数,a1*q+b1*p = 第n项的长度,解二元一次方程,算出q、p算第k项解方原创 2017-09-01 16:38:37 · 232 阅读 · 0 评论