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

            欧美日韩成人精品久久久免费看| 国产三级久久久精品麻豆三级| 色88久久久久高潮综合影院| 精品久久久久久无码人妻蜜桃| 久久99国产精品久久99| 国内精品久久久久影院日本| 国产精品99久久久精品无码| 国产精品一区二区久久精品涩爱| 久久久这里有精品中文字幕| 久久久精品日本一区二区三区 | 亚洲国产另类久久久精品 | 久久无码AV一区二区三区| 深夜久久AAAAA级毛片免费看| 看全色黄大色大片免费久久久| 日日狠狠久久偷偷色综合0| 亚洲国产精品无码久久久久久曰| 久久久久婷婷| 免费精品久久天干天干| 国产成人久久精品一区二区三区| 久久夜色精品国产噜噜噜亚洲AV| 久久久婷婷五月亚洲97号色| 精品久久一区二区三区| 国产成人精品久久综合| 久久伊人中文无码| 久久婷婷五月综合97色| 99久久人人爽亚洲精品美女| 合区精品久久久中文字幕一区| 精品多毛少妇人妻AV免费久久| 国产亚洲欧美精品久久久| 日本久久久久久中文字幕| 久久人人爽人人爽人人片AV麻豆| 精品国产乱码久久久久久人妻| 精品久久久久久久无码| 久久久久久国产a免费观看不卡| 久久精品国产免费观看| 热久久国产精品| 午夜久久久久久禁播电影| 91久久香蕉国产熟女线看| 久久人人爽人人爽人人片AV不| 国产精品久久免费| 伊人久久大香线蕉AV色婷婷色|