我们看到了及其可怕的数据范围

这个样子都没有办法直接读入的数据范围应该怎么算

我们观察一下递推式\(f[i][j]=a*f[i][j]+b(j!=1)\)

\(f[i][1]=c*f[i-1][m]+d\)

转移非常简单,于是可以考虑一下矩阵乘法

如果我们将这个矩阵破坏成一个链,那么就会有这种形式的递推

连续推\(m\)次第一个柿子,之后再推一次第二个柿子,之后反复

重复上面的过程\(n\)次就好了

于是我们可以将连续转移\(m\)次一式的到的矩阵和第二个式子的转移矩阵乘起来,之后将这个矩阵再转移\(n\)次就是答案了

由于\(n,m\)是在太大了,我们又发现模数是质数,于是我们可以利用费马小定理来降幂

则有

\[A^{mod-1}\%mod=1
\]

所以

\[A^m\equiv A^{m\%(mod-1)}(\%\ mod )
\]

至于为什么要特判\(a=1\)

我怎么知道啊

#include<iostream>
#include<cstring>
#include<cstdio>
#define re register
#define maxn 5
#define LL long long
const LL mod=1000000007;
char N_[1000005],M_[1000005];
LL n,m,a,b,c,d;
inline LL read()
{
char c=getchar();
LL x=0;
while(c<'0'||c>'9') c=getchar();
while(c>='0'&&c<='9')
x=x*10%(mod-1)+c-48,c=getchar();
return x;
}
struct Mat
{
LL a[4][4],ans[4][4];
inline void did_ans()
{
LL mid[4][4];
for(re int i=1;i<=2;i++)
for(re int j=1;j<=2;j++)
mid[i][j]=ans[i][j],ans[i][j]=0;
for(re int i=1;i<=2;i++)
for(re int j=1;j<=2;j++)
for(re int p=1;p<=2;p++)
ans[i][j]=(ans[i][j]+(a[i][p]*mid[p][j])%mod)%mod;
}
inline void did_a()
{
LL mid[4][4];
for(re int i=1;i<=2;i++)
for(re int j=1;j<=2;j++)
mid[i][j]=a[i][j],a[i][j]=0;
for(re int i=1;i<=2;i++)
for(re int j=1;j<=2;j++)
for(re int p=1;p<=2;p++)
a[i][j]=(a[i][j]+(mid[i][p]*mid[p][j])%mod)%mod;
}
}R,L,K;
inline void mul(Mat &A,Mat &B)
{
LL mid[4][4];
for(re int i=1;i<=2;i++)
for(re int j=1;j<=2;j++)
mid[i][j]=A.ans[i][j],A.ans[i][j]=0;
for(re int i=1;i<=2;i++)
for(re int j=1;j<=2;j++)
for(re int p=1;p<=2;p++)
A.ans[i][j]=(A.ans[i][j]+(mid[i][p]*B.ans[p][j])%mod)%mod;
for(re int i=1;i<=2;i++)
for(re int j=1;j<=2;j++)
A.a[i][j]=A.ans[i][j];
}
int main()
{
scanf("%s%s",N_+1,M_+1);
a=read(),b=read(),c=read(),d=read();
int lenn=strlen(N_+1);
if(a==1)
{
for(re int i=1;i<=lenn;i++) n=n*10%mod+N_[i]-48;
n=(n-2+mod)%mod;
}
else
{
for(re int i=1;i<=lenn;i++) n=n*10%(mod-1)+N_[i]-48;
n=(n-2+mod-1)%(mod-1);
}
lenn=strlen(M_+1);
if(c==1)
{
for(re int i=1;i<=lenn;i++) m=m*10%mod+M_[i]-48;
m=(m-2+mod)%mod;
}
else
{
for(re int i=1;i<=lenn;i++) m=m*10%(mod-1)+M_[i]-48;
m=(m-2+mod-1)%(mod-1);
}
R.a[1][1]=R.ans[1][1]=1;
R.a[2][1]=R.ans[2][1]=b;
R.ans[2][2]=R.a[2][2]=a;
R.ans[1][2]=R.a[1][2]=0;
K=R;
while(m)
{
if(m&1) R.did_ans();
m>>=1;
R.did_a();
}
L.a[1][1]=L.ans[1][1]=1;
L.a[2][1]=L.ans[2][1]=d;
L.a[2][2]=L.ans[2][2]=c;
L.ans[1][2]=L.a[1][2]=0;
mul(L,R);
while(n)
{
if(n&1) L.did_ans();
n>>=1;
L.did_a();
}
printf("%lld",((L.ans[2][1]+L.ans[2][2])%mod*R.ans[2][2]%mod+R.ans[2][1])%mod);
return 0;
}

