【模板】快速幂

普通的快速幂

#include<stdio.h>
#define lll long long
lll ksm(int,int);
int x,n;
int main()
{
    scanf("%d%d",&x,&n);
    printf("%lld",ksm(x,n));
    return 0;
}
lll ksm(int x,int n)
{
    lll ret = x; n--;
    while(n)
    {
        if(n & 1) ret *= x;
        x *= x;
        n >>= 1;
    }
    return ret;
}

 矩阵快速幂

#include<stdio.h> 
#define maxn 105
#define mo 1000000007
#define lll long long
int n,m;
struct mat
{
    lll ori[maxn][maxn];
    mat operator * (const mat &y) const
    {
        mat ret;
        int i,j,k;
        for(i=1;i<=n;i++)
            for(j=1;j<=n;j++)
            {
                ret.ori[i][j] = 0;
                for(k=1;k<=n;k++)
                    ret.ori[i][j] = (ret.ori[i][j]+ori[i][k]*y.ori[k][j])%mo;
            }
        return ret;
    }
};
mat A,ans;
mat ksm (int k,mat a)
{
    mat ret = a; k--;
    while(k)
    {
        if(k&1) ret = ret * a;
        a = a*a;
        k >>= 1; 
    }
    return ret;
}
int main()
{
    scanf("%d%d",&n,&m);
    int i,j;
    for(i=1;i<=n;i++)
        for(j=1;j<=n;j++)
            scanf("%d",&A.ori[i][j]);
    ans = ksm(m,A);
    for(i=1;i<=n;i++)
    {
        for(j=1;j<=n;j++) printf("%lld ",ans.ori[i][j]);
        printf("
");
    }
    return 0;
}