J - One-Dimensional HYSBZ - 4688

J - One-Dimensional HYSBZ - 4688

Problem Description

考虑一个含有 N 个细胞的一维细胞自动机。细胞从 0 到 N-1 标号。每个细胞有一个被表示成一个小于 M 的非负
整数的状态。细胞的状态会在每个整数时刻发生骤变。我们定义 S(i,t) 表示第 i 个细胞在时刻 t 的状态。在
时刻 t+1 的状态被表示为 S(i,t+1)=(A×S(i-1,t)+B×S(i,t)+C×S(i+1,t) ) mod M ，其中 A,B,C 是给定的非
负整数。对于 i<0 或 N≤i ，我们定义 S(i,t)=0 。给定一个自动机的定义和其细胞在时刻 0 的初始状态，你的
任务是计算时刻 T 时每个细胞的状态。

Input

输入包含多组测试数据。每组数据的第一行包含六个整数 N,M,A,B,C,T ，满足 0<N≤50,0<M≤1000,0≤A,B,C<M,0
≤T≤〖10〗^9 。第二行包含 N 个小于 M 的非负整数，依次表示每个细胞在时刻 0 的状态。输入以六个零作为
结束。

Output

对于每组数据，输出N个小于M的非负整数，每两个相邻的数字之间用一个空格隔开，表示每个细胞在时刻T的状态。

Sample Input

5 4 1 3 2 0

0 1 2 0 1

5 7 1 3 2 1

0 1 2 0 1

5 13 1 3 2 11

0 1 2 0 1

5 5 2 0 1 100

0 1 2 0 1

6 6 0 2 3 1000

0 1 2 0 1 4

20 1000 0 2 3 1000000000

0 1 2 0 1 0 1 2 0 1 0 1 2 0 1 0 1 2 0 1

30 2 1 0 1 1000000000

0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0

30 2 1 1 1 1000000000

1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

30 5 2 3 1 1000000000

0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0

Sample Output

0 1 2 0 1

2 0 0 4 3

212 10 9 11

3 0 4 2 1

0 4 2 0 4 4

0 376 752 0 376 0 376 752 0 376 0 376 752 0 376 0 376 752 0 376

1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0

1 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0

1 1 3 2 2 2 3 3 1 4 3 1 2 3 0 4 3 3 0 4 2 2 2 2 1 1 2 1 3 0

题解

矩阵快速幂，关键是找到矩阵的通项矩阵

BC0000
ABC000
0ABC00
00ABC0
000ABC

矩阵存储：

[i][j], i 为横行，j为纵列，下述程序中表示的矩阵为原矩阵的转置

#include<bits\stdc++.h>
using namespace std;

int N, M, A, B, C, T;
typedef struct node {
	int e[55][55];
	int Init(){
		memset(e, 0, sizeof(e));
		for(int i = 1; i <= N; i++)
			e[i][i] = 1;
	}
}Mar;

Mar mp, ans;

Mar Mut(Mar a, Mar b) { //  转置后矩阵的乘法
	Mar c;
	memset(c.e, 0, sizeof(c.e));
	for(int i = 1; i <= N; i++)
		for(int j = 1; j <= N; j++)
			for(int k = 1; k <= N; k++){
				c.e[i][j] += a.e[i][k] * b.e[k][j];
				c.e[i][j] %= M;
			}
	return c;	
}

Mar quick_mod(Mar a, int power) {
	Mar re;
	re.Init();
	while(power) {
		if(power&1) re = Mut(re, a);
		a = Mut(a, a);
		power >>= 1;
	}
	return re;
}

int main() {
	freopen("test.in", "r", stdin);
	while(cin >> N >> M >> A >> B >> C >> T && (N+M+A+B+C+T)) {
		for(int i = 1; i <= N; i++){
			mp.e[i][i+1] = A;
			mp.e[i][i] = B;
			mp.e[i][i-1] = C;
		}
		for(int i = 1; i <= N; i++) cin >> ans.e[1][i];
		ans = Mut(ans, quick_mod(mp, T));
		for(int i = 1; i < N; i++) cout << ans.e[1][i] << " ";
		cout << ans.e[1][N] << endl;
			
	}
	return 0;
}