矩阵快速幂模板

#include<stdio.h>
#include<iostream>
#include<string.h>
#include<ctype.h>
#include<math.h>
#include<map>
#include<set>
#include<vector>
#include<queue>
#include<string>
#include<algorithm>
#include<time.h>
#include<bitset>
using namespace std;
void fre(){freopen("c://test//input.in","r",stdin);freopen("c://test//output.out","w",stdout);}
#define MS(x,y) memset(x,y,sizeof(x))
#define MC(x,y) memcpy(x,y,sizeof(x))
#define ls o<<1
#define rs o<<1|1
typedef long long LL;
typedef unsigned long long UL;
typedef unsigned int UI;
typedef pair<int,int>PII;
template <class T> inline void gmax(T &a,T b){if(b>a)a=b;}
template <class T> inline void gmin(T &a,T b){if(b<a)a=b;}
const int N=0,M=0,maxint=2147483647,ms31=522133279,ms63=1061109567,ms127=2139062143;const double eps=1e-8,PI=acos(-1.0);//.0
int casenum,casei;
const int G=10;//矩阵大小
const int Z=1e9+7;//取模数
struct MX
{
    int v[G][G];
    void O(){MS(v,0);}
    void E(){MS(v,0);for(int i=0;i<G;i++)v[i][i]=1;}
    MX operator * (const MX &b) const
    {
        MX c;c.O();
        for(int i=0;i<G;i++)
        {
            for(int j=0;j<G;j++)
            {
                for(int k=0;k<G;k++)
                {
                    c.v[i][j]=(c.v[i][j]+v[i][k]*b.v[k][j])%Z;
                }
            }
        }
        return c;
    }
    MX operator + (const MX &b) const
    {
        MX c;c.O();
        for(int i=0;i<G;i++)
        {
            for(int j=0;j<G;j++)
            {
                c.v[i][j]=(v[i][j]+b.v[i][j])%Z;
            }
        }
        return c;
    }
    MX operator ^ (int p) const
    {
        MX y;y.E();
        MX x;MC(x.v,v);
        while(p)
        {
            if(p&1)y=y*x;
            x=x*x;
            p>>=1;
        }
        return y;
    }
};
MX mul(MX a,MX b,int I,int K,int J)
{
    MX c;c.O();
    for(int i=0;i<I;i++)
    {
        for(int j=0;j<J;j++)
        {
            for(int k=0;k<K;k++)
            {
                c.v[i][j]=(c.v[i][j]+a.v[i][k]*b.v[k][j])%Z;
            }
        }
    }
    return c;
}
int main()
{
    return 0;
}