原题可以转化为将n个人按顺序放在若干个区间中,可以不放
将区间离散化到4n个区间(一直不懂泥萌是怎么离散化到2n个区间的我可能只能弄到4n个qwq
f[i][j]表示前i个区间,分布j个人的方案数
为了去重,转移时枚举第i个区间有k个人,强制第k个人一定在这个区间,其他人可以不选,那么方案是
CkL+k−1
,可重组合,选重复的代表不选
dp时可以原地倒着dp省掉一维也不需要滚动
code:
#include<set>
#include<map>
#include<deque>
#include<queue>
#include<stack>
#include<cmath>
#include<ctime>
#include<bitset>
#include<string>
#include<vector>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<climits>
#include<complex>
#include<iostream>
#include<algorithm>
#define ll long long
using namespace std;
const int maxn = 2100;
const int mod = 1e9+7;
int add(int &a,const int &b){a+=b;if(a>=mod) a-=mod;}
int s[maxn],inv[maxn];
int pw(int x,int k)
{
int re=1;
for(;k;k>>=1,x=(ll)x*x%mod) if(k&1)
re=(ll)re*x%mod;
return re;
}
void pre()
{
s[0]=1; for(int i=1;i<maxn;i++) s[i]=(ll)s[i-1]*i%mod;
inv[maxn-1]=pw(s[maxn-1],mod-2);
for(int i=maxn-2;i>=0;i--) inv[i]=(ll)inv[i+1]*(i+1)%mod;
for(int i=1;i<maxn;i++) inv[i]=(ll)inv[i]*s[i-1]%mod;
}
int n,m;
struct inl{int p,i;}t[maxn]; int tp;
inline bool cmp(const inl x,const inl y){return x.p==y.p?x.i<y.i:x.p<y.p;}
struct node{int l,r;}a[maxn],c[maxn]; int cn;
int ci[maxn][maxn];
int f[maxn];
int main()
{
pre();
scanf("%d",&n);
for(int i=1;i<=n;i++)
{
scanf("%d%d",&a[i].l,&a[i].r);
t[++tp]=(inl){a[i].l,0},t[++tp]=(inl){a[i].r,1};
}
sort(t+1,t+tp+1,cmp); t[0].p=t[1].p-1;
for(int i=1;i<=tp;i++) if(t[i].p!=t[i-1].p)
{
bool flag=false;
if(t[i-1].p+1<t[i].p)
{
if(t[i-1].i==0) c[cn].r=t[i].p-1;
else if(t[i].i==1) flag=true,c[++cn]=(node){t[i-1].p+1,t[i].p};
else c[++cn]=(node){t[i-1].p+1,t[i].p-1};
}
if(!flag) c[++cn]=(node){t[i].p,t[i].p};
}
for(int i=1;i<=cn;i++)
{
int l=c[i].r-c[i].l+1,u=2*n;
ci[i][0]=1;
for(int j=1;j<=u;j++)
ci[i][j]=(ll)ci[i][j-1]*inv[j]%mod*(l+j-1)%mod;
}
for(int i=1;i<=n;i++)
{
int j=1;for(;c[j].r<a[i].l;j++); a[i].l=j;
j=cn;for(;a[i].r<c[j].l;j--); a[i].r=j;
}
f[0]=1;
for(int i=1;i<=cn;i++)
{
for(int j=n;j>=1;j--) if(a[j].l<=i&&i<=a[j].r)
{
int nowi=0;
for(int k=1;k<=j;k++)
{
if(a[j-k+1].l<=i&&i<=a[j-k+1].r) nowi++;
add(f[j],(ll)f[j-k]*ci[i][nowi]%mod);
}
}
}
int re=0;
for(int i=1;i<=n;i++) add(re,f[i]);
printf("%d\n",re);
return 0;
}