Hdu 2855 Fibonacci Check-up （线性代数_矩阵）

题目链接：http://acm.hdu.edu.cn/showproblem.php?pid=2855

题目大意：有为组合公式，F(k)为第k个Fibonacci数，问给定一个小于10亿的n和一个m（矩阵里的数对m取余），求最后结果是多少。

解题思路：刚看到这个公式觉得无比蛋疼，这个并不是线性的，用矩阵解不了吧，暴力又不可以...后来想YY下看有没线性的规律，找着找着果然发现规律了，其实就等于F(2n)，那瞬间变简单。

测试数据：

4

1 3000

2 3000

5 40

1000000000 30000

代码：

#include <stdio.h>
#include <math.h>
#include <string.h>


int n,m;

struct Mat{

	int mat[2][2],size;
	Mat() { 
		
		size = 2;
		mat[0][0] = 1,mat[0][1] = 1;
		mat[1][0] = 1,mat[1][1] = 0;
	}


	friend Mat operator *(Mat a,Mat b);
	friend Mat operator +(Mat a,Mat b);
	friend Mat operator ^(Mat a,int k);
}E;
Mat operator *(Mat a,Mat b){

	Mat c;
	memset(c.mat,0,sizeof(c.mat));
	for (int i = 0; i < c.size; ++i)
		for (int j = 0; j < c.size; ++j)
			for (int k = 0; k < c.size; ++k)
				if (a.mat[i][k] && b.mat[k][j])
					c.mat[i][j] = (c.mat[i][j] + a.mat[i][k] * b.mat[k][j]) % m;
	return c;
}
Mat operator +(Mat a,Mat b){

	Mat c;
	for (int i = 0; i < c.size; ++i)
		for (int j = 0; j < c.size; ++j)
			c.mat[i][j] = (a.mat[i][j] + b.mat[i][j]) % m;
	return c;
}
Mat operator ^(Mat a,int k) {

	Mat c = E;
	while (k) {

		if (k & 1) c = c * a;
		a = a * a,k >>= 1;
	}
	return c;
}


int main()
{
	int i,j,k,t;


	for (i = 0; i < E.size; ++i)
		E.mat[i][i] = 1,E.mat[i][1-i] = 0;


	scanf("%d",&t);
	while (t--){

			scanf("%d%d",&n,&m);
			Mat c; c = c^(2 * n);
			printf("%d\n",c.mat[0][1] % m);
	}
}

本文ZeroClock原创，但可以转载，因为我们是兄弟。