LA3704矩阵变换

/*LA3704：
细胞自动机(矩阵)：
关键是辨别和变换矩阵的构造
变换矩阵Mat：
第i列的构造，for(j=i-d;j<=i+d;j++) {k=(j+n)%n,M.mat[k][i]=1;},其余为零
取模的作用比较巧妙
*/
#include<iostream>
#include<stdio.h>
#include<string.h>
#include<algorithm>
#include<stdlib.h>
#include<math.h>
#include<queue>
#include<vector>
#include<map>
#define ULL unsigned long long
#define LL long long
#define Max_size 17
#define size n
#define MOD m
using namespace std;

int n,m,d,k;

struct Mat
{
    LL mat[Max_size][Max_size];

};
Mat operator*(Mat a,Mat b)
{
    Mat c;
    memset(c.mat,0,sizeof(c.mat));
    for (int i=0;i<size;i++)
        for (int j=0;j<size;j++)
            for (int k=0;k<size;k++)
                if (a.mat[i][k] && b.mat[k][j])
                   c.mat[i][j]=(c.mat[i][j]+a.mat[i][k]*b.mat[k][j])%MOD;
    return c;
}
Mat operator+(Mat a,Mat b)
{
    Mat c;
    memset(c.mat,0,sizeof(c.mat));
    for (int i=0;i<size;i++)
        for (int j=0;j<size;j++)
            c.mat[i][j]=(a.mat[i][j]+b.mat[i][j])%MOD;
    return c;
}
Mat E;

void builtE()
{
    memset(E.mat,0,sizeof(E.mat));
    for (int i=0;i<Max_size;i++)
    {
        E.mat[i][i]=1;
    }
}
Mat operator^(Mat A,int x)
{
    Mat c;
    builtE();
    c=E;
    for ( ; x; x>>=1)
    {
        if ( x & 1)
            c = c * A;
        A = A * A;
    }
    return c;
}
Mat K,F;

int main()
{
    while(cin>>n>>m>>d>>k)
    {
        memset(K.mat,0,sizeof(K.mat));
        memset(F.mat,0,sizeof(F.mat));
        for(int i=0;i<n;i++) cin>>F.mat[0][i];

        for(int i=0;i<n;i++)
        for(int j=i-d;j<=i+d;j++)
        {
            int p=(j+n)%n;K.mat[p][i]=1;
        }

        K=K^k;
        F=F*K;
        for(int i=0;i<n;i++) if (i<n-1) cout<<F.mat[0][i]<<" ";else cout<<F.mat[0][i]<<endl;

    }
    return 0;
}