【以前的空间】bzoj 1227 [SDOI2009]虔诚的墓主人

题解：hzw大神的博客说的很清楚嘛

http://hzwer.com/1941.html

朴素的做法就是每个点如果它不是墓地那么就可形成十字架的数量就是这个c（点左边的树的数量，k）*c（点右边的树的数量，k）*c（点上边的树的数量，k）*c（点下边的树的数量，k）。这样的话要枚举每个点，复杂度很明显爆表。

再仔细分析一下，其实完全可以不用枚举每个点，只需要枚举那些树的位置就可以，（下面的点指的都是树而非点？？蒟蒻语言不好）在同一列的两棵相邻的树之间的点的c（点左边的树的数量，k）*c（点右边的树的数量，k）是一样的，但是还是需要Σc（点上边的树的数量，k）*c（点下边的树的数量，k）（乘法分配率，自己yy一下），朴素枚举这些还是要复杂度爆表，于是问题就是如果在logn时间内搞定这个东西，于是想到如果维护“Σc（点上边的树的数量，k）*c（点下边的树的数量，k）”的信息（询问区间和，想到树状数组或者线段树），（有点类似于那个悬梁法）就是这个点的上一个点和这个点之间的“c（点上边的树的数量，k）*c（点下边的树的数量，k）”是有一点联系的，yy一下就是这个点在下一行中是作为下面的点上面的树的数量……于是修改量就是c（原来上边的树的数量+1，k）*c（原来下边的树（其实包括了这个点也就是这棵树）的数量-1，k）- c（原来上边的树的数量，k）*c（原来下边的树（其实包括了这个点也就是这棵树）的数量，k）。

然后只需要在一开始离散化一下树，再给树排个序，就差不多了，还用到了前缀和小小的搞一下。

type
  arr=record
    x,y:longint;
  end;

const
  mm=;

var
  i,j,kk,l,n,m,w,now,tot:longint;
  a:array[..]of arr;
  bit:array[..]of int64;
  b,d,totx,toty,hash:array[..]of longint;
  c:array[..,..]of int64;
  ans,ans1,ans2,ans3:int64;

procedure swap(var x,y:longint);
var
  i:longint;
begin
  i:=x;
  x:=y;
  y:=i
end;

function min(x,y:longint):longint;
begin
  if x>y then exit(y);
  exit(x)
end;

function lowbit(x:longint):longint;
begin
  exit(x and (-x))
end;

function find(x:longint):longint;
var
  l,r,mid:longint;
begin
  l:=;r:=tot;
  while l<=r do begin
    mid:=(l+r) >> ;
    if hash[mid]<x then l:=mid+
      else
        if hash[mid]>x then r:=mid-
          else exit(mid);
  end
end;

procedure work;
var
  i,j:longint;
begin
  c[][]:=;
  for i:= to w do begin
    c[i,]:=;
    for j:= to min(kk,i) do
      c[i,j]:=(c[i-,j]+c[i-,j-]) mod mm;
  end
end;

procedure qsort1(l,r:longint);
var
  i,j,k,mid:longint;
begin
  i:=l;
  j:=r;
  mid:=b[(l+r) >>];
  repeat
    while b[i]<mid do inc(i);
    while b[j]>mid do dec(j);
    if i<=j then begin
      swap(b[i],b[j]);
      inc(i);
      dec(j);
    end;
  until i>j;
  if l<j then qsort1(l,j);
  if i<r then qsort1(i,r)
end;

procedure qsort2(l,r:longint);
var
  i,j,k,mid1,mid2:longint;
begin
  i:=l;
  j:=r;
  mid1:=a[(l+r) >>].y;
  mid2:=a[(l+r) >>].x;
  repeat
    while (a[i].y<mid1) or (a[i].y=mid1) and (a[i].x<mid2) do inc(i);
    while (a[j].y>mid1) or (a[j].y=mid1) and (a[j].x>mid2) do dec(j);
    if i<=j then begin
      swap(a[i].x,a[j].x);
      swap(a[i].y,a[j].y);
      inc(i);
      dec(j);
    end;
  until i>j;
  if l<j then qsort2(l,j);
  if i<r then qsort2(i,r)
end;

procedure add(x:longint;y:int64);
begin
  while x<=tot do begin
    bit[x]:=(bit[x]+y) mod mm;
    inc(x,lowbit(x));
  end
end;

function ask(x:longint):int64;
var
  sum:int64;
begin
  sum:=;
  while x> do begin
    sum:=(sum+bit[x]) mod mm;
    dec(x,lowbit(x));
  end;
  exit(sum)
end;

begin
  readln(n,m);
  readln(w);
  fillchar(totx,sizeof(totx),);
  fillchar(toty,sizeof(toty),);
  fillchar(d,sizeof(d),);
  fillchar(bit,sizeof(bit),);
  fillchar(c,sizeof(c),);
  fillchar(hash,sizeof(hash),);
  for i:= to w do begin
    read(a[i].x,a[i].y);
    b[*i-]:=a[i].x;
    b[*i]:=a[i].y;
  end;
  qsort1(,*w);
  hash[]:=b[];
  tot:=;
  for i:= to *w do
    if b[i]<>b[i-] then begin
      inc(tot);
      hash[tot]:=b[i];
    end;
  read(kk);
  work;
  for i:= to w do begin
    inc(totx[find(a[i].y)]);
    inc(toty[find(a[i].x)]);
  end;
  qsort2(,w);
  now:=;
  ans:=;
  for i:= to w do begin
    j:=find(a[i].x);
    if (i>) and (a[i].y=a[i-].y) then begin
      inc(now);
      ans1:=ask(j-)-ask(find(a[i-].x)) mod mm;
      ans2:=c[now,kk]*c[totx[find(a[i].y)]-now,kk] mod mm;
      ans:=(ans+ans1*ans2) mod mm;
    end else now:=;
    inc(d[j]);
    ans3:=(c[d[j],kk]*c[toty[j]-d[j],kk]-c[d[j]-,kk]*c[toty[j]-d[j]+,kk])mod mm;
    if ans3<> then add(j,ans3);
  end;
  if ans< then ans:=ans+mm;
  writeln(ans);
end.