• <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)交流

               :: 首頁(yè) :: 新隨筆 :: 聯(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)  編輯 收藏 引用 所屬分類: 我的開源庫(kù)

            評(píng)論

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

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

            # re: findAnWithDegreeOfTwo(計(jì)算度數(shù)為2的齊次遞歸數(shù)列的第n項(xiàng)) 2011-04-11 17:47 Singa
            @今晚打老虎
            就這個(gè)問(wèn)題的話
            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)論
              

            精品亚洲综合久久中文字幕| 久久精品青青草原伊人| 久久久久久精品久久久久| 国产精自产拍久久久久久蜜| 99久久精品午夜一区二区| 男女久久久国产一区二区三区| 2021国产精品久久精品| 亚洲人成网站999久久久综合 | 久久精品一区二区影院| 亚洲国产成人久久精品影视| 青青国产成人久久91网| 91精品国产91久久| 国产无套内射久久久国产| 久久免费观看视频| 中文字幕精品无码久久久久久3D日动漫| 国产精品综合久久第一页| 欧美久久综合九色综合| 国产亚洲美女精品久久久2020| 新狼窝色AV性久久久久久| 国内精品伊人久久久久av一坑| 亚洲伊人久久大香线蕉苏妲己| 日本久久久精品中文字幕| 久久中文精品无码中文字幕| 7777精品伊人久久久大香线蕉| 国产69精品久久久久久人妻精品| 久久国产精品77777| 精品国产91久久久久久久a| 性高湖久久久久久久久AAAAA| 亚洲乱码中文字幕久久孕妇黑人| 久久久久久免费一区二区三区| 久久亚洲2019中文字幕| 久久久久久亚洲Av无码精品专口| 亚洲国产成人久久精品动漫| 中文国产成人精品久久亚洲精品AⅤ无码精品 | 精品久久久久久无码中文野结衣| 开心久久婷婷综合中文字幕| 无码久久精品国产亚洲Av影片| 91精品国产91久久久久久蜜臀 | 亚洲综合伊人久久大杳蕉| 久久免费高清视频| 久久精品人妻中文系列|