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

M.J的blog

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

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

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

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

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

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

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

 

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

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

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

這樣就計算了從i開始,長度為2^t次的區(qū)間和從r-2^i+1開始長度為2^t的區(qū)間的最大值(表達式比較煩瑣,細(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 閱讀(1562) 評論(2)  編輯 收藏 引用

評論

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

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

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

謝謝指點~@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>
            欧美激情中文字幕乱码免费| 亚洲一区区二区| 久久看片网站| 亚洲第一色中文字幕| 免费人成网站在线观看欧美高清| 欧美亚洲三级| 国内在线观看一区二区三区| 久久综合久久久| 久久综合99re88久久爱| 亚洲人成在线观看一区二区| 亚洲人成人一区二区在线观看 | 亚洲欧美文学| 亚洲一区影院| 一区二区亚洲精品国产| 欧美激情一二区| 欧美久久久久免费| 亚洲自拍三区| 久久九九精品| 一区二区高清视频| 亚洲欧美欧美一区二区三区| 一区二区在线视频播放| 亚洲欧洲精品一区二区精品久久久| 欧美精品www在线观看| 亚洲免费中文| 久久久噜噜噜久久狠狠50岁| 99精品99| 久久久av毛片精品| 一本色道久久| 久久精品国产亚洲一区二区| 日韩视频一区二区三区| 亚洲欧美高清| 亚洲黄色三级| 欧美一区二区三区在线播放| 99国产精品一区| 性欧美xxxx大乳国产app| 91久久精品国产91性色| 亚洲尤物影院| 一本色道久久加勒比88综合| 久久免费视频观看| 亚洲一区二区三区在线| 久久视频一区| 欧美一级精品大片| 欧美日韩国产成人精品| 欧美+日本+国产+在线a∨观看| 国产精品久久久久久久久久免费看 | 红桃视频成人| 一本色道久久加勒比88综合| 亚洲国产第一页| 亚欧成人在线| 亚洲欧美99| 欧美理论视频| 欧美成人精品一区| 国产有码在线一区二区视频| 国产精品99久久不卡二区| 亚洲日本久久| 久久综合久久美利坚合众国| 欧美一区二区三区婷婷月色| 欧美日韩精品高清| 欧美激情影音先锋| 禁断一区二区三区在线| 欧美一区二区三区婷婷月色| 午夜伦理片一区| 国产精品久久久99| 一本大道久久a久久精品综合| 一本色道久久88综合亚洲精品ⅰ| 免费在线观看成人av| 美国十次成人| 在线成人h网| 老司机午夜精品视频在线观看| 久久躁狠狠躁夜夜爽| 狠狠色狠狠色综合日日91app| 亚洲男人的天堂在线观看| 亚洲欧美日韩国产中文 | 亚洲免费中文字幕| 香蕉免费一区二区三区在线观看| 欧美午夜理伦三级在线观看| 9i看片成人免费高清| 在线午夜精品| 国产精品地址| 亚洲欧美中文字幕| 久久免费精品视频| 伊人婷婷欧美激情| 欧美成人午夜激情视频| 最新国产精品拍自在线播放| 亚洲视频福利| 国产精品日韩在线播放| 欧美在线观看视频在线| 麻豆久久久9性大片| 亚洲国产精品一区二区久| 欧美激情1区2区| 夜夜嗨av一区二区三区| 欧美一区永久视频免费观看| 伊人天天综合| 欧美日韩国产页| 先锋影音国产一区| 蜜臀av在线播放一区二区三区| 亚洲伦理精品| 国产美女诱惑一区二区| 久久久久网址| 亚洲精品在线一区二区| 先锋影院在线亚洲| 在线成人欧美| 国产精品v亚洲精品v日韩精品| 欧美有码在线观看视频| 亚洲国产日韩欧美综合久久| 亚洲专区一二三| 在线看国产日韩| 欧美视频在线观看免费网址| 久久精品人人做人人爽| 日韩亚洲欧美中文三级| 久久久国产精彩视频美女艺术照福利| 亚洲国产国产亚洲一二三| 国产精品www| 久久人人97超碰精品888| 9l视频自拍蝌蚪9l视频成人| 乱码第一页成人| 亚洲综合不卡| 亚洲精品免费在线| 国产精品永久入口久久久| 欧美成人小视频| 欧美一区二区在线播放| 9人人澡人人爽人人精品| 蜜桃久久av一区| 亚洲免费视频观看| 亚洲人成网站在线播| 国内精品久久久久久影视8| 国产精品久久久久久久久免费桃花| 免费国产自线拍一欧美视频| 午夜视频一区二区| 亚洲视频成人| 亚洲最新视频在线播放| 亚洲人成在线播放网站岛国| 欧美国产日韩一区| 久久久精彩视频| 欧美一区二区高清在线观看| 亚洲图片欧洲图片日韩av| 亚洲美女视频在线观看| 亚洲国产日韩欧美| 悠悠资源网久久精品| 激情久久综合| 国产在线拍偷自揄拍精品| 国产日本欧美一区二区三区| 国产精品乱子久久久久| 国产精品高潮呻吟视频| 国产精品成人久久久久| 欧美精品啪啪| 欧美精品不卡| 欧美另类久久久品| 欧美激情一区二区三级高清视频| 免费成人高清视频| 老司机一区二区三区| 蜜臀a∨国产成人精品| 欧美a级大片| 欧美日韩国产影院| 欧美肉体xxxx裸体137大胆| 欧美日韩一区综合| 欧美日韩国产综合视频在线观看 | 午夜精品久久久久久久久| 亚洲一区二三| 欧美一区二视频| 久久久久久网址| 美国十次成人| 欧美啪啪成人vr| 国产精品乱码妇女bbbb| 国内精品久久久久久| 亚洲国产日韩在线一区模特| 亚洲国产天堂久久综合网| 99综合精品| 欧美一区国产在线| 麻豆精品精品国产自在97香蕉| 亚洲国产成人av好男人在线观看| 亚洲精品在线观| 亚洲天堂网在线观看| 久久精品国产一区二区三| 欧美成人午夜剧场免费观看| 欧美午夜精彩| 国模吧视频一区| 亚洲精品视频免费观看| 亚洲欧美日韩国产| 麻豆成人在线观看| 亚洲精品乱码久久久久久| 亚洲欧美精品中文字幕在线| 久热精品视频| 国产精品vip| 亚洲电影免费观看高清| 亚洲婷婷综合色高清在线| 久久精品毛片| 日韩视频免费观看| 久久久xxx| 国产精品美女视频网站| 亚洲黄色成人| 欧美一区二区三区四区在线观看地址| 欧美刺激性大交免费视频| 在线视频免费在线观看一区二区| 久久久久久夜精品精品免费| 国产精品啊啊啊| 亚洲精品一区二区三区在线观看| 久久黄色级2电影| 一区二区激情| 欧美激情精品久久久六区热门 |