CodeForces 28 C.Bath Queue(概率DP+组合数学)

探讨了n个孩子在m个洗浴室洗手池排队洗漱的问题,通过动态规划算法求解最长队伍长度的期望值。

摘要生成于 C知道 ,由 DeepSeek-R1 满血版支持, 前往体验 >

Description

n个孩子去m个洗浴室洗漱,第i个洗浴室有ai个洗手池,每次只能给一个孩子用,其他孩子如果选择这个洗手池则需排队,每次孩子会选择最短的队排,求所有队伍中最长的队伍长度的期望值

Input

首先输入两个整数n,m分别表示孩子个数和洗浴室个数,之后m个整数ai表示第i个洗浴室里洗手池数量(1n,m,ai50)

Output

对于每组用例,输出最长队伍长度的期望值,误差不超过109

Sample Input

1 1
2

Sample Output

1.00000000000000000000

Solution

dp[i][j][k]表示还剩i个洗浴室和j个孩子没选择,最长队伍长度为k的期望值,初始化dp[0][0][i]=i

枚举第i个洗浴室中孩子数量l,那么这个洗浴室最长队伍长度为lai,本来这j个孩子每人有i种选择,现在选出l个孩子有Clj种情况,剩下的jl个孩子每人有i1种选择,故有转移方程

dp[i][j][k]+=dp[i1][jl][max(k,lai)](i1)jlCljij

答案即为dp[m][n][0]

Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
const int INF=0x3f3f3f3f,maxn=55;
double P[maxn][maxn],C[maxn][maxn],dp[maxn][maxn][maxn];
void init(int n=50)
{
    C[0][0]=1;
    for(int i=1;i<=n;i++)
    {
        C[i][0]=C[i][i]=1;
        for(int j=1;j<i;j++)C[i][j]=C[i-1][j-1]+C[i-1][j];
    }
    for(int i=0;i<=n;i++)
    {
        P[i][0]=1;
        for(int j=1;j<=n;j++)P[i][j]=P[i][j-1]*i;
    }
}
int main()
{
    init();
    int n,m,a[maxn];
    while(~scanf("%d%d",&n,&m))
    {
        for(int i=1;i<=m;i++)scanf("%d",&a[i]);
        memset(dp,0,sizeof(dp));
        for(int i=0;i<=n;i++)dp[0][0][i]=i;
        for(int i=1;i<=m;i++)
            for(int j=0;j<=n;j++)
                for(int k=0;k<=n;k++)
                    for(int l=0;l<=j;l++)
                        dp[i][j][k]+=dp[i-1][j-l][max(k,(l+a[i]-1)/a[i])]*P[i-1][j-l]/P[i][j]*C[j][l];
        printf("%.12f\n",dp[m][n][0]);

    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值