漫山遍野都是fake的光影。

題目

[LGP4859] 已經(jīng)沒(méi)有什么好害怕的了

給定兩個(gè)長(zhǎng)度為n的數(shù)組a和b，將a中元素與b中元素配對(duì)，求滿足ai>bj的配對(duì)(i,j)個(gè)數(shù)減去滿足ai<bi的配對(duì)(i,j)個(gè)數(shù)恰好為k的方案數(shù)，保證ab中無(wú)重復(fù)元素。

[某年NOI歡樂(lè)賽] 決斗

給定兩個(gè)長(zhǎng)度為n的數(shù)組a和b，將a中元素與b中元素隨機(jī)配對(duì)，求滿足ai≥bj的配對(duì)(i,j)個(gè)數(shù)k次方的期望。

題解

對(duì)于前一個(gè)問(wèn)題，我們轉(zhuǎn)換為求滿足ai>(≥)bj的配對(duì)(i,j)恰好為k=(n+k)/2的方案數(shù)。這樣就能與第二個(gè)問(wèn)題形式上保持一致。稱這樣配對(duì)的配對(duì)為“配對(duì)”（霧）。

其次將ab從小到大排序，然后依次為a數(shù)組配對(duì)，設(shè)f[i,j]表示前i個(gè)位置上確定了j個(gè)“配對(duì)”的方案數(shù)（跳過(guò)剩下的i-j對(duì)不為“配對(duì)”的配對(duì)的轉(zhuǎn)移），w[i]表示滿足bj≤ai的最大的j，有轉(zhuǎn)移 f[i,j]=f[i-1,j]+f[i-1,j-1]*(w[i]-j+1)。后邊那個(gè)系數(shù)其實(shí)是(w[i]-w[i-1])+(w[i-1]-(j-1))得來(lái)的。

如果你有興趣嘗試dp前i個(gè)位置上恰好有j個(gè)配對(duì)，會(huì)發(fā)現(xiàn)不為“配對(duì)”的情況根本dp不動(dòng)。

考慮對(duì)f[n,i]統(tǒng)一確定剩下的(n-i)個(gè)配對(duì)，記g[i]=f[n,i]*(n-i)!。顯然g[i]的統(tǒng)計(jì)是有重復(fù)的。具體的，設(shè)h[i]為恰好有i個(gè)“配對(duì)”的方案數(shù)，h[i]在g[j]中被統(tǒng)計(jì)C(i,j)次，其中i≥j。

即g[i]=Σ[j≥i] h[j]*C(j,i)，移項(xiàng)得h[i]=g[i] Σ[j>i] h[j]*C(j,i)，可以遞推求解了。

后一個(gè)問(wèn)題的后續(xù)操作已經(jīng)不重要了你說(shuō)是吧

參考實(shí)現(xiàn)

#include <bits/stdc++.h> #define ll long long using namespace std; const int N=2019; const int mod=1e9+9;int n,K,a[N],b[N],w[N],f[N][N],c[N][N];int main() {scanf("%d%d",&n,&K);for(int i=1; i<=n; ++i) scanf("%d",a+i);for(int i=1; i<=n; ++i) scanf("%d",b+i);if((n+K)&1) {puts("0");return 0;}K=(n+K)/2;sort(a+1,a+n+1);sort(b+1,b+n+1);for(int i=1,j=0; i<=n; ++i) {while(j<n&&b[j+1]<=a[i]) ++j;w[i]=j;}for(int i=0; i<=n; ++i) {c[i][0]=1;for(int j=1; j<=i; ++j)c[i][j]=(c[i-1][j-1]+c[i-1][j])%mod;}f[0][0]=1;for(int i=1; i<=n; ++i) {f[i][0]=f[i-1][0];int J=min(i,w[i]);for(int j=1; j<=J; ++j) f[i][j]=(f[i-1][j]+(ll)f[i-1][j-1]*(w[i]-j+1)%mod)%mod;for(int j=J+1; j<=i; ++j)f[i][j]=f[i-1][j];}int fc=1;for(int i=n; i>=K; --i) {w[i]=(ll)f[n][i]*fc%mod;for(int j=i+1; j<=n; ++j) w[i]=(w[i]+mod-(ll)w[j]*c[j][i]%mod)%mod;fc=(ll)fc*(n-i+1)%mod;}printf("%d\n",w[K]);return 0; }

轉(zhuǎn)載于:https://www.cnblogs.com/nosta/p/11174462.html

總結(jié)

以上是生活随笔為你收集整理的[LGP4859,...] 一类奇怪的容斥套DP的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。