矩阵快速幂
文章平均质量分 61
a7f650ebd327889c
这个作者很懒,什么都没留下…
展开
-
HDU 4686 Arc of Dream
题意:给出N A0 AX AY B0 BX BY ,求mod 1,000,000,007a 0 = A0 a i = a i-1*AX+AY b 0 = B0 b i = b i-1*BX+BY 思路:原创 2017-08-16 13:02:38 · 181 阅读 · 0 评论 -
uva 10655 Contemplation! Algebra(矩阵快速幂)
题意:给出p=a+b ,q=a*b ,n求 a^n+b^n思路:由于在矩阵快速幂的专题里一开始就想a2+b2 =(a+b)(a+b)-2ab又往下推发现a3+b3=(a2+b2)(a+b)-a2b-ab2 =(a2+b2)(a+b)-(a+b)ab又继续往下推原创 2017-08-16 13:54:26 · 219 阅读 · 0 评论 -
uva 10518 How Many Calls?
题意:递归求斐波那契数列,求递归次数,即节点个数思路:我们分析可以得到递归次数同样是前两个答案的和加1ans[i] =ans[i-1]+ans[i-2]+11 , 1, 3,5,9,15,25。。。。。同时又发现 ans[i]是第i个斐波那契数的两倍减1#includeusing namespace std;long long N;struct node原创 2017-08-24 11:45:48 · 240 阅读 · 0 评论 -
uva11551experienced endeavour
题意:给你个长为n的数列 a给你一个变换n行,第一个数为,后跟着几个数,表示经过一次变换后的数由原序列哪几个数组成#includeusing namespace std;const long long N=1000;int n;struct node{ long long a[55][55];};node cheng(node a,node b){ no原创 2017-09-04 16:12:08 · 208 阅读 · 0 评论 -
C - Matrix Power Series POJ - 3233
不多说贴代码牛顿迭代#include#includeint n,m,k;struct matrix{ int a[33][33]; }a,b,c;//s13=s6+a7+a7*s6;//s12=s6+a6*s6;//A+A^2+A^3+A^4+A^5+A^6=A+A^2+A^3 +A^3*(A+A^2+A^3)//A+A^2+A^3+A^4= A+A^2原创 2017-08-15 10:49:37 · 165 阅读 · 0 评论 -
UVA 11149-Power of Matrix
题意:给一个矩阵A和常数K,求一个C=A+A2+A3+。。。+AK思路:用矩阵的矩阵快速幂|0 A| |1 0| |A A^2||0 0| |1 A| |0 0 |还有牛顿迭代A+A2+A3+A4+A5+A6 = A+A2+A3 + A3(A+A2+A3))A+A2+A3+A4+A5+A6 +A7= A+A2+A3 + A4 +A4(A+原创 2017-08-15 10:33:08 · 185 阅读 · 0 评论 -
Codeforces 392C Yet Another Number Sequence
水题,不多说#include#includelong long N;#define M 2struct node{ long long a[M][M];};node cheng(node a,node b){ node c; memset(c.a,0,sizeof(c.a)); for(int i=0;i<M;i++){ for(int j=0;j<M;j++)原创 2017-08-14 15:01:14 · 258 阅读 · 0 评论 -
hdu4965---Fast Matrix Calculation(矩阵)
题意:给一个 NXK的矩阵 A 和 一个 KXN 的矩阵 B 求一个C=(A*B)^(N*N)得出 矩阵C中每个元素mod6的和K思路:(A*B)为一个NXN的矩阵我们知道 矩阵乘法是 O(N^3)的复杂度 一次乘法 1e9就快要超时了但是(B*A)是一个 KXK的矩阵 显然算(B*A)的快速幂是一个正确的方向(A*B)(N*N)= ABABABAB原创 2017-08-14 14:30:18 · 187 阅读 · 0 评论 -
hdu-4990 Reading comprehension
题意:Read the program below carefully then answer the question. #pragma comment(linker, "/STACK:1024000000,1024000000") #include #include #include #include #include #include co原创 2017-08-14 12:32:29 · 253 阅读 · 0 评论 -
HDU 5015 233 Matrix --矩阵快速幂
题意:ai,j = a i-1,j +a i,j-1( i,j > 0) , a0,i=23333...., i+1个3,ai,0为输入的数值给出n,m,给出(a0,1 ), ( a0,n)求an,m思路:我们知道这是NXM的矩阵试着画出该矩阵,想推出ai,j发现推不出来,但是可以推出下一列应该一目了然了吧#include#define N 100000原创 2017-08-14 11:05:20 · 143 阅读 · 0 评论 -
Codeforces 450B Jzzhu and Sequences
题意:给出x,y,n,求f[n]mod1000000007 (109 + 7).思路:基础的矩阵快速幂| fi fi-1 | | 1 1| |fi+1 fi|| 0 0 | | -1 0| | 0 0|#include#define N 1000000007using nam原创 2017-08-13 12:22:08 · 153 阅读 · 0 评论 -
HDU6050-Funny Function
多校好难啊补题吧。。。。看了题解发现有些东西不清楚啊。。。。对于任意i>=1,当j>=3时,有 通过归纳法可以得到来推推吧:1. 当 i=1时显然成立,2. 假设当 i= k-1时成立3. 当 i=k时 F i-1,j = F i-1,j-1 + 2*F i-1,j-2 = F i-1,j-1 + F i-1,j-2 +原创 2017-07-28 12:26:05 · 174 阅读 · 0 评论 -
HDU - 4565 So Easy!
题意:给出a,b,n,m;0求思路:不会写,看来别人的博客补的。。。。。。题目重点在于(a-1) 2即 a-1(a−√b) (a−√b)^n 令 cn= (a+√b)^n + (a−√b)^n 这里提到 (a+√b)^n 与 (a−√b)^n 共轭,所以二者相加的cn为整数 又因为 (a-sqrt(b原创 2017-08-18 14:33:15 · 238 阅读 · 0 评论 -
uva 10870 Recurrences
题意:给你d,n,m,a1,a2,a3,......f(1),f(2),f(3),f(4).......定义:f(n) = a1f(n − 1) + a2f(n − 2) + a3f(n − 3) + . . . + ad f(n − d), for n > d求 f(n)%m思路:不多说水题#includeusing namespace std;struct原创 2017-08-18 13:43:36 · 158 阅读 · 0 评论 -
hdu 4549 M斐波那契数列
题意:f[0]=a;f[1]=b;f[i]=f[i-1]*f[i-2];给出a,b,n;求f[n];思路:f[0]=a^1 b^0f[1]=a^0 b^1f[2]=a^1 b^1f[3]=a^1 b^2f[4]=a^2 b^3f[5]=a^3 b^5可以看到a,b的指数是斐波那契数列也就是先求得a,b的指数再求 (a^n * b^m )原创 2017-08-16 13:16:53 · 207 阅读 · 0 评论 -
Scout YYF I POJ - 3744
题意:有一个无限长的路,路上布满地雷,初始在位置1,每次 p 的概率到下一格,(1- p)的概率到下二格给出 p 以及地雷坐标a[i],问安全走过的概率是多少?思路:递推,用矩阵快速幂推,p[ a[i] +1 ] 推到 p[ a[i+1] -1 ] ,于是p[ a[i+1] +1 ] = p[ a[i+1] -1 ] *(1-p) #include原创 2017-09-14 15:09:08 · 208 阅读 · 0 评论