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

pku_1753_Flip Game--Gauss+Dfs

   
 

//Name: pku_1753_Flip Game

//Author: longxiaozhi

http://hi.baidu.com/xiehuixb/blog/item/9ce25f10ee8a2e77ca80c4d1.

//Root: 和前面的pku1681一個意思,但是,測試數據加強了!這個題目需要考慮自由元!!也就是我以前的gauss的所在,

//這個問題一直困擾了我兩天天。以前的題目的測試數據好像并沒有考慮這一點,所以可以水果去,但是這個題目就不行啦。

//但是,還是可以用gauss的算法解決的。如果不存在自由變元,那么說明原方程組只有一組解,也就是我們要求的;

//如果存在自由變元,則說明我們求道的方程的解只是其中的一種情況,但不一定是最優的。

//一次我們還要考慮自由變元因為我們求出來的知識滿足條件的一種解,并不是這里要求的最優解。

//我們要考慮自由變元的取值:

//每個自由元我們都枚舉她的值(或)讓后進行一次深搜就可以把自由元的情況加進去,正如xiehui博客里寫的一樣。

#include <iostream>

#include <string.h>

using namespace std;

#define eps 1e-10

#define fabs(x) ((x)>0?(x):-(x))

#define zero(x) (fabs(x) < eps)

const int maxn = 16;

int a[maxn][maxn+1], b[maxn][maxn+1];

int n, m, ans;

const int dir[4][2]={{0,1},{1,0},{0,-1},{-1,0}};

inline int isBound(int a, int b) {

     return a < 0 || a >= m || b < 0 || b >= m;

}

void pmat(int a[][maxn+1]) {

     for (int i = 0, j; i < n; i++) {

         for (j = 0; j < n; j++)

              printf("%d ", a[i][j]);

         printf("%d\n", a[i][j]);

     }

}

void pans(int a[][maxn+1]) {

     for (int i = 0; i < n; i++)

         printf((i+1)%m ? "%d ":"%d\n", a[i][n]);

}

int gans(int a[][maxn+1]) {

     int i, j, ret = a[n-1][n];

     for (i = n-2; i >= 0; i--) {

         for (j = i+1; j < n; j++)

              a[i][n] ^= a[i][j] && a[j][n];

         ret += b[i][n];

     }

     return ret;

}

void dfs(int p, int k) {

     if (p == k) {

         memcpy(b, a, sizeof(b));

         int ret = gans(b);

         if (ret < ans) ans = ret;

         return;

     }

     a[p][n] = 1; dfs(p-1, k);

     a[p][n] = 0; dfs(p-1, k);

}

 

int getRow(int p, int q, int &row) {

        for (int i = p; i < n; i++)

               if (!zero(a[i][q])) return a[row=i][q];

        return row = 0;

}

void swapRow(int i, int row, int q) {

        for (int k = q; k <= n; k++)

               swap(a[i][k], a[row][k]);

}

int gauss() {

        int i = 0, j = -1, k, p, q, ret, row;

        while(i < n && ++j < n) {

               ret = getRow(i, j, row);

               if (zero(ret)) continue;

               if (row != i) swapRow(i, row, j);

               for (p = i + 1; p < n; p++) if (a[p][j])

                       for (q = j; q <= n; q++)

                               a[p][q] ^= a[i][q];

               ++i;

        }pmat(a);

        for (k = i; k < n; k++) if (a[k][n]) return -1;

        dfs(n-1, i-1);
        return ans;

}

