序列(DP)(组合数)

Posted fengxunling

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了序列(DP)(组合数)相关的知识,希望对你有一定的参考价值。

这是一个DP题。

我们设(f[i][j][k])表示(i)序列长度中放入了(j)个元素,其中(k)是限定的众数的个数;状态转移方程是
[f[k][i][j]=f[k][i-1][j-1]+f[k][i-1][j]-f[k][i-k-1][j-1]]

前面的两个相加应该比较好理解,也就是我们考虑从上一个长度转移过来,新加入的那个数是原先已经加入过的元素还是没有加入过的元素。

后面减去是因为要去掉不合法的情况——如果新加入的这个数出现了k+1次那么就显然是不合法情况。(注:(i-(i-k-1)=k+1))

因为是多组数据,但是数据范围并不大,而且我们的答案是一步一步推出来的的,所以针对每次询问都重新算一遍不如预处理方便。那么我们就在询问前进行预处理。

之后就是开一个g数组,(g[i])表示的是众数个数小于等于i的时候的个数。之后我们就可以用原先算过的sum数组来累加g的值了。

但是要注意的是,最后一维我们不确定是那些数,而我们又要计算序列种类个数+需要的是不下降的序列,所以我们可以想到是组合数。然后也是同样的,考虑进行预处理。我们预处理从n个里面选出来i个数就行了qwq,但是因为n的范围很大,但我们需要的范围只在m以内,所以预处理到m即可。

最后我们用出现次数乘上它的种类个数再除以总和就是期望了。但是注意因为我们g数组相当于是在计算前缀和(因为要算总数个数),所以最后要差分。

最后注意精度问题,我们最好使用long double。(但是因为long double占用16个字节,是int 类型的4倍啊qwq),一定要注意空间是否会炸。。。。。)

#include<cmath>
#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
using namespace std;
int T,m,n;
long double C[255],sum[255][255][255],f[255][255][255],g[255];

int main(){
    C[0]=true;
    for(register int k=1;k<=251;++k){
        sum[k][0][0]=1;
        for(register int i=1;i<=251;++i){
            for(register int j=1;j<=i;++j){
                sum[k][i][j]=sum[k][i-1][j-1]+sum[k][i-1][j];
                if(k<i) sum[k][i][j]-=sum[k][i-k-1][j-1]; 
            }               
        }
    }
    while(cin>>m>>n){
        memset(g,0,sizeof(g)); 
        for(register int i=1;i<=m;++i)
            C[i]=C[i-1]*(n-i+1)/i;
        for(register int k=1;k<=m;++k)
            for(register int j=1;j<=m;++j)
                g[k]+=sum[k][m][j]*C[j];
        long double ans=0;
        for(register int k=1;k<=m;++k)
            ans+=k*(g[k]-g[k-1])/g[m];
        printf("%.4Lf
",(double)ans);
    }
    return 0;
}

以上是关于序列(DP)(组合数)的主要内容,如果未能解决你的问题,请参考以下文章

bzoj4517[Sdoi2016]排列计数 组合数+乘法逆元+dp

CF1024E Natasha, Sasha and the Prefix Sums——DP/数学(组合数)

P3643 [APIO2016]划艇 dp+组合数

C. On the Bench(dp + 组合数)

POJ 2229 计数DP

一道dp题目