bzoj4584: [Apio2016]赛艇

根据题意可想到离散化后dp：dp[i][j][k]表示前i个数在区间j内放k个的方案数，开滚动复杂度O(n^3)。然而被bzoj卡常。。。（约1.8s）优化之后变得更慢了。。

看了AC代码，发现可以换种思路：

#include
#include
#include
#include
#define N 505
#define ll long long
#define MO 1000000007
using namespace std;
int n,m,a[N],b[N],S[N*2],sum[N],dp[N*2][N],inv[N],C[N*2][N];
ll Ans;
int main()
{
	scanf("%d",&n);
	for (int i=1;i<=n;i++)
		scanf("%d%d",&a[i],&b[i]),S[i*2-1]=a[i],S[i*2]=b[i]+1;
	sort(S+1,S+n*2+1);m=1;
	for (int i=2;i<=n*2;i++)
		if (S[i]!=S[i-1]) S[++m]=S[i];
	inv[1]=1;
	for (int i=2;i<=n;i++)
		inv[i]=(ll)(MO-MO/i)*inv[MO%i]%MO;
	for (int i=1;i

复杂度O(n^3)，约0.6s。