青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

M.J的blog

algorithm,ACM-ICPC
隨筆 - 39, 文章 - 11, 評論 - 20, 引用 - 0
數(shù)據(jù)加載中……

POJ 2762. Balanced Lineup (區(qū)間求最值sparsetable算法)

啥也不說了,sparsetable算法的高效性,我真佩服那個(gè)發(fā)明它的人。Orz...

以下是關(guān)于這個(gè)算法的講解(我也沒理解透,回頭細(xì)細(xì)品味吧)

RMQ(Range Minimum/Maximum Query)問題是求區(qū)間最值問題。你當(dāng)然可以寫個(gè)O(n)的(怎么寫都可以吧=_=),但是萬一要詢問最值1000000遍,估計(jì)你就要掛了。這時(shí)候你可以放心地寫一個(gè)線段樹(前提是不寫錯(cuò))O(logn)的復(fù)雜度應(yīng)該不會掛。但是,這里有更牛的算法,就是ST算法,它可以做到O(nlogn)的預(yù)處理,O(1)!!!地回答每個(gè)詢問。

       來看一下ST算法是怎么實(shí)現(xiàn)的(以最大值為例):

       首先是預(yù)處理,用一個(gè)DP解決。設(shè)a[i]是要求區(qū)間最值的數(shù)列,f[i, j]表示從第i個(gè)數(shù)起連續(xù)2^j個(gè)數(shù)中的最大值。例如數(shù)列3 2 4 5 6 8 1 2 9 7 ,f[1,0]表示第1個(gè)數(shù)起,長度為2^0=1的最大值,其實(shí)就是3這個(gè)數(shù)。f[1,2]=5,f[1,3]=8,f[2,0]=2,f[2,1]=4……從這里可以看出f[i,0]其實(shí)就等于a[i]。這樣,Dp的狀態(tài)、初值都已經(jīng)有了,剩下的就是狀態(tài)轉(zhuǎn)移方程。我們把f[i,j]平均分成兩段(因?yàn)閒[i,j]一定是偶數(shù)個(gè)數(shù)字),從i到i+2^(j-1)-1為一段,i+2^(j-1)到i+2^j-1為一段(長度都為2^(j-1))。用上例說明,當(dāng)i=1,j=3時(shí)就是3,2,4,5 和 6,8,1,2這兩段。f[i,j]就是這兩段的最大值中的最大值。于是我們得到了動規(guī)方程F[i, j]=max(F[i,j-1], F[i + 2^(j-1),j-1]).

 

接下來是得出最值,也許你想不到計(jì)算出f[i,j]有什么用處,一般毛想想計(jì)算max還是要O(logn),甚至O(n)。但有一個(gè)很好的辦法,做到了O(1)。還是分開來。如在上例中我們要求區(qū)間[2,8]的最大值,就要把它分成[2,5]和[5,8]兩個(gè)區(qū)間,因?yàn)檫@兩個(gè)區(qū)間的最大值我們可以直接由f[2,2]和f[5,2]得到。擴(kuò)展到一般情況,就是把區(qū)間[L,R]分成兩個(gè)長度為2^n的區(qū)間(保證有f[i,j]對應(yīng))。直接給出表達(dá)式:

k := ln(R-L+1) / ln(2);

ans := max(F[L,k], F[R - 2^k+1, k]);

