我也不知道该咋分类--莫队算法
生活随笔
收集整理的這篇文章主要介紹了
我也不知道该咋分类--莫队算法
小編覺得挺不錯的,現在分享給大家,幫大家做個參考.
莫隊算法可以一個可高效解決絕大多數離線+無修改+區間查詢問題的算法。這類問題具體是指:如果知道[L,R]的答案時,可以在O(g(n))的時間下得到[L,R?1],[L,R+1],[L?1,R],[L+1,R]的答案的話,就可以O(nn?g(n))O(n\sqrt n · \mathrm{g}(n))O(nn??g(n))的時間內求出所有查詢。
假設我們算完[L,R]的答案后現在要算[L′,R′]的答案。由于可以在O(1)的時間下得到[L,R?1],[L,R+1],[L?1,R],[L+1,R]的答案,所以計算[L′,R′]的答案耗時|L?L′|+|R?R′|。如果把詢問[L,R]看做平面上的點a(L,R),詢問[L′,R′]看做點b(L′,R′)的話,那么時間開銷就為兩點的曼哈頓距離。
具體實現
我們先對序列分塊,每塊長度為sqrt(n),然后以詢問左端點所在的塊為第一關鍵字,右端點的大小為第二關鍵字進行排序
然后每個區間由它的上一個區間推出來,雖然單次詢問仍可能是O(n),
但是在塊中每次詢問l的移動最大為sqrt(n),r的移動總和最大為n,塊與塊之間的移動最大為n
所以總復雜度為O((n+m)sqrt(n))
我這里給出莫隊的模板,
#include <iostream> #include <cstdio> #include <cstring> #include <algorithm> #include <cmath> #define N 50010 #define LL long long using namespace std; int n, m, k, unit; struct Mo {int l;int r;int id;int pos;bool operator<(const Mo &rhs) const{if (pos == rhs.pos)return r < rhs.r;elsereturn l < rhs.l;} } a[N]; int b[N], cnt[N]; LL Ans[N]; int main() {scanf("%d%d%d", &n, &m, &k);unit = sqrt(n + 1);for (int i = 1; i <= n; i++)scanf("%d", &b[i]);for (int i = 1; i <= m; i++){scanf("%d%d", &a[i].l, &a[i].r);a[i].id = i;a[i].pos = a[i].l / unit;}sort(a + 1, a + 1 + m);int l = 1, r = 0;LL ans = 0;for (int i = 1; i <= m; i++){while (l > a[i].l){l--;//操作ans += }while (r < a[i].r){r++;//操作ans += }while (l < a[i].l){ans -= //操作l++;}while (r > a[i].r){ans -= //操作r--;}Ans[a[i].id] = ans;}for (int i = 1; i <= m; i++)printf("%lld\n", Ans[i]);return 0; }總結
以上是生活随笔為你收集整理的我也不知道该咋分类--莫队算法的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: 数学--数论--HDU 4675 GCD
- 下一篇: 感染infection攻略有哪些