Hdu 2604 Queuing (线性代数_矩阵)

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

题目大意：求2个字母f和m构成的长度为m的序列中不含fmf以及fff的种数。

解题思路：找规律，前几个数值分别为2,4,6,9,15,25，一看就看出规律F[n] = F[n-1] + F[n-3] + F[n-4],那么这个规律是怎么找出来的呢？假设F[n-1]之前的合法方案都出来了，那么在F[n]处加个m一样合理，也就是1 * F[n-1]，那如果加的是f呢？加f的话对前面两个的值有要求，如果F[n-3]是f的话，肯定会导致那么后三个序列就为f?f，都是不合法的，如果F[n-3]是m，则为m？f,中间的F[n-2]为m的话就可以，方案数可增加F[n-3] * 1，如果？是f，那么F[n-4]就要是m才行，要保证它合法啊，所以方法数增加F[n-4]。有人会问，F[n-3]、F[n-4]不会包含在F[n-1]里吗？之前就算进去了啊？是的，算进去了，但是前面增加的1 * F[n-1]，末尾是m，后面两种末尾是f,明显不一样，soga...

    本题暴力也可过，L太小，让我出题L就出1000万，然后卡下时限，这样就难倒一大片暴力哥。

测试数据：

3 8

4 7

1000000 30

代码：

#include <stdio.h>
#include <string.h>
#define MAX 5


int n,tot,m,arr[MAX][MAX];
struct Mat{

	int mat[MAX][MAX],size;
	Mat() {size = MAX;}
	Mat(int n) {size = n,memset(mat,0,sizeof(mat));}
	Mat(int arr[][MAX],int n) {size = n,memcpy(mat,arr,sizeof(mat));}


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

	Mat c(n);
	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(n);
	for (int i = 0; i < c.size; ++i)
		for (int j = 0; j < c.size; ++j)
			for (int k = 0; k < c.size; ++k)
				c.mat[i][j] = (a.mat[i][j] + b.mat[i][j]) % m;
	return c;
}
Mat operator ^(Mat a,int k) {

	Mat c(n);
	for (int i = 0; i < c.size; ++i)
		c.mat[i][i] = 1;
	while (k) {

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


int main()
{
	int i,j,sum,num[10] = {0,2,4,6,9};
	arr[0][0] = arr[0][2] = arr[0][3] = 1;
	arr[1][0] = arr[2][1] = arr[3][2] = 1;


	
	while (scanf("%d%d",&tot,&m) != EOF) {

		n = 4;
		Mat c(arr,n);
		if (tot < 4) printf("%d\n",num[tot]);
		else {

			c = c ^ (tot - 4);
			for (sum = i = 0; i < 4; ++i)
				sum = (sum + c.mat[0][i] * num[n-i]) % m;
			printf("%d\n",sum);
		}
	}
}

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