如果直接dp,状态里肯定要带上已走过的点的集合,感觉上不太好做。

  考虑一种对期望的minmax容斥:其中Max(S)为遍历完S集合的期望步数,Min(S)为遍历到S集合中一个点的期望步数。当然才不管怎么证,反正看上去非常优美。

  设f[i][S]为由i节点出发的Min(S),显然有f[i][S]=Σf[j][S]/di+1。暴力高斯消元复杂度就炸掉了。

  注意到给出的是一棵树,现在连这个性质都没用到当然没法做。根据一个我没见过的套路,可以考虑把f[i]表示成a·f[fa]+b的形式,大力推一波式子就可以了。

  求出f后,暴力枚举子集容斥进行预处理是O(3n)的,类似高维前缀和直接递推就是O(2nn)。然后就可以O(1)回答每个询问了。

#include<iostream>
#include<cstdio>
#include<cmath>
#include<cstdlib>
#include<cstring>
#include<algorithm>
using namespace std;
#define ll long long
#define N 18
#define P 998244353
char getc(){char c=getchar();while ((c<'A'||c>'Z')&&(c<'a'||c>'z')&&(c<'0'||c>'9')) c=getchar();return c;}
int gcd(int n,int m){return m==0?n:gcd(m,n%m);}
int read()
{int x=0,f=1;char c=getchar();while (c<'0'||c>'9') {if (c=='-') f=-1;c=getchar();}while (c>='0'&&c<='9') x=(x<<1)+(x<<3)+(c^48),c=getchar();return x*f;
}
int n,m,S,p[N],f[N][1<<N],size[1<<N],d[N],a[N],b[N],t;
struct data{int to,nxt;
}edge[N<<1];
void addedge(int x,int y){t++;edge[t].to=y,edge[t].nxt=p[x],p[x]=t;}
int ksm(int a,int k)
{int s=1;for (;k;k>>=1,a=1ll*a*a%P) if (k&1) s=1ll*s*a%P;return s;
}
int inv(int a){return ksm(a,P-2);}
void dfs(int k,int from,int S)
{for (int i=p[k];i;i=edge[i].nxt)if (edge[i].to!=from) dfs(edge[i].to,k,S);if (S&(1<<k)) a[k]=b[k]=0;else{int A=0,B=0;for (int i=p[k];i;i=edge[i].nxt)if (edge[i].to!=from) A=(A+a[edge[i].to])%P,B=(B+b[edge[i].to])%P;a[k]=inv((d[k]-A+P)%P),b[k]=1ll*(B+d[k])*a[k]%P;}
}
void dfs2(int k,int from,int S)
{for (int i=p[k];i;i=edge[i].nxt)if (edge[i].to!=from){f[edge[i].to][S]=(1ll*a[edge[i].to]*f[k][S]+b[edge[i].to])%P;dfs2(edge[i].to,k,S);}
}
int main()
{
#ifndef ONLINE_JUDGEfreopen("loj2542.in","r",stdin);freopen("loj2542.out","w",stdout);const char LL[]="%I64d\n";
#elseconst char LL[]="%lld\n";
#endifn=read(),m=read(),S=read()-1;for (int i=1;i<n;i++){int x=read()-1,y=read()-1;addedge(x,y),addedge(y,x);d[x]++,d[y]++;}for (int i=1;i<(1<<n);i++)for (int j=0;j<n;j++)if (i&(1<<j)) {dfs(j,j,i);dfs2(j,j,i);break;}for (int i=1;i<(1<<n);i++){size[i]=size[i^(i&-i)]+1;if (!(size[i]&1)) f[S][i]=(P-f[S][i])%P;}for (int i=0;i<n;i++)for (int j=1;j<(1<<n);j++)if (j&(1<<i)) f[S][j]=(f[S][j]+f[S][j^(1<<i)])%P;/*for (int i=(1<<n)-1;i;i--)for (int j=i^(i&-i);j;j=j-1&i)f[S][i]=(f[S][i]+f[S][j])%P;*/while (m--){int k=read(),x=0;for (int i=1;i<=k;i++) x|=1<<read()-1;printf("%d\n",f[S][x]);}return 0;
}

转载于:https://www.cnblogs.com/Gloid/p/10151978.html

