• <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 今晚打老虎 閱讀(1034) 評(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)論
              

            欧美亚洲另类久久综合婷婷| 精品久久久久久无码国产| 久久精品国产色蜜蜜麻豆| 久久免费看黄a级毛片| 久久亚洲春色中文字幕久久久| 国产精品久久久久久影院| 久久久久亚洲精品无码网址 | 欧美久久久久久精选9999| 无码人妻久久一区二区三区蜜桃| 新狼窝色AV性久久久久久| 国产精品永久久久久久久久久| 区久久AAA片69亚洲| 99久久精品无码一区二区毛片 | 欧美大战日韩91综合一区婷婷久久青草 | 99久久免费国产精精品| 亚洲国产成人久久综合区| 91精品国产9l久久久久| 久久婷婷五月综合97色直播| 国产激情久久久久影院老熟女免费| 午夜精品久久久久9999高清| 久久青草国产手机看片福利盒子| 99精品久久久久久久婷婷| 久久久久久久综合日本| 亚洲精品国产成人99久久| 99久久久国产精品免费无卡顿| 2021久久精品免费观看| 四虎久久影院| 亚洲国产精品成人久久蜜臀| 国产99久久九九精品无码| 粉嫩小泬无遮挡久久久久久| 热99RE久久精品这里都是精品免费| 精品国产91久久久久久久| 亚洲精品乱码久久久久久自慰 | 精品人妻久久久久久888| 色8久久人人97超碰香蕉987| 久久综合偷偷噜噜噜色| 伊人精品久久久久7777| 99久久国产宗和精品1上映| 欧美亚洲色综久久精品国产| 欧美牲交A欧牲交aⅴ久久 | 国产精品美女久久福利网站|