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

            5D空間

            學(xué)習(xí)總結(jié)與經(jīng)驗(yàn)交流

               :: 首頁 :: 新隨筆 :: 聯(lián)系 :: 聚合  :: 管理 ::
              遞歸數(shù)列a(n) = r1*a(n-1) + r2*a(n-2)
              其中r1 r2是常數(shù),a1 a2已知,按照順序a1 a2 r1 r2 n輸入?yún)?shù),則返回這樣一個(gè)遞歸數(shù)列的第n項(xiàng)值。類型均為double

            #ifndef FINDAN_H
            #define FINDAN_H

            #include 
            <cmath>

            // a(n) = r1*a(n-1) + r2*a(n-2), give the a1, a2, r1, r2 and n
            double findAnWithDegreeOfTwo( double a1, double a2, double r1, double r2, int n )
            {
                
            double x1;
                
            double x2;
                
            double u1;
                
            double u2;
                
            double det = r1*r1 + 4*r2;
                
            double result;

                
            if ( det < 0 )
                    
            return 0;
                
            else if ( det > 0 )
                
            {
                    det 
            = sqrt( det );
                    x1 
            = ( r1 + det ) / 2;
                    x2 
            = ( r1 - det ) / 2;
                    u1 
            = ( a1*x2 - a2 ) / ( x1*( x2 - x1 ) );
                    u2 
            = ( a2*x1 - a1 ) / ( x2*( x1 - x2 ) );
                    result 
            = u1*pow( x1, n ) + u2*pow( x2, n );
                    
            return result;
                }

                
            else
                
            {
                    x1 
            = r1 / 2;
                    u2 
            = ( a2 - x1*a1 ) / x1*x1;
                    u1 
            = a1 / x1 - u2;
                    result 
            = ( u1 + u2*n ) * pow( x1, n );
                    
            return result;
                }

            }


            #endif
            posted on 2011-04-03 19:53 今晚打老虎 閱讀(1026) 評(píng)論(4)  編輯 收藏 引用 所屬分類: 我的開源庫

            評(píng)論

            # re: findAnWithDegreeOfTwo(計(jì)算度數(shù)為2的齊次遞歸數(shù)列的第n項(xiàng)) 2011-04-05 15:24 Singa
            這樣做的時(shí)間復(fù)雜度是O(n),可以利用矩陣乘法來做,把時(shí)間復(fù)雜度降低到O(logn)。  回復(fù)  更多評(píng)論
              

            # re: findAnWithDegreeOfTwo(計(jì)算度數(shù)為2的齊次遞歸數(shù)列的第n項(xiàng)) 2011-04-08 10:00 今晚打老虎
            @Singa
            矩陣乘法么,沒學(xué)過耶,能指點(diǎn)一下么  回復(fù)  更多評(píng)論
              

            # re: findAnWithDegreeOfTwo(計(jì)算度數(shù)為2的齊次遞歸數(shù)列的第n項(xiàng)) 2011-04-11 17:47 Singa
            @今晚打老虎
            就這個(gè)問題的話
            a[n] a[n-1] r1 1 a[n+1] a[n]
            a[n-1] a[n-2] r2 0 a[n] a[n-1]
            上面分別是三個(gè)2*2的矩陣,前兩個(gè)相乘得到第三個(gè)
            這樣子就得到一個(gè)類似等比數(shù)列的遞推式了
            r1 1
            r2 0 類似于公比

            a[2] a[1]
            a[1] a[0] 類似于數(shù)列首項(xiàng)
            然后算公比的n次方的時(shí)候用下快速冪取模,就實(shí)現(xiàn)O(logn)的復(fù)雜度了  回復(fù)  更多評(píng)論
              

            # re: findAnWithDegreeOfTwo(計(jì)算度數(shù)為2的齊次遞歸數(shù)列的第n項(xiàng)) 2011-04-12 00:24 今晚打老虎
            @Singa
            感謝指點(diǎn)!  回復(fù)  更多評(píng)論
              

            久久人人爽人人爽人人片AV高清 | 久久水蜜桃亚洲av无码精品麻豆| 青青草原综合久久大伊人导航| 欧美一级久久久久久久大片| 久久国产成人午夜AV影院| 国产精品免费久久久久久久久 | 亚洲精品乱码久久久久久久久久久久 | 久久婷婷国产麻豆91天堂| 久久久一本精品99久久精品88| 亚洲人成无码www久久久| 狠狠狠色丁香婷婷综合久久五月| 四虎国产精品免费久久5151 | 久久经典免费视频| 国产精品久久波多野结衣| 青青青青久久精品国产h久久精品五福影院1421 | 亚洲欧洲中文日韩久久AV乱码| 99久久精品免费看国产一区二区三区| 2021少妇久久久久久久久久| 四虎国产精品成人免费久久| 精品久久8x国产免费观看| 久久精品国产亚洲αv忘忧草 | 久久发布国产伦子伦精品| 久久影视综合亚洲| 亚洲国产精品热久久| 色88久久久久高潮综合影院| 亚洲国产天堂久久综合| 亚洲国产成人久久综合一| 久久综合综合久久综合| 久久久久久久女国产乱让韩| 久久免费视频一区| 国产A级毛片久久久精品毛片| 无码人妻久久一区二区三区| 久久这里有精品视频| 久久精品这里只有精99品| 伊人久久综合热线大杳蕉下载| 人妻丰满AV无码久久不卡| 国产精品久久久久免费a∨| 久久久久久久亚洲精品| 国产精品99久久不卡| 国产激情久久久久影院小草 | 久久一区二区三区免费|