bzoj1305

让我们继续来练网络流;

很明显是一个最大流的问题;

二分枚举最多次数m,然后最大流判定;

具体就是男生女生都拆成两个点i1,i2,之间连一条流量为k的边(男生i1-->i2,女生i2-->i1);

i2连不喜欢的人,i1连喜欢的人

最后,男生i1连源点流量为m,女生i1连汇点流量为m

最后判断最大流是否等于n*m即可

但做着做着,我发现好像好像二分+最大流不是很优,因为进行了很多重复操作

但我也没管,先A了再说;

后来看status发现很多人代码很短,用时0ms(我的最大流180ms)

肯定有更简单的方法:

贪心!……其实是错的……

 code(using maxflow):
type node=record
       next,point,flow:longint;
     end; var edge:array[..] of node;
    a:array[..,..] of boolean;
    cur,pre,p,numh,h:array[..] of longint;
    j,m,n,k,i,len,t,ans,l,r:longint;
    c:string; function min(a,b:longint):longint;
  begin
    if a>b then exit(b) else exit(a);
  end; procedure add(x,y,f:longint);
  begin
    inc(len);
    edge[len].flow:=f;
    edge[len].point:=y;
    edge[len].next:=p[x];
    p[x]:=len;
  end; function sap(m:longint):boolean;
  var s,u,tmp,i,j,q:longint;
  begin
    len:=-;
    fillchar(p,sizeof(p),);
    for i:= to n do
    begin
      add(,i,m);
      add(i,,m);
      add(i+*n,t,m);
      add(t,i+*n,);
      add(i,i+n,k);
      add(i+n,i,);
      add(i+*n,i+*n,k);
      add(i+*n,i+*n,);
    end;
    for i:= to n do
      for j:= to n do
        if a[i,j] then
        begin
          add(i,j+*n,);
          add(j+*n,i,);
        end
        else begin
          add(i+n,j+*n,);
          add(j+*n,i+n,);
        end;
    fillchar(numh,sizeof(numh),);
    fillchar(h,sizeof(h),);
    numh[]:=t+;
    u:=;
    s:=;
    while h[]<t+ do
    begin
      if u=t then
      begin
        i:=;
        while i<>t do
        begin
          j:=cur[i];
          dec(edge[j].flow);
          inc(edge[j xor ].flow);
          i:=edge[j].point;
        end;
        u:=;
        inc(s);
        if s=n*m then exit(true);
      end;
      q:=-;
      i:=p[u];
      while i<>- do
      begin
        j:=edge[i].point;
        if (edge[i].flow>) and (h[u]=h[j]+) then
        begin
          q:=i;
          break;
        end;
        i:=edge[i].next;
      end;
      if q<>- then
      begin
        cur[u]:=q;
        pre[j]:=u;
        u:=j;
      end
      else begin
        dec(numh[h[u]]);
        if numh[h[u]]= then exit(false);
        tmp:=t+;
        i:=p[u];
        while i<>- do
        begin
          j:=edge[i].point;
          if edge[i].flow> then tmp:=min(tmp,h[j]);
          i:=edge[i].next;
        end;
        h[u]:=tmp+;
        inc(numh[h[u]]);
        if u<> then u:=pre[u];
      end;
    end;
    exit(false);
  end; begin
  readln(n,k);
  for i:= to n do
  begin
    readln(c);
    for j:= to n do
    begin
      if c[j]='Y' then a[i,j]:=true
      else a[i,j]:=false;
    end;
  end;
  t:=n*+;
  l:=;
  r:=n;
  ans:=;
  while l<=r do
  begin
    m:=(l+r) shr ;
    if sap(m) then
    begin
      ans:=m;
      l:=m+;
    end
    else r:=m-;
  end;
  writeln(ans);
end.
上一篇:vs2015中升级Nuget后,找不到Umbraco安装包问题


下一篇:linux chgrp命令【转载】