【BZOJ】4002: [JLOI2015]有意义的字符串
题意
求$\left \lfloor \left( \frac{b+\sqrt{d}}{2} \right)^n \right \rfloor \pmod {7528443412579576937} \(,\)\left( 0 \le n \le 10^{18}, 0 < b^2 \le d < (b+1)^2 \le 10^{18}, b \mbox{ mod } 2 = 1, d \mbox{ mod } 4=1 \right) $
分析
发现这个并不好算,而如果是\(\left( \frac{b-\sqrt{d}}{2} \right)^n\)那么就好算了。于是又想到数列的特征方程得到的解\(a_n = c_1 x_1^n + c_2 x_2^n\),于是我们搞搞。直接将\(c_1 = c_2 = 1\),则变成\(a_n = x_1^n + x_2^n\),而我们知道\(x_1、x_2\)是特征方程的两个解,和上面那个形式极为相似,于是我们继续假设。即\(x_1 = \left( \frac{b+\sqrt{d}}{2} \right), x_2 = \left( \frac{b-\sqrt{d}}{2} \right)\)。则\(a_{n+2} = pa_{n+1} + qa_{n}\)中,\(p = x_1 + x_2, q = - x_1 x_2\),因此得到\(a_{n+2} = ba_{n+1} - \frac{b^2-d}{4}a_{n}\)
题解
根据上面这个递推式,我们容易算出其中两项,容易得到\(a_1 = b, a_2 = \frac{b^2+d}{2}\)。而发现通项求出来的是整数,因此我们用矩阵乘法求出\(a_n\)即可。最后再根据条件特判一下\(\left( \frac{b-\sqrt{d}}{2} \right)^n\)即可,即\(ans = a_n - [b^2 \neq d \land n是偶数]\)
注意n=0要特判...
#include <bits/stdc++.h>
using namespace std;
typedef unsigned long long ll;
typedef ll mtx[2][2];
const ll mo=7528443412579576937ull, Lim=1e9;
inline void CK(ll &c) {
if(c>=mo)
c-=mo;
}
inline ll mul(ll a, ll b) {
if(a<=Lim && b<=Lim) {
return a*b;
}
if(a<b) {
swap(a, b);
}
ll c=0;
for(; b; b>>=1, CK(a<<=1)) {
if(b&1) {
CK(c+=a);
}
}
return c;
}
void mul(mtx a, mtx b, mtx c, int la, int lb, int lc) {
static mtx t;
memset(t, 0, sizeof t);
for(int i=0; i<la; ++i) {
for(int j=0; j<lc; ++j) {
for(int k=0; k<lb; ++k) {
CK(t[i][j]+=mul(a[i][k], b[k][j]));
}
}
}
memcpy(c, t, sizeof t);
}
ll b, d, n;
bool spj(ll n) {
if(n==1) {
printf("%lld\n", (ll)((((double)b+sqrt(d))/2.0)));
}
else if(n==2) {
printf("%lld\n", (b*b+d)/2);
}
return n<=2;
}
mtx a, c;
int main() {
scanf("%lld%lld%lld", &b, &d, &n);
if(spj(n)) {
return 0;
}
ll t1=b, t2=(d-b*b)/4;
CK(t1), CK(t2);
a[0][0]=t1, a[0][1]=1;
a[1][0]=t2, a[1][1]=0;
c[0][0]=c[1][1]=1;
for(ll tt=n-2; tt; tt>>=1, mul(a, a, a, 2, 2, 2)) {
if(tt&1) {
mul(c, a, c, 2, 2, 2);
}
}
ll a2=(b*b+d)/2, a1=b;
CK(a1), CK(a2);
ll ans;
CK(ans=mul(a2, c[0][0])+mul(a1, c[1][0]));
if(b*b!=d && (n&1)==0) {
if(ans==0) {
ans=mo-1;
}
else {
ans--;
}
}
printf("%llu\n", ans);
return 0;
}
【BZOJ】4002: [JLOI2015]有意义的字符串的更多相关文章
- bzoj 4002: [JLOI2015]有意义的字符串
这个题... #include <bits/stdc++.h> #define rep(i, a, b) for (int i = a; i <= b; i++) #define d ...
- [JLOI2015]有意义的字符串
4002: [JLOI2015]有意义的字符串 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 1000 Solved: 436[Submit][St ...
- 【BZOJ4002】[JLOI2015]有意义的字符串(数论,矩阵快速幂)
[BZOJ4002][JLOI2015]有意义的字符串(数论,矩阵快速幂) 题面 BZOJ 洛谷 题解 发现我这种题总是做不动... 令\(A=\frac{b+\sqrt d}{2},B=\frac{ ...
- BZOJ_4002_[JLOI2015]有意义的字符串_矩阵乘法
BZOJ_4002_[JLOI2015]有意义的字符串_矩阵乘法 Description B 君有两个好朋友,他们叫宁宁和冉冉.有一天,冉冉遇到了一个有趣的题目:输入 b;d;n,求 Input 一行 ...
- 【BZOJ4002】[JLOI2015]有意义的字符串 数学
[BZOJ4002][JLOI2015]有意义的字符串 Description B 君有两个好朋友,他们叫宁宁和冉冉.有一天,冉冉遇到了一个有趣的题目:输入 b;d;n,求 Input 一行三个整数 ...
- BZOJ4002 [JLOI2015]有意义的字符串
据说这两场加起来只要170= =而这是最简单的题目了QAQ 看到$(\frac {b + \sqrt {d} } {2} )^n$,第一反应是共轭根式$(\frac {b - \sqrt {d} } ...
- 【BZOJ4002】[JLOI2015]有意义的字符串 - 矩阵乘法
题意: 给出b,d,n,求$\lfloor(\frac{b+\sqrt{d}}{2})^n\rfloor \mod 999999999999999989$(原题是7528443412579576937 ...
- bzoj4002 [JLOI2015]有意义的字符串 特征根+矩阵快速幂
题目传送门 https://lydsy.com/JudgeOnline/problem.php?id=4002 题解 神仙题. 根据下面的一个提示: \[ b^2 \leq d \leq (b+1)^ ...
- [BZOJ4002][JLOI2015]有意义的字符串-[快速乘法+矩阵乘法]
Description 传送门 Solution 由于这里带了小数,直接计算显然会爆掉,我们要想办法去掉小数. 而由于原题给了暗示:b2<=d<=(b+1)2,我们猜测可以利用$(\fra ...
随机推荐
- client offset screen 的区别
clientX 设置或获取鼠标指针位置相对于窗口客户区域的 x 坐标,其中客户区域不包括窗口自身的控件和滚动条. clientY 设置或获取鼠标指针位置相对于窗口客户区域的 y 坐标,其中客户区域不包 ...
- BZOJ 3784: 树上的路径
Description 问一棵树上前 \(k\) 大路径的边权. Sol 边分治. 非常感谢数据没有菊花图. 为了写写边分治试试然后就开了这道题. 边分治非常好想,选一条重边,分成两部分,然后分别求最 ...
- sublime text3 常用插件安装
1.Package Control 按Ctrl+~调出console(或者view>show console) 粘贴以下代码到底部命令行并回车: import urllib.request,os ...
- MySql: 常见sql语句
1. show create table mysql> show create table t \G*************************** 1. row ************ ...
- 将数据导入Excel
/** * 查询未打印订单 * @param req * @param sort * @param order * @param rows * @param page * @return */ pub ...
- 多功能弹窗控件layer
开发网站的时候,如何合理运用好各种插件对开发的帮助是很大的. 免去了我们调试各种交互效果, 比如常用的弹窗.气泡.提示.加载.焦点.标签.导航.折叠等等 这里会推荐几个常用的js插件,丰富多样简单易移 ...
- JS中的if和else的用法以及基础语法
正常里的变量方式.var a = 10; 针对整数.var b = 3.14; 针对的小数点.var c = "你好":双引号或者单引号引起来的是定义字符串. 一.类型转换(强制转 ...
- PHP之session与cookie
1.session与cookie的关系 众所周知,session是存储在服务器端,cookie是存储在客户端,如果禁用了浏览器的cookie功能,很多时候(除非进行了特殊配置)服务器端就无法再读取se ...
- spring加载过程,源码带你理解从初始化到bean注入
spring在容器启动时,容器正式初始化入口refresh()如下图 ①包括初始化FactoryBean.解析XML注册所有BeanDefinition信息 ②包括注册scope管理类 ③初始化单 ...
- ViewPager适配器FragmentStatePagerAdapter 与FragmentPagerAdapter
使用FragmentPagerAdapter存在删除dataSet顺序错乱的问题 改用FragmentStatePagerAdapter