题目

我暴力过啦

看到这样的东西我们先搬出来\(min-max\)容斥

我们设\(max(S)\)表示\(x\)到达点集\(S\)的期望最晚时间,也就是我们要求的答案了

显然我们也很难求出这个东西,但是我们有\(min-max\)容斥

设\(min(S)\)表示\(x\)第一次到达\(S\)的期望时间,我们就有

\[max(S)=\sum_{T\subseteq S}(-1)^{|T|}min(T)
\]

我们现在只需要求出所有\(min(S)\)之后用\(fwt\)做一个子集和就好了

尽管这是一棵树,但是我并没有推出什么优美的转移方程,我们考虑暴力高消

设\(dp_{x,s}\)表示从\(x\)到集合\(s\)的期望步数

显然如果有\(x\in s\),那么\(dp_{x,s}=0\)

否则

\[dp_{x,s}=1+\sum_{(x,v)\in e}\frac{dp_{v,s}}{d_x}
\]

于是我们对于每一种\(s\)分别列方程转移就好了

复杂度是\(O(2^{n}n^3)\)

但是我们注意到没有包含\(x\)点的集合只有\(2^{n-1}\)个,同时高消的常数小至\(\frac{1}{8}\),同时很多高消都没有跑满,于是还是跑的挺快的

代码

#include<cstdio>
#include<cstring>
#include<iostream>
#include<algorithm>
#define re register
#define LL long long
inline int read() {
char c=getchar();int x=0;while(c<'0'||c>'9') c=getchar();
while(c>='0'&&c<='9') x=(x<<3)+(x<<1)+c-48,c=getchar();return x;
}
const int mod=998244353;
struct E{int v,nxt;}e[40];
int num,n,Q,X,len;
int dp[262145];
int head[19],d[19];
int cnt[262145];
int a[20][20],id[20],vis[20],to[20],inv[20],ans[20];
inline int ksm(int a,int b) {
int S=1;
while(b) {if(b&1) S=1ll*S*a%mod;b>>=1;a=1ll*a*a%mod;}
return S;
}
inline void add(int x,int y) {
e[++num].v=y;e[num].nxt=head[x];head[x]=num;
}
inline void solve(int S) {
int t=0;memset(vis,0,sizeof(vis));
for(re int i=1;i<=n;i++) {
if(S&(1<<(i-1))) continue;
vis[i]=1;id[++t]=i;to[i]=t;
}
if(!vis[X]) return;
memset(a,0,sizeof(a));
for(re int i=1;i<=t;i++) {
int x=id[i];
for(re int j=head[x];j;j=e[j].nxt)
if(vis[e[j].v]) a[i][to[e[j].v]]=inv[d[x]];
a[i][i]=mod-1;a[i][t+1]=mod-1;
}
for(re int i=1;i<=t;i++) {
int p=i;
for(p=i;p<=t;p++) if(a[p][i]) break;
if(p!=i) std::swap(a[i],a[p]);
int now=ksm(a[i][i],mod-2);
for(re int j=n+1;j>=i;--j) a[i][j]=1ll*a[i][j]*now%mod;
for(re int j=i+1;j<=t;j++)
for(re int k=t+1;k>=i;--k) {
a[j][k]=(a[j][k]-1ll*a[j][i]*a[i][k]%mod);
if(a[j][k]<0) a[j][k]=(a[j][k]+mod)%mod;
}
}
ans[t]=a[t][t+1];
for(re int i=t-1;i>=0;--i) {
ans[i]=a[i][t+1];
for(re int j=i+1;j<=t;j++) {
ans[i]-=1ll*a[i][j]*ans[j]%mod;
if(ans[i]<0) ans[i]=(ans[i]+mod)%mod;
}
}
dp[S]=ans[to[X]];
if(cnt[S]&1) return;
dp[S]=mod-dp[S];
}
inline void Fwt(int *f) {
for(re int i=2;i<=len;i<<=1)
for(re int ln=i>>1,l=0;l<len;l+=i)
for(re int x=l;x<l+ln;++x) {
f[x+ln]+=f[x];
if(f[x+ln]>=mod) f[x+ln]%=mod;
}
}
int main() {
n=read(),Q=read();X=read();
inv[1]=1;
for(re int i=2;i<=n;i++) inv[i]=1ll*(mod-mod/i)*inv[mod%i]%mod;
for(re int x,y,i=1;i<n;i++) {
x=read(),y=read(),d[x]++,d[y]++;
add(x,y),add(y,x);
}
len=(1<<n);
for(re int i=1;i<len;i++) cnt[i]=cnt[i>>1]+(i&1);
for(re int i=1;i<len;i++) solve(i);
Fwt(dp);
while(Q--) {
int k=read(),S=0;
for(re int i=1;i<=k;i++) S|=(1<<(read()-1));
printf("%d\n",dp[S]);
}
return 0;
}

