HDU - 2276 Kiki & Little Kiki 2（矩阵快速幂）

题目大意：本题所有的灯构成一个环，某时刻如果某盏灯的左边为on，则改变该灯的状态，问一段时间过后所有灯的状态。

题解：

   假设题目给的01串的长度为n，保存在a数组，下标从0开始

   根据上面的规则我们发现可以得出一秒过后的状态即为a[i] = (a[i]+a[i-1])%2 , 对于a[0] = (a[0]+a[n-1])%2

   那么我们就可以就能够找到递推的式子

   1 1 0 0....     a0        a1

   0 1 1 0...  *  a1   =   a2

   ..........1 1     .....      .....

   1 0 0.....1     an-1    a0

----->

即

那么我们要求at就直接矩阵的

注意乘的时候a0,a1,a2…….并不是只是一列

而是乘

#include<cstring>
#include<iostream>
#include<cstdio>
using namespace std;
#define ll long long
#define mod 2
int k;
int len;
struct matirx//注意矩阵的大小
{
    ll ma[110][110];
};
matirx a;
string s;
matirx mutimatirx(matirx a,matirx b)//矩阵乘
{
    matirx t;
    memset(t.ma,0,sizeof(t.ma));
    for(int i=0; i<len; i++)
        for(int k=0; k<len; k++)
            for(int j=0; j<len; j++)
            {
                t.ma[i][j]+=a.ma[i][k]*b.ma[k][j];
                t.ma[i][j]%=2;//注意要不要取模
            }
    return t;
}
matirx powmatirx(matirx a,long long b)//矩阵快速幂
//返回矩阵a的b次方
{
    matirx t;
    memset(t.ma,0,sizeof(t.ma));
    for(int i=0; i<len; i++)
        t.ma[i][0]=s[i]-'0';
    while(b)
    {
        if(b&1)
        {
            t=mutimatirx(a,t);
        }
        a=mutimatirx(a,a);
        b>>=1;
    }
    return t;
}
void init(int len)
{
    memset(a.ma,0,sizeof(a.ma));
    for(int i=0;i<len;++i)
        for(int j=0;j<len;++j)
        {
            if(i==0)
            {
                if(j==0 || j==len-1)
                    a.ma[i][j]=1;
                else a.ma[i][j]=0;
            }
            else
            {
                if(j==i||j==i-1)
                    a.ma[i][j]=1;
                else a.ma[i][j]=0;
            }
        }
}
int main()
{
    //freopen("input.txt","r",stdin);
    ll t;
    while(cin>>t)
    {
        cin>>s;
        len=s.length();
        init(len);
        a=powmatirx(a,t);
        for(int i=0;i<len;++i)
            cout<<a.ma[i][0];
        cout<<endl;
    }
    return 0;
}