文章分类 -  矩阵专辑

摘要：f[n]=2*f[n-2]+f[n-1]+1,f[1]=1,f[2]=2#include <stdio.h>#define CIRCLE 3186#define MOD 200907int main(){ int a[CIRCLE]={1,1,2}; int i;/* find the circle for(i=3;i<10000;i++) { a[i] = (2*a[i-2]... 阅读全文

摘要：#include <stdio.h>#define MAXN 20#define MOD 9973struct Matrix{ Matrix(int,int); Matrix(int,int,int arr[MAXN][MAXN]); friend Matrix operator*(const Matrix &,const Matrix &); int r,c,a[MA... 阅读全文

摘要：#include <stdio.h>#include <string.h>#define MAXN 10int MOD,mat[MAXN][MAXN];struct Matrix{ Matrix(); Matrix(int arr[MAXN][MAXN]); friend Matrix operator*(const Matrix &,const Matrix &a... 阅读全文

摘要：把对应的坐标，换成要走的步数F（n）表示从原点走n步的概率。易根据题意得到F(n)=p*F(n-1)+q*F(n-2) 其中F(0)=0,F(1)=1，q=1-p用矩阵二分加速：[F(n) F(n-1)] = [1 0] [p q]　　　　　　　 [0 1] * [1 0]^(n-1)　对于第一个要走n步的雷，特工要存活就必须跨过这个雷，也就是走n-1步。在跨2步。这个概率就是F(n-1)*q.... 阅读全文

摘要：http://acm.hust.edu.cn/thx/problem.php?id=1384#include <stdio.h>#define MOD 9901#define N 3struct Matrix{ Matrix(); Matrix(int arr[N][N]); int row,col,a[N][N];};Matrix::Matrix(){ row=col=N; for(... 阅读全文

摘要：#include <iostream>using namespace std;struct Matrix{ Matrix(int a,int b,int c,int d):a11(a),a12(b),a21(c),a22(d) {} friend Matrix operator*(const Matrix &M1,const Matrix &M2); long a11,... 阅读全文