• <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>

            為生存而奔跑

               :: 首頁 :: 聯(lián)系 :: 聚合  :: 管理
              271 Posts :: 0 Stories :: 58 Comments :: 0 Trackbacks

            留言簿(5)

            我參與的團隊

            搜索

            •  

            積分與排名

            • 積分 - 328403
            • 排名 - 74

            最新評論

            閱讀排行榜

            評論排行榜

                RMQ(Range Minimum/Maximum Query)問題是求區(qū)間最值問題。可以寫一個線段樹,但是預處理和查詢的復雜度都是O(logn)。這里有更牛的算法,就是ST算法,它可以做到O(nlogn)的預處理,O(1)!!!地回答每個詢問。
                來看一下ST算法是怎么實現(xiàn)的(以最大值為例):
                  
                首先是預處理,用一個DP解決。設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)、初值都已經有了,剩下的就是狀態(tài)轉移方程。我們把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))。用上例說明,當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-i),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]對應)。直接給出表達式:

             f[i,j] 表示 從第 i 個數(shù)數(shù) 2^j 中最小的數(shù)
            那么:
                    最大值 f[i,j]=max(f[i,j-1],f[i+2^(j-1),j-1]);
                    最小值 f[i,j]=min(f[i,j-1],f[i+2^(j-1),j-1]);


            模板:
            #include <iostream>
            #include 
            <math.h>
            #define max(a,b) ((a)>(b)?(a):(b))
            #define min(a,b) ((a)<(b)?(a):(b))
             
            using namespace std;
             
            const int maxn=50001;
             
            int h[maxn];
             
            int mx[maxn][16],mn[maxn][16];
            int n,q;
             
             
            void rmq_init()
             
            {
                 
            int i,j;
                 
            for(j=1;j<=n;j++) mx[j][0]=mn[j][0]=h[j];
                 
            int m=floor(log((double)n)/log(2.0));
                 
            for(i=1;i<=m;i++){
                     
            for(j=n;j>0;j--){
                         mx[j][i]
            =mx[j][i-1];
                         
            if(j+(1<<(i-1))<=n) mx[j][i]=max(mx[j][i],mx[j+(1<<(i-1))][i-1]);
                     }

                }

                
            for(i=1;i<=m;i++){
                     
            for(j=n;j>0;j--){
                        mn[j][i]
            =mn[j][i-1];
                         
            if(j+(1<<(i-1))<=n) mn[j][i]=min(mn[j][i],mn[j+(1<<(i-1))][i-1]);
                     }

                 }

             }

             
             
            int rmq(int l,int r)
             
            {
                 
            int m=floor(log((double)(r-l+1))/log(2.0));
                
            int a=max(mx[l][m],mx[r-(1<<m)+1][m]);
                 
            int b=min(mn[l][m],mn[r-(1<<m)+1][m]);
                
            return a-b;   
             }

             
             
            int main()
             
            {
                 
            int i,l,r;
                 scanf(
            "%d%d",&n,&q);
                 
            for(i=1;i<=n;i++) scanf("%d",&h[i]);
                 rmq_init();
                 
            for(i=0;i<q;i++){
                     scanf(
            "%d%d",&l,&r);
                     printf(
            "%d\n",rmq(l,r));
                 }


             }


            posted on 2009-08-06 09:59 baby-fly 閱讀(1617) 評論(2)  編輯 收藏 引用 所屬分類: Algorithm

            Feedback

            # re: ST算法求解RMQ問題 2010-05-01 11:24 Perhack
            寫的很好。。。  回復  更多評論
              

            # re: ST算法求解RMQ問題 2010-08-09 09:10 koma
            大牛 , 貌似這個版本沒線段樹快喲  回復  更多評論
              

            99久久精品国产免看国产一区| 午夜欧美精品久久久久久久| 高清免费久久午夜精品| 国产精品99久久免费观看| 久久99精品国产| 天堂无码久久综合东京热| 欧美噜噜久久久XXX| 久久久久久久国产免费看| 99久久免费国产精品特黄| 国产国产成人精品久久| 伊色综合久久之综合久久| 国产精品九九九久久九九 | 囯产极品美女高潮无套久久久| 久久无码人妻一区二区三区 | 性做久久久久久久| 国产精品99久久不卡| 久久久噜噜噜www成人网| 久久人人爽人人爽AV片| 99久久免费国产精品热| 伊人久久综合精品无码AV专区| 久久97久久97精品免视看秋霞| 色偷偷88888欧美精品久久久| 久久中文精品无码中文字幕| 久久久国产精品福利免费| 久久人人爽人人爽人人片av高请| 久久夜色撩人精品国产小说| 久久亚洲国产午夜精品理论片| 亚洲色大成网站www久久九| 思思久久好好热精品国产| 欧美亚洲另类久久综合婷婷| 九九久久99综合一区二区| av午夜福利一片免费看久久| 亚洲综合精品香蕉久久网| 亚洲香蕉网久久综合影视| 久久天天躁狠狠躁夜夜2020一| 亚洲国产精品一区二区三区久久 | 国产欧美一区二区久久| 久久久久亚洲AV无码麻豆| 人妻精品久久无码专区精东影业 | 欧美噜噜久久久XXX| 久久精品人成免费|