P5643-[PKUWC2018]随机游走【min-max容斥,dp】

正题

题目链接:https://www.luogu.com.cn/problem/P5643


题目大意

给出nnn个点的一棵树,一个人从点xxx开始随机游走,然后QQQ次询问给出一个点集SSS,求期望多少步这个人会经过这个点集中的所有点。

1≤n≤18,1≤Q≤50001\leq n\leq 18,1\leq Q\leq 50001n18,1Q5000


解题思路

整个点集都走完比较难统计,我们可以考虑用min−maxmin-maxminmax容斥转为求走到其中一个点的期望步数。

我们设我们目前枚举的集合是SSS,那么首先有fx=0(x∈S)f_x=0(x\in S)fx=0(xS)

然后有转移方程:
fx=1degx(ffax+∑x→yfy)f_{x}=\frac{1}{deg_x}(f_{fa_x}+\sum_{x\rightarrow y}f_{y})fx=degx1(ffax+xyfy)
惯例的我们设fx=Axffax+Bxf_x=A_xf_{fa_x}+B_xfx=Axffax+Bx
fx=1degx(ffax+∑x→y(Ayfx+By))f_{x}=\frac{1}{deg_x}\left(f_{fa_x}+\sum_{x\rightarrow y}(A_yf_x+B_y)\right)fx=degx1(ffax+xy(Ayfx+By))
fx=1degxffax+sumAfxdeg+1degxsumB+1f_{x}=\frac{1}{deg_x}f_{fa_x}+\frac{sumAf_x}{deg}+\frac{1}{deg_x}sumB+1fx=degx1ffax+degsumAfx+degx1sumB+1
degx−sumAdegxfx=1degxffax+1degxsumB+1\frac{deg_x-sumA}{deg_x}f_x=\frac{1}{deg_x}f_{fa_x}+\frac{1}{deg_x}sumB+1degxdegxsumAfx=degx1ffax+degx1sumB+1
fx=1degx−sumAffax+degx+sumBdegx−sumAf_x=\frac{1}{deg_x-sumA}f_{fa_x}+\frac{deg_x+sumB}{deg_x-sum_A}fx=degxsumA1ffax+degxsumAdegx+sumB

这样我们就可以推出AAABBB,而BxB_xBx就是节点xxxfff值,记gS=Bxg_S=B_xgS=Bx

那么根据min-max容斥如果我们询问集合SSS时答案就是
∑T⊆S(−1)∣T∣+1gT\sum_{T\sube S}(-1)^{|T|+1}g_TTS(1)T+1gT

用个高维前缀和就可以预处理所有集合的答案了。

时间复杂度:O(n2nlog⁡P)O(n2^n\log P)O(n2nlogP)


code

#include<cstdio>
#include<cstring>
#include<algorithm>
#define ll long long
using namespace std;
const ll N=18,P=998244353;
struct node{ll to,next;
}a[N<<1];
ll n,Q,rt,tot,ls[N],deg[N];
ll A[N],B[N],c[1<<N],f[1<<N];
ll power(ll x,ll b){ll ans=1;while(b){if(b&1)ans=ans*x%P;x=x*x%P;b>>=1;}return ans;
}
void addl(ll x,ll y){a[++tot].to=y;a[tot].next=ls[x];ls[x]=tot;deg[y]++;return;
}
void dfs(ll x,ll fa,ll S){ll sumA=0,sumB=0;if((S>>x)&1){A[x]=B[x]=0;return;}for(ll i=ls[x];i;i=a[i].next){ll y=a[i].to;if(y==fa)continue;dfs(y,x,S);sumA=(sumA+A[y])%P;sumB=(sumB+B[y])%P;}ll inv=power((deg[x]-sumA+P)%P,P-2);A[x]=inv;B[x]=(deg[x]+sumB)*inv%P;return;
}
signed main()
{scanf("%lld%lld%lld",&n,&Q,&rt);rt--;for(ll i=1,x,y;i<n;i++){scanf("%lld%lld",&x,&y);x--;y--;addl(x,y);addl(y,x);}ll MS=(1<<n);for(ll s=1;s<MS;s++){memset(A,0,sizeof(A));memset(B,0,sizeof(B));dfs(rt,n,s);f[s]=B[rt];}for(ll s=1;s<MS;s++){c[s]=c[s-(s&-s)]+1;f[s]=((c[s]&1)?f[s]:(P-f[s]));}for(ll i=0;i<n;i++)for(ll s=MS-1;s>=0;s--)if((s>>i)&1)(f[s]+=f[s-(1<<i)])%=P;while(Q--){ll k,s=0;scanf("%lld",&k);for(ll i=0,x;i<k;i++)scanf("%lld",&x),s|=(1<<x-1);printf("%lld\n",f[s]);}return 0;
}

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/317758.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

