bzoj1856: [Scoi2010]字符串

传送门
巧妙的思想。
首先我们可以将0看成(1,-1)的向量，1看成（1,1)的向量
然后题目就变成了从(0,0)走到（n+m,n-m)且不经过直线y==-1的方案数。
我们首先去掉不经过直线y==-1的方案数，发现答案=C(n+m,n）
然后我们考虑计算GG的方案数。
首先，GG的方案数和y=-1一定有交。
于是我们找到最后一次和y=-1的交，然后将之前的部分以y=-1为对称轴翻转。
发现现在我们要统计从（0，-2）走到（n+m,n-m)的方案数。
显然，这里每一种方案都是GG的。
发现这一部分答案是C（n+m,n+1）
答案=C（n+m,n)-C(n+m,n+1)

var
  a,b,p:array [0..2000005] of longint;
  n,m:longint;
procedure prime(x:longint);
  var i,j:longint;
  begin
    for i:=2 to x do begin
      if (a[i]=0) then begin inc(p[0]); p[p[0]]:=i; a[i]:=p[0]; end;
      for j:=1 to p[0] do if (i*p[j]<=x) then a[i*p[j]]:=j else break;
    end;
  end;
function C(x,y:longint):int64;
  var i,j,t:longint; ans:int64;
  begin
    fillchar(b,sizeof(b),0);
    for i:=y+1 to x do begin
      t:=i;
      while t<>1 do begin inc(b[a[t]]); t:=t div p[a[t]]; end;
    end;
    for i:=2 to x-y do begin
      t:=i;
      while t<>1 do begin dec(b[a[t]]); t:=t div p[a[t]]; end;
    end;
    ans:=1;
    for i:=1 to p[0] do
      for j:=1 to b[i] do ans:=(ans*p[i]) mod 20100403;
    exit(ans);
  end;
begin
  read(n,m);
  prime(n+m);
  write((C(m+n,n)-C(m+n,n+1)+20100403) mod 20100403);
end.