bzoj2466,poj1222

都是简单的异或高斯消元

由于bzoj2466要求解得最小和，所以我们最后还要穷举自由元取最优解

 type node=record
        po,next:longint;
      end;

 var f:array[..] of longint;
     a:array[..,..] of longint;
     i,x,y,ans,n:longint;

 procedure swap(var a,b:longint);
   var c:longint;
   begin
     c:=a;
     a:=b;
     b:=c;
   end;

 procedure dfs(k,sum:longint);
   var j,p:longint;
   begin
     if sum>=ans then exit;
     if k= then
     begin
       if ans>sum then ans:=sum;
       exit;
     end;
     if a[k,k]= then
     begin
       p:=a[k,n+];
       for j:=k+ to n do
         if a[k,j]= then p:=p xor f[j];   //回代
       f[k]:=p;
       dfs(k-,sum+f[k]);
     end
     else begin
       f[k]:=;
       dfs(k-,sum);
       f[k]:=;
       dfs(k-,sum+);
       f[k]:=;
     end;
   end;

 procedure work;
   var i,j,k,p:longint;
   begin
     for i:= to n do
     begin
       p:=i;
       for j:=i to n do
         if a[j,i]<> then
         begin
           p:=j;
           break;
         end;

      if p<>i then
        for j:= to n+ do
          swap(a[i,j],a[p,j]);
      if a[i,i]= then continue;
      for k:=i+ to n do
        if a[k,i]<> then
        begin
          for j:=i to n+ do
            a[k,j]:=a[k,j] xor a[i,j];
        end;
     end;
     ans:=;
     dfs(n,);
   end;

 begin
   readln(n);
   while n<> do
   begin
     fillchar(a,sizeof(a),);
     fillchar(f,sizeof(f),);
     for i:= to n- do
     begin
       readln(x,y);
       a[x,y]:=;
       a[y,x]:=;
     end;
     for i:= to n do
     begin
       a[i,i]:=;
       a[i,n+]:=;
     end;
     work;
     writeln(ans);
     readln(n);
   end;
 end.

bzoj2466