pku1815 割边

这题wa了N次，就是改某个图不改某个图的问题。。

题目大意：

一群人，1表示两个人能联系0表示两个人不能联系，一个起点s一个终点t，问最少删除多少个点能使s到t收不到信息，并且让你给出一组解，要求字典序最小。

简略题解：

显然就是最大流。

把每个点拆成两个点i，i'，从i向i'连一条容量为1的边，如果i和j能联系，那么从i'向j连一条容量为无穷的边。从s向s'连一条容量为无穷的边，从t向t'连一条容量为无穷的边。那么求这个网络的最大流就是第一问。

那么如何给出一组解呢，直接dfs是不行的，因为要求字典序了。我用了一个比较shax的办法就是按字典序枚举，然后删点，看最大流是否减小，若减小了那么就记录一下并且把这个点真正删掉。。这样会很慢，然后纠结的跑了500+ms。

其实还可以用dfs会很快，就是每删掉一个点，就从t到s开始dfs一条路径，这条路径包含i’--->i，如果到了说明有解。。我是后来才搜到这种方法的，所以没改。。

附程序：

program pku_1815;
const inf=10000000;
var a:array[0..201,0..201] of longint;
c,tc,av:array[0..401,0..401] of longint;
d,num:array[0..401] of longint;
n,st,ed,ans,now:longint;
{------------------------}
procedure init;
var i,j:longint;
begin
readln(n,st,ed);
for i:=1 to n do
for j:=1 to n do
read(a[i,j]);
if a[st,ed]=1 then
begin
writeln('NO ANSWER!');
halt;
end;
end;
{------------------------}
function min(x,y:longint):longint;
begin
if x<y then exit(x)
else exit(y);
end;
{------------------------}
function dfs(u,flow:longint):longint;
var i,tmp:longint;
begin
if u=ed+n then exit(flow);
dfs:=0;
for i:=1 to n*2 do
if (c[u,i]>0) and (d[u]=d[i]+1) then
begin
tmp:=dfs(i,min(flow-dfs,c[u,i]));
dec(c[u,i],tmp);
inc(c[i,u],tmp);
inc(dfs,tmp);
if dfs=flow then exit(dfs);
end;
if d[st]>=n*2 then exit;
dec(num[d[u]]);
if num[d[u]]=0 then d[st]:=n*2;
inc(d[u]);
inc(num[d[u]]);
end;
{------------------------}
procedure maxflow;
var flow:longint;
begin
ans:=0;
fillchar(d,sizeof(d),0);
fillchar(num,sizeof(num),0);
num[0]:=n*2;
while d[st]<n*2 do begin
flow:=dfs(st,maxlongint);
inc(ans,flow);
end;
end;
{------------------------}
function check(x:longint):boolean;
begin
if c[x,x+n]=1 then exit(false);
if (x=st) or (x=ed) then exit(false);
tc:=c;
c:=av;
c[x,x+n]:=0;
maxflow;
if now>ans then begin
now:=ans;
av[x,x+n]:=0;
exit(true);
end else begin
c:=tc;
exit(false);
end;
end;
{------------------------}
procedure main;
var i,j:longint;
begin
for i:=1 to n do
c[i,i+n]:=1;
c[st,st+n]:=inf;
c[ed,ed+n]:=inf;
for i:=1 to n do
for j:=1 to n do
if i<>j then
if a[i,j]=1 then c[i+n,j]:=inf;
av:=c;
tc:=c;
maxflow;
writeln(ans);
j:=ans;
now:=ans;
for i:=1 to n do
begin
if j=0 then break;
if check(i) then
begin
dec(j);
write(i,' ');
end;
end;
writeln;
end;
{------------------------}
begin
init;
main;
end.