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

            隨筆分類

            隨筆檔案

            搜索

            積分與排名

            最新評論

            閱讀排行榜

            中文字幕无码久久久| 91秦先生久久久久久久| 久久久久人妻一区精品色 | 久久不见久久见免费影院www日本| 久久久久亚洲av成人无码电影| 亚洲国产精品无码久久久不卡| 日本精品久久久久中文字幕8 | 狠狠色丁香久久综合五月| 一本色道久久综合狠狠躁篇| 国产91色综合久久免费| 亚洲乱码精品久久久久..| 久久精品国产第一区二区| 国内精品伊人久久久久av一坑 | 国产精品美女久久久m| 国产欧美久久久精品影院| 精品久久国产一区二区三区香蕉| 色妞色综合久久夜夜| 久久国内免费视频| 久久精品国产亚洲Aⅴ香蕉| 久久被窝电影亚洲爽爽爽| 精品久久久噜噜噜久久久| 欧洲精品久久久av无码电影| 久久人人爽人人爽人人片AV高清 | 五月丁香综合激情六月久久| 亚洲欧洲精品成人久久奇米网| 国产成人精品久久| 久久久中文字幕| 成人亚洲欧美久久久久| 日本久久久久久中文字幕| 香港aa三级久久三级| 亚洲国产精品久久久久婷婷软件| 精品一区二区久久| 大蕉久久伊人中文字幕| 欧美性猛交xxxx免费看久久久| 久久伊人影视| 伊人久久大香线蕉av一区| 久久99精品国产麻豆| 香蕉久久一区二区不卡无毒影院| 观看 国产综合久久久久鬼色 欧美 亚洲 一区二区 | 亚洲综合久久久| 欧美熟妇另类久久久久久不卡|