• <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 107 - The Cat in the Hat

               這是一個很神的數(shù)學題吧。基本上過這個題的很多都會wa10多次,而且這個題好像簡單的枚舉其中的一個指數(shù)值都能過,可能是
            數(shù)據(jù)量比較小。
               但是,這個題還是有數(shù)學的解法的。但是,即使找到了這個正確的解法,過題的話,也是一件很困難的事情。題意大致如下:一只貓,
            高度為H,戴了一個帽子,帽子里面有N只貓(N是常數(shù),且未知),同樣帽子里面的貓也戴了帽子,但是這些貓的高度變成了H / (N + 1),
            會向下取整。以此遞歸下去,直到最后的貓的高度都為1為止。現(xiàn)在,給出H和高度為1的貓的數(shù)量。要求的是高度大于1的貓的數(shù)量,
            以及所有貓的高度之和。
               很別扭吧。通過上面的信息,得出2個式子。假設one代表為高度為1的貓的數(shù)量。one = N的n次。H >= (N + 1)的n次。注意第
            二個式子不一定取等號,因為很多時候都是不能整除的。現(xiàn)在要求N和n。2個方程解2個未知數(shù),應該能解出來。但是,注意的是其中
            一個還是不等式。。。
               指數(shù)關系很多時候會轉(zhuǎn)換為對數(shù)的關系。所以,繼續(xù)求對數(shù),有l(wèi)gH >= n * lg(N + 1)。其中,由第一個式子可以得到n = lg(one)
            / lg(N)。那么最終轉(zhuǎn)換為:lgH >= (lg(one) / lgN) * lg(N + 1)。換個形式就是lgH / lg(One) >= lg(N + 1) / lgN。現(xiàn)在,已經(jīng)很
            清晰了。因為,函數(shù)lg(N + 1) / lg(N) 是單調(diào)遞減的。看到單調(diào)的函數(shù),馬上就會知道可以二分了。意思是,我們可以二分出一個N讓
             lg(N + 1) / lgN 最接近lgH / lg(One),而且是小于lgH / lg(One)的。剩下的工作就只是求和而已了。
               寫二分的時候,有一個地方可以注意一下。因為 lg(N + 1) / lgN 可能會出現(xiàn)除數(shù)為0的情況,所以可以進一步轉(zhuǎn)換為lgH * lgN >=
            lg(N + 1) * lg(one)
            。 也是求一個N讓上面那個不等式2邊的值最接近,而且右邊小于左邊。
               能很快寫對這個題真不是件容易的事情。。。

               代碼如下:
            #include <stdio.h>
            #include <math.h>

            int main()
            {
                int nInitH, nOnes;
                int nN, n;

                while (scanf("%d%d", &nInitH, &nOnes), nInitH + nOnes)
                {
                    int nBeg = 1;
                    int nEnd = nOnes;
                    int nMid;
                
                    while (nBeg <= nEnd)
                    {
                        nMid = (nBeg + nEnd) / 2;
                        
                        double fRes = log10(nInitH) * log10(nMid);
                        double fTemp = log10(nMid + 1) * log10(nOnes);
                        if (fabs(fRes - fTemp) < 1e-10)
                        {
                            //printf("Find nN:%d\n", nMid);
                            nN = nMid;
                            break;
                        }
                        else if (fTemp > fRes)
                        {
                            nBeg = nMid + 1;
                        }
                        else
                        {
                            nEnd = nMid - 1;
                        }
                    }
                    
                    n = floor(log10(nInitH) / log10(nN + 1) + 1e-9);
                    //printf("nN:%d, n:%d\n", nN, n);

                    int nSum = 0;
                    int nLazy = 0;
                    int nNum = 1;
                    for (int i = 0; i <= n; ++i)
                    {
                        nSum += nNum * nInitH;
                        nLazy += nNum;
                        nNum *= nN;
                        nInitH /= (nN + 1);
                    }
                    
                    printf("%d %d\n", nLazy - nOnes, nSum);
                }

                return 0;
            }

               

            posted on 2012-05-07 16:54 yx 閱讀(1685) 評論(0)  編輯 收藏 引用 所屬分類: 數(shù)學題

            <2013年6月>
            2627282930311
            2345678
            9101112131415
            16171819202122
            23242526272829
            30123456

            導航

            統(tǒng)計

            公告

            常用鏈接

            留言簿(3)

            隨筆分類

            隨筆檔案

            me

            好友

            同學

            網(wǎng)友

            搜索

            最新評論

            閱讀排行榜

            評論排行榜

            久久久网中文字幕| 亚洲国产精品嫩草影院久久| 久久国内免费视频| 久久亚洲国产成人影院网站| 久久九九久精品国产| 伊人久久五月天| 久久久国产打桩机| 亚洲AV日韩精品久久久久久| 久久久久99精品成人片试看| 香蕉久久夜色精品升级完成| 国产精品久久久香蕉| 99久久精品这里只有精品| 日韩精品国产自在久久现线拍 | 午夜精品久久久久久久久| 人妻无码αv中文字幕久久| 久久久精品视频免费观看| 日本福利片国产午夜久久| 久久精品国产精品亚洲下载| 亚洲色欲久久久久综合网| 久久精品人人做人人爽电影| 亚洲午夜久久久影院| 青青青青久久精品国产| 亚洲午夜福利精品久久| 久久久久亚洲精品天堂| 亚洲精品tv久久久久| 亚洲午夜无码久久久久| 日本国产精品久久| 亚洲国产成人精品91久久久| 久久免费国产精品| 久久99精品九九九久久婷婷| 久久精品国产91久久麻豆自制| 久久夜色精品国产噜噜噜亚洲AV| 中文字幕乱码久久午夜| AV无码久久久久不卡蜜桃| 亚洲精品无码久久毛片| 一本色道久久综合| 国产成人综合久久精品红| 天天综合久久一二三区| 一级女性全黄久久生活片免费 | 亚洲国产精品综合久久一线| 国产叼嘿久久精品久久|