BZOJ1770:[USACO]lights 燈(高斯消元,DFS)

Description

貝希和她的閨密們在她們的牛棚中玩遊戲。但是天不從人願，突然，牛棚的電源跳閘了，所有的燈都被關閉了。貝希是一個很膽小的女生，在伸手不見拇指的無盡的黑暗中，她感到驚恐，痛苦與絕望。她希望您能夠幫幫她，把所有的燈都給重新開起來！她才能繼續快樂地跟她的閨密們繼續玩遊戲！牛棚中一共有N（1 <= N <= 35）盞燈，編號為1到N。這些燈被置於一個非常複雜的網絡之中。有M（1 <= M <= 595）條很神奇的無向邊，每條邊連接兩盞燈。每盞燈上面都帶有一個開關。當按下某一盞燈的開關的時候，這盞燈本身，還有所有有邊連向這盞燈的燈的狀態都會被改變。狀態改變指的是：當一盞燈是開著的時候，這盞燈被關掉；當一盞燈是關著的時候，這盞燈被打開。問最少要按下多少個開關，才能把所有的燈都給重新打開。數據保證至少有一種按開關的方案，使得所有的燈都被重新打開。

Input

＊第一行：兩個空格隔開的整數：N和M。

＊第二到第M+1行：每一行有兩個由空格隔開的整數，表示兩盞燈被一條無向邊連接在一起。沒有一條邊會出現兩次。

Output

第一行：一個單獨的整數，表示要把所有的燈都打開時，最少需要按下的開關的數目。

Sample Input

5 6
1 2
1 3
4 2
3 4
2 5
5 3

輸入細節：

一共有五盞燈。燈1、燈4和燈5都連接著燈2和燈3。

Sample Output

輸出細節：

按下在燈1、燈4和燈5上面的開關。

Solution

不能再这么颓了不然真的就要凉透了
还有这题目为什么是繁体的

搞死小圆高斯消元好强啊QAQ这个题的做法也好神奇啊
很容易可以发现，开关的顺序不会影响最终结果

那么到底应该怎么搞呢？举个栗子
有3盏灯，其中1和2连接，1和3不连接，那么我们可以列出一个方程
(1*ans[1])^(1*ans[2])^(0*ans[3])=1
ans[]表示这个灯有没有按。灯1本身或者和灯1相连的灯系数为1，否则为0
为什么这么列呢？很容易发现，如果为1的项的异或和为1，那么最后结果肯定是1
感性理解一下，应该挺好懂的
以此类推我们可以列出n个这样的方程，然后我们就可以用高斯消元来把这个矩阵消了
消异或矩阵和消普通矩阵是一样的，只不过是把加减的操作换成异或

消了之后不要记着求答案，因为答案里面可能有自由变元。
怎么办呢？可以用爆搜替代普通高斯消元的答案回带，
若当前行的答案固定就算出来，否则的话就枚举0/1记入答案。继续搜下一层。

Code

 #include<iostream>

 #include<cstring>

 #include<cstdio>

 #define N (50)

 using namespace std;

 int n,m,minn=0x7fffffff,u,v;

 int f[N][N],ans[N],head[N],num_edge;

 void Gauss()

 {

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

     {

         int num=i;

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

             if (f[j][i]>f[i][i])

                 num=j;

         if (num!=i)

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

                 swap(f[i][j],f[num][j]);

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

             if (f[j][i])

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

                     f[j][k]^=f[i][k];

     }

 }

 void Dfs(int x,int now)

 {

     if (now>=minn) return;

     if (x==) {minn=now; return;}

     if (f[x][x])

     {

         int t=f[x][n+];

         for (int i=x+; i<=n; ++i) t^=f[x][i]*ans[i];

         ans[x]=t;

         Dfs(x-,now+(t==));

     }

     else

     {

         ans[x]=; Dfs(x-,now);

         ans[x]=; Dfs(x-,now+);

     }

 }

 int main()

 {

     scanf("%d%d",&n,&m);

     for (int i=; i<=m; ++i)

     {

         scanf("%d%d",&u,&v);

         f[u][v]=f[v][u]=;

     }

     for (int i=; i<=n; ++i) f[i][n+]=f[i][i]=;

     Gauss();

     Dfs(n,);

     printf("%d",minn);

 }