• <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 tutorial on computer science

              C++博客 :: 首頁 :: 新隨筆 :: 聯系 :: 聚合  :: 管理 ::
              21 隨筆 :: 0 文章 :: 17 評論 :: 0 Trackbacks

                  這個題糾結了好久。首先想到的是暴搜,必超時。
                然后注意到題目中,天平左邊的物體和右邊的物體交替取,可能會使更快得到最優解,所以就根據他們的位置是否>0把他們放到兩個數組里,這里有個優化就是按照到達-1.5和+1.5的力矩排序,這樣如果取了前面的物體天平失去了平衡,則后面的物體也不用取了,一定會失去平衡。這樣繼續暴搜,還是超時。我估計這個優化沒什么用,因為就算你把剩下的都取完,也不需要浪費多少時間嘛。。。平均意義上能減少一般的時間。但是A(20,20)。。這個有點大。
                 好吧,又繼續陷入僵局了。看別人的結題報告,發現主流版本是一堆亂糟糟的代碼,不懂在搞什么。。。后來偶然發現了一個版本,在暴搜的時候,寫的很奇怪,參見http://blog.csdn.net/wuli2496/article/details/7284641dfs(l,r),l,r表示前面l-1,r-1個已經取過了,這讓我煞是費解。(他的解法是把物體一個個的放到糾結的天平上),為什么不能先力矩放大的,后放力矩小的?然后我改了下他的dfs(),變成無參的,每次都從0開始找,交,超時了。作者沒有說明他為什么要這么做,或許是他不想說?或者他沒留意?而且這明明就不是個dfs了,如果可以證明,先放小的比先放大的好,那么這個dfs就是個DP了。和有限背包問題一樣的。好吧。然后我就想著證明為何先放力矩小的比先放力矩大的好,想了很久,發現,沒辦法。
                后來,偶然想到了一種方法,證明出
            如果在當前的情況中,拿掉左邊力矩最大的,天平向右歪,并且拿掉右邊力矩最大的,天平向左歪,那么無論以后怎么拿,都不可能成功。
                這樣反推一下,剛剛好是一個dp。r[i][j]表示放上去了左邊前i個,右邊前j個的天平狀況(這種狀況剛剛好依賴于0-i,0-j的質量)。它只能由r[i][j-1],r[i-1][j]推出。
                額,這樣就是一個完美的背包問題了,只要20*20個狀態,交了,A掉了,跑的足夠快,理論上應該是最快的解法了吧?肯定比搜索快,也要比狀態壓縮的DP快,因為狀態壓縮的DP需要2^20個狀態。。。可能是我的代碼寫的太亂了,所以才跑到15名。。
                好了,最重要的的問題就是證明上面一段中紅色字的結論。
                現在我們把天平分成3段:第一段是最左邊到-1.5,中間是-1.5到1.5,右邊是1.5-最右邊。
                左邊的物體的力矩我們記為:WL(左邊最大那個力矩) + Lrest(左邊剩余力矩) ,中間的物體的力矩我們記為 C+Mrest(C表示中間那一段的模板質量,Mrest表示中間物體質量),右邊的力矩我們記為
             WR(右邊最大那個力矩)+Rrest(右邊剩余力矩)。
            我們可以簡單地在左邊和右邊都剪掉Rrest。那么左邊,中間,右邊的力矩為WL+rest, C+Mrest,WR。如果拿掉左邊最大的右偏,我們可以得到式子 : rest + C + Mrest < WR 。那么很顯然,rest會越來越小,所以這個式子永遠成立,所以紅字得證。
            好了,這個證明想了很久,而且中間用到了一個結論:拿掉天平左右兩邊等重的物體,對天平沒影響。對于有兩個平衡點的這只是個臆測,額好吧,誰有更好更簡單的解法歡迎交流。LRJ出的題果真很神,不過也很讓人糾結。代碼寫的比較亂,莫看。
            #include <cstdio>
            #include <cstdlib>
            #include <cstring>

            class node
            {
              public:
                int pos,w;
            };

            class set
            {
             public:
               int ll,lr,rl,rr;
               int prev,ok;
            };

            set res[30][30];

            node lefts[110],rights[110];
            int lcount,rcount;
            int find;
            int L,W,C;

            int abs(int i)
            {
              if(i>0) return i;
              return -i;
            }

            void pans(int l,int r)
            {
              if(l == 0 && r== 0)
                return;
              if(res[l][r].prev == 1)
              {
                printf("%d %d\n",rights[r-1].pos/2,rights[r-1].w);
                pans(l,r-1);
              }
              else
              {
                printf("%d %d\n",lefts[l-1].pos/2,lefts[l-1].w);
                pans(l-1,r);
              }
            }
            int cmp1(const void* a,const void* b)
            {
              node i = *(node*)a;
              node j = *(node*)b;
              return (-i.pos-3)*i.w - (-j.pos-3)*j.w;
            }

            int cmp2(const void* a,const void* b)
            {
              node i = *(node*)a;
              node j = *(node*)b;
              return (i.pos-3)*i.w - (j.pos-3)*j.w;
            }


            int main()
            {
               //freopen("in.txt","r",stdin);
               
            //freopen("out.txt","w",stdout);
               int i,c=1,j; 
               int ll,lr,rl,rr;
               while(scanf("%d%d%d",&L,&W,&C)!=EOF && L>0)
               {
                 memset(res,0,sizeof(res));
                 L *= 2;
                 ll = rr =0; rl = lr = 3*W;
                 node tmp; 
                 lcount = rcount = 0;
                 for(i=0;i<C;i++)
                 {
                   scanf("%d%d",&tmp.pos,&tmp.w);
                   tmp.pos*=2;
                  if(tmp.pos<0)
                     lefts[lcount++] = tmp;
                   else
                     rights[rcount++] = tmp;
                 }
                 qsort(lefts,lcount,sizeof(node),cmp1);
                 qsort(rights,rcount,sizeof(node),cmp2);
                 
                 set tmpset = {ll,lr,rl,rr,0,1};
                 res[0][0] = tmpset;
                 for(i=0;i<=lcount;i++)
                  for(j=0;j<=rcount;j++)
                  {
                    if(i-1>=0 && res[i-1][j].ok == 1)
                        {
                              tmpset = res[i-1][j];
                              //放了第i個物體
                             if(lefts[i-1].pos < -3)
                               tmpset.ll += (-3-lefts[i-1].pos)*lefts[i-1].w; 
                             else
                               tmpset.lr += (lefts[i-1].pos+3)*lefts[i-1].w;
                             if(lefts[i-1].pos < 3)
                               tmpset.rl += (3 - lefts[i-1].pos)*lefts[i-1].w;
                             else
                               tmpset.rr += (lefts[i-1].pos - 3)*lefts[i-1].w;
                             if(tmpset.ll < tmpset.lr && tmpset.rl>tmpset.rr)
                             {
                                tmpset.ok = 1;
                                tmpset.prev = -1;
                                res[i][j] = tmpset;
                             }
                             else
                               {
                                 res[i][j].ok = 0;
                             }
                        }
                          
                          
                     if(res[i][j].ok) 
                          continue;
                     if(j-1>=0 && res[i][j-1].ok == 1)
                        {
                          tmpset = res[i][j-1];
                           if(rights[j-1].pos < -3)
                             tmpset.ll += (-3-rights[j-1].pos)*rights[j-1].w; 
                           else
                             tmpset.lr += (rights[j-1].pos+3)*rights[j-1].w;
                           if(rights[j-1].pos < 3)
                             tmpset.rl += (3 - rights[j-1].pos)*rights[j-1].w;
                           else
                             tmpset.rr += (rights[j-1].pos - 3)*rights[j-1].w;
                          if(tmpset.ll < tmpset.lr && tmpset.rl>tmpset.rr)
                             {
                                tmpset.ok = 1;
                                tmpset.prev = 1;
                                res[i][j] = tmpset;
                             }
                             else
                               {
                                 res[i][j].ok = 0;
                             }
                        }                      


                  }
                  printf("Case %d:\n",c++);
                   if(res[lcount][rcount].ok == 1)
                     pans(lcount,rcount);
                   else
                    printf("Impossible\n");
               }
               return 0;
            }

            posted on 2012-04-04 14:12 bigrabbit 閱讀(1800) 評論(2)  編輯 收藏 引用

            評論

            # re: uva 10123 - No Tipping 結題報告 2012-09-13 02:49 Jxy
            暴搜加上二進制記錄以及搜過的狀態可以到0.184s  回復  更多評論
              

            # re: uva 10123 - No Tipping 結題報告[未登錄] 2012-09-13 10:50 bigrabbit
            @Jxy
            我這個DP好像可以搞到0.02,求主頁  回復  更多評論
              

            一极黄色视频久久网站| 国产99久久九九精品无码| 久久久精品国产免大香伊| 精品蜜臀久久久久99网站| 久久久久亚洲精品无码网址| 少妇人妻综合久久中文字幕| 久久91综合国产91久久精品| 欧美无乱码久久久免费午夜一区二区三区中文字幕 | 久久久久久亚洲精品无码| 成人久久免费网站| 韩国三级中文字幕hd久久精品 | 久久精品无码一区二区无码| 久久这里只有精品视频99| 人妻精品久久无码专区精东影业| 久久se这里只有精品| 久久99精品久久久久久久久久| 久久久久99精品成人片| 久久精品国产影库免费看| 人妻少妇久久中文字幕| 超级碰碰碰碰97久久久久| 欧美一级久久久久久久大片| 欧美亚洲另类久久综合| 国产69精品久久久久777| 一本一本久久A久久综合精品| 青春久久| 亚洲国产天堂久久久久久| 国产精品美女久久久免费| 亚洲天堂久久精品| 久久精品国产亚洲网站| 97久久超碰国产精品2021| 久久99亚洲网美利坚合众国| 久久久精品国产sm调教网站| 久久久亚洲欧洲日产国码aⅴ| 亚洲中文字幕无码久久2017| 伊人久久五月天| 性做久久久久久久| 久久亚洲AV成人无码电影| 久久99精品久久久久久久久久| 97久久超碰国产精品2021| 99久久成人国产精品免费| 热久久这里只有精品|