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

            uva 10061 - How many zero's and how many digits ?

               這又是一個數學題,不過我還是比較喜歡做這類數學雜題的。題目意思很簡單,給2個十進制數,n和b。如果用b進制表示n!,
            需要多少位數,這個表示末尾會有多少個0。這個題并不需要什么高深的知識,這一點也是我喜歡做這類題的一個方法。
            大家顯然都知道求n!用10進制表示末尾會有多少個0的方法,就是求2*5最多有多少對。那么,b進制了。方法類似,發散一下想法而已。
               我還是先說求多少位數的方法吧。 b的m-1次 <= n!<= b的m次(PS,這個不等式如果把b換成10大家一定會明白的),
            看到這個不等式應該有想法了吧。兩邊同時取logb,就可以得到Σlogi(1<=i<=n) <= m,m直接就求出來了。m即是位數。
               再說怎么求末尾0的,發散下想法,我們也可以對n!中的每個因子試著求b的因子對,一共有多少對。但是,后面發現這樣不行,
            因為比如b是16,1和16是一對因子,2和8是一對因子,4和4是一對因子,也就是因為2也是4的因子,這樣計算因子對就會重復了。
            但是對于b等于10的情況,可以例外而已。
               呵呵,考慮其它的方法。求素數因子。任何數字都可以被分解為一些素數因子的乘積,這是毋容置疑的。那么,我們去分解n!中的
            小于等于b的素數因子,并將其個數存儲在數組里面。然后死循環的去分解b的素數因子,能夠完全分解一次
            (具體看下代碼,不好描述),ans就加1。
            否則,已經沒有末尾0了。
               雖然提交了16次才過。不過最后還算不錯吧,只用了0.508s。相比20s的時間界限,很小了。網上有些過的代碼,跑一個數據都要
            幾分鐘。。。PS:uva上那些神人,怎么用0.0s算出來的,很難想象啊。
               這個題目還有個很大的需要注意的地方,就是浮點數的精度問題。前面講到求位數需要用到log函數,log函數的計算精度就出問題了。
            最后需要對和加一個1e-9再floor才能過。特別需要注意這一點,因為開始我的程序過了所有的
            http://online-judge.uva.es/board/viewtopic.php?f=9&t=7137&start=30上說的數據還是wa了。而且我還發現log10計算精度高很多,如果log(這個是自然對數)去計算,這個網站上的數據都過不了。

            #include <stdio.h>
            #include <string.h>
            #include <math.h>

            int nN, nB;

            int nDivisor[1000];

            int GetDigit(int nN, int nB)
            {
                double fSum = 0.0;
                for (int i = 2; i <= nN; ++i)
                {
                    fSum += log10(i);
                }
                
                fSum /= log10(nB);
               
                return floor(fSum + 1e-9) + 1;
            }

            int GetZero(int nN, int nB)
            {
                memset(nDivisor, 0, sizeof(nDivisor));

                for (int i = 2; i <= nN; ++i)
                {
                    int nTemp = i;
                    
                    for (int j = 2; j <= nTemp && j <= nB; ++j)//這樣循環就可以進行素數因子分解了
                    {
                        while (nTemp % j == 0)
                        {
                            nDivisor[j]++;
                            nTemp /= j;
                        }
                    }
                }
                
                int nAns = 0;
                
                while (1)
                {
                    int nTemp = nB;
                    
                    for (int j = 2; j <= nTemp; ++j)//分解nB
                    {
                        while (nTemp % j == 0)
                        {
                            if (nDivisor[j] > 0)//如果還可以繼續分解
                            {
                                --nDivisor[j];
                            }
                            else //直接可以goto跳出多重循環了
                            {
                                goto out;
                            }
                            nTemp /= j;
                        }
                    }
                    ++nAns;
                }
                
            out:
                return nAns;
            }

            int main()
            {
                while (scanf("%d%d", &nN, &nB) == 2)
                {
                    int nDigit = GetDigit(nN, nB);
                    int nZero = GetZero(nN, nB);
                    printf("%d %d\n", nZero, nDigit);
                }

                return 0;
            }

            posted on 2012-05-02 20:05 yx 閱讀(1846) 評論(0)  編輯 收藏 引用 所屬分類: 數學題

            <2012年5月>
            293012345
            6789101112
            13141516171819
            20212223242526
            272829303112
            3456789

            導航

            統計

            公告

            常用鏈接

            留言簿(3)

            隨筆分類

            隨筆檔案

            me

            好友

            同學

            網友

            搜索

            最新評論

            閱讀排行榜

            評論排行榜

            99久久精品免费看国产一区二区三区 | 亚洲国产另类久久久精品| 91久久福利国产成人精品| 久久国产成人| 少妇被又大又粗又爽毛片久久黑人| 亚州日韩精品专区久久久| 久久人人爽人人爽人人片av麻烦| 欧美丰满熟妇BBB久久久| 亚洲国产精品久久久久婷婷软件| 久久久久亚洲AV成人网| 伊人久久大香线蕉综合影院首页| 久久精品欧美日韩精品| 久久久久一级精品亚洲国产成人综合AV区| 中文精品99久久国产 | 国产精品成人99久久久久91gav| 久久精品无码免费不卡| 无码国内精品久久人妻蜜桃 | 久久综合亚洲色HEZYO社区| 精品免费久久久久久久| 无码人妻久久一区二区三区蜜桃| 久久精品中文騷妇女内射| 亚洲国产成人精品女人久久久| 九九久久自然熟的香蕉图片| 久久婷婷是五月综合色狠狠| 久久综合九色综合精品| 久久精品蜜芽亚洲国产AV| 久久久久亚洲精品日久生情| 国内精品伊人久久久久网站| 国产精品美女久久久久| 久久精品www人人爽人人| 亚洲国产精品无码久久98| 亚洲午夜精品久久久久久app| 精品无码人妻久久久久久| 91精品国产高清久久久久久国产嫩草 | 中文字幕乱码人妻无码久久| 久久亚洲2019中文字幕| 国产成人精品久久| 岛国搬运www久久| 久久精品无码一区二区三区日韩| 91麻精品国产91久久久久 | 久久九九青青国产精品|