USACO JAN08 Bronze Costume Party 化装晚会

This post is written in Chinese. Please consider using Google Translate

由于N较大,O(N^2)的朴素算法会超时。所以需要用至少O(N*LogN)的算法。

方法如下

  1. 对数列A[i]升幂排序;
  2. 对于每一位i,在第i位之后找到小于等于S-A[i]的最大的数的最后一个的位置j
  3. 则能与i配对的数的个数为i-j

运用二分查找可以使时间复杂度降为O(N*LogN)

/*
ID: cmykrgb1
PROG: costume
LANG: C++
*/

#include <iostream>
#define MAX 20001
using namespace std;

int N,S,Ans;
int A[MAX];

inline int cmp(const void *a,const void *b)
{
    return *(int *)a-*(int *)b;
}

void init()
{
    int i;
    freopen("costume.in","r",stdin);
    freopen("costume.out","w",stdout);
    cin >> N >> S;
    for (i=1;i<=N;i++)
        cin >> A[i];
    qsort(A+1,N,sizeof(A[0]),cmp);
}

int bs(int p,int *S,int a,int b) //在[a,b]内找小于等于p的最大的数的最后一个
{
    while (b-a>1)
    {
        int m=(a+b)/2;
        if (p<S[m])
            b=m-1;
        else
            a=m;
    }
    if (p>=S[b])
        return b;
    return a;
}

int main()
{
    init();
    for (int i=1;i<=N;i++)
        Ans+=bs(S-A[i],A,i,N)-i;
    cout << Ans << endl;
    return 0;
}

Related posts