正題

題目鏈接:https://www.luogu.com.cn/problem/P5644

---

題目大意

$n$個人，每個人被選中的權重是$a_i$。每次按照權重選擇一個沒有死掉的人殺死，求第$1$個人最后死的概率。輸出答案對$998244353$取模。

$w_i>0,{\sum }_{i=1}^n{w}_i\le 10^5$

---

解題思路

這個死掉之后概率的分母會變所以挺麻煩的，考慮一下變成每次隨便選擇一個人，如果沒有死就殺掉，這樣每個人被選擇的概率就不變了。

然后考慮到計算恰好最后一個死很麻煩，可以假設第$1$個人死之后至少還剩下集合$T$的人，然后容斥這樣就不需要考慮到剩下的人必須在前面都選過一次了。

設$P(T)$表示$1$死之后剩下集合$T$的人的概率，怎么求這個東西，我們可以枚舉一下殺到$1$之前的輪數（記$S$為全集，$W(S)={\sum }_{x\in S}{w}_x$）
$$P(T)=\sum _{i=0}^{\infty }(\frac{W(S)?W(T)?w_1}{W(S)}{)}^i\frac{w_1}{W(S)}$$
等比數列求和展開一下就是
$$P(T)=\frac{(\frac{W(S)?W(T)?w_1}{W(S)}{)}^{\infty }?1}{\frac{W(S)?W(T)?w_1}{W(S)}?1}\frac{w_1}{W(S)}$$
然后因為那個$\infty$的東西是收斂（也就是等于$0$）的所以
$$P(T)=\frac{W(S)}{W(T)+w_1}\frac{w_1}{W(S)}=\frac{w_1}{w_1+W(T)}$$
就好了

然后答案就是
$$\sum _{T\in S}(?1{)}^{|T|}P(T)=\sum _{T\in S}(?1{)}^{|T|}\frac{w_1}{w_1+W(T)}$$

但是這樣的復雜度是$O(2^n)$的顯然不可能過。

但是我們不難發現的是因為$W(S)\le 10^5$，所以我們可以設$f(i)$表示對于所有集合$T$使得$W(T)=i$的容斥系數和那么答案就變成了
$$\sum _{i=0}^{W(S)}f(i)\frac{w_1}{w_1+i}$$

但是這個$f$怎么求，其實看上去就很生成函數，$f(i)$相等于${\prod }_{i=2}^n(1?x^{w_i})$的第$i$次項系數。

這個東西我們分治+$NTT$求就好了（因為這個做法實際上和分治$NTT$有區別）

時間復雜度$O(m{\log }^{2}m)$（$m=W(S)$

---

code

#include<cstdio> #include<cstring> #include<algorithm> #define ll long long using namespace std; const ll N=4e5+10,P=998244353; ll n,w[N],r[N],x[N],y[N]; struct Poly{ll n,a[N]; }F[20];bool v[20]; 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 NTT(ll *f,ll n,ll op){for(ll i=0;i<n;i++)if(i<r[i])swap(f[i],f[r[i]]);for(ll p=2;p<=n;p<<=1){ll len=p>>1,tmp=power(3,(P-1)/p);if(op==-1)tmp=power(tmp,P-2);for(ll k=0;k<n;k+=p){ll buf=1;for(ll i=k;i<k+len;i++){ll tt=f[i+len]*buf%P;f[i+len]=(f[i]-tt+P)%P;f[i]=(f[i]+tt)%P;buf=buf*tmp%P;}}}if(op==-1){ll inv=power(n,P-2);for(ll i=0;i<n;i++)f[i]=f[i]*inv%P;}return; } void Mul(Poly &a,Poly &b){for(ll i=0;i<a.n;i++)x[i]=a.a[i];for(ll i=0;i<b.n;i++)y[i]=b.a[i];ll l=1;while(l<a.n+b.n)l<<=1;for(ll i=0;i<l;i++)r[i]=(r[i>>1]>>1)|((i&1)?(l>>1):0);for(ll i=a.n;i<l;i++)x[i]=0;for(ll i=b.n;i<l;i++)y[i]=0;NTT(x,l,1);NTT(y,l,1);for(ll i=0;i<l;i++)x[i]=x[i]*y[i]%P;NTT(x,l,-1);for(ll i=0;i<l;i++)a.a[i]=x[i];a.n=a.n+b.n-1;return; } ll findq(){for(ll i=0;i<20;i++)if(!v[i]){v[i]=1;return i;} } ll Solve(ll l,ll r){if(l==r){ll p=findq();for(ll i=0;i<=w[l];i++)F[p].a[i]=0;F[p].a[0]=1;F[p].a[w[l]]=P-1;F[p].n=w[l]+1;return p;}ll mid=(l+r)>>1;ll ls=Solve(l,mid),rs=Solve(mid+1,r);Mul(F[ls],F[rs]);v[rs]=0;return ls; } signed main() {scanf("%lld",&n);for(ll i=1;i<=n;i++)scanf("%lld",&w[i]);ll p=Solve(2,n),ans=0;for(ll i=0;i<F[p].n;i++)(ans+=F[p].a[i]*w[1]%P*power(w[1]+i,P-2)%P)%=P;printf("%lld\n",(ans+P)%P);return 0; }

總結

以上是生活随笔為你收集整理的P5644-[PKUWC2018]猎人杀【NTT,分治】的全部內容，希望文章能夠幫你解決所遇到的問題。

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