青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

KISS(Keep It Simple, Standard)

  C++博客 :: 首頁 :: 聯系 :: 聚合  :: 管理
  10 Posts :: 0 Stories :: 24 Comments :: 0 Trackbacks

常用鏈接

留言簿(10)

我參與的團隊

搜索

  •  

最新評論

閱讀排行榜

評論排行榜

本算法只采用移位、加減法、判斷和循環實現,因為它不需要浮點運算,也不需要乘除運算,因此可以很方便地運用到各種芯片上去。

我們先來看看10進制下是如何手工計算開方的。
先看下面兩個算式,
x = 10*p + q  (1)
公式(1)左右平方之后得:
x^2 = 100*p^2 + 20pq + q^2 (2)
現在假設我們知道x^2和p,希望求出q來,求出了q也就求出了x^2的開方x了。
我們把公式(2)改寫為如下格式:
q = (x^2 - 100*p^2)/(20*p+q) (3)

這個算式左右都有q,因此無法直接計算出q來,因此手工的開方算法和手工除法算法一樣有一步需要猜值。

我們來一個手工計算的例子:計算1234567890的開方

首先我們把這個數兩位兩位一組分開,計算出最高位為3。也就是(3)中的p,最下面一行的334為余數,也就是公式(3)中的(x^2 - 100*p^2)近似值
    3
  ---------------
 / 12 34 56 78 90
    9
  ---------------
 /  3 34

下面我們要找到一個0-9的數q使它最接近滿足公式(3)。我們先把p乘以20寫在334左邊:
                           3  q
                         ---------------
                        / 12 34 56 78 90
                           9
                         ---------------
(20*3+q)*q      /  3 34

我們看到q為5時(60+q)*q的值最接近334,而且不超過334。于是我們得到:
      3  5
    ---------------
   / 12 34 56 78 90
      9
    ---------------
65 /  3 34
      3 25
    ---------------
         9 56

接下來就是重復上面的步驟了,這里就不再啰嗦了。

這個手工算法其實和10進制關系不大,因此我們可以很容易的把它改為二進制,改為二進制之后,公式(3)就變成了:
q = (x^2 - 4*p^2)/(4*p+q) (4)

我們來看一個例子,計算100(二進制1100100)的開方:
       1  0  1  0
      -----------
     / 1 10 01 00
       1
      -----------
 100 / 0 10
       0 00
      -----------
1001 /   10 01
         10 01
      -----------
          0 00

這里每一步不再是把p乘以20了,而是把p乘以4,也就是把p右移兩位,而由于q的值只能為0或者1,所以我們只需要判斷余數(x^2 - 4*p^2)和(4*p+1)的大小關系,如果余數大于等于(4*p+q)那么該上一個1,否則該上一個0。

下面給出完成的C語言程序,其中root表示p,rem表示每步計算之后的余數,divisor表示(4*p+1),通過a>>30取a的最高 2位,通過a<<=2將計算后的最高2位剔除。其中root的兩次<<1相當于4*p。程序完全是按照手工計算改寫的,應該不難理解。
unsigned short sqrt(unsigned long a){
  unsigned long rem = 0;
  unsigned long root = 0;
  unsigned long divisor = 0;
  for(int i=0; i<16; ++i){
    root <<= 1;
    rem = ((rem << 2) + (a >> 30));
    a <<= 2;
    divisor = (root<<1) + 1;
    if(divisor <= rem){
      rem -= divisor;
      root++;
    }
  }
  return (unsigned short)(root);
}

posted on 2008-01-23 14:21 QUIRE-0216 閱讀(5320) 評論(1)  編輯 收藏 引用 所屬分類: Arithmetic(算法)

Feedback

# re: 利用移位、加減法實現整數開平方算法的方法(轉) 2008-01-23 14:46 QUIRE-0216
為了大家能理解我把上面 1234567890 給做完!
              3 5 q
               ---------------
               / 12 34 56 78 90
               9
               ---------------
               65 / 3 34
               3 25
               ---------------
