【BZOJ2326】【HNOI2011】数学作业
2016-04-03 00:04
344 查看
【Description】
【Solution】
矩阵快速幂。递推式:⎡⎣⎢10k00110111⎤⎦⎥⎡⎣⎢f[i−1]i−11⎤⎦⎥=⎡⎣⎢f[i]i1⎤⎦⎥
代码如下:
/************************************************************** Problem: 2326 User: llgyc Language: Pascal Result: Accepted Time:56 ms Memory:224 kb ****************************************************************/ type matrix = array[1..3,1..3] of qword; var n,m,base:qword; ans,t:matrix; i:longint; function mul(x,y:qword):qword; var tmp:qword; begin tmp:=0; while (y>0) do begin if (y mod 2=1) then tmp:=(tmp+x) mod m; x:=(x*2) mod m; y:=y div 2; end; exit(tmp); end; procedure multi(var a,b,c:matrix); var i,j,k:longint; tmp:matrix; begin fillchar(tmp,sizeof(tmp),0); for i:=1 to 3 do for j:=1 to 3 do for k:=1 to 3 do tmp[i,j]:=(tmp[i,j]+mul(a[i,k],b[k,j])) mod m; c:=tmp; end; procedure cal(base,last:qword); var tmp:qword; i,j:longint; begin fillchar(t,sizeof(t),0); t[1,1]:=base; t[1,2]:=1; t[1,3]:=1; t[2,2]:=1; t[2,3]:=1; t[3,3]:=1; tmp:=last-base div 10+1; while (tmp>0) do begin if (tmp mod 2=1) then multi(t,ans,ans); multi(t,t,t); tmp:=tmp div 2; end; end; begin readln(n,m); for i:=1 to 3 do ans[i,i]:=1; base:=10; while (n>=base) do begin cal(base,base-1); base:=base*10; end; cal(base,n); writeln(ans[1,3]); end.
相关文章推荐
- Tcp的通信的简单程序
- TCP服务器和多个客户端传输数据
- 75道程序员面试逻辑思维题
- ORA-20000的解决过程
- 2013 黄金连分数
- 计算一个矩阵从左上角到右下角和最大
- 《Linux内核分析》 第六节 分析Linux内核创建一个新进程的过程
- PHP错误处理函数
- PHP中的Trait 特性及作用
- PHP命名空间和自动加载类
- PHP preg match正则表达式函数的操作实例
- js一维数组、多维数组和对象的混合使用方法
- 基于JavaScript实现 网页切出 网站title变化代码
- 深入浅析JSON.parse()、JSON.stringify()和eval()的作用详解
- AngularJs 60分钟入门基础教程
- Node.js 应用跑得更快 10 个技巧
- jQuery qrcode生成二维码的方法
- 再次谈论React.js实现原生js拖拽效果引起的一系列问题
- jQuery中的Deferred和promise 的区别
- Angular 根据 service 的状态更新 directive