int main() {

     //freopen("in.in", "r", stdin);

     int i, j, k, x, y, p, q; char s[5][5];

     n = 16; m = 4;

     for (i = 0; i < m; i++) {

         scanf("%s", s[i]);

         for (j = 0; j < m; j++) {

              a[i*m+j][n] = s[i][j] == 'b';

              a[i*m+j][i*m+j] = 1;

              for (k = 0; k < 4; k++) {

                   x = i + dir[k][0];

                   y = j + dir[k][1];

                   if (isBound(x, y)) continue;

                   a[i*m+j][x*m+y] = 1;

              }

         }

     }

     ans = 1 << 20;

     p = gauss();   // printf("p = %d\n", p);

     memset(a, 0, sizeof(a));

     for (i = 0; i < m; i++) {

         for (j = 0; j < m; j++) {

              a[i*m+j][n] = s[i][j] == 'w';

              a[i*m+j][i*m+j] = 1;

              for (k = 0; k < 4; k++) {

                   x = i + dir[k][0];

                   y = j + dir[k][1];

                   if (isBound(x, y)) continue;

                   a[i*m+j][x*m+y] = 1;

              }

         }

     }

     ans = 1 << 20;

     q = gauss();   // printf("q = %d\n", q);

     if (p == -1 && q == -1) puts("Impossible");

     else if (p == -1) printf("%d\n", q);ac

     else if (q == -1) printf("%d\n", p);

     else printf("%d\n", p <= q ? p:q);

     return 0;

}

 

posted on 2010-05-24 23:55 小志 閱讀(1479) 評論(1)  編輯 收藏 引用

評論

# re: pku_1753_Flip Game--Gauss+Dfs 2010-05-26 14:00 小志

聲明:第一次的Gauss模板存在Bug,現已修復  回復  更多評論   


只有注冊用戶登錄后才能發表評論。
網站導航: 博客園   IT新聞   BlogJava   博問   Chat2DB   管理


<2025年11月>
2627282930311
2345678
9101112131415
16171819202122
23242526272829
30123456

導航

統計

常用鏈接

留言簿

隨筆檔案(8)

文章檔案(1)

相冊

收藏夾

ACM --- Online Judge

小志

最新隨筆

積分與排名

最新隨筆

最新評論

閱讀排行榜

評論排行榜