LOJ2542 PKUWC2018随机游走(概率期望+容斥原理)相关推荐

  1. loj#2542 [PKUWC2018]随机游走 (概率期望、组合数学、子集和变换、Min-Max容斥)

    loj#2542 [PKUWC2018]随机游走 (概率期望.组合数学.子集和变换.Min-Max容斥) 很好很有趣很神仙的题! 题目链接: https://loj.ac/problem/2542 题 ...

  2. LOJ #2542 [PKUWC2018]随机游走 (概率期望、组合数学、子集和变换、Min-Max容斥)

    很好很有趣很神仙的题! 题目链接: https://loj.ac/problem/2542 题意: 请自行阅读 题解首先我们显然要求的是几个随机变量的最大值的期望(不是期望的最大值),然后这玩意很难求 ...

  3. 树上随机游走的期望距离

    无聊写写(顺便背板) d[i]表示i的度数,每条边的长度均为1(虽然不是1也能做,但是好写啊) 我们设FiFi表示i走到它父亲的期望距离,则 Fi=1+(∑j是i的儿子Fj+Fi+1)d[i] Fi= ...

  4. 第2章【思考与练习3】随机游走轨迹模拟。将随机游走的步数增加到100步,计算物体最终与原点的距离。重复多次随机游走过程,观察物体与原点距离的变化趋势。

    P33思考练习题3 1.将随机游走的步数增加到100步,计算物体最终与原点的距离. #第1题 #模拟每步生成方向 import numpy as np steps = 100 rndwlk = np. ...

  5. 最值反演[PKUWC2018][loj2542]随机游走

    前言 这是学习容斥过程中的一个比较裸的题了 题意简介 题目链接 题目大意 给出一棵nnn个点的树,给出树上的一个点xxx 现在进行QQQ次询问,每次询问一个点集,求从xxx点开始进行随机随机游走,第一 ...

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

    LOJ#2542. 「PKUWC2018」随机游走 题目描述 Solution 去过一个点集中所有节点的期望时间不好求,考虑min−maxmin-maxmin−max容斥,转化为求第一次到达某一个点集 ...

  7. 高斯消元配合概率dp-图上随机游走模型

    题目1:[HNOI2013]游走 题目大意 给你一张无向连通图,从第111点随机游走到nnn点.花费为经过的边的编号的和.问你如何安排边的编号使得期望花费最小化. n≤500,m≤n2n \leq 5 ...

  8. 「PKUWC2018」随机游走

    Min-Max容斥真好用 LOJ #2542 题意:给一棵不超过1818个节点的树,50005000次询问,每次问从根随机游走走遍一个集合的期望步数 Solution: 考虑Min-Max容斥 有Ma ...

  9. 概率统计笔记:白噪声与随机游走

    1 白噪声 1.1 定义         白噪声是一个理想中的时序模型,它有一个重要的特性,即序列不相关:一个白噪声序列中的每一个点都独立的来自某个未知的分布,即它们满足独立同分布. 考虑时间序列 . ...

  10. 【数学】稀疏图的随机游走问题

    Description 给出一张n个点,m条边的平面图,从1号点开始随机游走,抵达n号点则结束,问期望步数? n<=5000 Solution 这题在wxh的IOI2019国家候选队论文中也提到 ...

最新文章

  1. 数据蒋堂 | Hadoop - 一把杀鸡用的牛刀
  2. 计算机网络---非归零码、曼彻斯特编码和差分曼彻斯特编码
  3. 裴蜀定理(note)
  4. 利用Python自动发送电子邮件
  5. 谷歌浏览器如何将繁体字设置成中文?
  6. 深度学习编译:MLIR初步
  7. 微信小程序运行的底层逻辑
  8. vue 组件中的钩子函数 不能直接写this
  9. 第五章 多个消费者监听同一个队列
  10. mina mysql_Mina源码阅读笔记(四)—Mina的连接IoConnector2
  11. PPT设置自动保存时间 mac_第17期分享:如何控制PPT演讲汇报时间?
  12. 高德地图 删除marker(指定marker)
  13. 投影仪怎么安装才能得到最大的屏幕?学会这个投影距离公式轻松拿捏
  14. openstack 云计算平台 mitaka 搭建
  15. 明厨亮灶监控系统解决方案,看得见的食品安全
  16. Springboot实现过滤器
  17. 知名科幻作家称AI让富者更富、穷者更穷 | 思科推出全球首款人工智能会议语音助手
  18. 韦东山 IMX6ULL和正点原子_「正点原子NANO STM32开发板资料连载」第三章 MDK5 软件入门1...
  19. 【汇正财经】大盘震荡收涨
  20. oracle crs 4535,11gRAC报错CRS-4535,CRS-4000解决

热门文章

  1. Objective-C 协议(protocol)
  2. Leetcode刷题笔记(部分非原创)(1-20题)
  3. sql server 附加数据库的时候出现“无法为此请求检索数据”
  4. Win7配置SVN详细步骤(服务器和客户端)
  5. arguments的类型是Array吗?
  6. 职业方向网络词汇(不定时更新)
  7. 【图神经网络】从源头探讨 GCN 的行文思路
  8. 毕业一年后我转行NLP 这几点宝贵经验分享给大家
  9. 【ACL2019】最佳长论文阅读笔记,降低机器翻译中的exposure bias
  10. ICLR最佳论文“彩票假设”:如何通过彩票假设构建轻量化模型(上)