AC日记——摆花

思路：

　　矩阵加速dp；

代码：

#include <cstdio>
#include <cstring>
#include <iostream>
#include <algorithm>

using namespace std;

#define ll long long
#define mod (1000000007)

struct MatrixType {
    int n,m;

    ll ai[][];

    void mem()
    {
        for(int i=;i<=n;i++)
        {
            for(int j=;j<=m;j++) ai[i][j]=;
        }
    }

    void debug()
    {
        printf("\n");
        for(int i=;i<=n;i++)
        {
            for(int j=;j<=m;j++) printf("%d ",ai[i][j]);
            printf("\n");
        }
        printf("\n");
    }
};

int n,m;

MatrixType mul(MatrixType aa,MatrixType bb)
{
    MatrixType res;
    res.n=aa.n,res.m=bb.m,res.mem();
    for(int i=;i<=res.n;i++)
    {
        for(int j=;j<=res.m;j++)
        {
            for(int v=;v<=res.n;v++)
            {
                res.ai[i][j]=(res.ai[i][j]+aa.ai[i][v]*bb.ai[v][j])%mod;
            }
        }
    }
    return res;
}

MatrixType poww(int mi,MatrixType pos)
{
    MatrixType res=pos;mi--;
    while(mi)
    {
        if(mi&) res=mul(res,pos);
        pos=mul(pos,pos),mi=mi>>;
//        res.debug();
    }
    return res;
}

int main()
{
    freopen("harem.in","r",stdin);
    freopen("harem.out","w",stdout);
    char ch[];
    MatrixType pos;
    scanf("%d%d",&n,&m);
    pos.n=m,pos.m=m;
    for(int i=;i<=m;i++) pos.ai[][i]=;
    for(int i=;i<=m;i++)
    {
        scanf("%s",ch+),pos.ai[i][]=;
        for(int j=;j<=m;j++) pos.ai[i][j]=(ch[j]-'')^;
    }
//    pos.debug();
    MatrixType p;
    p.mem(),p.n=m,p.m=;
    for(int i=;i<=m;i++) p.ai[i][]=;
//    p.debug();
    MatrixType ans=mul(poww(n-,pos),p);
//    ans.debug();
    ll ans_=;
    for(int i=;i<=m;i++) ans_=(ans_+ans.ai[i][])%mod;
    cout<<ans_;
    return ;
}