JZOJ 4212. 【五校联考1day2】我想大声告诉你

题目

解析

设 (f_{i,j}) 表示 (i+1..n) 个人能受到 (j) 次攻击的概率
因为选人出局的顺序是无所谓的，所以我们设从 (1..n) 依次选人出局
那么转移时需要分类，分类即加

• 选第 (i) 个人出局，那么 (i+1..n) 会受到 (1) 次攻击。先前攻击了 (j-1) 次，(i) 必须要活下来。(f_{i,j} = f_{i-1,j-1} imes (1-p)^{j-1})
• 第 (i) 个人是遭受攻击后被杀，那么它必然在前 (j) 次攻击的某一次后出局。(f_{i,j} = f_{i-1,j} imes (1-(1-p)^j))

然后考虑答案是什么
对于 (k) 次攻击，小 (R) 可以是在所有的 (f_{i,k}) 后在被选中出局，那么他必然在受到攻击后仍活下来
所以 (ans_k = frac{1}{n} imes (1-p)^k imes sum_{i=0}^n f_{i,k})
乘上 (frac{1}{n}) 是因为小 (R) 有 (frac{1}{n}) 的概率在这个位置上

(Code)

#include<cstdio>
#include<cstring>
#define LL long long
using namespace std;

const int N = 2e3 + 5;
const LL P = 258280327;
LL f[N][N] , pow[N];

inline LL fpow(LL x , LL y)
{
	if (x == 0 && y == 0) return 0;
	x = (x + P) % P;
	LL res = 1;
	while (y)
	{
		if (y & 1) res = res * x % P;
		y >>= 1 , x = x * x % P;
	}
	return res;
}

int main()
{
	int T;
	scanf("%d" , &T);
	for(; T; T--)
	{
		int n; LL x , y , p;
		scanf("%d%lld%lld" , &n , &x , &y);
		p = x * fpow(y , P - 2) % P;
		pow[0] = 1;
		for(register int i = 1; i <= n; i++) pow[i] = pow[i - 1] * (1 - p + P) % P;
		memset(f , 0 , sizeof f);
		f[0][0] = 1;
		for(register int i = 1; i < n; i++)
		{
			f[i][0] = f[i - 1][0] * (1 - pow[0] + P) % P;
			for(register int j = 1; j < n; j++)
				f[i][j] = f[i - 1][j - 1] * pow[j - 1] % P ,
				f[i][j] += f[i - 1][j] * (1 - pow[j] + P) % P;
		}
		for(register int k = 0; k < n; k++)
		{
			LL ans = 0;
			for(register int i = 0; i < n; i++) ans = (ans + f[i][k]) % P;
			printf("%lld " , ans * pow[k] % P * fpow(n , P - 2) % P);
		}
		printf("
");
	}
}