[BJOI2019]光线

传送门

题目大意

有 $n$ 层玻璃，每一层都有 $a\%$ 的透光率和 $b\%$ 的反光率（剩下的会被吸收），现有 $1$ 单位光打向第一层玻璃，求有多少单位的光透过所有玻璃，对 $10^9+7$ 取模。

Solution

不要 nt 了，这题模拟会吐的。看样例的图就知道了。

那很自然会想到线性递推。首先很容易想到设 $d{p}_i$ 表示 $1$ 单位光线透过前 $i$ 层玻璃后的光线数。则：

$$d{p}_i=d{p}_{i-1}\times a_i+sth.$$

发现这个 $sth.$ 非常鬼畜，于是尝试再设一个东西方便于表示它。

设 $re{s}_k$ 是穿过前 $i-1$ 层后在第 $i$ 层反射后从前 $i-1$ 层反射回来，并如此 $k$ 次后透过第 $i$ 层的光线数。则：

$$re{s}_k=d{p}_{i-1}\times (b_i\times sth.{)}^k\times a_i$$

可以想到，这里的 $sth.$ 是光线从第 $i-1$ 层往回射时再反射回来的光线，不妨设为 $r{f}_{i-1}$。则：

$$re{s}_k=d{p}_{i-1}\times (b_i\times r{f}_{i-1}{)}^k\times a_i$$

然后可以得出 $dp$ 的递推式：

$$d{p}_i=\sum _{k=0}^{+\infty }re{s}_k=\sum _{k=0}^{+\infty }(d{p}_{i-1}\times (b_i\times r{f}_{i-1}{)}^k\times a_i)=d{p}_{i-1}\times a_i\times \sum _{k=0}^{+\infty }(b_i\times r{f}_{i-1}{)}^k$$

好我们假装这个东西可以线性递推了……个鬼啊。

好的我们发现这玩意儿有个无限等比数列，好我们把它求出来：

$$S=1+(b_i\times r{f}_{i-1})+(b_i\times r{f}_{i-1}{)}^2+(b_i\times r{f}_{i-1}{)}^3+\cdots$$

所以有：

$$(b_i\times r{f}_{i-1})\times S=(b_i\times r{f}_{i-1})+(b_i\times r{f}_{i-1}{)}^2+(b_i\times r{f}_{i-1}{)}^3+\cdots$$

所以有：

$$S-(b_i\times r{f}_{i-1})\times S=1$$

好直接把 $S$ 解出来，$S=\frac{1}{1-(b_i\times r{f}_{i-1})}$。

好了我们发现这就好了……个鬼啊！

$rf$ 怎么求呢？

好问题。

好的其实很简单。

$$r{f}_i=b_i+a_i\times \sum _{k=0}^{+\infty }(b_i\times r{f}_{i-1}{)}^k\times a_i$$

对吧，就是要么直接被 $i$ 直接反射，要么透过去后又是那个反复横跳 $k(k\ge 0)$ 次后透过 $i$ 的光加起来就是力！

所以直接把 $S$ 带进去就好力！

综上，有：
$$d{p}_i=\frac{d{p}_{i-1}\times a_i}{1-(b_i\times r{f}_{i-1})}$$

$$r{f}_i=b_i+\frac{a_i^2}{1-(b_i\times r{f}_{i-1})}$$

注意细节。

Code

#include<bits/stdc++.h>
#define ll long long
#define inf (1<<30)
#define INF (1ll<<60)
using namespace std;
const int MAXN=5e5+10;
const ll MOD=1e9+7;
ll ksm(ll a,ll p){
	ll ret=1;
	while(p){
		if(p&1) ret=ret*a%MOD;
		a=a*a%MOD; p>>=1;
	}return ret;
}
ll inv(ll x){return ksm(x,MOD-2);}//哦逆元搞心态
ll a[MAXN],b[MAXN],dp[MAXN],rf[MAXN];
int main()
{
	int n;
	scanf("%d",&n);
	ll inv100=inv(100);
	for(int i=1;i<=n;i++){
		scanf("%lld%lld",&a[i],&b[i]);
		a[i]=a[i]*inv100%MOD;
		b[i]=b[i]*inv100%MOD;
	}
	dp[1]=a[1];rf[1]=b[1];
	for(int i=2;i<=n;i++){
		rf[i]=(b[i]+a[i]*a[i]%MOD*rf[i-1]%MOD*inv(((1-b[i]*rf[i-1]%MOD)%MOD+MOD)%MOD)%MOD)%MOD;
		dp[i]=(dp[i-1]*a[i]%MOD*inv(((1-b[i]*rf[i-1]%MOD)%MOD+MOD)%MOD)%MOD)%MOD;//减法搞心态
	}
	printf("%lld\n",dp[n]);
}

End

好的校内模拟赛 T1。