BZOJ4584 & 洛谷3643 & UOJ204：[APIO2016]划艇——题解

https://www.lydsy.com/JudgeOnline/problem.php?id=4584

https://www.luogu.org/problemnew/show/P3643

http://uoj.ac/problem/204

dp神题，方程很难想……

（一瞬间我以为我已经傻到无可救药了直到我看了题解发现的确很难很神……）

参考1：https://www.luogu.org/blog/bestFy0731/solution-p3643

参考2：https://www.geek-share.com/detail/2694790527.html

参考3：http://45.76.49.80:8000/f/7bd0386ed1b24b8bb390/

看完这三个人还没看懂的话那接下来看我的理解吧。

不难想到对a和b离散化，并且令f[i][j]表示前i个学校，第i个学校必须取且第i个学校的划艇数量在j段区间内。

然后就有一个问题了，这样的话肯定会有几个学校的划艇数量区间重合，这时如何统计个数。

那么根据参考2的思路，考虑另一个问题：给m个长度为len的区间，选k个区间的方案数个数。

显然地为C(m,k)*C(len,k)。

对1~m的k求和即为sigma(C(m,m-k)*C(len,k))=C(m+len,m)。

现在我们就能求出当有m个学校在同一个数量段时的方案数了。

根据参考3，我们可以列出以下的式子：

令k+1~i中可以取j段的划艇数量的学校个数为m，且k+1~i的学校要么取0要么取j段，则：

f[i][j]=sigma(f[k][l]*C(m+len[j],m))(k<i,l<j)

（PS：k+1~i的学校取别的段（比如x）的情况早就被f[i][x]考虑过了。）

但是特殊的，当m=0时，C(m+len[j],m)=len[j]。

完后前缀和优化即可。

#include<map>
#include<cmath>
#include<stack>
#include<queue>
#include<cstdio>
#include<cctype>
#include<vector>
#include<cstdlib>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
typedef long long ll;
const int p=1e9+7;
const int N=505;
inline int read(){
int X=0,w=0;char ch=0;
while(!isdigit(ch)){w|=ch=='-';ch=getchar();}
while(isdigit(ch))X=(X<<3)+(X<<1)+(ch^48),ch=getchar();
return w?-X:X;
}
int n,m,a
,b
,c[N*2];
int s
,inv
,C
;
inline void LSH(){
sort(c+1,c+m+1);
m=unique(c+1,c+m+1)-c-1;
for(int i=1;i<=n;i++){
a[i]=lower_bound(c+1,c+m+1,a[i])-c;
b[i]=lower_bound(c+1,c+m+1,b[i]+1)-c;
}
}
inline void init(){
inv[1]=1;
for(int i=2;i<=n;i++)inv[i]=(ll)(p-p/i)*inv[p%i]%p;
}
int main(){
n=read();init();
for(int i=1;i<=n;i++){
a[i]=read(),b[i]=read();
c[++m]=a[i];c[++m]=b[i]+1;
}
LSH();
s[0]=1;C[0]=1;
for(int i=1;i<m;i++){
int len=c[i+1]-c[i];
for(int j=1;j<=n;j++)C[j]=(ll)C[j-1]*(len+j-1)%p*inv[j]%p;
for(int j=n;j>=1;j--){
if(a[j]<=i&&i+1<=b[j]){
int f=0,mul=len,cnt=1;
for(int k=j-1;k>=0;k--){
(f+=(ll)s[k]*mul%p)%=p;
if(a[k]<=i&&i+1<=b[k])mul=C[++cnt];
}
(s[j]+=f)%=p;
}
}
}
int ans=0;
for(int i=1;i<=n;i++)(ans+=s[i])%=p;
printf("%d\n",ans);
return 0;
}

+++++++++++++++++++++++++++++++++++++++++++

 +本文作者：luyouqi233。 　　　　　　　　　　　　　　+

 +欢迎访问我的博客：http://www.cnblogs.com/luyouqi233/+

+++++++++++++++++++++++++++++++++++++++++++