Fibonacci-ish II

题目链接：luogu CF633H

题目大意

给你一个序列，每次问你一个区间，把里面的数拿出来去重排序，第 i 个位置乘上斐波那契数列第 i 项之后所有数的和。

思路

这题卡常。
（而且好像能暴力优化草过去但是写的是标算）

首先看着数据范围会主观思考 n\sqrt{n}n 有关的，思考完分块不太行之后。
我们考虑莫队（有人忘了这个东西我不说是谁），因为发现可以离线。

那你就考虑每次插入或者删除一个数，然后因为要去重，我们就只需要考虑新多出这个数和这个数被全部删完的时候。
那以放新的数为例，那就是这个数的贡献出现，然后给后面的数的位置后移一个，那就是斐波那契各自后一位。

那这个可以用什么维护呢，其实就是矩阵乘法（有人又忘了我不说是谁），那就是乘上一个转移矩阵嘛。
那至于去掉一个数，就是把它的贡献弄掉，然后把后面的数乘上一个转移矩阵的逆矩阵。
那这个区间乘矩阵，而且最后你要的是每个乘上的一个系数（就它自己的值），所以我们考虑用线段树来维护，除了矩阵那些再弄一个记录数组记录着系数，这样子我们就可以很快的让你当前位置的数的贡献出现（=ai=a_i=ai）或消失（000）

然后就是一些卡常说的：
因为你线段树嘛，所以要离散化，你就好离散化好了之后把编号对应一下就直接给权值取模了。
而且过程中也是能先不取模就别取模，而且别开 long long。
然后发现莫队那个奇偶优化真的有快到。（由于是最后加了这个过了心里的想法，说不定没加多少）
然后调调块长，加点 inline 和 register 应该就差不多了。

代码

#include<cmath>
#include<cstdio>
#include<algorithm>using namespace std;const int N = 3e4 + 100;
int n, m, a[N], b[N], bn, dy[N], q, num[N];
int B, bla[N], ans[N];
struct Quest {int id, l, r;
}qs[N];int re, zf; char c;
int read() {re = 0; zf = 1; c = getchar();while (c < '0' || c > '9') {if (c == '-') zf = -zf; c = getchar();}while (c >= '0' && c <= '9') {re = (re << 3) + (re << 1) + c - '0';c = getchar();}return re * zf;
}void writen(int x) {if (x > 9) writen(x / 10);putchar(x % 10 + '0');
}
void write(int x) {if (x < 0) putchar('-'), x = -x;writen(x);
}bool cmp(Quest x, Quest y) {if (bla[x.l] == bla[y.l]) return bla[x.l] & 1 ? x.r < y.r : x.r > y.r;return bla[x.l] < bla[y.l];
}struct matrix {int a[2][2];matrix() {a[0][0] = 1; a[0][1] = 0;a[1][0] = 0; a[1][1] = 1;}matrix(int op) {if (op == 1) {a[0][0] = 1; a[1][0] = 1;a[0][1] = 1; a[1][1] = 0;}if (op == -1) {a[0][0] = 0; a[1][0] = 1;a[0][1] = 1; a[1][1] = m - 1;}if (op == 0) {a[0][0] = 0; a[1][0] = 0;a[0][1] = 0; a[1][1] = 0;}}inline int* operator[](int x) {return a[x];}
}E, G, Gv;inline matrix operator +(matrix x, matrix y) {matrix z = matrix(0);for (int i = 0; i < 2; i++) for (int j = 0; j < 2; j++) z[i][j] = (x[i][j] + y[i][j]) % m;return z;
}
inline matrix operator *(matrix x, int y) {if (y == 1) return x;matrix z = matrix(0);for (int i = 0; i < 2; i++)for (int j = 0; j < 2; j++)z[i][j] = x[i][j] * y % m;return z;
}
inline matrix operator *(matrix x, matrix y) {matrix z = matrix(0);for (int k = 0; k < 2; k++)for (int i = 0; i < 2; i++)for (int j = 0; j < 2; j++)z[i][j] += x[i][k] * y[k][j];for (int i = 0; i < 2; i++) for (int j = 0; j < 2; j++) z[i][j] %= m;return z;
}
inline bool operator !=(matrix x, matrix y) {for (int i = 0; i < 2; i++) for (int j = 0; j < 2; j++)if (x[i][j] != y[i][j]) return 1;return 0;
}struct XD_tree {int fb[N << 2];matrix lzy[N << 2], val[N << 2];inline void up(int now) {val[now] = (val[now << 1] * fb[now << 1]) + (val[now << 1 | 1] * fb[now << 1 | 1]);}inline void downm(int now, matrix va) {val[now] = val[now] * va; lzy[now] = lzy[now] * va;}inline void down(int now) {if (lzy[now] != E) {downm(now << 1, lzy[now]); downm(now << 1 | 1, lzy[now]);lzy[now] = E;}}inline void build(int now, int l, int r) {if (l == r) {val[now] = G;return ;}fb[now] = 1;int mid = (l + r) >> 1;build(now << 1, l, mid); build(now << 1 | 1, mid + 1, r);up(now);}inline void insert(int now, int l, int r, int pl) {if (l == r) {fb[now] = dy[pl]; return ;}down(now); int mid = (l + r) >> 1;if (pl <= mid) insert(now << 1, l, mid, pl), downm(now << 1 | 1, G);else insert(now << 1 | 1, mid + 1, r, pl);up(now);}inline void erase(int now, int l, int r, int pl) {if (l == r) {fb[now] = 0; return ;}down(now); int mid = (l + r) >> 1;if (pl <= mid) erase(now << 1, l, mid, pl), downm(now << 1 | 1, Gv);else erase(now << 1 | 1, mid + 1, r, pl);up(now);}
}T;inline void add(int x) {if (!num[x]++) T.insert(1, 1, bn, x);
}inline void del(int x) {if (!--num[x]) T.erase(1, 1, bn, x);
}int main() {n = read(); m = read();for (int i = 1; i <= n; i++) {b[i] = a[i] = read();}E = matrix(); G = matrix(1); Gv = matrix(-1);sort(b + 1, b + n + 1); bn = unique(b + 1, b + n + 1) - b - 1;for (int i = 1; i <= bn; i++) dy[i] = b[i] % m;for (int i = 1; i <= n; i++) a[i] = lower_bound(b + 1, b + bn + 1, a[i]) - b;B = 300;
//  B = sqrt(n);for (int i = 1; i <= n; i++)bla[i] = (i - 1) / B + 1;q = read();for (int i = 1; i <= q; i++) {qs[i] = (Quest){i, read(), read()};}sort(qs + 1, qs + q + 1, cmp);T.build(1, 1, bn);int l = 1, r = 0;for (int i = 1; i <= q; i++) {while (l < qs[i].l) del(a[l]), l++;while (l > qs[i].l) l--, add(a[l]);while (r < qs[i].r) r++, add(a[r]);while (r > qs[i].r) del(a[r]), r--;ans[qs[i].id] = T.val[1][0][1];}for (int i = 1; i <= q; i++) {write((ans[i] % m + m) % m); putchar('\n');}return 0;
}

