UVA 1152 4 Values Whose Sum is Zero 和为0的4个值 （中途相遇）

摘要：中途相遇。对比map，快排+二分查找，Hash效率。

n是4000的级别，直接O（n^4）肯定超，所以中途相遇法，O（n^2）的时间枚举其中两个的和，O（n^2）的时间枚举其他两个的和的相反数，然后O（logN）的时间查询是否存在。

首先试了下map，果断TLE

//TLE
#include<cstdio>
#include<algorithm>
#include<map>
using namespace std;

const int maxn = ;
int data[][maxn];

map<int,int> cnt;

int main()
{

    int T ; scanf("%d",&T);
    int *A = data[], *B = data[], *C = data[],*D = data[];
    map<int,int>::iterator it;
    while(T--){
        int n; scanf("%d",&n);
        for(int i = ; i < n; i++){
            scanf("%d%d%d%d",A+i,B+i,C+i,D+i);
        }

        for(int i = ; i < n; i++)
        for(int j = ; j < n; j++){
            cnt[A[i]+B[j]]++;
        }

        int ans = ;
        for(int i = ; i < n; i++)
        for(int j = ; j < n; j++){
            int tmp = -C[i]-D[j];
            it = cnt.find(tmp);
            if(it!=cnt.end()) ans += it->second;
        }
        printf("%d\n",ans);
        if(T) putchar('\n');
    }

    return ;
}

map，TLE

然后改成了快排+二分查找，4920ms

// runtime 4920ms
#include<cstdio>
#include<algorithm>
using namespace std;

const int maxn = ;
int data[][maxn];
int vec[maxn*maxn];

int _lower_bound(int *A,int L,int R,int v)
{
    int m;
    while(L<R){
        m = (L+R)>>;
        if(A[m]>=v) R = m;
        else L = m+;
    }
    return L;
}

int _upper_bound(int *A,int L,int R,int v)
{
    int m;
    while(L<R){
        m = (L+R)>>;
        if(A[m]>v) R = m;
        else L = m+;
    }
    return L;
}

int main()
{

    int T ; scanf("%d",&T);
    int *A = data[], *B = data[], *C = data[],*D = data[];
    while(T--){
        int n; scanf("%d",&n);
        for(int i = ; i < n; i++){
            scanf("%d%d%d%d",A+i,B+i,C+i,D+i);
        }

        int sz = ;
        for(int i = ; i < n; i++)
        for(int j = ; j < n; j++){
            vec[sz++] = A[i]+B[j];
        }
        sort(vec,vec+sz);

        int ans = ;
        for(int i = ; i < n; i++)
        for(int j = ; j < n; j++){
            int tmp = -(C[i]+D[j]);
            ans +=  _upper_bound(vec,,sz,tmp) - _lower_bound(vec,,sz,tmp);
        }
        printf("%d\n",ans);
        if(T) putchar('\n');
    }

    return ;
}

快拍+二分，4920ms

实际上没有必要每次二分查找，用两个指针，两边都从最小的数开始比较。

快排+计数，2832ms

#include<cstdio>
#include<algorithm>
#include<map>
using namespace std;
typedef pair<int,int> pii;
#define fi first
#define se second
const int maxn = ;
int data[][maxn];

pii cnt1[maxn*maxn];
pii cnt2[maxn*maxn];
int vec[maxn*maxn];

int main()
{
    int T ; scanf("%d",&T);
    int *A = data[], *B = data[], *C = data[],*D = data[];
    while(T--){
        int n; scanf("%d",&n);
        for(int i = ; i < n; i++){
            scanf("%d%d%d%d",A+i,B+i,C+i,D+i);
        }
        int sz = ;
        for(int i = ; i < n; i++)
        for(int j = ; j < n; j++){
            vec[sz++] = A[i]+B[j];
        }
        sort(vec,vec+sz);
        int len1 = ;
        cnt1[len1] = pii(vec[len1],);
        for(int i = ; i < sz; i++){
            if(vec[i] == cnt1[len1].fi) cnt1[len1].se++;
            else { cnt1[++len1].fi = vec[i]; cnt1[len1].se = ; }
        }
        sz = ;
        for(int i = ; i < n; i++)
        for(int j = ; j < n; j++){
            vec[sz++] = -C[i]-D[j];
        }
        sort(vec,vec+sz);
        int len2 = ;
        cnt2[len2] = pii(vec[len2],);
        for(int i = ; i < sz; i++){
            if(vec[i] == cnt2[len2].fi) cnt2[len2].se++;
            else { cnt2[++len2].fi = vec[i]; cnt2[len2].se = ; }
        }

        int p = ,q = ,ans = ;
        while(p<=len1&&q<=len2){
            if(cnt1[p].fi == cnt2[q].fi){
                ans += cnt1[p++].se*cnt2[q++].se;
            }else if(cnt1[p].fi>cnt2[q].fi) q++;
                  else p++;
        }
        printf("%d\n",ans);
        if(T) putchar('\n');
    }

    return ;
}

快排+计数

还有Hash表，写挂了，待补。。。