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

            pku 1222 EXTENDED LIGHTS OUT 高斯消元

            題意:
            一個(gè)5*6的開(kāi)關(guān)矩陣,撥動(dòng)每個(gè)開(kāi)關(guān),都會(huì)使得它本身以及前后左右四個(gè)開(kāi)關(guān)反轉(zhuǎn)。現(xiàn)在給出所有開(kāi)關(guān)的初始狀態(tài),問(wèn)使得所有開(kāi)關(guān)處于關(guān)狀態(tài),需要撥動(dòng)的開(kāi)關(guān)。

            解法:
            首先看到這題就應(yīng)該想到方程組。或者說(shuō),是一個(gè)模2方程組。
            但是,這個(gè)會(huì)給我們帶來(lái)求解的麻煩。求解方程組有經(jīng)典的高斯消元法,但是出現(xiàn)模運(yùn)算,確實(shí)讓人頭疼。這里,我們想到了異或運(yùn)算,聯(lián)系高斯消元的本質(zhì),是用一個(gè)方程來(lái)代換另外一個(gè)方程,十進(jìn)制模二運(yùn)算中的加、減運(yùn)算與二進(jìn)制中的異或運(yùn)算正好對(duì)應(yīng)!然后下面的事情就簡(jiǎn)單多了

            有一個(gè)小技巧,使用位運(yùn)算能夠大大簡(jiǎn)化編程復(fù)雜度。
            可以將原來(lái)矩陣的每一行壓縮成一個(gè)32位整數(shù),這樣每次消元的過(guò)程中選擇列主元的過(guò)程可以用排序輕松解決~,然后消去的過(guò)程和回代的過(guò)程也就非常好實(shí)現(xiàn)了。用這種方法,這題的代碼量可以控制在60行以?xún)?nèi)。
            忽然想起了老隊(duì)長(zhǎng)說(shuō)過(guò)的話(huà):100行以?xún)?nèi)的程序才是正解,做了這么多題,越來(lái)越發(fā)現(xiàn)這句話(huà)是多么的經(jīng)典。話(huà)說(shuō)現(xiàn)在郭老大在馬化騰那應(yīng)該算個(gè)紅人了~

            還有件很?chē)宓氖拢褂肧TL里的greater仿函數(shù)竟然要包括一個(gè)叫functional的頭文件,甚是詭異。。

            代碼:
             1 # include <cstdio>
             2 # include <algorithm>
             3 # include <functional>
             4 # include <cstring>
             5 using namespace std;
             6 inline void setbit(int &num,int bit)
             7 {
             8     num|=1<<(30-bit);
             9 }
            10 inline bool getbit(int &num,int bit)
            11 {
            12     if(num&(1<<(30-bit))) return true;
            13     else return false;
            14 }
            15 int main()
            16 {
            17     int test;
            18     scanf("%d",&test);
            19     for(int t=1;t<=test;t++)
            20     {
            21         int e[30];
            22         memset(e,0,sizeof(e));
            23         for(int i=0;i<30;i++)
            24         {
            25             int tmp;
            26             scanf("%d",&tmp);
            27             setbit(e[i],i);
            28             if(i%6!=0) setbit(e[i],i-1);
            29             if(i%6!=5) setbit(e[i],i+1);
            30             if(i/6!=0) setbit(e[i],i-6);
            31             if(i/6!=4) setbit(e[i],i+6);
            32             if(tmp) setbit(e[i],30);
            33         }
            34         for(int i=0;i<30;i++)//消元
            35         {
            36             sort(e+i,e+30,greater<int>());
            37             if(getbit(e[i],i))
            38                 for(int j=i+1;j<30;j++)
            39                     if(getbit(e[j],i))
            40                         e[j]^=e[i];
            41             
            42         }
            43         for(int i=29;i>=0;i--)//回代
            44             if(getbit(e[i],i))
            45               for(int j=i-1;j>=0;j--)
            46                   if(getbit(e[j],i))
            47                       e[j]^=e[i];
            48         printf("PUZZLE #%d\n",t);
            49         for(int i=0;i<30;i++)
            50         {
            51            if(e[i])
            52              printf("%d",getbit(e[i],30));
            53            else printf("0");
            54            if(i%6==5) printf("\n");
            55            else printf(" ");
            56         }
            57     }
            58     return 0;
            59 }



            posted on 2011-01-16 01:35 yzhw 閱讀(344) 評(píng)論(0)  編輯 收藏 引用 所屬分類(lèi): combination math

            <2010年10月>
            262728293012
            3456789
            10111213141516
            17181920212223
            24252627282930
            31123456

            導(dǎo)航

            統(tǒng)計(jì)

            公告

            統(tǒng)計(jì)系統(tǒng)

            留言簿(1)

            隨筆分類(lèi)(227)

            文章分類(lèi)(2)

            OJ

            最新隨筆

            搜索

            積分與排名

            最新評(píng)論

            閱讀排行榜

            91久久精品国产成人久久| 久久久久国产| 国产成人综合久久精品红| 国产精品久久亚洲不卡动漫| 久久亚洲精品成人无码网站 | 久久99精品国产一区二区三区| 亚洲欧美日韩精品久久亚洲区| 久久精品一区二区影院| 99久久国产热无码精品免费久久久久 | 亚洲伊人久久大香线蕉苏妲己| 国产成人无码久久久精品一| 久久综合狠狠综合久久综合88| 狠狠色综合网站久久久久久久高清| 久久国产AVJUST麻豆| 2021国产精品午夜久久| 亚洲国产精品无码久久SM| 久久综合九色综合网站| 国产精品久久久久久吹潮| 亚洲欧美精品伊人久久| 久久精品无码一区二区三区日韩| 久久久久综合国产欧美一区二区| 一本久久综合亚洲鲁鲁五月天| 久久精品国产亚洲av麻豆图片| 嫩草伊人久久精品少妇AV| 99久久99这里只有免费的精品| 91亚洲国产成人久久精品| 久久精品成人| 久久综合给合久久狠狠狠97色69| 97久久超碰成人精品网站| 久久国产乱子伦精品免费午夜| 色妞色综合久久夜夜| 国产精品九九九久久九九| 久久久久国色AV免费看图片| 久久人妻少妇嫩草AV蜜桃| 精品国产一区二区三区久久久狼 | 久久毛片免费看一区二区三区| 亚洲国产精品成人AV无码久久综合影院| 久久久久亚洲av综合波多野结衣| 国产精品久久网| 亚洲成色www久久网站夜月| 国产伊人久久|