bzoj2466,poj1222

都是简单的异或高斯消元

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

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

var f:array[0..110] of longint;
    a:array[0..110,0..110] 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=0 then
    begin
      if ans>sum then ans:=sum;
      exit;
    end;
    if a[k,k]=1 then
    begin
      p:=a[k,n+1];
      for j:=k+1 to n do
        if a[k,j]=1 then p:=p xor f[j];   //回代
      f[k]:=p;
      dfs(k-1,sum+f[k]);
    end
    else begin
      f[k]:=0;
      dfs(k-1,sum);
      f[k]:=1;
      dfs(k-1,sum+1);
      f[k]:=0;
    end;
  end;

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

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

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