小B的夏令营（前缀和优化DP）

嗯，普及组好题。

#include<bits/stdc++.h>
#define maxn 1505
#define mod 1000000007
using namespace std;

int n,m,a,b,p,K;
int f[2][maxn],g[2][maxn],gq[2][maxn],h[maxn],sh[maxn],inv[maxn];
void add(int &a,int b){ if((a+=b)>=mod) a-=mod; }
int Pow(int b,int k){ int r=1;for(;k;k>>=1,b=1ll*b*b%mod) if(k&1) r=1ll*r*b%mod; return r; }

int main(){
	scanf("%d%d%d%d%d",&n,&m,&a,&b,&K);
	p = a * 1ll * Pow(b,mod-2) % mod;
	inv[0] = inv[1] = 1;
	for(int i=2;i<=m;i++) inv[i] = 1ll * (mod - mod / i) * inv[mod % i] % mod;
	int _ivp = Pow(mod+1-p , mod-2); 
	for(int i=0,C=1,P=Pow(mod+1-p,K);i<m;i++)
		h[i] = 1ll * C * P % mod,
		C = 1ll * C * (K-i) % mod * inv[i+1] % mod,
		P = 1ll * P * p % mod * _ivp % mod;
	sh[0] = h[0];
	for(int i=1;i<=m;i++)
		sh[i] = (h[i] + sh[i-1]) % mod;
	int now = 1 , pre = 0;
	g[now][m] = 1;
	for(int i=1;i<=n;i++){
		swap(now,pre);
		memset(f[now],0,sizeof f[now]);
		memset(g[now],0,sizeof g[now]);
		memset(gq[now],0,sizeof gq[now]);
		for(int j=1;j<=m;j++)
			f[now][j] = h[m-j] * 1ll * (sh[j-1] * 1ll * (g[pre][m] - g[pre][m-j]) % mod - gq[pre][j-1]) % mod;
		for(int j=1;j<=m;j++)
			g[now][j] = (g[now][j-1] + f[now][j]) % mod,
			gq[now][j] = (gq[now][j-1] + 1ll * h[j] * g[now][j]) % mod;
	}
	printf("%d\n",(g[now][m]+mod)%mod);
}