USACO JAN08 Bronze Costume Party 化装晚会
This post is written in Chinese. Please consider using Google Translate
由于N较大,O(N^2)的朴素算法会超时。所以需要用至少O(N*LogN)的算法。
方法如下
- 对数列A[i]升幂排序;
- 对于每一位i,在第i位之后找到小于等于S-A[i]的最大的数的最后一个的位置j
- 则能与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;
}