点的统计

给定一个点阵，求五元组(a,b,c,d,e)满足xa=xb小于xc小于xd=xe且ya=yd小于yc小于yb=ye 的个数对10^9+7 取模的结果
【输入】
第一行一个整数 n，表示点的个数。
第二行到第 n+1 行，每行一个数对(x,y)，表示一个点的坐标
【输出】
一个整数为五元组的个数。
【样例】
输入文件 w.in 输出文件 w.out
7
1 1
1 3
2 2
3 1
3 3
4 1
4 3
2
【样例解释】
可以看出，总共有 2 个”w”形
【数据范围】
对于 30%的数据，1<=n<=10
对于 100%的数据，1<=n<=1000,1<=x,y<=1000,保证无重点

这个是青岛二中的模拟题
你看条件，如果你可以看出这个是一个矩形的条件，那这个题就好做了；想象一下，再看看条件。
题解是先找一个矩形，再找这个矩形的对角线，如果两条对角线都存在的话，再找矩形的内部点，这个矩形内部有几个点，这个矩形（只是这个矩形）就又几个解，最后把所有矩形的所有解加起来就是答案啦
这是大概的代码，明白思想；

include cstdio

using namespace std;
int d[1005][1005],v[1005][1005];
struct node
{
int x,y;
} p[1005];
int n,ans=0,x1,x2,y1,y2;
int main()
{
scanf(“%d”,&n);
for(int i=1;i<=n;i++)
{
scanf(“%d %d”,&p[i].x,&p[i].y);
v[p[i].x][p[i].y]=1;
d[p[i].x][p[i].y]=1;
}
for(int i=1;i<=1000;i++)
for(int j=1;j<=1000;j++)
d[i][j]+=d[i-1][j]+d[i][j-1]-d[i-1][j-1];

for(int j=1;j<=n;j++)
{
    for(int i=1;i<=n;i++)
    printf("%d ",d[j][i]);
    printf("\n");
}
    for(int j=1;j<=n;j++)
{
    for(int i=1;i<=n;i++)
    printf("%d ",v[j][i]);
    printf("\n");
}
for(int i=1;i<=n;i++)
for(int j=1;j<=n;j++)
{
    if(p[i].x<p[j].x&&p[i].y<p[j].y)
    {
        x1=p[i].x;x2=p[j].x;
        y1=p[i].y;y2=p[j].y;
        if(v[x1][y2]&&v[x2][y1])
            ans+=d[x2-1][y2-1]-d[x1][y2-1]-d[x2-1][y1]+d[x1][y1];
    }
} 
printf("%d",ans);

}