【CF1788D】Moving Dots【组合数学】

分析

看到题目枚举所有点集都不行，应该是统计类问题（全情况），思考组合数学，DP。DP貌似不太行，状态太多。

先考虑特殊情况：

1. 点集中只有两个点：最后剩下一个点
2. 点集中有三个点：最后剩一个点
3. 点集中有三个点以上，情况有点复杂

我们想到虽然点集有很多，但是发生碰撞的点的组合只有 n² ,于是考虑枚举最先碰撞的两个点，当这两点最先碰撞的时候，满足条件的点集有多少种。

首先看$(i,j)$最先碰撞的的条件，看哪些点能否选。

1. $(i,j)$中间的点不能选
2. 在 $[j+1,j+d_{i,j}-1]$的点以及 $[i-d_{i,j},i-1]$ 的点不可以。

于是剩下的点都是可选可不选，假设有cnt个，构成好多点集，共有$2^{cnt}$ 个，查找段就用二分，最终是 $O(n^{2}logn)$

升华

这个题的问法是全情况计数，也是计数题，想到排列组合或者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;
}