Rolling The Polygon Gym - 102222B

Rolling The Polygon Gym - 102222B 题意&#xff1a; 给你一个多边形&#xff0c;给你内部一个点Q&#xff0c;多边形在平面上滚动一周&#xff08;当有一个边第二次触地滚动停止&#xff09;&#xff0c;问Q的轨迹长度 题解&#xff1a; 计算几何题目 自己一直不是很擅长…

[杂题训练]CF1228E Another Filling the Grid(容斥),CF936C Lock Puzzle(构造)

文章目录T1&#xff1a;CF1228E Another Filling the GridsolutioncodeT2&#xff1a;CF936C Lock PuzzlesolutioncodeT1&#xff1a;CF1228E Another Filling the Grid 点我 solution 反过来思考&#xff0c;用所有方案数➖不合法方案数 很容易想到的是——容斥&#xff01…

Asp.Net Core 轻松学-经常使用异步的你,可能需要看看这个文章

前言事情的起因是由于一段简单的数据库连接代码引起&#xff0c;这段代码从语法上看&#xff0c;是没有任何问题&#xff1b;但是就是莫名其妙的报错了&#xff0c;这段代码极其简单&#xff0c;就是打开数据库连接&#xff0c;读取一条记录&#xff0c;然后立即更新到数据库中…

CF1146F: Leaf Partition(树形dp)

解析 阴间dp题qwq 不难设计dp&#xff1a; dpx,0:x节点没有被包含、子树内的方案数dp_{x,0}:x节点没有被包含、子树内的方案数dpx,0​:x节点没有被包含、子树内的方案数 dpx,1:x节点被包含、子树内的方案数dp_{x,1}:x节点被包含、子树内的方案数dpx,1​:x节点被包含、子树内的…

Take Your Seat Gym - 102222D

Take Your Seat Gym - 102222D 题意&#xff1a; 第一次是n个人坐飞机&#xff0c;按照1到n的顺序登机&#xff0c;第一个人登机牌丢了&#xff0c;他随机做一个座位&#xff0c;2到n个人上来按照自己的登机牌坐座位&#xff0c;如果他的座位被坐了&#xff0c;就在随机找一个…

CF446D-DZY Loves Games【高斯消元,矩阵乘法】

正题 题目链接:https://www.luogu.com.cn/problem/CF446D 题目大意 给出nnn个点mmm条边的一张无向图&#xff0c;一些点有陷阱&#xff0c;走到时会损失一条生命&#xff0c;总共有kkk条生命&#xff0c;求从111出发随机游走到nnn没有死亡且到终点时仅剩一条命的概率。 1≤n≤…

「BJOI2019」奥术神杖(AC自动机+DP)

