ACM-ICPC Live Archive:6184 --- One-Dimensional Cellular Automaton（矩阵快速幂）

思路：

如果将题目中原先的序列，（假设n为5）我们可以通过一个构造转换矩阵来实现对原先序列的一个转换，转换如下：

S(0,T-1) 　 B C 0 0 0 　 S(0,T)

S(1,T-1) 　 A B C 0 0 　 S(1,T)

S(2,T-1) X 0 A B C 0 = 　S(2,T)

S(3,T-1) 　 0 0 A B C 　 S(3,T)

S(4,T-1) 　 0 0 0 A B 　 S(4,T)

这样转化后就变成一个求某矩阵的T次幂，然后再做一次矩阵乘法就可以了，T很大，但是可以通过快速幂来实现。

AC代码：

#include <iostream>
#include <cstring>
#include <cstdlib>
#include <algorithm>
#include <queue>
#include <cmath>
#include <cstdio>
#define mem(a, b) memset(a, b, sizeof(a))
#define PI 3.1415926535
using namespace std;
typedef long long ll;
const int MAXn = 55;

struct matirix
{
    int a[MAXn][MAXn];
} r, o;  //o为转换矩阵
int sum[MAXn], n, m, a, b, c, t; //sum[MAXn]为初始矩阵
matirix cheng(matirix a, matirix b)
{//矩阵相乘
    int x, y, z;
    matirix t;
    mem(t.a, 0);
    for (x = 0; x < n; x++)
        for (y = 0; y < n; y++)
        {
            for (z = 0; z < n; z++)
            {
                t.a[x][y] += (a.a[x][z] * b.a[z][y]) % m;
                t.a[x][y] %= m;
            }
        }
    return t;
}
void sort_pow(int t)
{//矩阵快速幂
    while (t)
    {
        if (t & 1)
            r = cheng(r, o);
        o = cheng(o, o);
        t >>= 1;
    }
}
void init()
{
    mem(r.a,0); //清零
    mem(o.a,0);
    int i, j;
    scanf("%d%d%d%d%d", &m, &a, &b, &c, &t);
    r.a[0][0] = r.a[n - 1][n - 1] = 1; //构建矩阵
    o.a[0][0] = b;
    o.a[0][1] = c;
    o.a[n - 1][n - 2] = a;
    o.a[n - 1][n - 1] = b;
    for (i = 1; i < n - 1; i++)
    {
        o.a[i][i - 1] = a;
        o.a[i][i] = b;
        o.a[i][i + 1] = c;
        r.a[i][i] = 1;
    }
    sort_pow(t); 
    for (i = 0; i < n; i++)
        scanf("%d", sum + i);
    for (i = 0; i < n; i++)
    {
        int ans = 0;
        for (j = 0; j < n; j++)
        {
            ans += sum[j] * r.a[i][j];
            ans %= m;
        }
        if (i != 0)
            printf(" ");
        printf("%d", ans);
    }
    printf("\n");
}
int main()
{
    while (1)
    {
        scanf("%d", &n);
        if (n == 0)
            break;
        init();
    }
}