【2019 江苏省队集训】Day1 解题报告

【T1】 光影交错

【思路要点】

• $N$ 轮后仪式仍然进行的概率为 $(1-p{)}^N$ ，当 $p=10^{-5}$ 时，取 $N>10^7$ 可以保证该概率在 $10^{-20}$ 以内，因此忽略这部分情况不会导致精度要求不能接受的误差产生。
• 记 $f(i)$ 表示出现 $i$ 次非中性灵力的情况出现的期望次数， $g(i)$ 表示出现 $i$ 点非中性灵力后阳性的灵力严格多于阴性灵力的概率，那么有 $$Ans=\sum _{i=1}^{+\infty }f(i)g(i)\approx \sum _{i=1}^{N}f(i)g(i)$$
• 关于 $g(x)$ 的计算显然有
  $$g(x)=\sum _{i=\lfloor \frac{x}{2}\rfloor +1}^x\left(\genfrac{}{}{0ex}{}{x}{i}\right)p_L^i{p}_D^{x-i}$$
  $$g(x)=\{\begin{array}{rcl}0& & x=0\\ g(x-1)-p_D\times \left(\genfrac{}{}{0ex}{}{x-1}{\frac{x}{2}}\right)p_L^{\frac{x}{2}}{p}_D^{\frac{x-2}{2}}& & x\equiv 0(mod2)\\ g(x-1)+p_L\times \left(\genfrac{}{}{0ex}{}{x-1}{\frac{x-1}{2}}\right)p_L^{\frac{x-1}{2}}{p}_D^{\frac{x-1}{2}}& & x\equiv 1(mod2)\end{array}$$
• 递推需要用到的组合数，该部分时间复杂度为 $O(N)$ 。
• 令 $F(x)={\sum }_{i=0}^{+\infty }f(i)x^i$ ，则有
  $$F(x)=\sum _{i=1}^{+\infty }(1-p{)}^{i-1}[(p_L+p_D)x+(1-p_L-p_D){]}^i=\frac{1-p}{1-(1-p)[(p_L+p_D)x+(1-p_L-p_D)]}$$
• 因此对于 $x\ge 2$ ， $f(x)$ 存在递推式 $$f(x)=\frac{f(x-1)\times (1-p)(p_L+p_D)}{1-(1-p)(1-p_L-p_D)}$$
• 暴力计算 $f(0),f(1)$ ，递推即可，该部分时间复杂度同样为 $O(N)$ 。
• 时间复杂度 $O(N)$ ，其中$N=10^7$。

【代码】

#include<bits/stdc++.h>
using namespace std;
const int MAXN = 1e7 + 5;
const double eps = 1e-12;
typedef long long ll;
typedef long double ld;
typedef unsigned long long ull;
template <typename T> void chkmax(T &x, T y) {
    x = max(x, y); }
template <typename T> void chkmin(T &x, T y) {
    x = min(x, y); } 
template <typename T> void read(T &x) {
    
	x = 0; int f = 1;
	char c = getchar();
	for (; !isdigit(c); c = getchar()) if (c == '-') f = -f;
	for (; isdigit(c); c = getchar()) x = x * 10 + c - '0';
	x *= f;
}
template <typename T> void write(T x