亚洲香蕉成人av网站在线观看_欧美精品成人91久久久久久久_久久久久久久久久久亚洲_热久久视久久精品18亚洲精品_国产精自产拍久久久久久_亚洲色图国产精品_91精品国产网站_中文字幕欧美日韩精品_国产精品久久久久久亚洲调教_国产精品久久一区_性夜试看影院91社区_97在线观看视频国产_68精品久久久久久欧美_欧美精品在线观看_国产精品一区二区久久精品_欧美老女人bb

首頁 > 學院 > 開發設計 > 正文

莫隊算法——解決序列上詢問的利器

2019-11-11 04:12:05
字體:
來源:轉載
供稿:網友

問題: 有一個長為N序列,有M個詢問:在區間[L,R]內,出現了多少個不同的數字。(序列中所有數字均小于K)。題目會給出K。

莫隊算法就是滋磁解決這類問題的離線算法。(其實很簡單

首先來看看暴力: 由于暴力還是比較水的,所以直接上:

#include <bits/stdc++.h>using namespace std ;const int maxn = 50010 ;int n, m, a[maxn] ;bool vis[maxn] ;int main() { int i, j, k, query_time, L, R ; cin >> n >> query_time >> k ; for ( i = 1 ; i <= n ; i ++ ) cin >> a[i] ; while ( query_time -- ) { cin >> L >> R ; memset ( vis, 0, sizeof(vis) ) ; for ( i = L ; i <= R ; i ++ ) vis[a[i]] = true ; int ans = 0 ; for ( i = 0 ; i <= k ; i ++ ) ans += vis[i] ; cout << ans << endl ; } return 0 ;}

這個復雜度顯然是 O(N2) 的。有些題目的范圍可能到100000或者更大,那么顯然就不能了。

這里還有一種稍作改進的方法,比上述做法大多數情況要快些。

void add ( int pos ) { ++cnt[a[pos]] ; if ( cnt[a[pos]] == 1 ) ++ answer ;}void remove ( int pos ) { -- cnt[a[pos]] ; if ( cnt[a[pos]] == 0 ) -- answer ;}void solve() { int curL = 1, curR = 0 ; // current L R for ( each query [L,R] ) { while ( curL < L ) remove ( curL++ ) ; while ( curL > L ) add ( --curL ) ; while ( curR < R ) add ( ++curR ) ; while ( curR > R ) remove ( curR-- ) ; cout << answer << endl ; // Warning : please notice the order "--","++" and "cur" ; }}

其實還算好理解的,如果不明白,隨便搞組數據手玩一下就明白了

手玩數據:6 4 31 3 2 1 1 31 42 63 55 6輸出答案:322

不幸的是,雖然這個東西比我們的暴力要快,但是它的時間復雜度仍然是 O(N2) 的。 但好消息是,這個東西可以算是莫隊算法的核心了。(你在逗我笑????) 其實是真的 :)

莫隊算法是怎么做的呢? 考慮到上述改進的辦法的效率低下主要是因為curL和curR兩個指針前前后后跑來跑去太多次。于是,我們的做法就是不要讓他們跑太多沒用的距離。 我們可以通過離線下所有的詢問,然后通過某種排序,讓兩個指針跑動的距離盡量變少。具體的做法是把N劃分成N??√段,每段長度都是N??√,然后在把所有詢問按照L端點排序,看各個詢問被劃分到哪一塊里。接著,對于各個劃分出的段,在各自的段里,將它包含的所有區間再按照R端點排序。 舉個例子:假設我們有3個長度為3的段(0-2,3-5,6-8): {0, 3} {1, 7} {2, 8} {7, 8} {4, 8} {4, 4} {1, 2} 先根據所在段落的編號重排 {0, 3} {1, 7} {2, 8} {1, 2} | {4, 8} {4, 4} | {7, 8} 現在按R的值重排 {1, 2} {0, 3} {1, 7} {2, 8} | {4, 4} {4, 8} | {7, 8}

然后?還要然后嗎?就用剛剛那個算法來玩就好了。畢竟我們只是交換了一下詢問的順序而已,并沒有對算法做什么改動。

對于這個的復雜度嘛,其實是比較鬼畜的,就這么改了下順序,然后就變成了O(N?N??√)

上面代碼所有查詢的復雜性是由4個while循環決定的。前2個while循環是curL的移動總量”,后2個while循環是curR的移動總量”。這兩者的和將是總復雜度。 有趣的是。先考慮右指針:對于每個塊,查詢是遞增的順序排序,所以右指針curR按照遞增的順序移動。在下一個塊的開始時,指針是盡量靠右的 ,將移動到下一個塊中的最小的R處。這意味著對于一個給定的塊,右指針移動的量是 O(N)。我們有O(N??√)個塊。所以總共是O(N?N??√)。 關于左指針:所有查詢的左指針都在同一段中,當我們完成一個查詢到下一個查詢時,左指針會移動,但由于兩次詢問的L在同一塊中,此移動是 O(N??√)的。所以,左指針的移動總量是O(M?N??√)。 所以,總復雜度就是O((N+M)?N??√),就當做是O(N?N??√)吧。 是不是很簡單的吶2333

接下來給幾個例題:

例題1:BZOJ3781 小B的詢問

Description

小B有一個序列,包含N個1~K之間的整數。他一共有M個詢問,每個詢問給定一個區間[L..R],求Sigma(c(i)^2)的值,其中i的值從1到K,其中c(i)表示數字i在[L..R]中的重復次數。小B請你幫助他回答詢問。

Input

第一行,三個整數N、M、K。 第二行,N個整數,表示小B的序列。 接下來的M行,每行兩個整數L、R。

Output

M行,每行一個整數,其中第i行的整數表示第i個詢問的答案。

Sample Input

6 4 3 1 3 2 1 1 3 1 4 2 6 3 5 5 6

Sample Output

6 9 5 2

HINT

對于全部的數據,1<=N、M、K<=50000

很裸的吧~

/************************************************************** Source Code : GoAway Date : 2017-02-06****************************************************************/#include <iostream>#include <cstdio>#include <cstring>#include <cstdlib>#include <vector>#include <map>#include <stack>#include <queue>#include <set>#include <cmath>#include <algorithm>#include <ctime>using namespace std ;const int zhf = 1<<30 ;const int maxn = 50010, tim = 250 ;bool Read ( int &x ) { bool f = 0 ; x = 0 ; char c = getchar() ; while ( !isdigit(c) ) { if ( c == '-' ) f = 1 ; if ( c == EOF ) return false ; c = getchar() ; } while ( isdigit(c) ) { x = 10 * x + c - '0' ; c = getchar() ; } if ( f ) x = -x ; return true ;}struct query { int L, R, id ; friend bool Operator < ( query a, query b ) { return (a.L/tim) == (b.L/tim) ? a.R < b.R : a.L < b.L ; }} e[maxn] ;int n, m, a[maxn], cnt[maxn], ans[maxn], answer ;void add ( int pos ) { answer += (cnt[a[pos]]++)<<1|1 ;}void remove ( int pos ) { answer -= (--cnt[a[pos]])<<1|1 ;}int main() { int i, j, k, curL = 1, curR = 0 ; Read(n) ; Read(m) ; Read(j) ; for ( i = 1 ; i <= n ; i ++ ) Read(a[i]) ; for ( i = 1 ; i <= m ; i ++ ) { Read(e[i].L) ; Read(e[i].R) ; e[i].id = i ; } sort ( e+1, e+m+1 ) ; for ( i = 1 ; i <= m ; i ++ ) { int L = e[i].L, R = e[i].R ; while ( curL < L ) remove ( curL++ ) ; while ( curL > L ) add ( --curL ) ; while ( curR < R ) add ( ++curR ) ; while ( curR > R ) remove ( curR-- ) ; ans[e[i].id] = answer ; } for ( i = 1 ; i <= m ; i ++ ) PRintf ( "%d/n", ans[i] ) ; return 0 ;}

例題2:SDOI2009 HH的項鏈 洛谷1972

Description

HH有一串由各種漂亮的貝殼組成的項鏈。HH相信不同的貝殼會帶來好運,所以每次散步 完后,他都會隨意取出一段貝殼,思考它們所表達的含義。HH不斷地收集新的貝殼,因此, 他的項鏈變得越來越長。有一天,他突然提出了一個問題:某一段貝殼中,包含了多少種不同 的貝殼?這個問題很難回答。。。因為項鏈實在是太長了。于是,他只好求助睿智的你,來解 決這個問題。

Input

第一行:一個整數N,表示項鏈的長度。 第二行:N個整數,表示依次表示項鏈中貝殼的編號(編號為0到1000000之間的整數)。 第三行:一個整數M,表示HH詢問的個數。 接下來M行:每行兩個整數,L和R(1 ≤ L ≤ R ≤ N),表示詢問的區間。

Output

M行,每行一個整數,依次表示詢問對應的答案。

Sample Input

6 1 2 3 4 3 5 3 1 2 3 5 2 6

Sample Output

2 2 4

HINT

對于20%的數據,N ≤ 100,M ≤ 1000; 對于40%的數據,N ≤ 3000,M ≤ 200000; 對于100%的數據,N ≤ 50000,M ≤ 200000。

還是很裸的2333

/************************************************************** Source Code : GoAway Date : 2017-02-06****************************************************************/#include <iostream>#include <cstdio>#include <cstring>#include <cstdlib>#include <vector>#include <map>#include <stack>#include <queue>#include <set>#include <cmath>#include <algorithm>#include <ctime>using namespace std ;const int zhf = 1<<30 ;const int maxn = 1000010 ;bool Read ( int &x ) { bool f = 0 ; x = 0 ; char c = getchar() ; while ( !isdigit(c) ) { if ( c == '-' ) f = 1 ; if ( c == EOF ) return false ; c = getchar() ; } while ( isdigit(c) ) { x = 10 * x + c - '0' ; c = getchar() ; } if ( f ) x = -x ; return true ;}int n, m, cnt[maxn], a[maxn], answer, ans[maxn], tim ;struct query { int l, r, id ; friend bool operator < ( query a, query b ) { return (a.l/tim) == (b.l/tim) ? a.r < b.r : a.l<b.l ; }} e[maxn] ;void add ( int pos ) { if ( (++cnt[a[pos]]) == 1 ) ++ answer ;}void remove ( int pos ) { if ( (--cnt[a[pos]]) == 0 ) -- answer ;}int main() { int i, j, k, curL = 1, curR = 0 ; Read(n) ; for ( i = 1 ; i <= n ; i ++ ) Read(a[i]) ; Read(m) ; tim = sqrt(m) ; for ( i = 1 ; i <= m ; i ++ ) { Read(e[i].l) ; Read(e[i].r) ; e[i].id = i ; } sort(e+1,e+m+1) ; for ( i = 1 ; i <= m ; i ++ ) { int L = e[i].l, R = e[i].r ; while ( curL < L ) remove(curL++) ; while ( curL > L ) add(--curL) ; while ( curR < R ) add(++curR) ; while ( curR > R ) remove(curR--) ; ans[e[i].id] = answer ; } for ( i = 1 ; i <= m ; i ++ ) printf ( "%d/n", ans[i] ) ; return 0 ;}

例題3:2009國家集訓隊 小Z的襪子 清橙OJ1206

Description

HH有一串由各種漂亮的貝殼組成的項鏈。HH相信不同的貝殼會帶來好運,所以每次散步 完后,他都會隨意取出一段貝殼,思考它們所表達的含義。HH不斷地收集新的貝殼,因此, 他的項鏈變得越來越長。有一天,他突然提出了一個問題:某一段貝殼中,包含了多少種不同 的貝殼?這個問題很難回答。。。因為項鏈實在是太長了。于是,他只好求助睿智的你,來解 決這個問題。

Input

第一行:一個整數N,表示項鏈的長度。 第二行:N個整數,表示依次表示項鏈中貝殼的編號(編號為0到1000000之間的整數)。 第三行:一個整數M,表示HH詢問的個數。 接下來M行:每行兩個整數,L和R(1 ≤ L ≤ R ≤ N),表示詢問的區間。

Output

M行,每行一個整數,依次表示詢問對應的答案。

Sample Input

6 1 2 3 4 3 5 3 1 2 3 5 2 6

Sample Output

2 2 4

HINT

對于20%的數據,N ≤ 100,M ≤ 1000; 對于40%的數據,N ≤ 3000,M ≤ 200000; 對于100%的數據,N ≤ 50000,M ≤ 200000。

樣例解釋

詢問1:共C(5,2)=10種可能,其中抽出兩個2有1種可能,抽出兩個3有3種可能,概率為(1+3)/10=4/10=2/5。 詢問2:共C(3,2)=3種可能,無法抽到顏色相同的襪子,概率為0/3=0/1。 詢問3:共C(3,2)=3種可能,均為抽出兩個3,概率為3/3=1/1。 注:上述C(a, b)表示組合數,組合數C(a, b)等價于在a個不同的物品中選取b個的選取方案數。

這道題倒是有些需要想想。 但是有一個神奇的事情:C(N,M)=N!M!?(N?M)! 那么C(N,2)呢? 這不就是一個∑N?1i=1i嘛? (有沒有感覺自己被續了一秒

/************************************************************** Source Code : GoAway Date : 2017-02-06****************************************************************/#include <iostream>#include <cstdio>#include <cstring>#include <cstdlib>#include <vector>#include <map>#include <stack>#include <queue>#include <set>#include <cmath>#include <algorithm>#include <ctime>#define ll long longusing namespace std ;const ll zhf = 1<<30 ;const ll maxn = 50010 ;bool Read ( ll &x ) { bool f = 0 ; x = 0 ; char c = getchar() ; while ( !isdigit(c) ) { if ( c == '-' ) f = 1 ; if ( c == EOF ) return false ; c = getchar() ; } while ( isdigit(c) ) { x = 10 * x + c - '0' ; c = getchar() ; } if ( f ) x = -x ; return true ;}ll tim ;struct query { ll L, R, id ; friend bool operator < ( query a, query b ) { return (a.L/tim) == (b.L/tim) ? a.R < b.R : a.L < b.L ; }} e[maxn] ;ll gcd ( ll x, ll y ) { return y ? gcd ( y, x%y ) : x ;}struct Answer { ll x, y ; void out() { if ( !x ) puts("0/1") ; else { ll d = gcd(x, y) ; x /= d ; y /= d ; printf ( "%lld/%lld/n", x, y ) ; } }} ans[maxn] ;ll n, m, a[maxn], cnt[maxn], answer ;void add ( ll pos ) { ++ cnt[a[pos]] ; if ( cnt[a[pos]] > 1 ) answer += cnt[a[pos]] - 1 ;}void remove ( ll pos ) { -- cnt[a[pos]] ; if ( cnt[a[pos]] > 0 ) answer -= cnt[a[pos]] ;}ll sum ( ll x ) { return x*(x+1)/2 ; }int main() { ll i, j, k, curL = 1, curR = 0 ; Read(n) ; Read(m) ; tim = sqrt(m) ; for ( i = 1 ; i <= n ; i ++ ) Read(a[i]) ; for ( i = 1 ; i <= m ; i ++ ) { Read(e[i].L) ; Read(e[i].R) ; e[i].id = i ; } sort ( e+1, e+m+1 ) ; for ( i = 1 ; i <= m ; i ++ ) { ll L = e[i].L, R = e[i].R ; while ( curL < L ) remove ( curL++ ) ; while ( curL > L ) add ( --curL ) ; while ( curR < R ) add ( ++curR ) ; while ( curR > R ) remove ( curR-- ) ; ans[e[i].id] = (Answer){answer,sum(R-L)} ; } for ( i = 1 ; i <= m ; i ++ ) ans[i].out() ; return 0 ;}

其實莫隊還可以套上樹狀數組或者在一棵樹上搞的,但是這篇文章就簡單先介紹下啦~ 要是俺有空就出莫隊的續集,嘿嘿~


發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
亚洲香蕉成人av网站在线观看_欧美精品成人91久久久久久久_久久久久久久久久久亚洲_热久久视久久精品18亚洲精品_国产精自产拍久久久久久_亚洲色图国产精品_91精品国产网站_中文字幕欧美日韩精品_国产精品久久久久久亚洲调教_国产精品久久一区_性夜试看影院91社区_97在线观看视频国产_68精品久久久久久欧美_欧美精品在线观看_国产精品一区二区久久精品_欧美老女人bb
国产精品久久中文| 欧美激情久久久| 国产精品美女在线观看| 亚洲国产成人精品一区二区| 久久精品国产96久久久香蕉| 久久久久久久久久久国产| 69国产精品成人在线播放| 日韩av手机在线观看| 91沈先生作品| 久久理论片午夜琪琪电影网| 亚洲a∨日韩av高清在线观看| 88xx成人精品| 国产成人精品视频在线观看| 国产日韩亚洲欧美| 亚洲无限乱码一二三四麻| 中文字幕亚洲一区| 91高清视频免费观看| 亚洲国产天堂久久综合| 欧美男插女视频| 国产不卡av在线免费观看| 国产成人精品综合久久久| 久久人体大胆视频| 亚洲欧美成人一区二区在线电影| 亚洲精品自在久久| 欧美性xxxxx| 精品一区二区三区四区| 91免费的视频在线播放| 欧美一乱一性一交一视频| 成人久久久久久| 国产色婷婷国产综合在线理论片a| 97免费视频在线| 精品二区三区线观看| 欧美性xxxx极品hd欧美风情| 欧美电影《睫毛膏》| 欧美肥老太性生活视频| 久久福利视频导航| 人九九综合九九宗合| 亚洲人午夜精品| 尤物yw午夜国产精品视频| 欧美成年人在线观看| 亚洲国产精品成人va在线观看| 国产剧情久久久久久| 国产成人极品视频| 欧美在线播放视频| 欧美性视频网站| 欧美一性一乱一交一视频| 国产精品美女免费看| 国外色69视频在线观看| 国产精品草莓在线免费观看| 久久久国产精品免费| 欧美国产日本在线| 精品国偷自产在线视频| 在线观看视频亚洲| 亚洲精品美女免费| 精品夜色国产国偷在线| 国产精品午夜国产小视频| 亚洲视频国产视频| 亚洲人成自拍网站| 国产精品久久久久久久久久新婚| 日韩久久午夜影院| www.日韩不卡电影av| 九九久久久久久久久激情| 亚洲另类xxxx| 久久久久久成人精品| 午夜精品久久久久久久99热| 亚洲精品720p| 91香蕉亚洲精品| 国产精品视频自拍| 亚洲精品久久久久中文字幕欢迎你| 久久久久成人网| www.亚洲成人| 亚洲国产成人精品久久久国产成人一区| 91精品国产综合久久久久久蜜臀| 欧美精品在线观看91| 国模极品一区二区三区| 一区二区三区四区精品| 色婷婷久久av| 5252色成人免费视频| 亚洲中国色老太| 日韩精品视频中文在线观看| 国产不卡在线观看| 精品少妇v888av| 国产美女久久精品香蕉69| 亚洲美女久久久| 日韩在线播放一区| 久久久久久久久久婷婷| 国产午夜精品免费一区二区三区| 亚洲欧美日韩视频一区| 国产精品中文字幕在线观看| 成人在线国产精品| 日韩av电影在线播放| 日韩精品在线观看视频| 尤物yw午夜国产精品视频| 日韩精品在线观看一区| 91av在线视频观看| 久久手机免费视频| 中文字幕欧美精品日韩中文字幕| 欧美性受xxxx白人性爽| 国产精品久久久久久久久免费| 91成人在线视频| 日韩美女中文字幕| 日韩中文av在线| 国产香蕉一区二区三区在线视频| 国产在线观看精品| 亚洲天堂第二页| 欧美电影《睫毛膏》| 国产一区二区视频在线观看| 精品无人区乱码1区2区3区在线| 亚洲精品av在线| 萌白酱国产一区二区| 日韩av免费在线| 国产在线拍揄自揄视频不卡99| 国产日本欧美一区二区三区在线| 中国日韩欧美久久久久久久久| 成人黄色短视频在线观看| 亚洲性猛交xxxxwww| 日韩在线视频二区| 日韩大片在线观看视频| 亚洲图片在区色| 国产在线观看一区二区三区| 亚洲欧美国产制服动漫| 国产日韩欧美成人| 国产99久久久欧美黑人| 久久成人这里只有精品| 亚洲第一区中文字幕| 国产精品久久久久77777| 成人网欧美在线视频| 亚洲人成在线播放| 欧美日韩在线视频观看| 亚洲黄页视频免费观看| 欧美怡春院一区二区三区| 一本大道亚洲视频| 国产亚洲aⅴaaaaaa毛片| 国产精品第8页| 美女av一区二区三区| 狠狠色香婷婷久久亚洲精品| 国产欧美久久久久久| 国产精品视频白浆免费视频| 高清亚洲成在人网站天堂| 日韩综合视频在线观看| 国产精品美女久久| 91精品久久久久久久久久入口| 日韩中文字幕免费| 久久99热精品这里久久精品| 亚洲无线码在线一区观看| 亚洲精品美女在线观看播放| 97成人精品视频在线观看| 亚洲精品国产拍免费91在线| 福利一区视频在线观看| 性欧美视频videos6一9| 欧美黄网免费在线观看| 亚洲成人动漫在线播放| 国产精品白丝av嫩草影院| 欧美日韩精品在线视频| 秋霞午夜一区二区| 性欧美亚洲xxxx乳在线观看| 日韩美女激情视频| 3344国产精品免费看| 综合av色偷偷网| 亚洲男人天堂2024| 久久精品成人欧美大片古装| 久久久久久久久国产| 57pao成人国产永久免费| 国产精品视频26uuu|