注册 登录  
 加关注
查看详情
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

n+e

NewWeb:http://trinkle.is-programmer.com/

 
 
 

日志

 
 

[BZOJ3240][NOI2013]矩阵游戏  

2014-07-20 21:29:59|  分类: BZOJ |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |

Description

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

F[1][1]=1
F[i,j]=a*F[i][j-1]+b (j!=1)
F[i,1]=c*F[i-1][m]+d (i!=1)
递推式中a,b,c,d都是给定的常数。

现在婷婷想知道F[n][m]的值是多少,请你帮助她。由于最终结果可能很大,你只需要输出F[n][m]除以1,000,000,007的余数。

Input

一行有六个整数n,m,a,b,c,d。意义如题所述

Output

包含一个整数,表示F[n][m]除以1,000,000,007的余数

Sample Input

3 4 1 3 2 6

Sample Output

85

HINT

样例中的矩阵为:

1 4 7 10

26 29 32 35

76 79 82 85


1<=n,m<=10^1000000,1<=a,b,c,d<=10^9


Solution

  还是矩阵快速幂,这回用到费马小定理,直接把n,m降成10^9级别。

  注意特判a=c=1

  我不知道为何别人的代码比我不缩行还长2倍以上= =

/**************************************************************
Problem: 3240
User: wjy1998
Language: C++
Result: Accepted
Time:332 ms
Memory:2760 kb
****************************************************************/

#include<cstdio>
typedef long long ll;struct M{ll x,y;}t1,t2,t3;ll P=1e9+7,a,b,c,d,n,m,phi;
char s1[1000010],s2[1000010];
void F(char*s,ll&aa){for(int t=0;s[t];t++)aa=(aa*10+s[t]-'0')%phi;}
M operator*(const M&a,const M&b){return (M){a.x*b.x%P,(a.x*b.y+a.y)%P};}
M operator^(M t,ll k){M f=t;for(--k;k;k>>=1,t=t*t)if(k&1)f=f*t;return f;}
int main(){
scanf("%s%s%lld%lld%lld%lld",s1,s2,&a,&b,&c,&d);
if(a==1&&c==1)phi=P;else phi=P-1;
F(s1,n),F(s2,m),t1=(M){a,b},t2=(M){c,d};t1=t1^(m-1),t3=t2*t1;t3=t3^(n-1);t1=t1*t3;
printf("%lld\n",(t1.x+t1.y)%P);
}


  评论这张
 
阅读(21)| 评论(0)
推荐 转载

历史上的今天

在LOFTER的更多文章

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2018