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;
}