POJ 2785 4 Values whose Sum is 0

sort+二分（求上下界）：qsort + 纯二分 会TE：

sort用法：

头文件：
#include <algorithm>
using namespace std;

1.默认的sort函数是按升序排。对应于1）
sort(a,a+n);   //两个参数分别为待排序数组的首地址和尾地址
2.可以自己写一个cmp函数，按特定意图进行排序。对应于2）
例如：
int cmp( const int &a, const int &b ){
    if( a > b )
       return 1;
    else
       return 0;
}
sort(a,a+n,cmp);
是对数组a降序排序
又如：
int cmp( const POINT &a, const POINT &b ){
    if( a.x < b.x )
       return 1;
    else
       if( a.x == b.x ){
          if( a.y < b.y )
             return 1;
          else
             return 0;
        }
       else
          return 0;
}
sort(a,a+n,cmp);
是先按x升序排序，若x值相等则按y升序排

代码如下：

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
#include<algorithm>
#define MAXN 4010

using namespace std;

long long int num;
int n, t2;
int a[MAXN], b[MAXN], c[MAXN], d[MAXN],f1[MAXN*MAXN],f2[MAXN*MAXN];
int cmp(const void *_p, const void *_q)
{
return *((int *)_p) - *((int *)_q);
}
int lower_bd(int* f, int t, int v)
{
/*
    int max, mid, min;
    min = 0;
    max = t;
    for(;;)
    {
        mid = (min + max)    / 2;
        if(mid == min)
            break;
        if(f[mid] < v)
            min = mid;
        else
            max = mid;
    }
    return f[mid] >= v ? mid : mid + 1;
*/
int m;
int first = 0;
while(first < t)
    {
        m = first + (t - first)/2;
if(f[m] >= v) t = m;
else first = m+1;
    }
return first;
}
int up_bd(int* f,int t, int v)
{
/*int max, min, mid;
    max = t;
    min = 0;
    for(;;)
    {
        mid = (max + min) / 2;
        if(mid == min)
            break;
        if(f[mid] <= v)
            min = mid;
        else
            max = mid;
    }
    return f[mid] > v ? mid : mid + 1;
*/
int m;
int first = 0;
while(first < t)
    {
        m = first + (t - first)/2;
if(f[m] <= v) first = m + 1;
else t = m;
    }
return first;
}
void input()
{
int p, ch, i, j, max, min;
while(scanf("%d",&n) == 1)
    {
for(i = 0; i < n; i ++)
            scanf("%d%d%d%d",&a[i],&b[i],&c[i],&d[i]);
        t2 = 0;
for(i = 0; i < n; i ++)
for(j = 0; j < n; j ++)
                {
                    f1[t2] = a[i] + b[j];
                    f2[t2] = c[i] + d[j];
                    ++ t2;
                }
        sort(f1, f1 + t2);
        sort(f2, f2 + t2);
//qsort(f1,t2,sizeof(f1[0]),cmp);
//    qsort(f2,t2,sizeof(f2[0]),cmp);
//for(int i = 0; i < t1; i ++)
//    printf("i=%d f1[i]=%d f2[i]=%d\n",i,f1[i],f2[i]);
        num = 0;
        p = -1;
for(i = 0; i < t2; i ++)
        {
if(f1[i] == p)
                num += ch;
else
            {
                max = up_bd(f2,t2,-f1[i]);
                min = lower_bd(f2,t2,-f1[i]);
                ch = max - min;
                num += ch;
            }
            p = f1[i];
        }
        printf("%lld\n",num);
    }    
}
int main()
{
    input();
return 0;
}

转载于:https://www.cnblogs.com/yuzhaoxin/archive/2012/03/30/2426138.html