【POJ3498】March of the Penguins（最大流，裂点）

题意：在靠近南极的某处，一些企鹅站在许多漂浮的冰块上。由于企鹅是群居动物，所以它们想要聚集到一起，在同一个冰块上。企鹅们不想把自己的身体弄湿，所以它们在冰块之间跳跃，但是它们的跳跃距离，有一个上限。 
随着气温的升高，冰块开始融化，并出现了裂痕。而企鹅跳跃的压力，使得冰块的破裂加速。幸运的是，企鹅对冰块十分有研究，它们能知道每块冰块最多能承受多少次跳跃。对冰块的损害只在跳起的时候产生，而落地时并不对其产生伤害。 
现在让你来帮助企鹅选择一个冰面使得它们可以聚集到一起。

(1≤N ≤100) (0 ≤D ≤100 000), 
接下来N行，xi, yi, ni and mi,分别表示冰块的X和Y坐标，该冰块上的企鹅数目，以及还能承受起跳的次数。 

输出所有可能的相聚冰块的编号，以0开始。如果不能相遇，输出-1。

思路：这道题是一年前上课最大流的例题，如果去年就拿了一等多好

        考虑强行限制起跳次数很难，尝试裂点

        将每块冰都裂成两个点(i,1),(i,2)

        (i,1)-->(i,2)连流量为m[i]的边

        对于原来在冰面上的企鹅，建立超级源S

         s-->(i,1)连流量为n[i]的边 

        对于平面距离小于D的两点i,j

        (i,2)-->(j,1)连流量为maxlongint的边

        枚举(j,1)作为汇点判断最大流量是否>=企鹅总数即可

1 var head,vet,next,gap,dis,len,c,fan,a,b,save:array[0..100000]of longint;
2     x,y:array[1..100000]of double;
3     num:array[1..100000,1..2]of longint;
4     n,m,qq,tot,i,j,ans,v,cas,s,source,src,st:longint;
5     d:double;
6
7 procedure add(a,b,c:longint);
8 begin
9  inc(tot);
10  next[tot]:=head[a];
11  vet[tot]:=b;
12  len[tot]:=c;
13  head[a]:=tot;
14 end;
15
16 function min(x,y:longint):longint;
17 begin
18  if x<y then exit(x);
19  exit(y);
20 end;
21
22 function dfs(u,aug:longint):longint;
23 var e,v,flow,t,val:longint;
24 begin
25  if u=src then exit(aug);
26  e:=head[u]; flow:=0; val:=s-1;
27  while e<>0 do
28  begin
29   v:=vet[e];
30   if len[e]>0 then
31   begin
32    if dis[u]=dis[v]+1 then
33    begin
34     t:=dfs(v,min(len[e],aug-flow));
35     len[e]:=len[e]-t;
36     len[fan[e]]:=len[fan[e]]+t;
37     flow:=flow+t;
38     if dis[source]>=s then exit(flow);
39     if aug=flow then break;
40    end;
41    val:=min(val,dis[v]);
42   end;
43   e:=next[e];
44  end;
45  if flow=0 then
46  begin
47   dec(gap[dis[u]]);
48   if gap[dis[u]]=0 then dis[source]:=s;
49   dis[u]:=val+1;
50   inc(gap[dis[u]]);
51  end;
52  exit(flow);
53 end;
54
55 function maxflow:longint;
56 var ans:longint;
57 begin
58  fillchar(dis,sizeof(dis),0);
59  fillchar(gap,sizeof(gap),0);
60  gap[0]:=s; ans:=0;
61  while dis[source]<s do ans:=ans+dfs(source,maxlongint);
62  exit(ans);
63 end;
64
65 begin
66  assign(input,'poj3498.in'); reset(input);
67  assign(output,'poj3498.out'); rewrite(output);
68  readln(cas);
69  for v:=1 to cas do
70  begin
71   fillchar(head,sizeof(head),0);
72   tot:=0; qq:=0; s:=0;
73   read(n,d);
74   for i:=1 to n do
75   begin
76    read(x[i],y[i],a[i],b[i]);
77    qq:=qq+a[i];
78   end;
79   for i:=1 to n do
80   begin
81    inc(s); num[i,1]:=s;
82    inc(s); num[i,2]:=s;
83   end;
84  inc(s); st:=s;
85  for i:=1 to n do
86   for j:=1 to n do
87    if (i<>j)and(sqrt(sqr(x[i]-x[j])+sqr(y[i]-y[j]))<=d) then
88    begin
89     fan[tot+1]:=tot+2;
90     fan[tot+2]:=tot+1;
91     add(num[i,2],num[j,1],maxlongint);
92     add(num[j,1],num[i,2],0);
93    end;
94  for i:=1 to n do
95  begin
96   fan[tot+1]:=tot+2;
97   fan[tot+2]:=tot+1;
98   add(num[i,1],num[i,2],b[i]);
99   add(num[i,2],num[i,1],0);
100  end;
101  for i:=1 to n do
102  begin
103   fan[tot+1]:=tot+2;
104   fan[tot+2]:=tot+1;
105   add(st,num[i,1],a[i]);
106   add(num[i,1],st,0);
107  end;
108  source:=st; ans:=0;
109  for i:=1 to tot do save[i]:=len[i];
110  for i:=1 to n do
111  begin
112   src:=num[i,1];
113   if maxflow>=qq then
114   begin
115    inc(ans); c[ans]:=i;
116   end;
117   for j:=1 to tot do len[j]:=save[j];
118  end;
119  if ans=0 then writeln(-1)
120   else
121   begin
122    for i:=1 to ans-1 do write(c[i]-1,' ');
123    write(c[ans]-1);
124    writeln;
125   end;
126  end;
127  close(input);
128  close(output);
129 end.