這樣就計(jì)算了從i開始,長度為2^t次的區(qū)間和從r-2^i+1開始長度為2^t的區(qū)間的最大值(表達(dá)式比較煩瑣,細(xì)節(jié)問題如加1減1需要仔細(xì)考慮.

Code:

 1 #include<iostream>
 2 #include<cmath>
 3 #define M 100005
 4 int rmq1[M][30];
 5 int rmq2[M][30];
 6 int a[M];
 7 int N,Q,R,L;
 8 using namespace std;
 9 int max(int a,int b)
10 {
11     if(a>b) return a;
12     else   return b;
13 }
14 int min(int a,int b)
15 {
16     if(a<b) return a;
17     else    return b;
18 }
19 void DP(int N)
20 {
21     int i,j,k;
22     for(i=1;i<=int(log(double(N))/log(2.0));i++)
23         for(j=1;j<=N;j++)
24         {
25             rmq1[j][i]=max(rmq1[j][i-1],rmq1[j+int(pow(2.0,i-1))][i-1]);
26             rmq2[j][i]=min(rmq2[j][i-1],rmq2[j+int(pow(2.0,i-1))][i-1]);
27         }
28 }
29 int search(int i,int j)
30 {
31      int k=(int)(log((double)(j-i+1))/log(2.0));
32      return max(rmq1[i][k],rmq1[j-(1<<k)+1][k])-min(rmq2[i][k],rmq2[j-(1<<k)+1][k]);
33 }
34 int main()
35 {
36     int i,j,k;
37     scanf("%d%d",&N,&Q);
38     for(i=1;i<=N;i++)
39     {
40         scanf("%d",&a[i]);
41         rmq1[i][0]=a[i]; 
42         rmq2[i][0]=a[i];                                                  //初始化
43     }
44     DP(N);
45     while(Q--)
46     {
47         scanf("%d%d",&i,&j);
48         printf("%d\n",search(i,j));
49     }
50 }
51 

posted on 2010-04-25 23:12 M.J 閱讀(1554) 評論(2)  編輯 收藏 引用

評論

# re: POJ 2762. Balanced Lineup (區(qū)間求最值sparsetable算法)  回復(fù)  更多評論   

算法基于倍增的思想。 以空間換取時(shí)間的典范(空間開銷O(nlogn))。
擴(kuò)展閱讀: LCA與RMQ的相互轉(zhuǎn)換等。
2010-04-28 22:45 | Ocean

# re: POJ 2762. Balanced Lineup (區(qū)間求最值sparsetable算法)  回復(fù)  更多評論   

謝謝指點(diǎn)~@Ocean
2010-04-29 16:41 | M.J

只有注冊用戶登錄后才能發(fā)表評論。
網(wǎng)站導(dǎo)航: 博客園   IT新聞   BlogJava   博問   Chat2DB   管理


青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <ins id="pjuwb"></ins>
    <blockquote id="pjuwb"><pre id="pjuwb"></pre></blockquote>
    <noscript id="pjuwb"></noscript>
          <sup id="pjuwb"><pre id="pjuwb"></pre></sup>
            <dd id="pjuwb"></dd>
            <abbr id="pjuwb"></abbr>
            久久久欧美精品| 中文高清一区| 久久精品日韩欧美| 久久久久中文| 亚洲一区二区三区在线视频| 亚洲美女视频在线观看| 国产精品视频午夜| 亚洲日本va午夜在线电影| 免费成人在线观看视频| 国产精品毛片a∨一区二区三区|国| 亚洲国产一区在线| 国产深夜精品| 欧美午夜精品电影| 亚洲一区免费网站| 欧美一区2区三区4区公司二百| 国产日韩亚洲欧美| 欧美一二三视频| 国产精品爽黄69| 欧美日韩高清免费| 99re6热只有精品免费观看| 亚洲欧美一区二区三区极速播放| 亚洲电影一级黄| 黑人极品videos精品欧美裸| 国产精品三级久久久久久电影| 欧美日韩午夜视频在线观看| 久久福利精品| 亚洲国产成人午夜在线一区| 久久av一区二区| 亚洲自拍都市欧美小说| 亚洲欧美一区二区三区极速播放| 欧美在线中文字幕| 欧美一区二区三区在线观看| 亚洲一区免费网站| 麻豆精品网站| 欧美国产先锋| 免费久久99精品国产自| 欧美日韩国产在线播放| 久久综合久久久| 欧美成人亚洲| 一区二区在线视频| 樱桃视频在线观看一区| 亚洲国产成人久久综合一区| 国产日韩欧美一区在线| 亚洲美女中文字幕| 亚洲免费在线看| 欧美激情国产高清| 玖玖综合伊人| 精品不卡一区二区三区| 午夜精品久久久久影视 | 欧美在线视频一区二区三区| 久久国产精品99国产| 久久综合久久综合久久综合| 99国产精品国产精品久久| 日韩一区二区高清| 快射av在线播放一区| 欧美午夜不卡视频| 亚洲美女色禁图| 欧美顶级少妇做爰| 久久久久久一区| 欧美母乳在线| 亚洲一区精品视频| 亚洲一区在线免费观看| 亚洲缚视频在线观看| 欧美黄色一区| 国产日韩精品在线| 久久国产精品一区二区| 久久成人久久爱| 在线日韩av永久免费观看| 一区二区三区国产在线| 久久―日本道色综合久久| 亚洲免费观看高清完整版在线观看熊| 亚洲免费在线视频一区 二区| 亚洲网站在线播放| 欧美成人乱码一区二区三区| 亚洲最新视频在线| 欧美成人在线免费视频| 在线亚洲一区二区| 亚洲国产影院| 久久精品国产精品| 亚洲丝袜av一区| 欧美大片va欧美在线播放| 一区二区三区波多野结衣在线观看| 欧美女同视频| 亚洲最新在线| 久久久久久夜| 久久视频一区| 欧美国产欧美亚州国产日韩mv天天看完整| 亚洲精选国产| 欧美中文在线观看| 国产精品一区一区| 一区二区亚洲| 久久精品国产第一区二区三区最新章节 | 狠狠综合久久| 久久综合九色综合久99| 欧美噜噜久久久xxx| 久久久夜精品| 国产欧美日韩伦理| 亚洲国产精品激情在线观看| 国产精品入口日韩视频大尺度| 日韩午夜精品| 亚洲国产精彩中文乱码av在线播放| 亚洲欧美色婷婷| 免费一级欧美片在线播放| 欧美区一区二| av不卡在线观看| 亚洲视频在线播放| 久久久91精品| 久久精品国产亚洲精品| 亚洲欧美日韩人成在线播放| 欧美精品一区二区三区一线天视频| 六月婷婷久久| 亚洲毛片在线观看| 亚洲欧美日韩一区二区三区在线| 欧美日韩在线看| 亚欧成人在线| 亚洲免费成人av电影| 久久婷婷一区| 夜夜爽www精品| 欧美成人69av| 久久成人久久爱| 99在线精品视频| 午夜电影亚洲| 亚洲人成77777在线观看网| 久久免费国产| 亚洲欧美在线另类| 一本大道久久a久久综合婷婷| 欧美日本在线| 久热精品在线视频| 欧美gay视频激情| 欧美精品免费视频| 欧美另类极品videosbest最新版本| 国产精品白丝jk黑袜喷水| 免费欧美在线| 美女在线一区二区| 欧美精品福利| 欧美日韩精选| 免费不卡在线观看av| 欧美韩国一区| 欧美韩日一区二区三区| 欧美精品色网| 欧美精品偷拍| 国产精品毛片大码女人| 激情懂色av一区av二区av| 在线成人小视频| 国产精品视频精品| 国产色综合久久| 亚洲第一页中文字幕| 中日韩美女免费视频网址在线观看 | 久久国产欧美精品| 国产精品老女人精品视频| 欧美精品一区二区精品网| 欧美在线影院在线视频| 欧美涩涩网站| 欧美日韩日本网| 久久久久久午夜| 亚洲福利视频一区| 欧美中在线观看| 久久精品国产99精品国产亚洲性色| 欧美国产日韩在线| 黄色工厂这里只有精品| 亚洲一级免费视频| 日韩一区二区精品视频| 另类专区欧美制服同性| 欧美精品系列| 国产亚洲美州欧州综合国| 亚洲影视在线| 美女免费视频一区| 欧美凹凸一区二区三区视频| 欧美影院久久久| 美女图片一区二区| 亚洲第一毛片| 欧美va天堂在线| 久久精品视频在线看| 国产日韩欧美在线视频观看| 国产精品99久久久久久久女警| 蜜乳av另类精品一区二区| 久久久久99| 在线欧美日韩国产| 久久精品一区二区三区不卡| 久久精品首页| 国产精品99久久久久久久vr| 亚洲一区亚洲二区| 国产精品国产福利国产秒拍| 亚洲国产一区二区三区青草影视| 另类国产ts人妖高潮视频| 欧美一区二区免费观在线| 欧美日韩视频不卡| 性色av一区二区三区在线观看| 羞羞答答国产精品www一本| 国产精品福利久久久| 久久久久久久一区| 久久综合999| 一区二区三区你懂的| 国产精品久久久久久久免费软件 | 久久中文字幕一区| 久久精品国产精品 | 一本久久综合亚洲鲁鲁五月天 | 久久久久久97三级| 欧美日韩一区二区三区视频 | 久久久精品国产免费观看同学|