[HDU4624]Endless Spin 最值反演+DP

本文介绍了最值反演(min-max容斥)算法的基本原理,并详细解释了如何利用该算法求解特定概率问题,包括如何转换问题为求最小值期望,以及如何使用动态规划进行高效计算。

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

首先最值反演(也有人叫min-max容斥)就是:

max{S}=TS(1)|T|+1min{T}

假设位置i被染黑的时间是Xi,那么要求的就是E[max{Xi}],于是转化为对于所有点集S,求E[miniS{Xi}]
这个就好求多了,就是第一次染到该集合的某个点的期望,就是n(n1)2numSnumS表示覆盖至少S中一个点的区间个数,可以用所有区间减去未覆盖的区间数得到,未覆盖的区间数就可以通过相邻两点距离得到。
然后就可以DP,设fi,j,0/1表示考虑前i个球,有j个区间是未覆盖的,被选的球的奇偶性,转移即可。
还要写高精度小数,我没写。。。我只过了样例。。。
代码:
#include<iostream>
#include<cstdio>
#include<cstring>
#include<algorithm>
#define G(x) (((x+1)*(x))>>1)
#define ll long long
using namespace std;
int n;
ll f[55][2555][2];
int main()
{
    int ca;
    scanf("%d",&ca);
    while(ca--)
    {
        scanf("%d",&n);
        memset(f,0,sizeof(f));
        f[0][0][0]=1;
        for(int i=1;i<=n;i++)
            for(int j=0;j<i;j++)
                for(int s=G(i-j-1);s<=G(i);s++)
                    for(int b=0;b<=1;b++)
                        f[i][s][b]+=f[j][s-G(i-j-1)][b^1];                
        double ans=0,tot=G(n);
        for(int i=1;i<=n;i++)
            for(int s=0;s<=G(i);s++)
                if(f[i][s][1]-f[i][s][0]!=0)
                ans+=tot*(f[i][s][1]-f[i][s][0])/(tot-s-G(n-i));
        printf("%.15lf\n",ans);                        
    }
    return 0;
}
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值