bzoj 4584: [Apio2016]赛艇 dp

最新推荐文章于 2022-03-22 22:09:28 发布

_beginend

最新推荐文章于 2022-03-22 22:09:28 发布

阅读量562

点赞数

分类专栏：动态规划

本文链接：https://blog.csdn.net/qq_33229466/article/details/71085087

版权

动态规划专栏收录该内容

127 篇文章 0 订阅

订阅专栏

题意

在首尔城中，汉江横贯东西。在汉江的北岸，从西向东星星点点地分布着个划艇学校，编号依次为到。每个学校都拥有若干艘划艇。同一所学校的所有划艇颜色相同，不同的学校的划艇颜色互不相同。颜色相同的划艇被认为是一样的。每个学校可以选择派出一些划艇参加节日的庆典，也可以选择不派出任何划艇参加。如果编号为的学校选择派出划艇参加庆典，那么，派出的划艇数量可以在Ai至Bi之间任意选择（Ai<=Bi）。值得注意的是，编号为i的学校如果选择派出划艇参加庆典，那么它派出的划艇数量必须大于任意一所编号小于它的学校派出的划艇数量。输入所有学校的Ai、Bi的值，求出参加庆典的划艇有多少种可能的情况，必须有至少一艘划艇参加庆典。两种情况不同当且仅当有参加庆典的某种颜色的划艇数量不同。
（1<=Ai<=Bi<=10^9），N<=500

分析

首先一个很显然的dp方程就是 $f[i,j]$ 表示前i个最后一个为j的方案数。转移自己yy一下就好了。
第二维显然会炸，但是我们可以注意到题目中的区间会把数轴分成最多n*2-1个部分，那么我们可以把dp状态设为f[i,j]表示前i个最后一个在区间j的方案数。
注意到不好转移，于是我们就多加一维状态：f[i,j,k]表示前i个，最后一个在区间j且区间j已经放了k个的方案数。
$f[i,j,k]=f[i-1,j,1]+\sum_{l=1}^{j-1}\sum_{w=1}^nf[i-1,l,w]*len[j]$

$f[i,j,k]=f[i-1,j,k]+f[i-1,j,k-1]*\frac{len[j]-k+1}{k}$

第一维滚动掉，第一个式子用前缀和优化一下，其余的加点小优化就能过了。

居然在bzoj上垫底的说2333

代码

#include<iostream>
#include<cstdio>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;

typedef long long LL;

const int N=505;
const int MOD=1000000007;

int n,g1,g[N*2],a[N],b[N],f[N*2][N],ny[N],len[N*2],s[N*2];

int ksm(int x,int y)
{
    int ans=1;
    while (y)
    {
        if (y&1) ans=(LL)ans*x%MOD;
        x=(LL)x*x%MOD;y>>=1;
    }
    return ans;
}

int main()
{
    scanf("%d",&n);
    for (int i=1;i<=n;i++) scanf("%d%d",&a[i],&b[i]),b[i]++,g[++g1]=a[i],g[++g1]=b[i];
    for (int i=1;i<=n;i++) ny[i]=ksm(i,MOD-2);
    sort(g+1,g+g1+1);
    g1=unique(g+1,g+g1+1)-g-1;
    for (int i=1;i<=n;i++) a[i]=lower_bound(g+1,g+g1+1,a[i])-g,b[i]=lower_bound(g+1,g+g1+1,b[i])-g;
    g[g1+1]=g[g1]+1;
    for (int i=1;i<=g1;i++) len[i]=g[i+1]-g[i];
    f[0][0]=1;
    for (int i=a[1];i<=b[1]-1;i++) f[i][1]=len[i];
    int mx=b[1];
    s[0]=1;
    for (int j=1;j<g1;j++) s[j]=(s[j-1]+f[j][1])%MOD;
    for (int i=2;i<=n;i++)
    {
        for (int j=a[i-1];j<=g1;j++)
        {
            s[j]=s[j-1];
            for (int k=1;k<=i;k++) s[j]=(s[j]+f[j][k])%MOD;
        }
        for (int j=mx;j<=b[i];j++) s[j]=s[j-1];
        for (int j=a[i];j<b[i];j++)
        {
            for (int k=min(i,len[j]);k>=2;k--) f[j][k]=(f[j][k]+(LL)f[j][k-1]*(len[j]-k+1)%MOD*ny[k]%MOD)%MOD;
            f[j][1]=(f[j][1]+(LL)s[j-1]*len[j]%MOD)%MOD;
        }
        mx=max(mx,b[i]);
    }
    int ans=0;
    for (int i=1;i<=g1;i++)
        for (int j=1;j<=n;j++)
            ans=(ans+f[i][j])%MOD;
    printf("%d",ans);
    return 0;
}