[日常摸鱼]bzoj2875[NOI2012]随机数生成器-矩阵快速幂

好裸的矩阵快速幂…然而我一开始居然构造不出矩阵…

平常两个的情况都是拿相邻两项放在矩阵里拿去递推的…然后我就一直构造不出来…其实把矩阵下面弄成1就好了啊orz

#include<cstdio>
#include<cstring>
#define rep(i,n) for(register lint i=1;i<=n;i++)
typedef long long lint;
lint m,n,g;
struct matrix
{
    lint m[3][3];
    matrix(){memset(m,0,sizeof(m));}
};
inline lint mul(lint a,lint b)
{
    lint res=0;
    for(;b;b>>=1,a=(a+a)%m)if(b&1)res=(res+a)%m;
    return res;
}
inline matrix operator *(matrix a,matrix b)
{
    matrix res;
    rep(i,2)rep(j,2)rep(k,2)res.m[i][j]=(res.m[i][j]+mul(a.m[i][k],b.m[k][j]))%m;
    return res;
}
inline matrix pow_mod(matrix a,lint b)
{
    matrix res;res.m[1][1]=res.m[2][2]=1;
    for(;b;b>>=1,a=a*a)if(b&1)res=a*res;
    return res;
}
int main()
{
    matrix t,r;
    scanf("%lld%lld%lld%lld%lld%lld",&m,&t.m[1][1],&t.m[1][2],&r.m[1][1],&n,&g);
    r.m[2][1]=t.m[2][2]=1;r=pow_mod(t,n)*r;printf("%lld",r.m[1][1]%g);
    return 0;
}

 

转载于:https://www.cnblogs.com/yoooshinow/p/8329340.html

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值