【NOIP提高组模拟】树塔狂想曲

Description

相信大家都在长训班学过树塔问题，题目很简单求最大化一个三角形数塔从上往下走的路径和。走的规则是：（i，j）号点只能走向（i+1，j）或者（i+1，j+1）。如下图是一个数塔，映射到该数塔上行走的规则为：从左上角的点开始，向下走或向右下走直到最底层结束。

   1
   3 8
   2 5 0
   1 4 3 8
   1 4 2 5 0

   路径最大和是1+8+5+4+4 = 22，1+8+5+3+5 = 22或者1+8+0+8+5 = 22。

   小S觉得这个问题so easy。于是他提高了点难度，他每次ban掉一个点（即规定哪个点不能经过），然后询问你不走该点的最大路径和。
   当然他上一个询问被ban掉的点过一个询问会恢复（即每次他在原图的基础上ban掉一个点，而不是永久化的修改）。

Solution

因为修改不是永久的，所以不能直接修改原图。那我们设两个数组$f_{i,j}$和$g_{i,j}$，$f_{i,j}$表示由（1，1）走到（i，j）的最大值，$g_{i,j}$表示从底层走到（i，j）的最大值，当询问某一个点被禁止时，在同行的点中求最大值。

Code

var
    f,g,a:array[0..1000,0..1000] of longint;
    t:array[1..1000,1..4000] of longint;
    n,m,i,j,ans,x,y:longint;
function max(x,y:longint):longint;
begin
    if x>y then exit(x) else exit(y);
end;
procedure get(x,l,r,wz,z,y:longint);
var mid:longint;
begin
    if l>r then exit;
    if (l=z)and(r=y) then
    begin
        ans:=max(ans,t[x,wz]);
        exit;
    end;
    mid:=(l+r) div 2;
    if y<=mid then get(x,l,mid,wz*2,z,y)
    else if mid<z then get(x,mid+1,r,wz*2+1,z,y)
    else
    begin
        get(x,l,mid,wz*2,z,mid);
        get(x,mid+1,r,wz*2+1,mid+1,y);
    end;
end;
procedure make(x,l,r,wz:longint);
var mid:longint;
begin
    if l=r then
    begin
        t[x,wz]:=g[x,l]+f[x,l]-a[x,l];
        exit;
    end;
    mid:=(l+r)div 2;
    make(x,l,mid,wz*2);
    make(x,mid+1,r,wz*2+1);
    t[x,wz]:=max(t[x,wz*2],t[x,wz*2+1]);
end;
procedure dp;
var i,j:longint;
begin
    fillchar(g,sizeof(g),255);
    fillchar(f,sizeof(f),255);
    f[1,1]:=a[1,1];
    for i:=1 to n-1 do
        for j:=1 to i do
        begin
            f[i+1,j]:=max(f[i+1,j],f[i,j]+a[i+1,j]);
            f[i+1,j+1]:=max(f[i+1,j+1],f[i,j]+a[i+1,j+1]);
        end;
    for i:=1 to n do g[n,i]:=a[n,i];
    make(n,1,n,1);
    for i:=n-1 downto 1 do
    begin
        for j:=i downto 1 do
            g[i,j]:=max(g[i+1,j],g[i+1,j+1])+a[i,j];
        make(i,1,i,1);
    end;
end;
begin
    readln(n,m);
    for i:=1 to n do
    begin
        for j:=1 to i do read(a[i,j]);
        readln;
    end;
    dp;
    for i:=1 to m do
    begin
        readln(x,y);
        ans:=-1;
        if 1<=y-1 then get(x,1,x,1,1,y-1);
        if y+1<=x then get(x,1,x,1,y+1,x);
        writeln(ans);
    end;
end.