在一块梯形田地上,一群蚯蚓在做收集食物游戏。蚯蚓们把梯形田地上的食物堆积整理如下:
a(1,1) a(1,2)…a(1,m)
a(2,1) a(2,2) a(2,3)…a(2,m) a(2,m+1)
a(3,1) a (3,2) a(3,3)…a(3,m+1) a(3,m+2)
……
a(n,1) a(n,2) a(n,3)… a(n,m+n-1)
它们把食物分成n行,第1行有m堆的食物,每堆的食物量分别是a(1,1),a(1,2),…,a(1,m);
第2行有m+1堆食物,每堆的食物量分别是a(2,1),a(2,2),…, a(2,m+1);以下依次有m+2堆、m+3堆、…m+n-1堆食物。
现在蚯蚓们选择了k条蚯蚓来测试它们的合作能力(1≤ k ≤m)。测试法如下:第1条蚯蚓从第1行选择一堆食物,然后往左下或右下爬,并收集1堆食物,例如从a(1,2)只能爬向a(2,2) 或a(2,3),而不能爬向其它地方。接下来再爬向下一行收集一堆食物,直到第n行收集一堆食物。第1条蚯蚓所收集到的食物量是它在每一行所收集的食物量之和;第2条蚯蚓也从第1行爬到第n行,每行收集一堆食物,爬的方法与第1条蚯蚓相类似,但不能碰到第1条蚯蚓所爬的轨迹;一般地,第i 条蚯蚓从第1行爬到第 n行,每行收集一堆食物,爬的方法与第1条蚯蚓类似,但不能碰到前 I-1 条蚯蚓所爬的轨迹。这k条蚯蚓应该如何合作,才能使它们所收集到的食物总量最多?收集到的食物总量可代表这k条蚯蚓的合作水平。
- Ø编程任务:
给定上述梯形m、n和k的值(1≤k≤m≤30;1≤n≤30)以及梯形中每堆食物的量(小于10的非整数),编程计算这k条蚯蚓所能收集到的食物的最多总量。
输入数据由文件名为INPUT1.*的文本文件提供,共有n+1行。每行的两个数据之间用一个空格隔开。
●第1行是n、m和k的值。
- 接下来的n行依次是梯形的每一行的食物量a(i,1),a(i,2),…,a(i,m+i-1),i=1,2,…,n。
将第一排与源点连接 最后一排与汇点连接
对于每个点拆点 中间连一条流量 1 费用 方格里的数 的边
注意连别的边的时候流量要置为 1 因为每个格子只能经过一次而不是k次
代码:
type
bi=record
next,cost,flow,go:longint;
end;
var
n,k,i,j,x,a,xb,t,w,ans,p,m,po:longint;
e:array[0..200000] of bi;
h,ll:array[0..60000] of longint;
b:array[0..60000] of boolean;
f,l,mi,r:array[0..60000] of longint;
q:array[0..10000000] of longint;
procedure adde(a,b,c,d:longint);
begin
inc(xb);
e[xb].cost:=d;
e[xb].flow:=c;
e[xb].go:=b;
if h[a]=0 then
begin
h[a]:=xb;
ll[a]:=xb
end else
begin
e[ll[a]].next:=xb;
ll[a]:=xb
end
end;
procedure cl(x:longint);
begin
if (x<=m)or(x=xb) then exit;
if odd(x-m) then inc(e[x+1].flow,mi[po+2])
else inc(e[x-1].flow,mi[po+2])
end;
begin
readln(n,m,k);
for i:=1 to m do adde(0,2*i-1,1,0);
po:=(m shl 1+n-1)*n;
for i:=1 to n do
for j:=1 to m+i-1 do
begin
read(x);
a:=(m+m+i-2)*(i-1) shr 1+j;
adde(2*a-1,2*a,1,x);
adde(2*a,2*a-1,0,-x);
if i<n then
begin
addE(2*a,2*a+(m+i-1)*2-1,1,0);
addE(2*a+(m+i-1)*2-1,2*a,0,0);
addE(2*a,2*a+(m+i-1)*2+1,1,0);
addE(2*a+(m+i-1)*2+1,2*a,0,0);
end else
begin
addE(2*a,po+1,1,0);
addE(po+1,2*a,0,0)
end
end;
adde(po+1,po+2,k,0);
repeat
fillchar(b,sizeof(b),false);
b[0]:=true;
t:=0;
w:=1;
for i:=0 to po+2 do
begin
f[i]:=-100000000;
mi[i]:=-f[i]
end;
for i:=0 to po+2 do l[i]:=-1;
fillchar(r,sizeof(r),0);
q[1]:=0;
f[0]:=0;
while t<w do
begin
inc(t);
b[q[t]]:=false;
p:=h[q[t]];
while p<>0 do
begin
if (f[e[p].go]<f[q[t]]+e[p].cost)and(e[p].flow>0) then
begin
f[e[p].go]:=f[q[t]]+e[p].cost;
l[e[p].go]:=q[t];
r[e[p].go]:=p;
if e[p].flow<mi[q[t]] then mi[e[p].go]:=e[p].flow
else mi[e[p].go]:=mi[q[t]];
if not b[e[p].go] then
begin
inc(w);
b[e[p].go]:=true;
q[w]:=e[p].go
end;
end;
p:=e[p].next
end
end;
j:=po+2;
if f[j]=-100000000 then break;
inc(ans,f[j]*mi[j]);
while l[j]<>-1 do
begin
dec(e[r[j]].flow,mi[po+2]);
cl(r[j]);
j:=l[j]
end
until false;
writeln(ans)
end.
type
bi=record
next,cost,flow,go:longint;
end;
var
n,k,i,j,x,a,xb,t,w,ans,p,m,po:longint;
e:array[0..200000] of bi;
h,ll:array[0..60000] of longint;
b:array[0..60000] of boolean;
f,l,mi,r:array[0..60000] of longint;
q:array[0..10000000] of longint;
procedure adde(a,b,c,d:longint);
begin
inc(xb);
e[xb].cost:=d;
e[xb].flow:=c;
e[xb].go:=b;
if h[a]=0 then
begin
h[a]:=xb;
ll[a]:=xb
end else
begin
e[ll[a]].next:=xb;
ll[a]:=xb
end
end;
procedure cl(x:longint);
begin
if (x<=m)or(x=xb) then exit;
if odd(x-m) then inc(e[x+1].flow,mi[po+2])
else inc(e[x-1].flow,mi[po+2])
end;
begin
readln(n,m,k);
for i:=1 to m do adde(0,2*i-1,1,0);
po:=(m shl 1+n-1)*n;
for i:=1 to n do
for j:=1 to m+i-1 do
begin
read(x);
a:=(m+m+i-2)*(i-1) shr 1+j;
adde(2*a-1,2*a,1,x);
adde(2*a,2*a-1,0,-x);
if i<n then
begin
addE(2*a,2*a+(m+i-1)*2-1,1,0);
addE(2*a+(m+i-1)*2-1,2*a,0,0);
addE(2*a,2*a+(m+i-1)*2+1,1,0);
addE(2*a+(m+i-1)*2+1,2*a,0,0);
end else
begin
addE(2*a,po+1,1,0);
addE(po+1,2*a,0,0)
end
end;
adde(po+1,po+2,k,0);
repeat
fillchar(b,sizeof(b),false);
b[0]:=true;
t:=0;
w:=1;
for i:=0 to po+2 do
begin
f[i]:=-100000000;
mi[i]:=-f[i]
end;
for i:=0 to po+2 do l[i]:=-1;
fillchar(r,sizeof(r),0);
q[1]:=0;
f[0]:=0;
while t<w do
begin
inc(t);
b[q[t]]:=false;
p:=h[q[t]];
while p<>0 do
begin
if (f[e[p].go]<f[q[t]]+e[p].cost)and(e[p].flow>0) then
begin
f[e[p].go]:=f[q[t]]+e[p].cost;
l[e[p].go]:=q[t];
r[e[p].go]:=p;
if e[p].flow<mi[q[t]] then mi[e[p].go]:=e[p].flow
else mi[e[p].go]:=mi[q[t]];
if not b[e[p].go] then
begin
inc(w);
b[e[p].go]:=true;
q[w]:=e[p].go
end;
end;
p:=e[p].next
end
end;
j:=po+2;
if f[j]=-100000000 then break;
inc(ans,f[j]*mi[j]);
while l[j]<>-1 do
begin
dec(e[r[j]].flow,mi[po+2]);
cl(r[j]);
j:=l[j]
end
until false;
writeln(ans)
end.