CF708E Student's Camp

CF708E Student's Camp [* easy]

有一个 \(n\times m\) 的矩形。

每天两种操作,持续 \(k\) 天:

  • 所有最左边的方块有 \(p\) 的概率消失。
  • 所有最右边的方块有 \(p\) 的概率消失。

然后我们在最顶部和最底部增加一层,如果某个时刻,此图分成了两个联通分量,那么就 GG。

求不 GG 的概率。答案对 \(10^9+7\) 取模。

\(n,m\le 1500,k\le 10^5\)

Solution

考虑给每一行确定一个区间,合法当且仅当每个区间都和上一个区间相交。

\(f_{i,l,r}\) 表示考虑到第 \(i\) 行,区间为 \([l,r]\) 的概率和,由于 \([l,r]\) 的长度大于等于 \(1\),所以这件事发生的概率可以视为左边走到 \(l\),右边走到 \(r\) 的概率的乘积,即 \(p_l\times p_r\)

考虑转移,形如:

\[f_{i,l,r}=p_lp_r\times \sum f_{i-1,l',r'}([l',r']\land [l,r]) \]

这样难以优化,考虑减去不相交的区间,那么后者只关乎区间的 \(r\) 端点,设 \(F_{i,R}=\sum f_{i,l,r}[r\le R],G_{i,L}=\sum f_{i,l,r}[l\ge L]\)

不难得到:

\[f_{i,l,r}=p_lp_r(\sum f-F_{i-1,l-1}-G_{i-1,r+1}) \]

然而复杂度还是太高了,考虑直接转移 \(g_{i,r'}=\sum f_{i,l,r}[r=r']\),那么就有:

\[g_{i,r}=p_r(\sum p_l(\sum f-F_{i-1,l-1}-G_{i-1,r+1})) \]

\[g_{i,r}=p_r(\sum p_l)(\sum f-G_{i-1,r+1})-p_r\sum p_lF_{i-1,l-1} \]

\[f_{i,l}=p_l(\sum p_r)(\sum f-F_{i-1,l-1})-p_l\sum p_rG_{i-1,r+1} \]

然后就可以算 \(G\) 了,类似处理 \(F\) 即可。

复杂度 \(\mathcal O(nm)\)

\(Code:\)

#include<bits/stdc++.h>
using namespace std ;
#define Next( i, x ) for( register int i = head[x]; i; i = e[i].next )
#define rep( i, s, t ) for( register int i = (s); i <= (t); ++ i )
#define drep( i, s, t ) for( register int i = (t); i >= (s); -- i )
#define re register
#define int long long
int gi() {
	char cc = getchar() ; int cn = 0, flus = 1 ;
	while( cc < '0' || cc > '9' ) {  if( cc == '-' ) flus = - flus ; cc = getchar() ; }
	while( cc >= '0' && cc <= '9' )  cn = cn * 10 + cc - '0', cc = getchar() ;
	return cn * flus ;
}
const int P = 1e9 + 7 ; 
const int N = 1500 + 5 ; 
const int M = 1e5 + 5 ; 
int n, m, p, ip, t, h[N], H[N], F[N][N], G[N][N] ;
int f[N][N], g[N][N], pl[N], pr[N], fac[M], inv[M] ; 
int fpow(int x, int k) {
	if( k < 0 ) return 0 ; 
	int ans = 1, base = x ;
	while(k) {
		if(k & 1) ans = 1ll * ans * base % P ;
		base = 1ll * base * base % P, k >>= 1 ;
	} return ans ;
}
int C(int x, int y) {
	if(y > x || x < 0 || y < 0) return 0 ;
	return fac[x] * inv[y] % P * inv[x - y] % P ; 
}
signed main()
{
	n = gi(), m = gi() ; int x, y ; 
	x = gi(), y = gi(), p = x * fpow(y, P - 2) % P,
	t = gi(), fac[0] = inv[0] = 1, ip = (1 - p + P) % P ; 
	rep( i, 1, t ) fac[i] = fac[i - 1] * i % P, inv[i] = fpow(fac[i], P - 2) ; 
	rep( i, 1, m ) 
		pl[i] = C(t, i - 1) * fpow(p, i - 1) % P * fpow(ip, t - (i - 1)) % P, 
		pr[i] = C(t, m - i) * fpow(p, m - i) % P * fpow(ip, t - (m - i)) % P ;
	
	H[0] = 1 ; 
	rep( i, 1, n ) {
		int k = H[i - 1], d = 0, df = 0 ;
		rep( r, 1, m ) {
			d = (d + pl[r]) % P, df = (df + pl[r] * G[i - 1][r - 1]) % P,
			g[i][r] = pr[r] * d % P * (k - F[i - 1][r + 1] + P) % P,
			g[i][r] = (g[i][r] - pr[r] * df % P + P) % P ; //1/8
			//1/8
		}
		d = 0, df = 0 ; 
		for(re int l = m; l >= 1; -- l ) {
			d = (d + pr[l]) % P, df = (df + pr[l] * F[i - 1][l + 1]) % P, 
			f[i][l] = pl[l] * d % P * (k - G[i - 1][l - 1] + P) % P, 
			f[i][l] = (f[i][l] - pl[l] * df % P + P) % P ;
		}
		rep( r, 1, m ) G[i][r] = (G[i][r - 1] + g[i][r]) % P ;
		drep( l, 1, m ) F[i][l] = (F[i][l + 1] + f[i][l]) % P ; 
		H[i] = G[i][m] ;
	}
	cout << H[n] << endl ; 
	return 0 ;
}
posted @ 2020-10-08 15:43  Soulist  阅读(140)  评论(0编辑  收藏  举报