2019牛客暑期多校训练营(第五场)- B generator 1 (齐次线性递推+矩阵快速幂)
时间:2019-08-17
本文章向大家介绍2019牛客暑期多校训练营(第五场)- B generator 1 (齐次线性递推+矩阵快速幂),主要包括2019牛客暑期多校训练营(第五场)- B generator 1 (齐次线性递推+矩阵快速幂)使用实例、应用技巧、基本知识点总结和需要注意事项,具有一定的参考价值,需要的朋友可以参考一下。
题目链接:https://ac.nowcoder.com/acm/contest/885/B
题意:已知齐次线性式xn=a*xn-1+b*xn-2,已知a,b,x0,x1,求xn,n很大,n<=1010^6.
思路:矩阵快速幂模板题,构造矩阵t:
a b 1 0
矩阵ans:
x1 0 x0 0
显然ans1=t×ans,ans1为:
x2 0 x1 0
那么ansn=t^n*ans,ansn为:
xn+1 0 xn 0
所以用矩阵快速幂计算t^n,n很大,快速幂要用十进制倍增,对每一位的计算不能直接乘,还要用二进制的快速幂,不然会TLE。
最近写代码像中邪了一样,照着别人的代码写总会有个小地方写错,然后找一天的错误,连续几天了。
AC代码:
#include<cstdio> #include<algorithm> #include<cstring> using namespace std; typedef long long LL; LL x0,x1,a,b; int MOD; char s[1000005]; struct Matrix{ LL a[2][2]; Matrix(){memset(a,0,sizeof(a));} Matrix operator * (const Matrix y){ Matrix ret; for(int i=0;i<2;++i) for(int j=0;j<2;++j) for(int k=0;k<2;++k) ret.a[i][j]=(ret.a[i][j]+a[i][k]*y.a[k][j])%MOD; return ret; } void operator = (const Matrix y){ for(int i=0;i<2;++i) for(int j=0;j<2;++j) a[i][j]=y.a[i][j]; } }tmp; Matrix qpow2(Matrix m,int b){ Matrix ret=tmp; ret.a[0][0]=ret.a[1][1]=1; while(b){ if(b&1) ret=ret*m; m=m*m; b>>=1; } return ret; } Matrix qpow10(Matrix m,char *s){ Matrix ret=tmp; ret.a[0][0]=ret.a[1][1]=1; for(int i=strlen(s)-1;i>=0;--i){ int num=s[i]-'0'; ret=ret*qpow2(m,num); m=qpow2(m,10); } return ret; } int main(){ scanf("%lld%lld%lld%lld",&x0,&x1,&a,&b); scanf("%s%d",s,&MOD); Matrix t,ans; t.a[0][0]=a,t.a[0][1]=b,t.a[1][0]=1,t.a[1][1]=0; ans.a[0][0]=x1,ans.a[0][1]=0,ans.a[1][0]=x0,ans.a[1][1]=1; t=qpow10(t,s); ans=t*ans; printf("%lld\n",ans.a[1][0]); return 0; }
原文地址:https://www.cnblogs.com/FrankChen831X/p/11370331.html
- JavaScript 教程
- JavaScript 编辑工具
- JavaScript 与HTML
- JavaScript 与Java
- JavaScript 数据结构
- JavaScript 基本数据类型
- JavaScript 特殊数据类型
- JavaScript 运算符
- JavaScript typeof 运算符
- JavaScript 表达式
- JavaScript 类型转换
- JavaScript 基本语法
- JavaScript 注释
- Javascript 基本处理流程
- Javascript 选择结构
- Javascript if 语句
- Javascript if 语句的嵌套
- Javascript switch 语句
- Javascript 循环结构
- Javascript 循环结构实例
- Javascript 跳转语句
- Javascript 控制语句总结
- Javascript 函数介绍
- Javascript 函数的定义
- Javascript 函数调用
- Javascript 几种特殊的函数
- JavaScript 内置函数简介
- Javascript eval() 函数
- Javascript isFinite() 函数
- Javascript isNaN() 函数
- parseInt() 与 parseFloat()
- escape() 与 unescape()
- Javascript 字符串介绍
- Javascript length属性
- javascript 字符串函数
- Javascript 日期对象简介
- Javascript 日期对象用途
- Date 对象属性和方法
- Javascript 数组是什么
- Javascript 创建数组
- Javascript 数组赋值与取值
- Javascript 数组属性和方法