题目
\(N\)个位置,每个位置要么不选,要么选\([ a_i, b_i ]\)中的一个数;
问最后的单调上升序列(mod 1e9+7)有多少种;
\(1 \le N \le 500\)
题解
orz abclzr
直接\(dp\)最后一位是什么数字的话只能得到31分
将数字离散化分段,第\(i\)段为\([l_i,r_i)\),设\(f_{i,j}\)表示第i个位置选的数字在第j段的方案数(第0段表示没有)
\[f_{i,j} \ = \sum_{k=0}^{i-1} \sum_{l=0}^{j-1} f_{k,l} \times cal(k+1,i,j) \\
ans = \sum_{i=1}^n \sum_{j=1}^m f_{i,j} \\
\]其中 $ cal(l,r,x) $ 表示 $ [l,r) $ 都不选或者选在第j段并且单调上升的方案数
设 $ [ l,r) $ 这里面有 $ S $ 个包含x区间,x区间的长度为 $ L $
\[cal(l,r,x) = \sum_{i=0}^{S}(^S_i)(^L_{i+1}) = \sum_{i=0}^{S}(^S_{S-i})(^L_{i+1}) \\
思考组合意义:左边选S-i个再在右边选i+1个相当与一起选S+1个\\
cal(l,r,x) = (^{S+L}_{S+1}) \\
\]前缀和优化dp即可:\(O(n^3)\)
#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=1010,mod=1e9+7;
int n,tot,sub[N],L[N],R[N],ny[N],l[N],f[N][N];
void inc(int&x,int y){x+=y;if(x>=mod)x-=mod;}
int main(){
// freopen("boat.in","r",stdin);
// freopen("boat.out","w",stdout);
scanf("%d",&n);
for(int i=1;i<=n;++i){
scanf("%d%d",&L[i],&R[i]);
sub[++tot]=L[i];
sub[++tot]=++R[i];
}
sort(sub+1,sub+tot+1);
tot=unique(sub+1,sub+tot+1)-sub-1;
for(int i=1;i<=n;++i){
L[i]=lower_bound(sub+1,sub+tot+1,L[i])-sub;
R[i]=lower_bound(sub+1,sub+tot+1,R[i])-sub;
}
ny[1]=1;for(int i=2;i<=n;++i)ny[i]=(ll)(mod-mod/i)*ny[mod%i]%mod;
for(int i=0;i<tot;++i)f[0][i]=1;
for(int i=1;i<tot;++i)l[i]=sub[i+1]-sub[i];
for(int i=1;i<=n;++i){
for(int j=L[i];j<R[i];++j){
int C=l[j],a=l[j],b=1;
for(int k=i-1;~k;--k){
inc(f[i][j],(ll)f[k][j-1]*C%mod);
if(L[k]<=j&&j<R[k])C=(ll)C*(++a)%mod*ny[++b]%mod;
}
}
for(int j=1;j<tot;++j)inc(f[i][j],f[i][j-1]);
}
int ans=0;for(int i=1;i<=n;++i)inc(ans,f[i][tot-1]);
cout<<ans<<endl;
return 0;
}
//菜兔兔写的部分分
#include<bits/stdc++.h>
#define pb push_back
using namespace std;
const int N=510,M=2000010,mod=1e9+7;
int n,a[N],b[N],len[N],sub[M],tot;
int f[M],g[M];
void inc(int&x,int y){x+=y;if(x>=mod)x-=mod;}
int main(){
// freopen("boat.in","r",stdin);
// freopen("boat.out","w",stdout);
scanf("%d",&n);
for(int i=1;i<=n;++i){
scanf("%d%d",&a[i],&b[i]);
for(int j=a[i];j<=b[i];++j)sub[++tot]=j;
}
sort(sub+1,sub+tot+1);
tot=unique(sub+1,sub+tot+1)-sub-1;
f[0]=1;for(int i=0;i<=tot;++i)g[i]=1;
for(int i=1;i<=n;++i){
a[i]=lower_bound(sub+1,sub+tot+1,a[i])-sub;
b[i]=lower_bound(sub+1,sub+tot+1,b[i])-sub;
for(int j=a[i];j<=b[i];++j)f[j]=g[j];
for(int j=a[i];j<=tot;++j)inc(g[j]=f[j],g[j-1]);
}
cout<<g[tot]-1<<endl;
return 0;
}