• <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_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 小志 閱讀(1444) 評論(1)  編輯 收藏 引用

            評論

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

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

            <2025年5月>
            27282930123
            45678910
            11121314151617
            18192021222324
            25262728293031
            1234567

            導航

            統計

            常用鏈接

            留言簿

            隨筆檔案(8)

            文章檔案(1)

            相冊

            收藏夾

            ACM --- Online Judge

            小志

            最新隨筆

            積分與排名

            最新隨筆

            最新評論

            閱讀排行榜

            評論排行榜

            日本高清无卡码一区二区久久| 亚洲国产精品无码久久| 久久久久无码中| 久久国产福利免费| 国产精品久久久天天影视香蕉 | 日本人妻丰满熟妇久久久久久| 伊人久久精品无码二区麻豆| 久久精品国产亚洲AV麻豆网站| 国产成人久久精品麻豆一区| 久久久午夜精品福利内容| 99久久国产综合精品麻豆| 精品久久久无码中文字幕天天| 久久久老熟女一区二区三区| 久久久久女教师免费一区| 亚洲欧洲日产国码无码久久99| 国产成人久久精品二区三区| 久久人爽人人爽人人片AV| 青青草原综合久久大伊人导航| 国产精品99久久99久久久| 久久丫忘忧草产品| 久久国产精品国语对白| 久久婷婷国产麻豆91天堂| 久久精品国产亚洲AV麻豆网站| 精品国产乱码久久久久软件| 国产精品久久久久乳精品爆| 2020久久精品国产免费| 无码精品久久久久久人妻中字| 久久综合色区| 久久精品国产99久久丝袜| 狠狠色丁香婷婷综合久久来来去| 国内精品伊人久久久久| 国产成人久久AV免费| 久久精品亚洲中文字幕无码麻豆| 99久久夜色精品国产网站| 中文字幕久久精品无码| 狠狠色噜噜色狠狠狠综合久久| 亚洲国产成人久久综合碰| 久久婷婷色综合一区二区| 性做久久久久久久久浪潮| 国内精品综合久久久40p| 久久人人爽人人爽人人爽|