bzoj1143 2718

最小可相交路径覆盖

先预处理可到达的点然后转化为最小不相交路径覆盖

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

var edge:array[0..101000] of node;
    p,cx,cy:array[0..200] of longint;
    v:array[0..200] of boolean;
    a,b:array[0..200,0..200] of boolean;
    j,len,n,m,x,y,ans,i:longint;

procedure add(x,y:longint);
  begin
    inc(len);
    edge[len].point:=y;
    edge[len].next:=p[x];
    p[x]:=len;
  end;

function find(x:longint):longint;
  var i,y:longint;
  begin
    i:=p[x];
    while i<>-1 do
    begin
      y:=edge[i].point;
      if not v[y] then
      begin
        v[y]:=true;
        if (cy[y]=-1) or (find(cy[y])>0) then
        begin
          cx[x]:=y;
          cy[y]:=x;
          exit(1);
        end;
      end;
      i:=edge[i].next;
    end;
    exit(0);
  end;

procedure dfs(x:longint);
  var i:longint;
  begin
    v[x]:=true;
    for i:=1 to n do
      if a[x,i] and not v[i] then dfs(i);
  end;

begin
  readln(n,m);
  fillchar(p,sizeof(p),255);
  for i:=1 to m do
  begin
    readln(x,y);
    a[x,y]:=true;
  end;
  for i:=1 to n do
  begin
    fillchar(v,sizeof(v),false);
    dfs(i);
    for j:=1 to n do
      if v[j] and (i<>j) then
        add(i,j);
  end;
  fillchar(cx,sizeof(cx),255);
  fillchar(cy,sizeof(cy),255);
  for i:=1 to n do
    if cx[i]=-1 then
    begin
      fillchar(v,sizeof(v),false);
      ans:=ans+find(i);
    end;
  writeln(n-ans);
end.