hdu 1590（Recursive sequence）矩阵快速幂

来源：互联网发布：拍照软件里有小动物编辑：程序博客网时间：2024/05/29 14:50

题目链接：点击打开链接

第一次做快速幂的题目，题目很简单，就当做个笔记，题解就难得写了~~~

#include <iostream>#include <algorithm>#include <cstdio>#include <cstring>using namespace std;typedef long long LL;const LL fast_mi[7][7]={{1,4,6,4,1,0,0},{0,1,3,3,1,0,0},{0,0,1,2,1,0,0},{0,0,0,1,1,0,0},{0,0,0,0,1,0,0},{0,0,0,0,0,0,1},{1,4,6,4,1,2,1}};struct Martix{    const static int SIZE=7;    const static LL MOD=2147493647;    LL m[SIZE][SIZE];    Martix() {memset(m,0,sizeof(m));}    Martix(const Martix& r) {memcpy(m,r.m,sizeof(r.m));}    Martix(const int);    Martix(const LL (*p)[SIZE]);    Martix operator*(const Martix&)const;    Martix fast_pow(int n)const;};Martix::Martix(const int r){    memset(m,0,sizeof(m));    for(int i=0;i<SIZE;i++)        m[i][i]=r;}Martix::Martix(const LL (*p)[SIZE]){    for(int i=0;i<SIZE;i++)        for(int j=0;j<SIZE;j++)            m[i][j]=p[i][j];}Martix Martix::operator*(const Martix& r)const{    Martix ret;    for(int i=0;i<SIZE;i++)        for(int j=0;j<SIZE;j++)            for(int k=0;k<SIZE;k++)                ret.m[i][j]=(ret.m[i][j]+m[i][k]*r.m[k][j])%MOD;    return ret;}Martix Martix::fast_pow(int n)const{    Martix ret(1),t(*this);    int m=n;    while(m)    {        if(m&1) ret=ret*t;        t=t*t;        m>>=1;    }    return ret;}int main(){    const Martix sub(fast_mi);    const LL MOD=2147493647;    int t;    scanf("%d",&t);    while(t--)    {        int n,a,b;        scanf("%d%d%d",&n,&a,&b);        const int seq[]={16,8,4,2,1,a,b};        Martix ret=sub.fast_pow(n-2);        LL ans=0;        for(int i=0;i<7;i++)            ans=(ans+ret.m[6][i]*seq[i])%MOD;        printf("%lld\n",ans);    }    return 0;}

阅读全文

0 0