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

            USACO 3.1 Humble Numbers



            每次找下一個最小的humble數直到第n個為止。
            對每一個prime,用二分法找出乘以已有的humble數且>=現有最大的humble數的最小的那一個。
            然后取最小的那一個作為新找到的humble數。
            開始的時候最后一個例子總是超時,后來用一個next_min數組緩存一下結果,效率一下提高了10倍。

            #include?<iostream>
            #include?
            <fstream>

            using?namespace?std;

            int?humbles[100002];
            int?primes[100];
            int?k,n;
            int?top;

            int?next_min[100];

            long?long?find_next(int?i)
            {
            ????
            int?l?=?0;
            ????
            int?h?=?top;
            ????
            int?mid;

            ????
            while(l<=h){
            ????????mid?
            =?(l+h)/2;
            ????????
            if((long?long)humbles[mid]*primes[i]<=humbles[top]){
            ????????????l?
            =?mid+1;
            ????????}
            else{
            ????????????h?
            =?mid-1;
            ????????}
            ????}

            ????
            return?humbles[h+1]*primes[i];
            }

            void?solve()
            {
            ????humbles[
            0]?=?1;
            ????top?
            =?0;

            ????freopen(
            "humble.in","r",stdin);
            ????freopen(
            "humble.out","w",stdout);

            ????scanf(
            "%d%d",&k,&n);

            ????
            for(int?i=0;i<k;++i)
            ????????scanf(
            "%d",&primes[i]);

            ????
            for(int?i=0;i<k;++i)
            ????????next_min[i]?
            =?find_next(i);

            ????
            while(top<n){
            ????????
            long?long?next?=?LONG_MAX;
            ????????
            long?long?tmp;
            ????????
            for(int?i=0;i<k;++i){
            ????????????
            if(next_min[i]<=humbles[top])
            ????????????????tmp?
            =?next_min[i]?=?find_next(i);
            ????????????
            else
            ????????????????tmp?
            =?next_min[i];

            ????????????
            if(tmp<next){
            ????????????????next?
            =?tmp;
            ????????????}
            ????????????
            if(next==humbles[top]+1)
            ????????????????
            break;
            ????????}
            ????????humbles[
            ++top]?=?next;
            ????}

            ????printf(
            "%d\n",humbles[n]);
            }

            int?main(int?argc,char?*argv[])
            {
            ????solve();?
            ????
            return?0;
            }

            posted on 2009-06-30 00:05 YZY 閱讀(1679) 評論(0)  編輯 收藏 引用 所屬分類: Algorithm 、USACO

            導航

            <2009年6月>
            31123456
            78910111213
            14151617181920
            21222324252627
            2829301234
            567891011

            統計

            常用鏈接

            留言簿(2)

            隨筆分類

            隨筆檔案

            搜索

            積分與排名

            最新評論

            閱讀排行榜

            看全色黄大色大片免费久久久| 久久精品国产亚洲AV香蕉| 成人a毛片久久免费播放| 无码乱码观看精品久久| 久久婷婷五月综合色奶水99啪| 久久精品国产精品亚洲人人 | 亚洲午夜无码AV毛片久久| 思思久久好好热精品国产| 亚洲综合久久综合激情久久| 亚洲欧洲日产国码无码久久99| 欧美亚洲国产精品久久高清| 亚洲国产另类久久久精品小说 | 午夜久久久久久禁播电影| 人妻无码精品久久亚瑟影视| 精品久久一区二区| 激情久久久久久久久久| 久久精品三级视频| 伊人久久大香线蕉精品不卡| 久久久久97国产精华液好用吗| 国内精品久久久久影院老司| 色天使久久综合网天天 | 一级做a爰片久久毛片毛片 | 手机看片久久高清国产日韩| 国产成人精品久久亚洲| 久久综合色老色| 亚洲va久久久噜噜噜久久天堂| 久久永久免费人妻精品下载| 亚洲欧美日韩中文久久 | 狠狠综合久久综合88亚洲| 漂亮人妻被中出中文字幕久久| 久久国产精品免费一区二区三区| 人人狠狠综合久久亚洲高清| 亚洲第一永久AV网站久久精品男人的天堂AV | 精品久久久久久无码专区不卡| 久久免费线看线看| 中文字幕久久久久人妻| 99久久免费国产特黄| 久久精品国产久精国产果冻传媒| 91精品国产9l久久久久| 99久久精品免费看国产一区二区三区| 久久国产精品久久精品国产|