我还是来补一下正解吧,据说这是树上随机游走的套路

我们设\(f_x\)表示从\(x\)到点集\(s\)的期望步数

据说树上路径唯一我们可以设\(f_x=A_xf_x+B_x\)

我们写出\(f_x\)的转移

\[f_x=\frac{f_t+\sum_{x\rightarrow c}f_c}{d_x}+1
\]

其中\(c\)是\(x\)的儿子

也就是

\[d_xf_x=f_t+\sum_{x\rightarrow c}f_c+d_x
\]

\[d_xf_x=f_t+\sum_{x\rightarrow c}A_cf_x+\sum_{x\rightarrow c}B_c+d_x
\]

\[(d_x-\sum_{x\rightarrow c}A_c)f_x=f_t+\sum_{x\rightarrow c}B_c+d_x
\]

于是我们现在解得

\[A_x=\frac{1}{d_x-\sum_{x\rightarrow c}A_c}
\]

\[B_x=\frac{\sum_{x\rightarrow c}B_c+d_x}{d_x-\sum_{x\rightarrow c}A_c}
\]

对于在\(S\)集合的点显然满足\(A=B=0\),叶子结点的\(A,B\)我们能直接算,我们一路推到根由于根没有父亲,所以\(f_{rt}=B_{rt}\),这样我们就能把所有的\(f_x\)都算出来了,复杂度是\(O(2^nn)\)

代码就不写了

「PKUWC2018」随机游走的更多相关文章

  1. Loj #2542. 「PKUWC2018」随机游走

    Loj #2542. 「PKUWC2018」随机游走 题目描述 给定一棵 \(n\) 个结点的树,你从点 \(x\) 出发,每次等概率随机选择一条与所在点相邻的边走过去. 有 \(Q\) 次询问,每次 ...

  2. 「PKUWC2018」随机游走(min-max容斥+FWT)

    「PKUWC2018」随机游走(min-max容斥+FWT) 以后题目都换成这种「」形式啦,我觉得好看. 做过重返现世的应该看到就想到 \(min-max\) 容斥了吧. 没错,我是先学扩展形式再学特 ...

  3. LOJ2542. 「PKUWC2018」随机游走

    LOJ2542. 「PKUWC2018」随机游走 https://loj.ac/problem/2542 分析: 为了学习最值反演而做的这道题~ \(max{S}=\sum\limits_{T\sub ...

  4. LOJ #2542「PKUWC2018」随机游走

    $ Min$-$Max$容斥真好用 $ PKUWC$滚粗后这题一直在$ todolist$里 今天才补掉..还要更加努力啊.. LOJ #2542 题意:给一棵不超过$ 18$个节点的树,$ 5000 ...

  5. loj2542「PKUWC2018」随机游走

    题目描述 给定一棵 nn 个结点的树,你从点 xx 出发,每次等概率随机选择一条与所在点相邻的边走过去. 有 QQ 次询问,每次询问给定一个集合 SS,求如果从 xx 出发一直随机游走,直到点集 SS ...

  6. 【LOJ2542】「PKUWC2018」随机游走

    题意 给定一棵 \(n\) 个结点的树,你从点 \(x\) 出发,每次等概率随机选择一条与所在点相邻的边走过去. 有 \(Q\) 次询问,每次询问给定一个集合 \(S\),求如果从 \(x\) 出发一 ...

  7. loj#2542. 「PKUWC2018」随机游走(树形dp+Min-Max容斥)

    传送门 首先,关于\(Min-Max\)容斥 设\(S\)为一个点的集合,每个点的权值为走到这个点的期望时间,则\(Max(S)\)即为走遍这个集合所有点的期望时间,\(Min(S)\)即为第一次走到 ...

  8. LOJ2542. 「PKUWC2018」随机游走【概率期望DP+Min-Max容斥(最值反演)】

    题面 思路 我们可以把到每个点的期望步数算出来取max?但是直接算显然是不行的 那就可以用Min-Max来容斥一下 设\(g_{s}\)是从x到s中任意一个点的最小步数 设\(f_{s}\)是从x到s ...

  9. LOJ 2542 「PKUWC2018」随机游走 ——树上高斯消元(期望DP)+最值反演+fmt

    题目:https://loj.ac/problem/2542 可以最值反演.注意 min 不是独立地算从根走到每个点的最小值,在点集里取 min ,而是整体来看,“从根开始走到点集中的任意一个点就停下 ...