(20*35+q)*q /  9 56
我們看到q為1時(700+q)*q的值最接近956,而且不超過956。于是我們得到:
              3 5 1 q
               ---------------
               / 12 34 56 78 90
               9
               ---------------
               65 / 3 34
               3 25
               ---------------
701 /   9 56
7 01
----------------
(20*351+q)*q / 2 55 78

我們看到q為3時(20*351+q)*q的值最接近25578,而且不超過25578。于是我們得到:

              3 5 1 3 q
               ---------------
               / 12 34 56 78 90
               9
               ---------------
               65 / 3 34
               3 25
               ---------------
701 /   9 56
7 01
----------------
7023 / 2 55 78
2 10 69
----------------
(20*3513+q)*q / 45 0990

我們看到q為6時(20*3513+q)*q的值最接近450990,而且不超過450990。于是我們得到:
              3 5 1 3 6
               ---------------
               / 12 34 56 78 90
               9
               ---------------
               65 / 3 34
               3 25
               ---------------
701 /   9 56
7 01
----------------
7023 / 2 55 78
2 10 69
----------------
70266 / 45 0990
42 1596
----------------
2 9394

至此1234567890的根為35136.我想能看明白吧!




  回復  更多評論
  

青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            欧美另类99xxxxx| 国产精品免费网站| 亚洲激精日韩激精欧美精品| 玖玖在线精品| 欧美77777| 99国产精品久久久久久久久久 | 一区二区三区视频在线观看| 91久久综合| 欧美片第一页| 先锋a资源在线看亚洲| 性欧美暴力猛交69hd| 狠久久av成人天堂| 欧美激情一区二区三区蜜桃视频| 欧美成人精品一区| 亚洲香蕉网站| 久久精品91| 一区二区黄色| 欧美一区午夜视频在线观看| 亚洲人午夜精品免费| 亚洲性感激情| 一区二区在线观看视频| 亚洲区在线播放| 国产亚洲永久域名| 亚洲欧洲一区二区三区在线观看| 欧美婷婷在线| 男女精品网站| 国产精品九九久久久久久久| 六十路精品视频| 欧美四级伦理在线| 免费观看国产成人| 国产精品伦一区| 欧美成人午夜77777| 国产精品素人视频| 亚洲国产高清在线| 国产视频欧美视频| 999亚洲国产精| 亚洲高清不卡在线观看| 亚洲自拍偷拍色片视频| 一道本一区二区| 另类av一区二区| 久久九九精品| 欧美性做爰毛片| 亚洲精品社区| 亚洲黄色尤物视频| 欧美一区日韩一区| 篠田优中文在线播放第一区| 欧美人与禽猛交乱配视频| 麻豆精品国产91久久久久久| 国产农村妇女精品| 一区二区三区 在线观看视| 亚洲日本va午夜在线影院| 久久国产精品久久精品国产 | 1204国产成人精品视频| 亚洲欧美精品伊人久久| 亚洲一区二区视频在线| 欧美精品一区二| 亚洲国产精品久久人人爱蜜臀| 国内精品久久久久久久97牛牛| 亚洲一区免费看| 亚洲欧美日韩在线| 国产精品久久久免费| 一本色道久久综合精品竹菊| 日韩午夜电影av| 欧美精品午夜| 亚洲三级国产| 亚洲视频一区二区| 欧美视频日韩| 中文精品视频一区二区在线观看| 亚洲视频国产视频| 欧美系列电影免费观看| 一区二区三区视频在线播放| 亚洲一区二区三区在线| 国产精品久久久久9999吃药| 亚洲一区二区三区视频播放| 亚洲欧美一区二区精品久久久| 国产精品高清一区二区三区| 中国女人久久久| 欧美中在线观看| 韩国女主播一区| 免费永久网站黄欧美| 91久久夜色精品国产网站| 亚洲天堂av在线免费观看| 国产精品乱人伦一区二区| 欧美在线视频导航| 麻豆精品在线观看| 日韩午夜激情电影| 在线一区欧美| 国产综合av| 国产精品日韩精品欧美在线| 亚洲尤物视频网| 亚洲影院一区| 久久av一区二区三区| 伊人久久亚洲影院| 欧美激情中文字幕一区二区| 亚洲丝袜av一区| 免费成人av资源网| 99精品国产高清一区二区| 国产精品久久久久久久久久ktv| 午夜精品国产更新| 亚洲二区视频| 午夜欧美不卡精品aaaaa| 一区二区三区在线免费观看| 欧美激情中文字幕乱码免费| 午夜精品影院| 亚洲精品国产精品国自产在线 | 免费看av成人| 亚洲一区二区免费| 亚洲第一网站| 久久精品国产99国产精品澳门| 亚洲国产成人精品女人久久久| 欧美手机在线视频| 久久亚洲春色中文字幕久久久| 一区二区欧美在线观看| 欧美成人一区二区在线| 午夜精品久久久久久久| 亚洲精品久久久久久久久久久久| 国产精品亚洲欧美| 欧美二区在线观看| 欧美一级网站| 亚洲亚洲精品在线观看 | 亚洲一区二区三区在线视频| 伊人久久亚洲美女图片| 国产欧美大片| 国产精品激情偷乱一区二区∴| 六月婷婷一区| 久久精品人人做人人爽电影蜜月| 这里是久久伊人| 亚洲精品久久久久久下一站| 欧美成人一区二区三区片免费| 久久精品亚洲精品| 午夜国产一区| 亚洲欧美日韩精品在线| 宅男精品视频| 日韩视频―中文字幕| 亚洲精品乱码久久久久久日本蜜臀 | 免费中文日韩| 久久人体大胆视频| 久久精品视频播放| 久久精品国产第一区二区三区最新章节| 亚洲精品国精品久久99热| 欧美国产一区视频在线观看| 狂野欧美性猛交xxxx巴西| 久久久久久**毛片大全| 久久精品国产欧美亚洲人人爽| 香蕉亚洲视频| 欧美在线亚洲一区| 久久国产视频网| 久久精品99久久香蕉国产色戒| 欧美在线关看| 久久亚洲一区| 欧美成人免费在线| 亚洲国产精品一区二区尤物区| 欧美激情va永久在线播放| 亚洲国产欧美精品| 亚洲精品在线免费观看视频| 99国产精品国产精品久久| 国产精品99久久久久久久vr| 亚洲一区二区三区高清 | 欧美成年人视频网站欧美| 免费成人美女女| 亚洲高清资源| 夜夜爽www精品| 校园激情久久| 免费成人网www| 欧美日韩国产探花| 国产乱肥老妇国产一区二| 国内一区二区三区在线视频| 亚洲激情在线播放| 中日韩高清电影网| 久久久www| 亚洲精品国精品久久99热| 制服丝袜亚洲播放| 久久丁香综合五月国产三级网站| 免费观看亚洲视频大全| 欧美午夜免费| 亚洲第一色在线| 亚洲小视频在线观看| 久久综合成人精品亚洲另类欧美| 亚洲国产小视频在线观看| 亚洲一区在线直播| 久久这里只精品最新地址| 欧美午夜精品久久久久久人妖 | 欧美日韩国产美女| 国产亚洲精品福利| 99视频一区二区三区| 欧美一区二区在线观看| 亚洲国产一区二区三区在线播| 亚洲在线播放电影| 欧美激情一区二区三区在线| 国产亚洲精品bt天堂精选| 在线亚洲一区| 欧美高清视频在线观看| 午夜欧美大尺度福利影院在线看| 欧美国产精品中文字幕| 国产一区视频观看| 欧美亚洲免费高清在线观看| 亚洲精品国产精品久久清纯直播| 欧美在线91| 国产乱码精品一区二区三区五月婷 | 国产精品影视天天线|