Processing math: 12%

luoguP3175 [HAOI2015]按位或 min-max容斥 + 高维前缀和


考虑min-max容斥

E[max(S)] = \sum \limits_{T \subset S} min(T)

min(T)是可以被表示出来

即所有与T有交集的数的概率的和的倒数

通过转化一下,可以考虑求所有与T没有交集的数的概率和

即求T的补集的子集的概率和

用FMT随意做下吧...

注意:概率为1的时候需要特判

复杂度O(2^n * n)


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

#define de double
#define ri register int
#define rep(io, st, ed) for(ri io = st; io <= ed; io ++)
#define drep(io, ed, st) for(ri io = ed; io >= st; io --)
    
const int sid = (1 << 20) + 25;

int n, show;
de Max, sub[sid];

int main() {
    scanf("%d", &n);
    rep(i, 0, (1 << n) - 1) {
        scanf("%lf", &sub[i]);
        show |= i * (sub[i] > 1e-8);
    }
    if(show != (1 << n) - 1) { puts("INF"); return 0; }
    
    rep(i, 1, n) rep(S, 0, (1 << n) - 1)
        if(!(S & (1 << i - 1))) 
            sub[S ^ (1 << i - 1)] += sub[S];
    
    int T = (1 << n) - 1;
    rep(S, 1, (1 << n) - 1) { // no 0
        if(__builtin_popcount(S) & 1) Max += 1.0 / (1.0 - sub[T ^ S]);
        else Max -= 1.0 / (1.0 - sub[T ^ S]);
    }
    printf("%.12lf\n", Max);
    return 0;
}
posted @   remoon  阅读(132)  评论(0编辑  收藏  举报
编辑推荐:
· 智能桌面机器人:用.NET IoT库控制舵机并多方法播放表情
· Linux glibc自带哈希表的用例及性能测试
· 深入理解 Mybatis 分库分表执行原理
· 如何打造一个高并发系统?
· .NET Core GC压缩(compact_phase)底层原理浅谈
阅读排行:
· 手把手教你在本地部署DeepSeek R1,搭建web-ui ,建议收藏!
· 新年开篇:在本地部署DeepSeek大模型实现联网增强的AI应用
· Janus Pro:DeepSeek 开源革新,多模态 AI 的未来
· 互联网不景气了那就玩玩嵌入式吧,用纯.NET开发并制作一个智能桌面机器人(三):用.NET IoT库
· 【非技术】说说2024年我都干了些啥
点击右上角即可分享
微信分享提示