【[NOI2013]矩阵游戏】的更多相关文章

  1. bzoj 3240: [Noi2013]矩阵游戏 矩阵乘法+十进制快速幂+常数优化

    3240: [Noi2013]矩阵游戏 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 613  Solved: 256[Submit][Status] ...

  2. BZOJ 3240: [Noi2013]矩阵游戏

    3240: [Noi2013]矩阵游戏 Time Limit: 10 Sec  Memory Limit: 256 MBSubmit: 1586  Solved: 698[Submit][Status ...

  3. BZOJ 3240([Noi2013]矩阵游戏-费马小定理【矩阵推论】-%*s-快速读入)

    3240: [Noi2013]矩阵游戏 Time Limit: 10 Sec   Memory Limit: 256 MB Submit: 123   Solved: 73 [ Submit][ St ...

  4. (十进制高速幂+矩阵优化)BZOJ 3240 3240: [Noi2013]矩阵游戏

    题目链接: http://www.lydsy.com/JudgeOnline/problem.php?id=3240 3240: [Noi2013]矩阵游戏 Time Limit: 10 Sec  M ...

  5. P1397 [NOI2013]矩阵游戏(递推)

    P1397 [NOI2013]矩阵游戏 一波化式子,$f[1][m]=a^{m-1}+b\sum_{i=0}^{m-2}a^i$,用快速幂+逆元求等比数列可以做到$logm$ 设$v=a^{m-1}, ...

  6. 【bzoj3240】 Noi2013—矩阵游戏

    http://www.lydsy.com/JudgeOnline/problem.php?id=3240 (题目链接) 题意$${F[1][1]=1}$$$${F[i][j]=a*F[i][j-1]+ ...

  7. BZOJ3240 [Noi2013]矩阵游戏

    本文版权归ljh2000和博客园共有,欢迎转载,但须保留此声明,并给出原文链接,谢谢合作. 本文作者:ljh2000作者博客:http://www.cnblogs.com/ljh2000-jump/转 ...

  8. 3240: [Noi2013]矩阵游戏

    Description 婷婷是个喜欢矩阵的小朋友,有一天她想用电脑生成一个巨大的n行m列的矩阵(你不用担心她如何存储).她生成的这个矩阵满足一个神奇的性质:若用F[i][j]来表示矩阵中第i行第j列的 ...

  9. NOI2013矩阵游戏

    Description 婷婷是个喜欢矩阵的小朋友,有一天她想用电脑生成一个巨大的n行m列的矩阵(你不用担心她如何存储).她生成的这个矩阵满足一个神奇的性质:若用F[i][j]来表示矩阵中第i行第j列的 ...

  10. 洛谷P1397 [NOI2013]矩阵游戏

    矩阵快速幂+费马小定理 矩阵也是可以跑费马小定理的,但是要注意这个: (图是盗来的QAQ) 就是说如果矩阵a[i][i]都是相等的,那么就是mod p 而不是mod p-1了 #include< ...

随机推荐

  1. vi命令——修改文件内容

    以下内容转载自http://blog.csdn.net/bruce0532/article/details/7842384 vi编辑器是所有Unix及Linux系统下标准的编辑器,介绍一下它的用法和一 ...

  2. python 生成嵌套字典

    import collections import json tree=lambda:collections.defaultdict(tree) some_dict=tree() some_dict[ ...

  3. IE中使用TerraExplorerPro ActiveX控件问题总结

    近段时间,由于工作需要,涉及到Skyline Web的二次开发,遇到一些问题: 问题1:IE中无法使用Skyline 提供的ActiveX Component 的API 解决方法:这个问题一般是没有正 ...

  4. 深入理解JavaScript系列(2):揭秘命名函数表达式

    前言 网上还没用发现有人对命名函数表达式进去重复深入的讨论,正因为如此,网上出现了各种各样的误解,本文将从原理和实践两个方面来探讨JavaScript关于命名函数表达式的优缺点. 简单的说,命名函数表 ...

  5. Win2D 官方文章系列翻译 - 预乘 Alpha

    本文为个人博客备份文章,原文地址: http://validvoid.net/win2d-premultiplied-alpha/ 在计算机绘图中有两种表示颜色值不透明度的方法.Win2D 中两种方法 ...

  6. FontSize sp 和 dp 的区别

    dp不会随着“设置->显示->字体大小”的改变而改变,sp会. sp会随着configeration的配置来scale, dp不会. 所以,什么时候用sp, 什么时候用dp需要斟酌.

  7. Implementation: Quick Sort 2014-08-19

    #include <stdio.h> void print(int *a, int start , int end); void quick_sort(int *a, int start, ...

  8. document.write()重写问题

    document.write(); 可用于重写给某个元素追加内容; 当document.write(); 用于JS文件中,会重写整个页面,解决这个问题有多种方法. 重写原因:当onload的时候执行 ...

  9. iOS中使用RNCryptor对资源文件加密

    原文:http://blog.csdn.net/chenpolu/article/details/46277587 RNCryptor源码https://github.com/RNCryptor/RN ...

  10. linux c 遍历目录及文件

    #include <dirent.h>void recovery_backend() { DIR * pdir ; struct dirent * pdirent; struct stat ...