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.