P1397 [NOI2013]矩阵游戏

2019-09-15 12:03:21 阅读：208 来源： 互联网

标签：ch P1397 矩阵 long while NOI2013 include getchar

首先显然可以矩乘快速幂然后 $T$ 飞

看一眼题解发现因为这一题矩阵的特殊性所以可以对矩阵的次数欧拉降幂

然而我并不懂证明，所以我选择暴力乱搞的做法

十进制快速幂，然后注意一下常数，还有矩阵乘的顺序，别反了

#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cstring>
#include<cmath>
using namespace std;
typedef long long ll;
inline int read()
{
    int x=0,f=1; char ch=getchar();
    while(ch<'0'||ch>'9') { if(ch=='-') f=-1; ch=getchar(); }
    while(ch>='0'&&ch<='9') { x=(x<<1)+(x<<3)+(ch^48); ch=getchar(); }
    return x*f;
}
const int N=1e6+7,mo=1e9+7;
char n[N],m[N];
struct Matrix {
    int a[3][3];
    Matrix () { memset(a,0,sizeof(a)); }
    inline Matrix operator * (const Matrix &tmp) const {
        Matrix res;
        res.a[1][1]=(1ll*a[1][1]*tmp.a[1][1]+1ll*a[1][2]*tmp.a[2][1])%mo;
        res.a[1][2]=(1ll*a[1][1]*tmp.a[1][2]+1ll*a[1][2]*tmp.a[2][2])%mo;
        res.a[2][1]=(1ll*a[2][1]*tmp.a[1][1]+1ll*a[2][2]*tmp.a[2][1])%mo;
        res.a[2][2]=(1ll*a[2][1]*tmp.a[1][2]+1ll*a[2][2]*tmp.a[2][2])%mo;
        //循环展开优化常数
        return res;
    }
}F,G,Ans;
Matrix ksm(Matrix x,char *y)//十进制快速幂
{
    Matrix res,t; res.a[1][1]=res.a[2][2]=1;
    for(int i=strlen(y+1);i;i--)
    {
        t=x;
        if(y[i]=='9') { t=t*t; t=t*t; t=t*t; t=t*x; res=res*t; }
        else if(y[i]=='8') { t=t*t; t=t*t; t=t*t; res=res*t; }
        else for(int j=1;j<=y[i]-'0';j++) res=res*x;
        t=x; x=x*x; x=x*x; x=x*x; x=x*t; x=x*t;
    }
    return res;
}
void Minus(char *s)//把数减一
{
    for(int i=strlen(s+1);i;i--)
        if(s[i]=='0') s[i]='9';
        else { s[i]=s[i]-1; break; }
}
int main()
{
    int a,b,c,d; scanf("%s",n+1); scanf("%s",m+1);
    a=read(),b=read(),c=read(),d=read();
    Ans.a[1][1]=Ans.a[1][2]=1;
    F.a[1][1]=a; F.a[2][1]=b; F.a[2][2]=1;
    G.a[1][1]=c; G.a[2][1]=d; G.a[2][2]=1;
    Minus(n); Minus(m);
    F=ksm(F,m); G=F*G; G=ksm(G,n);
    Ans=Ans*G*F;
    printf("%d\n",Ans.a[1][1]);
    return 0;
}

标签：ch,P1397,矩阵,long,while,NOI2013,include,getchar
来源： https://www.cnblogs.com/LLTYYC/p/11521638.html

本站声明： 1. iCode9 技术分享网（下文简称本站）提供的所有内容，仅供技术学习、探讨和分享；
2. 关于本站的所有留言、评论、转载及引用，纯属内容发起人的个人观点，与本站观点和立场无关；
3. 关于本站的所有言论和文字，纯属内容发起人的个人观点，与本站观点和立场无关；
4. 本站文章均是网友提供，不完全保证技术分享内容的完整性、准确性、时效性、风险性和版权归属；如您发现该文章侵犯了您的权益，可联系我们第一时间进行删除；
5. 本站为非盈利性的个人网站，所有内容不会用来进行牟利，也不会利用任何形式的广告来间接获益，纯粹是为了广大技术爱好者提供技术内容和技术思想的分享性交流网站。

ICode9

P1397 [NOI2013]矩阵游戏