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;
}