P2154 虔诚的墓主人(树状数组+排序降维)

 

  1.进行二维离散化。

  2.对坐标排序，先按x从小到大排序，如果x相同则按y从小到大排序。

  3.求出各行各列有几棵树，预处理组合数。

  4.从1开始到离散后x坐标最大值结束，遍历x坐标。

  5.遍历每个x坐标时y坐标从小到大，需要维护当前已经遍历到的棵数(down)，用该列的总数(col)减去已经遍历的(down)就是上面一部分的。

  6.对遍历到的坐标(x,y),检测x是否为此次循环枚举的x坐标，若不是则说明横坐标为x的所有点已遍历完毕，开始遍历下一个x坐标。否则检测这棵树和上面一棵树是否相邻，相邻则中间没有空位。若不相邻，更新答案。无论是否相邻，都要把该点对y轴方向的影响加入树状数组。

#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef long long ll;
const ll mod=2147483648;
pair<int,int>p[100010];
int disx[100010],disy[100010],k;
int r[100010],c[100010],l[100010],ex,ey;
ll com[100010][16],tree[100010];
int lowbit(int &x)
{return x&(-x);}
void update(int x,ll y)
{
    while(x<=ey)
    {
        tree[x]=(tree[x]+y)%mod;
        x+=lowbit(x);
    }
}
ll query(int x)
{
    ll res=0;
    while(x>0)
    {
        res=(res+tree[x])%mod;
        x-=lowbit(x);
    }
    return res;
}
int main()
{
   // freopen("1.txt","r",stdin);
    int n,m,w;
    scanf("%d%d%d",&n,&m,&w);
    for(int i=1;i<=w;++i)
    {
        scanf("%d%d",&p[i].first,&p[i].second);
        disx[i]=p[i].first;
        disy[i]=p[i].second;
    }
    scanf("%d",&k);
    sort(disx+1,disx+w+1);
    ex=unique(disx+1,disx+w+1)-disx-1;
    sort(disy+1,disy+w+1);
    ey=unique(disy+1,disy+w+1)-disy-1;
    sort(p+1,p+w+1);
    for(int i=1;i<=w;++i)
    {
        p[i].first=lower_bound(disx+1,disx+ex+1,p[i].first)-disx;
        p[i].second=lower_bound(disy+1,disy+ey+1,p[i].second)-disy;
    }
    for(int i=1;i<=w;++i)
    {
    r[p[i].second]++;
    c[p[i].first]++;
    }
    com[0][0]=1;
    for(int i=1;i<=w;++i)
    {
        com[i][0]=1;
        for(int j=1;j<=10;++j)
        {
        com[i][j]=(com[i-1][j-1]+com[i-1][j])%mod;
        }
    }
    int con=1;
    ll ans=0;
    p[w+1].first=-1;
    for(int i=1;i<=ex;++i)
    {
        int down=0;
        for(int &j=con;j<=w;++j)
        {
        if(p[j].first!=i)
        break;
        down++;
        l[p[j].second]++;
        if(p[j+1].second-p[j].second==1||p[j+1].first!=i)
        {
        ll mid=(com[l[p[j].second]][k]*com[r[p[j].second]-l[p[j].second]][k]%mod-(query(p[j].second)-query(p[j].second-1))%mod)%mod;
        update(p[j].second,mid);
        continue;
        }
        ans=(ans+com[down][k]*com[c[p[j].first]-down][k]%mod*((query(p[j+1].second-1)-query(p[j].second)+mod)%mod)%mod)%mod;
        ll mid=(com[l[p[j].second]][k]*com[r[p[j].second]-l[p[j].second]][k]%mod-(query(p[j].second)-query(p[j].second-1))%mod+mod)%mod;
        update(p[j].second,mid);
        }
    }
    printf("%lld\n",ans);
    return 0;
}