青青草原综合久久大伊人导航_色综合久久天天综合_日日噜噜夜夜狠狠久久丁香五月_热久久这里只有精品
  • <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>
            99re成人精品视频| 日韩视频免费观看| 极品中文字幕一区| 亚洲一区二区日本| 欧美激情国产高清| 亚洲伊人伊色伊影伊综合网 | 一本久久综合| 欧美国产在线电影| 久久久久99| 国产午夜一区二区三区| 欧美中文在线字幕| 亚洲资源av| 国产精品入口福利| 欧美一区二区三区日韩| 亚洲免费影视第一页| 国产精品久久久久永久免费观看| 亚洲一区国产精品| 亚洲一区二区三区免费观看| 国产精品久久久久久av下载红粉| 亚洲午夜激情网页| 亚洲午夜极品| 国产亚洲精品久久飘花 | 亚洲精品国产精品国自产在线| 午夜影院日韩| 激情视频一区二区| 欧美激情a∨在线视频播放| 老**午夜毛片一区二区三区| 在线精品视频一区二区三四| 免费在线成人av| 欧美成人一区二区在线 | 欧美激情四色| 欧美成人性网| 亚洲一区二区三区中文字幕在线 | 亚洲国产精品久久91精品| 欧美精品少妇一区二区三区| 亚洲婷婷免费| 欧美专区第一页| 亚洲国产婷婷| av成人免费在线观看| 国产乱码精品一区二区三区忘忧草| 欧美在线一区二区三区| 浪潮色综合久久天堂| av成人毛片| 午夜精品福利视频| 亚洲激情在线观看视频免费| 9人人澡人人爽人人精品| 国产精品日韩二区| 免费成人在线观看视频| 欧美日韩三区四区| 久久久久久高潮国产精品视| 欧美+日本+国产+在线a∨观看| 在线一区二区视频| 欧美在线免费视屏| 一卡二卡3卡四卡高清精品视频| 亚洲一级片在线观看| 激情欧美一区二区| 国产精品99久久久久久宅男| 国际精品欧美精品| 一区二区精品在线观看| 老牛国产精品一区的观看方式| 亚洲精品女av网站| 欧美一级二区| 欧美激情日韩| 亚洲精品综合| 欧美日韩综合一区| 99精品视频一区| 欧美一区二区私人影院日本| 国产免费观看久久黄| 久久久久久夜| 日韩午夜在线电影| 亚洲欧美中文日韩v在线观看| 欧美成人中文字幕在线| 午夜在线一区二区| 欧美日韩福利视频| 猛干欧美女孩| 国产乱码精品| 亚洲天堂第二页| 99re66热这里只有精品3直播| 久久精品欧洲| 久久成人精品视频| 国产精品一区二区在线| 一本一本久久| 夜夜嗨av一区二区三区| 久久中文在线| 你懂的一区二区| 国语自产在线不卡| 亚洲在线电影| 午夜精品在线观看| 国产精品免费看| 亚洲网站在线看| 亚洲一区二区在线| 欧美日韩一区二区在线视频| 亚洲国产精品视频一区| 亚洲黄色在线看| 欧美 日韩 国产一区二区在线视频| 久久综合久久美利坚合众国| 国产午夜精品视频| 欧美一级午夜免费电影| 欧美一区亚洲| 国产亚洲午夜| 久久精品伊人| 欧美va天堂| 亚洲国产视频一区| 欧美精品久久久久久| 亚洲精品久久久久久下一站 | 亚洲摸下面视频| 午夜欧美大片免费观看 | 久久久亚洲精品一区二区三区 | 亚洲电影免费观看高清| 久热精品视频在线免费观看| 狠狠干成人综合网| 久久久久久久网| 欧美成人嫩草网站| 欧美久色视频| 9久re热视频在线精品| 亚洲综合成人婷婷小说| 国产精品丝袜久久久久久app | 欧美视频在线观看一区| 在线一区二区三区四区五区| 午夜在线不卡| 在线精品国产欧美| 欧美激情在线有限公司| 亚洲香蕉成视频在线观看 | 欧美国产精品一区| 在线视频中文亚洲| 国产精品外国| 久久影院午夜论| 99国产一区二区三精品乱码| 午夜精品国产更新| 1024精品一区二区三区| 欧美三级日本三级少妇99| 欧美一级片在线播放| 亚洲日韩中文字幕在线播放| 亚洲欧美日韩天堂| 亚洲电影免费观看高清完整版在线| 欧美日韩国产123区| 欧美一区视频| 亚洲精品黄色| 欧美18av| 欧美一区激情| 亚洲网址在线| 亚洲日本中文| 一区二区三区中文在线观看 | 久久超碰97人人做人人爱| 日韩午夜av| 激情五月综合色婷婷一区二区| 欧美精品videossex性护士| 性亚洲最疯狂xxxx高清| 亚洲美女av网站| 欧美福利一区二区三区| 欧美在线在线| 亚洲影院色无极综合| 91久久中文字幕| 狠狠做深爱婷婷久久综合一区| 欧美日韩亚洲一区三区| 久久亚洲电影| 久久成人人人人精品欧| 亚洲专区欧美专区| 一区二区欧美视频| 最新亚洲电影| 欧美黄色aaaa| 美女诱惑黄网站一区| 久久久91精品国产一区二区三区| 亚洲图片欧洲图片日韩av| 日韩亚洲视频| 亚洲精品乱码久久久久久蜜桃麻豆 | 久久中文字幕一区| 欧美在线视频导航| 亚洲欧美中文另类| 亚洲欧美成人一区二区在线电影| 亚洲精品网站在线播放gif| 亚洲国产精品v| 伊人天天综合| 亚洲高清免费在线| 在线看片成人| 亚洲三级色网| 日韩一区二区久久| 99精品免费| 亚洲无线视频| 亚洲一区免费视频| 久久久之久亚州精品露出| 久久久久久久久久久一区| 亚洲欧美久久| 亚洲欧美日本另类| 亚洲欧美bt| 久久精品成人一区二区三区蜜臀 | 欧美国产欧美亚洲国产日韩mv天天看完整| 久久久久久久网站| 欧美成人一区在线| 欧美日韩一区二区三区免费| 欧美午夜精品久久久久久久| 欧美视频在线观看| 国产伦精品一区| 狠狠久久婷婷| 亚洲美女精品一区| 亚洲一区二区网站| 久久精品人人爽| 欧美成人一区二区三区片免费| 亚洲电影网站| 一区二区国产精品|