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

            A Za, A Za, Fighting...

            堅信:勤能補拙

            PKU 1190 生日蛋糕

            問題:
            http://acm.pku.edu.cn/JudgeOnline/problem?id=1190

            思路:
            DFS+減枝,好題

            代碼:
             1 /*
             2  * N = R[1]^2*H[1] + R[2]^2*H[2] +  + R[M]^2*H[M]
             3  * S = R[1]^2 + 2R[1]*H[1] + 2R[2]*H[2] +  + 2R[M]H[M]
             4  */
             5 #include<stdio.h>
             6 #include<stdlib.h>
             7 #include<string.h>
             8 #include<math.h>
             9 #define MAX_LEVEL 21
            10 #define INF 0x7FFFFFFF
            11 /* from top level to the i[th] level, the minimum total volumn and area */
            12 int min_volumn[MAX_LEVEL], min_area[MAX_LEVEL];
            13 int n, m;
            14 int rt;
            15 
            16 void
            17 init()
            18 {
            19     int i;
            20     rt = INF;
            21     min_volumn[0= min_area[0= 0;
            22     for(i=1; i<MAX_LEVEL; i++) {
            23         min_volumn[i] = min_volumn[i-1+ i*i*i;
            24         min_area[i] = min_area[i-1+ 2*i*i;
            25     }
            26 }
            27 
            28 /* from bottom(m[th] level) to the top */
            29 void
            30 dfs(int level, int last_r, int last_h, int cur_volumn, int cur_area)
            31 {
            32     int r, h, tmp, v, a;
            33     if(cur_volumn+min_volumn[level]>|| cur_area+min_area[level]>=rt)
            34         return;
            35     /* ADD this pruning according the volumn&area formula */
            36     if(2*(n-cur_volumn)/last_r+cur_area >= rt)
            37         return;
            38     if(level==0) {
            39         if(cur_volumn == n)
            40             rt = cur_area<rt ? cur_area : rt;
            41         return;
            42     }
            43     /* the minimal r in [level] would be level */
            44     for(r=last_r-1; r>=level; r--) {
            45         tmp = (int)((n-cur_volumn-min_volumn[level-1])/(double)(r*r));
            46         tmp = tmp>(last_h-1? (last_h-1) : tmp;
            47         for(h=tmp; h>=level; h--) {
            48             v = r*r*h;
            49             a = 2*r*h;
            50             if(level == m)
            51                 a += (r*r);
            52             dfs(level-1, r, h, cur_volumn+v, cur_area+a);
            53         }
            54     }
            55 }
            56 
            57 int
            58 main(int argc, char **argv)
            59 {
            60     int max_m_r, max_m_h;
            61     while(scanf("%d %d"&n, &m) != EOF) {
            62         init();
            63         max_m_r = (int)(sqrt((n-min_volumn[m-1])/(double)m)) + 1;
            64         max_m_h = (int)((n-min_volumn[m-1])/(double)(m*m)) + 1;
            65         dfs(m, max_m_r, max_m_h, 00);
            66         if(rt == INF)
            67             printf("0\n");
            68         else
            69             printf("%d\n", rt);
            70     }
            71 }

            posted on 2010-08-03 12:33 simplyzhao 閱讀(515) 評論(0)  編輯 收藏 引用 所屬分類: B_搜索

            導航

            <2010年8月>
            25262728293031
            1234567
            891011121314
            15161718192021
            22232425262728
            2930311234

            統計

            常用鏈接

            留言簿(1)

            隨筆分類

            隨筆檔案

            搜索

            最新評論

            閱讀排行榜

            評論排行榜

            久久一区二区三区99| 99久久人妻无码精品系列| 日本免费久久久久久久网站| 久久综合狠狠色综合伊人| 国产精品九九久久免费视频 | 成人国内精品久久久久一区| 韩国无遮挡三级久久| 国产午夜精品理论片久久| 中文字幕精品无码久久久久久3D日动漫 | 亚洲色婷婷综合久久| 欧美精品丝袜久久久中文字幕| 99久久精品国产一区二区| 中文字幕久久欲求不满| 久久久久无码精品国产app| 99国内精品久久久久久久| 性色欲网站人妻丰满中文久久不卡| 99久久精品国产一区二区三区 | 亚洲精品高清一二区久久| 久久亚洲精品中文字幕三区| 色婷婷综合久久久久中文一区二区 | 欧美精品一本久久男人的天堂 | 品成人欧美大片久久国产欧美...| 国产成人无码精品久久久性色| 久久久久久久国产免费看| 99久久精品免费看国产一区二区三区| 日本强好片久久久久久AAA| 亚洲一区精品伊人久久伊人| 国产精品VIDEOSSEX久久发布| 久久国产乱子精品免费女| 青青草原精品99久久精品66| 久久AV高潮AV无码AV| 精品久久久久久久国产潘金莲| 激情综合色综合久久综合| 国产精品欧美久久久久无广告| 99久久人人爽亚洲精品美女| 一级做a爰片久久毛片16| 99精品伊人久久久大香线蕉| 伊人色综合久久天天| 久久se精品一区精品二区国产| 国产亚洲精午夜久久久久久| 成人精品一区二区久久久|