分析
看到题目枚举所有点集都不行,应该是统计类问题(全情况),思考组合数学,DP。DP貌似不太行,状态太多。
先考虑特殊情况:
- 点集中只有两个点:最后剩下一个点
- 点集中有三个点:最后剩一个点
- 点集中有三个点以上,情况有点复杂
我们想到虽然点集有很多,但是发生碰撞的点的组合只有 n2 ,于是考虑枚举最先碰撞的两个点,当这两点最先碰撞的时候,满足条件的点集有多少种。
首先看 ( i , j ) (i,j) (i,j)最先碰撞的的条件,看哪些点能否选。
- ( i , j ) (i,j) (i,j)中间的点不能选
- 在 [ j + 1 , j + d i , j − 1 ] [j+1,j+d_{i,j}-1] [j+1,j+di,j−1]的点以及 [ i − d i , j , i − 1 ] [i-d_{i,j},i-1] [i−di,j,i−1] 的点不可以。
于是剩下的点都是可选可不选,假设有cnt个,构成好多点集,共有 2 c n t 2^{cnt} 2cnt 个,查找段就用二分,最终是 O ( n 2 l o g n ) O(n^2logn) O(n2logn)
升华
这个题的问法是全情况计数,也是计数题,想到排列组合或者DP。但是全情况的问法,不要傻乎乎去讨论所有可能的点集,而是要找出全情况中的“原子”。
上代码
#include<iostream>
#include<cstdio>
#include<algorithm>
#include<cmath>
using namespace std;
typedef long long ll;
ll ans;
const int mod=1e9+7;
int n,a[3010];
ll ksm(ll x,ll y)
{
ll s=1;
while(y)
{
if(y&1) s=s*x%mod;
x=x*x%mod;
y>>=1;
}
return s%mod;
}
int main()
{
cin>>n;
for(int i=1;i<=n;i++) scanf("%d",&a[i]);
sort(a+1,a+n+1);
for(int i=1;i<=n-1;i++)
{
for(int j=i+1;j<=n;j++)
{
ll cnt=0,d=a[j]-a[i];
int x=lower_bound(a+1,a+n+1,a[j]+d)-a;
int y=lower_bound(a+1,a+n+1,a[i]-d)-a;
cnt=n-x+1+y-1;
ans=(ans+ksm(2,cnt))%mod;
}
}
cout<<ans;
return 0;
}