【解题报告】牛客挑战赛70 maimai

news/2024/6/16 10:02:33 标签: 动态规划, 算法

题目链接

这个挑战赛的 F F F是我出的,最后 zhoukangyang 爆标了。。。orzorz

记所有有颜色的边的属性集合 S S S

首先在外层容斥,枚举 S ∈ [ 0 , 2 w ) S\in [0,2^w) S[0,2w),计算被覆盖的的边中不包含 S S S 中属性,并且没有被覆盖的边的数目恰好为 i i i 的配对方案数。

暴力的 DP 做法是记录子树内还没有被匹配的点的数目,复杂度 O ( n 5 ) O(n^5) O(n5),不能通过。出题人特意卡了这种做法。

如果一条边没有覆盖,那么所有点对间的路径都不能经过这条边,这样我们可以把一个连通块分成两个子联通块进行求解。但是这样就要记录连通块里面所有的点,无法通过

考虑二项式反演。记 g ( i ) g(i) g(i) 表示钦定断了 i i i 条边(即 i i i条边没有被覆盖)的方案数, f ( i ) f(i) f(i) 表示恰好断了 i i i 条边的方案数,注意这里的下标 i i i 不包含一定不被覆盖的边。那么有:

g ( i ) = ∑ j = i n − 1 ( j i ) f ( j ) ⇒ f ( i ) = ∑ j = i n − 1 ( − 1 ) j − i ( j i ) g j g(i)=\sum _{j=i}^{n-1}\binom{j}{i}f(j)\Rightarrow f(i)=\sum_{j=i}^{n-1}(-1)^{j-i}\binom{j}{i}g_j g(i)=j=in1(ij)f(j)f(i)=j=in1(1)ji(ij)gj

g ( i ) g(i) g(i) 是好算的,也就是剩下的每个连通块内部任意连边的方案数的乘积

h ( n ) h(n) h(n) 表示大小为 n n n 的连通块任意连边的方案数,如果 n n n 为奇数那么答案是 0 0 0,如果 n n n 为偶数那么答案是 ( n − 1 ) × ( n − 3 ) × . . . × 1 (n-1)\times (n-3)\times ...\times 1 (n1)×(n3)×...×1

考虑 DP。设 d p ( u , i , j ) dp(u, i, j) dp(u,i,j) 表示以 u u u 为根的子树,已经断了 i i i 条边,连通块大小为 j j j 的方案数。对于一条边 ( u , v , w ) (u,v,w) (u,v,w) 转移式子如下:

1.1 1.1 1.1 d p ( u , i , j ) × d p ( v , i 2 , j 2 ) × h ( j 2 ) → d p ( u , i + i 2 + 1 , j ) dp(u, i, j) \times dp(v, i_2, j_2) \times h(j_2) \to dp(u, i + i_2+1, j) dp(u,i,j)×dp(v,i2,j2)×h(j2)dp(u,i+i2+1,j)

1.2 1.2 1.2 如果 w ∉ S w\notin S w/S d p ( u , i , j ) × d p ( v , i 2 , j 2 ) → d p ( u , i + i 2 , j + j 2 ) dp(u,i,j)\times dp(v,i_2,j_2)\to dp(u,i+i_2,j+j_2) dp(u,i,j)×dp(v,i2,j2)dp(u,i+i2,j+j2)

这个 DP 的时间复杂度上界是 O ( n 4 ) O(n^4) O(n4) 的,因此总复杂度 O ( 2 w n 4 ) O(2^wn^4) O(2wn4)

但是注意到每个连通块大小都是必须偶数,因此常数极小,实测单次 DP 计算量在 1 0 6 10^6 106 左右,链的情况可以卡满。注意要把 DP 值为 0 的状态跳过,否则无法通过

数据里面造了一些几条链并起来的情况,暴力要跑 4s 以上,std 能稳定在 0.5s 内出解。随机数据下基本卡不了。如果有人暴力冲过去了或者正解被卡常了,出题人在这里谢罪:(

考虑到打这场比赛的大佬肯定还是比较多的,如果场切这道题的大佬们有更精确的分析复杂度的方式欢迎赛后分享。

upd:F 存在 O ( 2 w n 3 ) O(2^wn^3) O(2wn3) 的做法,具体是将DP状态的第一维看成多项式并用点值来维护。详见zhoukangyang的代码。


http://www.niftyadmin.cn/n/5088402.html

相关文章

Android JNI代码语法解释

文章目录 JNI中的JNIEXPORT、JNIIMPORT和JNICALLJVM如何查找native方法①按照JNI规范的命名规则②调用JNI提供的RegsterNatives函数,将本地函数注册到JVM中示例代码 JNI数据类型JNI字符串的处理①获取字符串②释放字符串③创建字符串④其他字符串处理API JNI中的JNI…

统信UOS 1060上通过Fail2Ban来Ban IP

原文链接:统信UOS 1060上通过Fail2Ban来Ban IP hello,大家好啊,今天给大家带来一篇在统信UOS 1060上安装Fail2Ban并且当ip被ban后通过邮件发送通知的文章。Fail2Ban 是一个用于防止暴力攻击的开源软件。它可以扫描日志文件(例如&a…

CleanMyMacX4.12.3最新免费版mac电脑管家

当我们收到一台崭新的mac电脑,第一步肯定是找到一款帮助我们管理电脑运行的“电脑管家”,监控内存运行、智能清理系统垃圾、清理Mac大文件旧文件、消除恶意软件、快速卸载更新软件、隐私保护、监控系统运行状况等。基本在上mac电脑防护一款CleanMyMac就够…

uniapp左滑列表删除

目录 效果index.vuedelSlideLeft.vuedelRightIconfont.css参考插件最后 效果 index.vue 使用页面 引入封装的页面delSlideLeft.vue <template><view class""><view v-for"(item, index) in busOrderList" :key"index" :data-in…

chatgpt GPT-4V是如何实现语音对话的

直接上代码 https://chat.openai.com/voice/get_token 1. 请求内容 Request:GET /voice/get_token HTTP/1.1 Host: ios.chat.openai.com Content-Type: application/json Cookie: _puiduser***Fc9T:16962276****Nph%2Fb**SU%3D; _uasid"Z0FBQUF***nPT0"; __cf_bmBUg…

PHP和JAVA AES加解密问题

1.问题 php加密的java不能正常机密&#xff0c;java加密的php不能正常解密 2.前置参数 key 32位 iv 16位 3.问题解决 因为key是32位的&#xff0c;所以参数需要$methodAES-256-CBC,$options1 4.PHP代码 <?php /**** Created by PhpStorm* User: Noah* Date: 2023/10/10*…

LINUX定时解压缩方案

需求背景 对接客户中某个上游为外包系统&#xff0c;外包系统每日推送压缩文件至指定文件夹下&#xff0c;文件格式为YYYYMMDD_RegReport.zip。由于每日采集文件&#xff0c;无法对接压缩包内文件&#xff0c;需要将推送的压缩文件每日解压为文件夹 需求分析 与客户沟通后&a…

设计模式~调停者(中介者)模式(Mediator)-21

调停者&#xff08;中介者&#xff09;模式(Mediator) &#xff08;1&#xff09;优点 &#xff08;2&#xff09;缺点 &#xff08;3&#xff09;使用场景 &#xff08;4&#xff09;注意事项&#xff1a; &#xff08;5&#xff09;应用实例&#xff1a; 代码 调停者&a…