[BZOJ4002][JLOI2015]有意义的字符串(结论+矩阵乘法)
来源:互联网 发布:彩虹秒赞最新源码 编辑:程序博客网 时间:2024/06/06 01:05
首先得出,
把
两边同乘以
设
此时就容易得出
最后考虑怎样通过
注意特判
代码:
#include <cmath>#include <cstdio>#include <cstring>#include <iostream>#include <algorithm>using namespace std;typedef long long ll;ll b, d, n, tm; const ll ZZQ = 7528443412579576937ll;ll add(ll a, ll b) { return (1ull * a + 1ull * b) % ZZQ;}ll prod(ll a, ll b) { ll res = 0; while (b) { if (b & 1) res = add(res, a); a = add(a, a); b >>= 1; } return res;}struct cyx { int n, m; ll v[4][4]; cyx() {} cyx(int _n, int _m) : n(_n), m(_m) {memset(v, 0, sizeof(v));} friend inline cyx operator * (cyx a, cyx b) { int i, j, k; cyx res = cyx(a.n, b.m); for (i = 1; i <= res.n; i++) for (j = 1; j <= res.m; j++) for (k = 1; k <= a.m; k++) res.v[i][j] = add(res.v[i][j], prod(a.v[i][k], b.v[k][j])); return res; } friend inline cyx operator ^ (cyx a, ll b) { int i; cyx res = cyx(a.n, a.m); for (i = 1; i <= res.n; i++) res.v[i][i] = 1; while (b) { if (b & 1) res = res * a; a = a * a; b >>= 1; } return res; }} P, Q;int main() { cin >> b >> d >> n; P = cyx(2, 2); Q = cyx(2, 1); if (!n) return printf("1\n"), 0; tm = (d >> 2) - prod(b + 1 >> 1, b - 1 >> 1); P.v[1][1] = b; P.v[1][2] = tm; P.v[2][1] = 1; Q.v[1][1] = b; Q.v[2][1] = 2; P = (P ^ n - 1) * Q; ll ans = P.v[1][1]; if (d != b * b && !(n & 1)) ans--; if (ans < 0) ans += ZZQ; cout << ans << endl; return 0;}
阅读全文
0 0
- [BZOJ4002][JLOI2015]有意义的字符串(结论+矩阵乘法)
- [BZOJ4002]JLOI2015有意义的字符串|矩阵乘法
- 【BZOJ4002】【JLOI2015】有意义的字符串 推公式+矩阵乘法
- [题解]bzoj4002(JLOI2015)有意义的字符串
- [数列通项 矩阵快速幂] BZOJ4002 [JLOI2015]有意义的字符串
- BZOJ4002 有意义的字符串
- JLOI2015有意义的字符串
- 4002: [JLOI2015]有意义的字符串
- bzoj 4002: [JLOI2015]有意义的字符串(特征根法+矩阵快速幂)
- bzoj 4002 [JLOI2015]有意义的字符串 数学
- 【结论&例题】【图(Floyd)&矩阵】图邻接矩阵上的乘法
- 矩阵的一些常用结论
- 【矩阵乘法】:矩阵乘法的基本实现
- 51nod 1197 字符串的数量 V2 dp+矩阵乘法
- 矩阵的思考-矩阵乘法
- 【codevs1287】矩阵乘法(矩阵乘法)
- 【codevs1287】矩阵乘法(矩阵乘法)
- bzoj 4180: 字符串计数 (后缀自动机+矩阵乘法)
- 写在人生特别的一天,烙在内心深处的感触
- 算法导论的C++实现(签到)
- Database—数据库基础
- Verilog序列检测器-两例
- Java 23种设计模式之原型模式
- [BZOJ4002][JLOI2015]有意义的字符串(结论+矩阵乘法)
- C#去除HTML标签
- linux elf加载过程
- 变分推断
- 网络判断
- [BZOJ2286][SDOI2011]消耗战(虚树+树形DP)
- 没有意大利炮的团长
- 文章标题
- 机器学习实战——岭回归、缩减法