HDU：2604 Queuing（发现似乎所有…

此题的题意是求出长度为L的 ，有f 和m 组成的窜中找出没有存在 fmf 和fff 的窜的个数

 

求解f[n]=f[n-1]+f[n-3]+f[n-4]；我们就可以构造出这样的 矩阵

这是搜的代码：

 

#include <iostream>

using namespace std;

const int N = 4;

struct Mat

{

int matrix[N][N];

};

Mat mat, mt;

int n, m;

Mat mul(Mat a, Mat b)

{

int i, j, k;

Mat c;

for (i = 0; i < N; i++)

for (j = 0; j < N; j++)

{

c.matrix[i][j] = 0;

for (k = 0; k < N; k++)

{

c.matrix[i][j] += a.matrix[i][k] * b.matrix[k][j];

c.matrix[i][j] %= m;

}

}

return c;

}

Mat solve(int m)

{

Mat mt;

if (m == 1)

return mat;

if (m & 1)

return mul(solve(m - 1), mat);

else

{

mt = solve(m / 2);

return mul(mt, mt);

}

}

void init()

{

int i, j;

for (i = 0; i < N; ++i)

{

mat.matrix[0][i] = 1;

}

mat.matrix[0][1] = 0;

for (i = 1; i < N; ++i)

{

for (j = 0; j < N; ++j)

{

if (i - j == 1)

mat.matrix[i][j] = 1;

else

mat.matrix[i][j] = 0;

}

}

}

int main()

{

Mat ans;

int f41[] =

{ 9, 6, 4, 2 }, sum, f04[] =

{ 0, 2, 4, 6, 9 };

//freopen("t","r",stdin);

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

{

if (n < 5)

{

printf("%d\n", f04[n] % m);

}

else

{

init();

ans = solve(n - 4);

sum = 0;

for (int k = 0; k < N; k++)

{

sum += ans.matrix[0][k] * f41[k];

sum %= m;

}

printf("%d\n", sum);

}

}

return 0;

}