兩個 soft 滿足 brush div 之間 == script opera

【BZOJ4688】One-Dimensional

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
2 12 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

題解：矩乘傻題。

#include <iostream>
#include <cstring>
#include <cstdio>
using namespace std;
int n,m,A,B,C,T;
struct M
{
	int a[60][60];
	int * operator [] (int b) {return a[b];}
	M() {memset(a,0,sizeof(a));}
	M operator * (M b)
	{
		M c;
		int i,j,k;
		for(i=0;i<n;i++)	for(j=0;j<n;j++)	for(k=0;k<n;k++)	c[i][j]=(c[i][j]+a[i][k]*b[k][j])%m;
		return c;
	}
}tr,ans;
void pm(int y)
{
	while(y)
	{
		if(y&1)	ans=ans*tr;
		tr=tr*tr,y>>=1;
	}
}
int main()
{
	while(1)
	{
		scanf("%d%d%d%d%d%d",&n,&m,&A,&B,&C,&T);
		if(!m)	return 0;
		int i;
		memset(tr.a,0,sizeof(tr.a)),memset(ans.a,0,sizeof(ans.a));
		for(i=0;i<n;i++)
		{
			scanf("%d",&ans[0][i]);
			if(i)	tr[i-1][i]=A;
			tr[i][i]=B;
			if(i<n-1)	tr[i+1][i]=C;
		}
		pm(T);
		for(i=0;i<n;i++)	printf("%d%c",ans[0][i],i==n-1?‘\n‘:‘ ‘);
	}
}

【BZOJ4688】One-Dimensional 矩陣乘法