青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品

a tutorial on computer science

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

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

評論

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

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

青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            夜夜嗨av一区二区三区四区 | 国内成人精品视频| 在线播放日韩| 亚洲免费中文| 欧美激情亚洲| 久久久亚洲高清| 欧美午夜在线观看| 日韩亚洲欧美在线观看| 美女成人午夜| 久久亚洲精品一区二区| 国产精品区免费视频| 亚洲私人影院| 一本久道久久综合婷婷鲸鱼| 欧美国产日韩二区| 亚洲精品系列| 亚洲精品一区二| 欧美日韩大片| 中文国产成人精品久久一| 亚洲精品美女在线观看| 久久人人97超碰精品888| 韩国三级电影一区二区| 美女网站久久| 嫩草国产精品入口| 99视频+国产日韩欧美| 亚洲日本久久| 国产精品网站在线观看| 欧美亚洲一区在线| 久久国产直播| 亚洲国产高清高潮精品美女| 亚洲成人在线网站| 欧美视频一区二区三区| 亚洲欧美日韩久久精品| 香蕉成人伊视频在线观看 | 免费人成网站在线观看欧美高清| 久久精品国产一区二区三| 影音先锋另类| 亚洲精品色图| 国产乱理伦片在线观看夜一区| 久久久91精品| 久久精品国产亚洲aⅴ| 亚洲精品免费一区二区三区| 在线免费日韩片| 亚洲伦理中文字幕| 韩日精品视频一区| 欧美成人亚洲成人| 欧美日本精品一区二区三区| 午夜精品久久久久久久99黑人| 一区二区三区高清在线| 午夜在线视频观看日韩17c| 欧美精品久久久久久久久久| 久久一区二区三区超碰国产精品| 欧美日韩无遮挡| 亚洲日韩欧美视频| 一区二区在线观看视频| 久久精品视频在线看| 久久国产视频网站| 国产精品三级视频| 欧美精品电影在线| 国产精品久久国产精品99gif | 老司机免费视频久久| 欧美激情视频一区二区三区免费| 欧美在线视频观看免费网站| 免费看黄裸体一级大秀欧美| 欧美一区1区三区3区公司| 美脚丝袜一区二区三区在线观看 | 国产精品久久久久国产a级| 蜜臀av一级做a爰片久久| 国产精品另类一区| 亚洲精品黄色| 亚洲日本欧美在线| 久久天堂av综合合色| 亚洲欧洲99久久| 欧美女同视频| 亚洲国产天堂久久国产91| 国产曰批免费观看久久久| 一区二区三区欧美成人| 亚洲精品欧美专区| 美女尤物久久精品| 欧美成人第一页| 在线观看亚洲视频| 久久激情五月激情| 久久久伊人欧美| 国产主播一区二区三区四区| 亚洲欧美大片| 久久精品一二三区| 国产亚洲精品aa午夜观看| 亚洲视频免费在线观看| 日韩午夜在线电影| 欧美激情综合色综合啪啪| 欧美成人黑人xx视频免费观看| 精品999在线观看| 久久精品人人爽| 久热精品视频在线观看一区| 国户精品久久久久久久久久久不卡| 亚洲欧美国产一区二区三区| 国产欧美日韩精品专区| 久久裸体艺术| 久久精品一区二区三区不卡| 欧美午夜激情小视频| 91久久精品国产91久久| 亚洲国产精品欧美一二99| 久久五月婷婷丁香社区| 久久一区中文字幕| 亚洲缚视频在线观看| 久久综合99re88久久爱| 欧美高清在线观看| 99在线精品视频| 欧美日韩国产高清视频| 99热这里只有精品8| 亚洲在线成人| 国产欧美一二三区| 久久精品国产欧美激情| 欧美激情一区二区三区| 一二三区精品| 国产欧美日韩三区| 久久精品国产2020观看福利| 欧美国产高潮xxxx1819| 亚洲视频电影图片偷拍一区| 国产欧美丝祙| 欧美第十八页| 亚洲无线视频| 欧美+日本+国产+在线a∨观看| 亚洲美女黄网| 国产精品爽黄69| 久久全球大尺度高清视频| 亚洲欧洲在线免费| 久久成人免费视频| 亚洲国产黄色片| 国产精品任我爽爆在线播放 | 国产精品综合色区在线观看| 久久久久久久久久久久久9999| 亚洲国产精品一区二区尤物区 | 欧美激情按摩在线| 亚洲在线观看免费| 亚洲成人在线观看视频| 国产精品久久国产愉拍 | 狠狠色狠狠色综合日日小说| 欧美精品在线免费| 久久久久中文| 亚洲专区在线视频| 亚洲欧洲一区二区三区| 久久午夜国产精品| 亚洲一区影院| 亚洲精品一区二区三区樱花| 国产香蕉久久精品综合网| 欧美日韩国产一区二区三区| 久久久久一区| 欧美伊人久久| 亚洲在线观看| 99国产精品国产精品久久| 欧美成人一区二区三区| 性欧美1819性猛交| 一区二区av| 亚洲人永久免费| 在线播放日韩欧美| 黑人巨大精品欧美黑白配亚洲| 国产精品视频在线观看| 欧美视频免费在线观看| 欧美va天堂va视频va在线| 欧美在线视频观看| 亚洲欧美日韩另类| 亚洲影院免费| 中文亚洲视频在线| 免费成人在线视频网站| 亚洲欧美精品在线| av成人免费| 亚洲精品影院| 亚洲日本成人| 91久久久久久久久久久久久| 国产亚洲成av人片在线观看桃| 欧美性猛交99久久久久99按摩| 欧美sm重口味系列视频在线观看| 久久久久9999亚洲精品| 欧美中文字幕久久| 欧美专区第一页| 久久精品99久久香蕉国产色戒| 亚洲欧美日韩精品在线| 午夜精品亚洲一区二区三区嫩草| 亚洲一级电影| 亚洲欧美一区二区原创| 亚欧成人精品| 久久久91精品国产一区二区三区 | **性色生活片久久毛片| 极品日韩久久| 亚洲第一精品久久忘忧草社区| 精品91视频| 在线观看精品| 亚洲精品国精品久久99热| 91久久一区二区| 正在播放欧美一区| 亚洲一区高清| 久久本道综合色狠狠五月| 欧美中文字幕在线观看| 久久久久久高潮国产精品视| 久久亚洲春色中文字幕| 亚洲成人在线视频网站| 一区二区三区高清| 久久国产精品色婷婷| 免播放器亚洲一区| 欧美色道久久88综合亚洲精品|