POJ 2785 4 Values whose Sum is 0

#include<iostream>
#include<cstring>
#include<algorithm>
#include<cstdio>
using namespace std;
#define maxn 4005
int a[maxn],b[maxn],c[maxn],d[maxn];
int ab[maxn*maxn],cd[maxn*maxn];
int main()
{
    int i,j,N,lo,hi,mid,cou=0,p=0;
    cin>>N;
    for(i=0;i<N;i++)
        cin>>a[i]>>b[i]>>c[i]>>d[i];
    for(i=0;i<N;i++)
        for(j=0;j<N;j++)
        {
            ab[p]=a[i]+b[j];
            p++;
        }
        sort(ab,ab+N*N);
    p=0;
    for(i=0;i<N;i++)
        for(j=0;j<N;j++)
        {
            cd[p]=c[i]+d[j];
            p++;
        }
        sort(cd,cd+N*N);
    for(i=0;i<p;i++)
    {
        lo=0;hi=p;//从p-1和p开始二分都行
        while(lo<hi)//这个循环会找到一个与要找的数相等的最小下标
        {
            mid=(lo+hi)/2;
            if(ab[mid]<-cd[i])//因为是排好序的，所以可以直接利用数组下标二分
                lo=mid+1;
            else hi=mid;
        }
        while(ab[lo]==-cd[i]&&lo<p)//判断lo和lo之后是否相等，因为可能会有数据相同，注意lo<p！！！
        {
            cou++;
            lo++;
        }
    }
    cout<<cou<<endl;
    return 0;
}