jzoj3505. 【NOIP2013模拟11.4A组】积木(brick) （正解）

好久没写题解了*2

声明

想看DP水法者请自觉点击右上角的红叉

引入

观察题目可以发现，最终结果一定是一个类似于山峰的图形。
因为每次更新区间时，首和尾都不会被增加( [A+1..B-1] )

所以很显然，我们可以把区间分成一段段来求解，求出的解之积就是最终答案。

推导

根据上面可以直到，每一段首尾已知的区间内无论怎样变化，在区间外的其它值都不会为之改变。

所以一段区间内的值只和首与尾有关，与其他值无关。

根据这个性质，我们可以把一个总的大问题，变成每一段区间的小问题求解。

预备

因为阶梯状的原因，设前一个位置积木高度为x，则当前位置为(x-1)或x或(x+1)。
理由同上。

所以可以把求解方案转化为在图形上求解两点间路径数。
如下：

比如我要求a~b区间的方案数，我们可以计算从a到b的路径总数。
同上，每次我们可以向→↗↘移动。

计算

枚举向右移动r步，可以求出向右上移动u步，向右下移动d步。

先假设a高度小于b。
比如我一路向前怼到黄圈处，这时如果我还继续往前走，就会发现无论如何也无法走到b点。

那么我们可以知道最少要向右上走多少步。
但是如果我们要向下走，我们必须要多走回等量的右上步，不然又会走不会b点。

所以说，多走的右下和右上相等。

公式

直接上公式。
C(r+u+d,r)*C(u+d,u)

解释一下。
第一个把所有的移动方向归纳成直走和斜走，总共要走(r+u+d)步，其中直走有r步，可以计算出方案数。
第二个是把斜走分成上斜和下斜，总共有(u+d)次斜走，其中上斜有u次，可以计算出方案数。

两个相称就是当前这一段的全部方案数。

组合数逆元不解释。

特殊

注意题目有一个性质，每条积木的高度必为非负。
所以有可能会出现这样一种情况(紫线为0

这TM就很尴尬了

答案

考虑：
不走下0的方案数=全部方案-一定走下0的方案

把b点根据-1行翻折

那么从a走到b’的方案就是从a走到b必经绿线的方案。
举个例子：

为什么这样能走到b呢？
因为b’是b的对称点，可以把它根据绿线翻折回来：

源码

计算部分

function calc(x,y,len:longint):int64;
var
        i,z:longint;
begin
        if x<y then
        begin
                z:=x;
                x:=y;
                y:=z;
        end;
        down:=0;
        up:=0;
        calc:=0;
        for i:=0 to len-(x-y) do
        if (len-(x-y)-i) mod 2=0 then
        begin
                down:=(len-(x-y)-i) div 2;
                up:=x-y+down;
                calc:=(calc+(c(i+up+down,i)*c(up+down,up))) mod md;
        end;
end;