随机推荐

  1. 来自于一个问题的回答对自己的反思 php怎么发送邮件?发送邮件插件PHPMailer

    前言: 昨天用手机无意点了一下博问,看见了一个朋友问了一个关于php发邮件不能添加发件人名称的问题,试着看了一下代码,觉得自己发现了问题所在,谁知道只是一知半解没有真正发现问题所在,看来有一段时间没有 ...

  2. ES6+ 开发 React 组件

    在这里简要的说一下这些语言新特性对 React 应用的开发有什么影响,这些 ES6+ 特性使得 React 开发更简单更有趣. 类 迄今为止,最能体现我们使用 ES6+ 来编写 React 组件的就是 ...

  3. WebFrom 小程序【条件查询与分页整合】

    将前面的条件查询功能与分页显示整合到一个页面中 <%@ Page Language="C#" AutoEventWireup="true" CodeFil ...

  4. Python入门到精通学习书籍推荐!

    1.Python基础教程(第2版 修订版)<Python基础教程(第2版修订版)>包括Python程序设计的方方面面,内容涉及的范围较广,既能为初学者夯实基础,又能帮助程序员提升技能,适合 ...

  5. [转] Mac下 快速写博客的软件 MarsEdit

    正文 marsedit最好. 这东西还是收费的,这里, 我给个注册码: Name: The Blade SN: RSME3-DA4KUN-3EL6Y-MXD2X-LYMT9-6KGX8-4 ~~~~~ ...

  6. Does the C standard guarantee buffers are not touched past their null terminator?

    Question: In the various cases that a buffer is provided to the standard library's many string funct ...

  7. Servlet—Cookie(显示用户上次访问时间、显示商品浏览历史)

    1 . 什么是会话? 会话可简单理解为:用户开一个浏览器,点击多个超链接,访问服务器多个web资源,然后关闭浏览器,整个过程称之为一个会话. 1.1 会话过程中要解决的一些问题? 每个用户在使用浏览器 ...

  8. expressjs项目里使用npm run dev,实现热更新

    只需要在package.json文件里添加一句代码 "dev": "supervisor -i ./views ./bin/www" 添加的位置如上图所示, 然 ...

  9. 微信小程序 数组索引 data-“”解释

    按照官方最新文档循环的方式,索引值是以  wx:for-index="index" 方式写的, 以   parseInt(event.currentTarget.dataset.i ...

  10. python txt文件数据转excel

    txt content: perf.txt 2018-11-12 16:48:58 time: 16:48:58 load average: 0.62, 0.54, 0.56 mosquitto CP ...