日韩性视频-久久久蜜桃-www中文字幕-在线中文字幕av-亚洲欧美一区二区三区四区-撸久久-香蕉视频一区-久久无码精品丰满人妻-国产高潮av-激情福利社-日韩av网址大全-国产精品久久999-日本五十路在线-性欧美在线-久久99精品波多结衣一区-男女午夜免费视频-黑人极品ⅴideos精品欧美棵-人人妻人人澡人人爽精品欧美一区-日韩一区在线看-欧美a级在线免费观看

歡迎訪問 生活随笔!

生活随笔

當前位置: 首頁 > 编程资源 > 编程问答 >内容正文

编程问答

【PKUWC2018】随机游走【Min-Max容斥】【树形dp】【FWT】

發布時間:2023/12/3 编程问答 33 豆豆
生活随笔 收集整理的這篇文章主要介紹了 【PKUWC2018】随机游走【Min-Max容斥】【树形dp】【FWT】 小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

題意:給一棵nnn個點的無權樹和xxxqqq次詢問,每次給定一個點集SSS,詢問從xxx開始每次隨機走一步,SSS中的每個點至少被經過一次的期望步數。

n≤18,q≤5000n\leq 18,q\leq 5000n18,q5000

題目求的相當于是SSS中的所有點 第一次被訪問的時間 的最大值 的期望,發現nnn很小,可以Min-Max容斥成第一次到SSS 中任意一個點的期望時間。

然后考慮 dp。設xxx為根,fuf_ufu?表示從uuu開始到SSS中任意一個點的期望時間。

fu=1degu(ffau+∑v∈sonufv)+1f_u=\frac{1}{deg_u}(f_{fa_u}+\sum _{v\in son_u}f_v)+1fu?=degu?1?(ffau??+vsonu??fv?)+1

degu?fu=ffau+∑v∈sonufv+degudeg_u\cdot f_u=f_{fa_u}+\sum _{v\in son_u}f_v+deg_udegu??fu?=ffau??+vsonu??fv?+degu?

一個套路,待定系數法設fu=kuffau+buf_u=k_uf_{fa_u}+b_ufu?=ku?ffau??+bu?

degu?fu=ffau+∑v∈sonu(kvfu+bv)+degudeg_u\cdot f_u=f_{fa_u}+\sum _{v\in son_u}(k_vf_u+b_v)+deg_udegu??fu?=ffau??+vsonu??(kv?fu?+bv?)+degu?

(degu?∑kv+1)?fu=ffau+∑bv+degu(deg_u-\sum k_v+1)\cdot f_u=f_{fa_u}+\sum b_v+deg_u(degu??kv?+1)?fu?=ffau??+bv?+degu?

解得

kv=1degu?∑kv+1k_v=\frac{1}{deg_u-\sum k_v+1}kv?=degu??kv?+11?

bv=∑bv+degudegu?∑kv+1b_v=\frac{\sum b_v+deg_u}{deg_u-\sum k_v+1}bv?=degu??kv?+1bv?+degu??

如果u∈Su\in SuS,有fu=0f_u=0fu?=0,對它的父親沒有貢獻,直接ku=bu=0k_u=b_u=0ku?=bu?=0

SSS的答案就是fx=bxf_x=b_xfx?=bx?,乘上容斥系數再FWT一下就可以快速回答詢問

#include <iostream> #include <cstdio> #include <cstring> #include <cctype> #include <vector> using namespace std; const int MOD=998244353; typedef long long ll; inline int qpow(int a,int p) {int ans=1;while (p){if (p&1) ans=(ll)ans*a%MOD;a=(ll)a*a%MOD,p>>=1;}return ans; } #define inv(x) qpow(x,MOD-2) inline int add(const int& x,const int& y){return x+y>=MOD? x+y-MOD:x+y;} inline int dec(const int& x,const int& y){return x<y? x-y+MOD:x-y;} int k[20],b[20]; vector<int> e[20]; void dfs(int u,int f,int S) {if (S&(1<<u)) return (void)(k[u]=b[u]=0);int sumk=0,sumb=0;for (int i=0;i<(int)e[u].size();i++)if (e[u][i]!=f)dfs(e[u][i],u,S),sumk=add(sumk,k[e[u][i]]),sumb=add(sumb,b[e[u][i]]);k[u]=inv(dec((int)e[u].size(),sumk)),b[u]=((ll)e[u].size()+sumb)*k[u]%MOD; } int f[1<<20]; int main() {int n,q,x;scanf("%d%d%d",&n,&q,&x);for (int i=1;i<n;i++){int u,v;scanf("%d%d",&u,&v);e[u].push_back(v),e[v].push_back(u);}for (int S=1;S<(1<<n);S++){dfs(x,0,S<<1);f[S]=b[x];int c=0;for (int i=0;i<n;i++) c^=(S>>i)&1;if (!c) f[S]=dec(0,f[S]);}for (int l=0;l<n;l++){int mid=1<<l,len=mid<<1;for (int s=0;s<(1<<n);s+=len)for (int k=0;k<mid;k++)f[s+mid+k]=add(f[s+mid+k],f[s+k]);}while (q--){int S=0;int k;scanf("%d",&k);for (int i=1;i<=k;i++){int x;scanf("%d",&x);S|=1<<x;}printf("%d\n",f[S>>1]);}return 0; }

總結

以上是生活随笔為你收集整理的【PKUWC2018】随机游走【Min-Max容斥】【树形dp】【FWT】的全部內容,希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯,歡迎將生活随笔推薦給好友。