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

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

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

2019-11-11 05:40:41
字體:
來源:轉載
供稿:網友

問題: 有一個長為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 ;}

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


上一篇:POJ 3059 Wormholes

下一篇:template類

發表評論 共有條評論
用戶名: 密碼:
驗證碼: 匿名發表
亚洲香蕉成人av网站在线观看_欧美精品成人91久久久久久久_久久久久久久久久久亚洲_热久久视久久精品18亚洲精品_国产精自产拍久久久久久_亚洲色图国产精品_91精品国产网站_中文字幕欧美日韩精品_国产精品久久久久久亚洲调教_国产精品久久一区_性夜试看影院91社区_97在线观看视频国产_68精品久久久久久欧美_欧美精品在线观看_国产精品一区二区久久精品_欧美老女人bb
国产综合香蕉五月婷在线| 国产精品99蜜臀久久不卡二区| 亚洲一区二区三区香蕉| 亚洲成人激情在线观看| 精品日韩中文字幕| 亚洲精品国产精品久久清纯直播| 狠狠躁18三区二区一区| 欧美www视频在线观看| 91人成网站www| 国产精品久久色| 热门国产精品亚洲第一区在线| 中文字幕久热精品在线视频| 久久国产精品电影| 国产一区二区三区直播精品电影| 久久精品国产2020观看福利| 久久久久国产精品免费网站| 欧美猛交免费看| 午夜剧场成人观在线视频免费观看| 国产香蕉一区二区三区在线视频| 日韩欧美一区二区三区| 91国内免费在线视频| 成人精品视频99在线观看免费| 91影视免费在线观看| 日韩欧美中文字幕在线播放| 亚洲成年人影院在线| 欧美高清自拍一区| 亚洲综合在线小说| 国产成人精品国内自产拍免费看| 国产日韩视频在线观看| 91探花福利精品国产自产在线| 日产精品久久久一区二区福利| 日韩亚洲综合在线| 91免费综合在线| 亚洲国产精品嫩草影院久久| 在线观看欧美成人| 国产精品久久久久久久久久久不卡| 美日韩精品免费视频| 国内精品视频久久| 精品久久久91| 日韩欧美高清视频| 国产亚洲欧美一区| 91精品国产综合久久香蕉的用户体验| 国内精品一区二区三区四区| 麻豆乱码国产一区二区三区| 久久精品视频亚洲| 国产精品视频最多的网站| 色噜噜亚洲精品中文字幕| 亚洲第一区第一页| 国产日韩欧美在线观看| 亚洲国模精品一区| 久久亚洲精品毛片| 国产精品专区h在线观看| 夜夜嗨av一区二区三区四区| 91豆花精品一区| 7777精品视频| 精品国产999| 国产主播喷水一区二区| 久久精品夜夜夜夜夜久久| 欧美巨猛xxxx猛交黑人97人| 亚洲а∨天堂久久精品喷水| www.久久色.com| 亚洲xxxx视频| 久久久精品一区二区| 日韩激情av在线免费观看| 亚洲精品在线观看www| 久久精品福利视频| 日韩欧美国产骚| 亚洲成人av片在线观看| 亚洲国产成人久久综合一区| 国产精品1234| 欧美日韩国产综合视频在线观看中文| 九九热精品视频在线播放| 亚洲一区二区中文| 国产91色在线|免| 中文字幕亚洲欧美日韩在线不卡| 国产91|九色| 欧美午夜xxx| 久久韩剧网电视剧| 日韩精品免费看| 在线日韩欧美视频| 国产精品一香蕉国产线看观看| 国产美女久久精品香蕉69| 亚洲性69xxxbbb| 97不卡在线视频| 欧美电影电视剧在线观看| 91精品视频免费| 欧美黄色片视频| 欧美一级淫片aaaaaaa视频| 亚洲天堂av在线免费观看| 亚洲欧美制服第一页| 亚洲欧美国产精品| 日韩欧美在线视频| 国内精品久久久久久久| 成人亚洲欧美一区二区三区| 中国人与牲禽动交精品| 国产一区二区三区在线播放免费观看| 亚洲电影免费观看高清完整版在线观看| 精品香蕉在线观看视频一| 黄色一区二区在线| 亚洲美女视频网| 日韩欧美在线第一页| 亚洲精品美女久久| 亚洲国产精品视频在线观看| 欧美电影免费在线观看| 91精品国产综合久久久久久久久| 国产精品日韩一区| 国产精品丝袜视频| 一本色道久久88亚洲综合88| 这里只有精品在线观看| 91免费的视频在线播放| 国产成人亚洲综合| 91亚洲精品视频| 狠狠操狠狠色综合网| 国产精品人人做人人爽| 亚洲一区二区日本| 亚洲第一精品久久忘忧草社区| 国产区精品视频| 亚洲日本欧美日韩高观看| 日韩中文在线不卡| 日韩av大片在线| 性欧美在线看片a免费观看| 亚洲老头同性xxxxx| 久久99亚洲热视| 亚洲综合色av| 亚洲女人初尝黑人巨大| 亚洲aⅴ男人的天堂在线观看| 欧美精品成人在线| 国语自产精品视频在线看一大j8| 精品福利免费观看| 欧美精品午夜视频| 668精品在线视频| 欧美国产视频一区二区| 欧美日韩在线一区| 精品福利樱桃av导航| 亚洲欧美成人精品| 亚洲女人天堂色在线7777| 亚洲精品国产精品久久清纯直播| 欧美电影在线观看| 亚洲男人天堂网| 国产精品日本精品| 国产激情999| 2019中文字幕在线| 成人精品在线观看| 日韩欧美在线第一页| 91中文字幕在线| 精品一区二区三区四区在线| 欧美噜噜久久久xxx| 中文字幕免费精品一区| 国产精品偷伦视频免费观看国产| 日韩精品高清在线观看| 亚洲第一区中文99精品| 国产午夜精品全部视频播放| 久久久视频免费观看| 国产在线精品播放| 亚洲欧美日韩高清| 亚洲女人被黑人巨大进入al| 成人精品久久av网站| 久久精品国产久精国产一老狼| 欧美性猛交xxxx黑人猛交| xxx一区二区| 国产精品视频网站| 久久影视电视剧凤归四时歌| 久久久99久久精品女同性| 国产精品678|