文章目录titlesolutioncodetitle solution 令MagicViVjVk...MagicV_i\times V_j\times V_k...MagicVi​Vj​Vk​... 这里对Magicc\sqrt[c]{Magic}cMagic​有一个很巧妙的转换——取对数 Magicc(Magic)1celoge(Magic)1c\sqrt[c]{Magic}(Magic)^{\frac{1}{c}}e^{log_e(Magic)^{…

try.dot.net 的正确使用姿势

来源&#xff1a;https://www.cnblogs.com/7tiny/p/10277600.html【简介】微软官方前不久发布了 try.dot.net 这个有趣的网址&#xff0c;开始只是图个新鲜看了一下&#xff0c;后面通过自身实践过后&#xff0c;发现这着实算是个“有趣”的站点&#xff01;首先我们大概地列举…

A - TOYS POJ - 2318

A - TOYS POJ - 2318 题意&#xff1a; 一个盒子中有n个隔板&#xff0c;分出n1个空间&#xff08;从左往右空间的编号分别是0…n&#xff09;&#xff0c;&#xff08;隔板之间不会相交&#xff0c;且按照从左往右的顺序给出&#xff09;&#xff0c;现在给你m个坐标的物品&…

CF1016F:Road Projects(树形dp)

解析 好题 意思就是我没做出来 稍微分析一下就可以发现加边的位置始终是一样的 换句话说询问完全可以O1 关键就是找到这条边加在哪里 一开始我完全把这道题看成了彻头彻尾的数据结构题 容易想到二分答案 然后上个树状树组搞一搞就行了 但是遇到一个关键的问题 它无法解决加边…

CF750E-New Year and Old Subsequence【动态dp】

正题 题目链接:https://www.luogu.com.cn/problem/CF750E 题目大意 给出一个长度为nnn的数字字符串&#xff0c;qqq次询问给出其的一个子串ttt&#xff0c;询问至少要删除多少个数字才能使得其中包含201720172017这个子序列却不包含201620162016这个子序列。 1≤n,q≤21051\l…

[BJOI2017]魔法咒语(AC自动机+DP+矩阵快速幂)

文章目录titlesolutioncodetitle solution 针对数据编程才是坠吊的&#xff01;&#xff01;&#xff01; 观察数据&#xff0c;发现分隔数据的LLL跨度过大&#xff0c;没有衔接——推测很有可能是分数据做法 ①&#xff1a;考虑L≤100L\le100L≤100的情况 可以暴力DPDPDP转移…

P6242-[模板]线段树3【吉司机线段树】

正题 题目链接:https://www.luogu.com.cn/problem/P6242 题目大意 给出一个长度为nnn的序列aaa&#xff0c;mmm次要求支持操作 区间加上一个值kkk区间所有aia_iai​变为min{ai,k}min\{a_i,k\}min{ai​,k}区间求和区间求最大值区间求历史最大值 1≤n,q≤51051\leq n,q\leq 5\…

洛谷P1912:诗人小G(二分栈、决策单调性)

二分栈&#xff0c;就是通过二分维护的栈 &#xff08;逃&#xff09; 解析 本题的决策单调性可以说是显然 但是本题是同维度&#xff08;其实只有一维&#xff09;自左向右转移&#xff0c;分治的写法是不能奏效的 所以我们使用决策点调性的另一种实现方法&#xff1a;二分栈…

[国家集训队]middle(二分+主席树[中位数思维题])

文章目录点击查看solutioncode点击查看 solution 简单口胡一下就跑 考虑二分答案ansansans 区间[x1,x2],x1∈[a,b]&#xff0c;x2∈[c,d][x1,x2],x1∈[a,b]&#xff0c;x2∈[c,d][x1,x2],x1∈[a,b]&#xff0c;x2∈[c,d] 大于等于ansansans的设为111&#xff0c;小于ansans…

使用Roslyn脚本化C#代码,C#动态脚本实现方案

来源&#xff1a;https://www.cnblogs.com/7tiny/p/10279349.html【前言】Roslyn 是微软公司开源的 .NET 编译器。编译器支持 C# 和 Visual Basic 代码编译&#xff0c;并提供丰富的代码分析 API。Roslyn不仅仅可以直接编译输出&#xff0c;难能可贵的就是上述描述中的开放了编…

Moving On Gym - 102222F

Moving On Gym - 102222F 题意&#xff1a; 有 n 个城市&#xff0c;q 次询问. 给出每个城市的危险度 r 和 城市的邻接矩阵. 每次询问给出 u、v、w&#xff0c;求从 u 到 v 且不经过其他危险度超过 w 的城市的最短路. 题解&#xff1a; floyd 变形 我队友一开始想的是每次…

10.27模拟 总结

前言 220pts 100200100 按照gg的建议从《暴力模式》转化为《切题模式》 但是感觉本次有点难阿… 进入石头门困境… 再四道题中反复横跳 说实话心态是炸的 T1&#xff08;伪&#xff09;正解出来后才安了一些心 由于这次数据的强度确实不高 T1和T4的做法其实都很假 过掉也有运气…

[构造训练]CF1227G Not Same,CF1375H Set Merging,CF1364E X-OR

文章目录T1&#xff1a;CF1227G Not SamesolutioncodeT2&#xff1a;CF1364E X-ORsolutioncodeT3&#xff1a;CF1375H Set Mergingsolutioncode~~脑子是个好东西&#xff0c;希望人人都有构造真的不是个东西&#xff0c;看了一天视频&#xff0c;没有一道题会做~~ T1&#xff…

asp.net core 环境(Development、Staging 、Production)

一.在asp.net core中使用多个环境ASP.NET Core 配置是基于运行时环境, 使用环境变量。ASP.NET Core 在应用启动时读取环境变量ASPNETCORE_ENVIRONMENT&#xff0c;并将该值存储在 IHostingEnvironment.EnvironmentName 中。ASPNETCORE_ENVIRONMENT 可设置为任意值&#xff0c;但…