poj 3150 Cellular Automaton

http://poj.org/problem?id=3150

不擅长数学类的题呀

网上很多此题的解题报告 都不错 但是没有注释

推荐解题报告：http://www.cppblog.com/varg-vikernes/archive/2011/02/08/139804.html

我的代码及其注释：

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
#define LL long long

using namespace std;

const int N=505;
LL a[N];
LL b[N];
LL temp[N];
int n,m,d,k;
void Multi(LL l1[],LL l2[])//把 l1*l2 赋值给 l1
{
    for(int i=0;i<n;++i)
    {
        temp[(n-i)%n]=0;
        for(int j=0;j<n;++j)
        {
            temp[(n-i)%n]+=(l2[(i+j)%n]*l1[j]);
        }
    }
    for(int i=0;i<n;++i)
    l1[i]=temp[i];
}
int main()
{
    //freopen("data.txt","r",stdin);
    while(scanf("%d %d %d %d",&n,&m,&d,&k)!=EOF)
    {
        for(int i=0;i<n;++i)
        scanf("%I64d",&a[i]);
        memset(b,0,sizeof(b));
        for(int i=0,j=0;i<=d;++i,--j)
        b[i]=b[j+n]=1;
        while(k)
        {
            if(k&1)//如果这时的k末位为1
            Multi(a,b);//则正好乘上这时的b
            Multi(b,b);//b每次循环都变成了原来的b*b
            k=k>>1;
        }
        for(int i=0;i<n;++i)
        {
            if(i)
            printf(" ");
            printf("%I64d",a[i]%m);
        }
        printf("\n");
    }
    return 0;
}