牛客多校第五场 B generator 1 矩阵快速幂
题意:
给定$x_0,x_1,a,b,n,mod, x_i=a*x_{i-1}+b*x_{i-2}$ ,求$x_n % mod$
n最大有1e6位
题解:
矩阵快速幂。
巨大的n并不是障碍,写一个十进制的矩阵快速幂就行了。
$ \begin{bmatrix}x_n \\ x_{n-1} \end{bmatrix}=\begin{bmatrix}a &b \\ 1 &0 \end{bmatrix} *\begin{bmatrix}x_{n-1} \\ x_{n-2} \end{bmatrix}=\begin{bmatrix}a &b \\ 1 &0 \end{bmatrix}^{n-1}\begin{bmatrix}x_1 \\ x_0 \end{bmatrix} $
#include<iostream>
#include<cstring>
#include<cassert>
#define LL long long
using namespace std;
LL mod;
struct Mtx{
LL x[][];
friend Mtx operator *(const Mtx &a,const Mtx &b){
Mtx c;
LL tmp00=a.x[][]*b.x[][] % mod+a.x[][]*b.x[][] % mod;
LL tmp01=a.x[][]*b.x[][] % mod+a.x[][]*b.x[][] % mod;
LL tmp10=a.x[][]*b.x[][] % mod+a.x[][]*b.x[][] % mod;
LL tmp11=a.x[][]*b.x[][] % mod+a.x[][]*b.x[][] % mod;
assert(tmp00>=);
assert(tmp01>=);
assert(tmp10>=);
assert(tmp11>=);
c.x[][]=tmp00 % mod;
c.x[][]=tmp01 % mod;
c.x[][]=tmp10 % mod;
c.x[][]=tmp11 % mod;
return c;
}
friend Mtx operator ^ (Mtx base,int n){
Mtx ans=Mtx();
while(n){
if(n%){
ans=ans*base;
}
base=base*base;
n>>=;
}
return ans;
}
Mtx(){}
Mtx(int a){
x[][]=;
x[][]=;
x[][]=;
x[][]=;
}
};
char n[];
int main(){
LL x0,x1,a,b;
scanf("%lld %lld %lld %lld",&x0,&x1,&a,&b);
scanf("%s",n);
scanf("%lld",&mod);
int len=strlen(n);
int kk=len-;
Mtx ans=Mtx();
Mtx base;
if(len==){
if(n[]==''){
printf("%lld\n",x0%mod);
goto E;
}
if(n[]==''){
printf("%lld\n",x1%mod);
goto E;
}
} while(){
n[kk]--;
if(n[kk]<''){
n[kk]+=;
kk--;
}else break;
}
// printf("%s\n",n);
base.x[][]=a;
base.x[][]=b;
base.x[][]=;
base.x[][]=;
for(int i=;i<len;i++){
ans=ans^;
ans=ans*(base^(n[i]-''));
}
printf("%lld\n",(x1*ans.x[][]+x0*ans.x[][])%mod); E:return ;
}
小贴士:矩阵快速幂,以及一些其他的,比较复杂的,比较套路的东西,一定要封装好,这样在不太损失效率的前提下最大限度的保证了代码的可读性,这是我某次打cf时wa到自闭的教训。
8月6日PS:这种概念叫做DRY原则,全称是Don't repeat yourself.今天刚学到的名词。
牛客多校第五场 B generator 1 矩阵快速幂的更多相关文章
- 2019牛客多校第五场 B - generator 1 矩阵快速幂+十倍增+二进制倍增优化
B - generator 1 题意 给你\(x_{0}.x_{1}.a.b.b.mod\),根据\(x_{i} = a*x_{i-1} + b*x_{i-2}\)求出\(x_{n}\) 思路 一般看 ...
- generator 1(2019年牛客多校第五场B题+十进制矩阵快速幂)
目录 题目链接 思路 代码 题目链接 传送门 思路 十进制矩阵快速幂. 代码 #include <set> #include <map> #include <deque& ...
- 十进制快速幂(牛客多校第五场)-- generator 1
思路: 十进制快速幂. #include <stdio.h>//sprintf #include <cstdlib>////malloc exit strcat itoa sy ...
- 2019 牛客多校第五场 B generator 1
题目链接:https://ac.nowcoder.com/acm/contest/885/B 题目大意 略. 分析 十进制矩阵快速幂. 代码如下 #include <bits/stdc++.h& ...
- 2019牛客多校第五场C generator 2 hash,bsgs模板
generator 2 题意 给出\(x_0,a,b,p\),有方程\(x_i\equiv (a*x_{i-1}+b)(\% p)\),求最小的i,使得\(x_i=v\),不存在输出-1 分析 经过公 ...
- 2019牛客多校第五场B generator 十进制快速幂
generator 1 题意 给出\(x_0,x_1,a,b\)已知递推式\(x_i=a*x_{i-1}+b*x_{i-2}\),出个n和mod,求\(x_n\) (n特别大) 分析 比赛的时候失了智 ...
- 2019牛客多校第五场C generator 2(BSGS)题解
题意: 传送门 已知递推公式\(x_i = a*x_{i - 1} + b\mod p\),\(p\)是素数,已知\(x_0,a,b,p\),给出一个\(n\)和\(v\),问你满足\(x_i = v ...
- 牛客多校第五场 F take
链接:https://www.nowcoder.com/acm/contest/143/F来源:牛客网 题目描述 Kanade has n boxes , the i-th box has p[i] ...
- 牛客多校第五场 J:Plan
链接:https://www.nowcoder.com/acm/contest/143/J 来源:牛客网 时间限制:C/C++ 1秒,其他语言2秒 空间限制:C/C++ 262144K,其他语言524 ...
随机推荐
- 2019-9-2-win10-uwp-布局
title author date CreateTime categories win10 uwp 布局 lindexi 2019-09-02 12:57:38 +0800 2018-2-13 17: ...
- linux安装lolcat实现彩色文字输出信息
[root@localhost ~]# mount /dev/sr0 /media/[root@localhost ~]# rpm -ivh epel-release-latest-7.noarch. ...
- python_django_在views模块中操作状态保持(session)
什么叫状态保持? 就比如说我们登陆一个网站,登陆之后的当前页面显示的是登陆状态,但是我们要再跳转同一网站的其他页面,则显示的未登录状态,状态保持就是:我们在当前页面登陆后,再访问其他页面时也显示为登陆 ...
- net core静态文件 访问除默认目录文件配置
在我们项目的实际应用中,不光是需要访问默认静态文件夹 wwwroot ,还有可能要要去访问除默认目录以外的文件夹,接下来我们进行配置 1.在根目录创建一个文件夹,继续创建它的子文件夹Images,在I ...
- vue+elementui 中 @keyup 键盘上下左右移动聚焦
<template> <el-table :data="CreditUnclearOutlineList" border style="width: 1 ...
- delphi常见的错误
******************************* * 编 译 错 误 信 息 * ******************************* ';' not allowed befo ...
- 莫比乌斯反演+二维前缀和——hdu4746二刷
第二次做这题,求前缀和的时候还是卡住了 fg函数的反演是可以直接用莫比乌斯基本代换式来代换的 #include<bits/stdc++.h> using namespace std; #d ...
- 自定义alert 确定、取消功能
以删除为例,首先新建html <table border="1" cellpadding="0" cellspacing="0" id ...
- 标准H.460公私网穿越视频解决方案
一.概述 H.460协议是一种网络通信协议,主要用于音视频的网络穿越,可以解决客户私网到公网,以及公网到私网的互相通信. 大连羽化集团是中国较大的零售业集团之一.目前羽化集团百货店已达20多家,营业面 ...
- ubantu 编译mysql++
sudo apt-get install mysql-server mysql-client 指定mysql-lib位置: ./configure --with-mysql-lib=/usr/lib/ ...