SSL 1558——科技庄园

Description

由于PFT不是机器人，所以他的体力并不是无限的，他不想摘很多的桃以至体力为0，而白白把桃给Life。同时PFT每次只能摘一棵桃树，摘桃花费时间体力，每棵桃树都可以摘K次（对于同一棵桃每次摘的桃数相同）。每次摘完后都要返回出发点（PFT一次拿不了很多）即Life的所在地（0，0）{试验田左上角的桃坐标是（1，1）}。
PFT每秒只能移动一个单位，每移动一个单位耗费体力1(摘取不花费时间和体力,但只限上下左右移动)。

Input

第一行：四个数为N，M，TI，A 分别表示试验田的长和宽，Life给PFT的时间，和PFT的体力。
下面一个N行M列的矩阵桃田。表示每次每棵桃树上能摘的桃数。
接下来N行M列的矩阵，表示每棵桃最多可以采摘的次数K。

Output

一个数：PFT可以获得的最大的桃个数。

Sample Input

4 4 13 20
10 0 0 0
0 0 10 0
0 0 10 0
0 0 0 0
1 0 0 0
0 0 2 0
0 0 4 0
0 0 0 0
【样例解释】
可以摘到1次（1，1）和1次（2，3），体力和时间不满足再摘桃了。

Sample Output

10
Hint

对于M N TI A
10<=30%<=50
10<=100%<=100
对于K
10<=100%<=100
保证结果在longint范围内

---

因为每摘一次都要回到起点，所以我们在预处理时，就将摘这棵桃树一次耗费的时间求出。
枚举桃树、时间和这棵桃树能摘的次数
状态转移方程：f[j]:=max(f[j],f[j-w[i]*k]+c[i]*k);

---

代码如下：

var
  a,b:array[1..1000,1..1000] of longint;
  w,c,u,f:array[0..10000]of longint;
  n,m,i,j,k,v,t,p:longint;
function max(x,y:longint):longint;
  begin
    if x>y then exit(x);
    exit(y);
end;
begin
  readln(n,m,t,p);
  k:=0;
  for i:=1 to n do
    for j:=1 to m do
      read(a[i,j]);
  for i:=1 to n do
    for j:=1 to m do
      begin
        read(b[i,j]);
        if a[i,j]>0 then
          begin
            inc(k);
            c[k]:=a[i,j];
            u[k]:=b[i,j];
            w[k]:=(i+j)*2;
          end;
      end;
  v:=k;
  for i:=1 to v do
    for j:=t downto 0 do
      for k:=0 to u[i] do
        if (j>=w[i]*k) then f[j]:=max(f[j],f[j-w[i]*k]+c[i]*k);
  if t>=p then t:=p-1;
  writeln(f[t]);
end.