【luogu CF633H】Fibonacci-ish II（莫队）（线段树）（矩阵乘法）相关推荐

牛客练习赛85 数学家的迷题（带修莫队/线段树）
题意: 1:将a[id]a[id]a[id]的值改为xxx. 2:令t=a[l]×a[l+1]×...×a[r−1]×a[r]t=a[l]×a[l+1]×...×a[r−1]×a[r]t=a[l]×a ...
HDU - 5381 The sum of gcd(莫队/线段树区间合并)
题目链接:点击查看题目大意:给出一个长度为 nnn 的序列,再给出 mmm 次询问,每次询问需要回答区间 [L,R][L,R][L,R] 内所有子区间的 gcdgcdgcd 之和.更具体的,对于询问 ...
SPOJ - DQUERY D-query(莫队/线段树+离线/主席树)
题目链接:点击查看题目大意:给出一个由n个数组成的序列,再给出m次查询,每次查询区间[l,r]中有多少个不同的数题目分析:莫队模板题,直接套板子就好了有点意思的是函数返回值为布尔类型,然后没有r ...
51nod-猴猴的比赛【莫队,线段树】
正题题目链接:https://www.51nod.com/Contest/Problem.html#contestProblemId=1150 题目大意给出两颗nnn个点的树,求有多少个点(i,j ...
jzoj6293-迷宫【ddp,线段树,矩阵乘法】
正题题目大意一个n∗mn*mn∗m的迷宫,不能往左走,有墙,每次修改一个点或询问两个点之间的最短距离. 解题思路考虑到nnn的值很小,所以我们可以用矩阵转移,然后要求支持修改和查询所以我们考虑d ...
LOJ2980 THUSC2017大魔法师（线段树+矩阵乘法）
线段树每个节点维护(A,B,C,len)向量,操作即是将其乘上一个矩阵. #include<iostream> #include<cstdio> #include<cma ...
牛客 - 求函数(线段树+区间合并/线段树+矩阵维护)
题目链接:点击查看题目大意:现在有 n 个函数,每个函数都是诸如 f( x ) = k * x + b 的形式,只是每个函数的 k 和 b 都是相互独立的,现在给出两个操作: 1 pos k b:将 ...
P4113 [HEOI2012]采花【树状数组 AC】【莫队/主席树 TLE】
传送门题意简述:给定一个长度为 n 的序列,有 m 次询问,每次询问一段区间,求区间中有多少个数出现次数超过 1 次这道题真的有点恶心...在洛谷上属于莫队,但是莫队会T,数据居然有2e6.莫队的 ...
Hdu 6534 Chika and Friendly Pairs 莫队算法+树状数组
题目链接题意求给区间[L,R]中有少对(i,j)满足i<j且abs(a[i]-a[j])<=k. 首先来说暴力的方法就是离散化,然后用树状数组来维护,但是m次询问,m很大,所以说一定会t ...
luogu P5142 区间方差（线段树、乘法逆元）
luogu P5142 区间方差本题要求维护模区间方差,很明显是一道数据结构题. 我们化简方差公式: 而平均数等于可以发现,我们只需要维护序列的区间和和区间平方和,就可以维护平均数和方差. 区间和 ...

【luogu CF633H】Fibonacci-ish II（莫队）（线段树）（矩阵乘法）

Fibonacci-ish II

题目链接：luogu CF633H

题目大意

思路

代码

【luogu CF633H】Fibonacci-ish II（莫队）（线段树）（矩阵乘法）相关推荐

最新文章

热门文章