• <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)  編輯 收藏 引用 所屬分類: AlgorithmUSACO

            導航

            <2025年8月>
            272829303112
            3456789
            10111213141516
            17181920212223
            24252627282930
            31123456

            統計

            常用鏈接

            留言簿(2)

            隨筆分類

            隨筆檔案

            搜索

            積分與排名

            最新評論

            閱讀排行榜

            99久久无码一区人妻| 久久久久久a亚洲欧洲aⅴ| 久久99精品久久久久久噜噜 | 国产午夜电影久久| 久久www免费人成看国产片| 欧美亚洲国产精品久久| 久久99精品国产自在现线小黄鸭| 久久夜色精品国产欧美乱| 久久se精品一区二区影院 | 亚洲&#228;v永久无码精品天堂久久 | 欧美久久亚洲精品| 久久ww精品w免费人成| 久久久久久青草大香综合精品| 777午夜精品久久av蜜臀| 91久久精品国产成人久久| 日韩av无码久久精品免费| 久久久久国产一级毛片高清板| 久久亚洲私人国产精品vA | 香蕉久久夜色精品升级完成| 久久WWW免费人成—看片| 国产精品无码久久久久久| 久久午夜夜伦鲁鲁片免费无码影视| 久久久青草久久久青草| 久久久久久国产精品免费无码| 日韩十八禁一区二区久久| 国产综合成人久久大片91| 99久久99这里只有免费费精品| 午夜精品久久久久久影视riav | 偷窥少妇久久久久久久久| 久久成人18免费网站| 99久久亚洲综合精品成人| 久久香蕉国产线看观看乱码| 国产精品美女久久久久| 久久亚洲欧美国产精品 | 成人国内精品久久久久影院VR| 日本五月天婷久久网站| 中文字幕无码av激情不卡久久| 久久久久成人精品无码| 亚洲精品视频久久久| 狠狠色噜噜色狠狠狠综合久久| 久久久久久久久久久久久久 |