POJ 2785 4 Values whose Sum is 0（折半枚举+二分）

4 Values whose Sum is 0

Time Limit: 15000MS		Memory Limit: 228000K
Total Submissions: 25675		Accepted: 7722
Case Time Limit: 5000MS

Description

The SUM problem can be formulated as follows: given four lists A, B, C, D of integer values, compute how many quadruplet (a, b, c, d ) ∈ A x B x C x D are such that a + b + c + d = 0 . In the following, we assume that all lists have the same size n .

Input

The first line of the input file contains the size of the lists n (this value can be as large as 4000). We then have n lines containing four integer values (with absolute value as large as 2²⁸ ) that belong respectively to A, B, C and D .

Output

For each input file, your program has to write the number quadruplets whose sum is zero.

Sample Input

6
-45 22 42 -16
-41 -27 56 30
-36 53 -37 77
-36 30 -75 -46
26 -38 -10 62
-32 -54 -6 45

Sample Output

5

Hint

Sample Explanation: Indeed, the sum of the five following quadruplets is zero: (-45, -27, 42, 30), (26, 30, -10, -46), (-32, 22, 56, -46),(-32, 30, -75, 77), (-32, -54, 56, 30).

Source

Southwestern Europe 2005

 

　　ForwardIter lower_bound(ForwardIter first, ForwardIter last,const _Tp& val)算法返回一个非递减序列[first, last)中的第一个大于等于值val的位置。

     ForwardIter upper_bound(ForwardIter first, ForwardIter last, const _Tp& val)算法返回一个非递减序列[first, last)中第一个大于val的位置。

     lower_bound和upper_bound如下图所示：

题意：

给定各有n个整数的四个数列A、B、C、D。要从每个数列中各取出1个数，使四个数的和为0。求出这样的组合的个数。当一个数列中有多个相同的数字时，把它们作为不同的数字看待。

分析：

所有全都判断一遍不可行。不过将它们对半分成AB和CD再考虑，就可以快速解决了。从2个数列中选择的话只有n2种组合，所以可以进行枚举。先从A、B中取出a、b后，为了使总和为0则需要从C、D中取出c + d = a - b。因此先将从C、D中取数字的n2种方法全部枚举出来，将这些和排好序，这样就可以运用二分搜索了。

 

#include <cstdio>  
#include <algorithm>  
using namespace std;
typedef long long ll;
const int maxn = 4000 + 10;

int n;
ll a[maxn], b[maxn], c[maxn], d[maxn];
ll cd[maxn * maxn];    //C和D中数字的组合方法  

void solve()
{
    //枚举从C和D中取出数字的所有方法  
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++) 
        {
            cd[i * n + j] = c[i] + d[j];
        }
    }
    sort(cd, cd + n * n);

    ll res = 0;
    for (int i = 0; i < n; i++) 
    {
        for (int j = 0; j < n; j++) 
        {
            int CD = -(a[i] + b[j]);
            //取出C和D中和为CD的部分  
            //二分搜索  
            res += upper_bound(cd, cd + n * n, CD) - lower_bound(cd, cd + n * n, CD);//    可能有多个答案
            //lower_bound(first, last, const _Tp& val)算法返回一个非递减序列[first, last)中的第一个大于等于值val的位置。
            //upper_bound(first, last, const _Tp& val)算法返回一个非递减序列[first, last)中第一个大于val的位置。
        }
    }
    printf("%lld
", res);
}

int main()
{
    scanf("%d", &n);
    for (int i = 0; i < n; i++) 
    {
        scanf("%lld%lld%lld%lld", &a[i], &b[i], &c[i], &d[i]);
    }
    solve();
    return 0;
}