「BJOI 2019」光线

最新推荐文章于 2021-02-04 16:39:45 发布

forever_dreams

最新推荐文章于 2021-02-04 16:39:45 发布

阅读量152

点赞数

分类专栏： # 递推 # 思维题文章标签：递推数学思维题

本文链接：https://blog.csdn.net/forever_dreams/article/details/102259808

版权

思维题同时被 2 个专栏收录

6 篇文章 0 订阅

订阅专栏

递推

3 篇文章 0 订阅

订阅专栏

传送门

problem

当一束光打到一层玻璃上时，有一定比例的光会穿过这层玻璃，一定比例的光会被反射回去，剩下的光被玻璃吸收。

设对于任意 $x$ ，有 $x\times a_i\%$ 单位的光会穿过它，有 $x\times b_i\%$ 的会被反射回去。

现在 $n$ 层玻璃叠在一起，有 $1$ 单位的光打到第 $1$ 层玻璃上，那么有多少单位的光能穿过所有 $n$ 层玻璃呢？

数据范围： $n≤5×10^5$ 。

solution

参考博客

由于第 $n$ 层的玻璃也有透光率， $n$ 层玻璃也可以看做 $1$ 层，所以我们试着以合并玻璃的角度思考。

如果我们将 $a_1,b_1)$ 和 $a_2,b_2)$ 两块玻璃合并，考虑得到的新玻璃 $(A, B)$ 。

注：先说明一点，前 $i$ 块复合玻璃的反射率不是从 $1$ 射入时的反射率，而是从 $i$ 射入的反射率。

$\begin{aligned}A&=a_1a_2\sum_{i=0}^{\infty}(b_1b_2)^i=\frac{a_1a_2}{1-b_1b_2}\\B&=b_2+a_2^{\,2}b_1\sum_{i=0}^{\infty}(b_1b_2)^i=b_2+\frac{a_2^{\,2}b_1}{1-b_1b_2}\end{aligned}$

推导过程可以找规律，我就不写了。

两两合并后，最后的答案就是 $A$ 。

时间复杂度 $O (n)$ 。

code

#include<bits/stdc++.h>
#define N 500005
#define P 1000000007
using namespace std;
int n,a[N],b[N];
int add(int x,int y)  {return x+y>=P?x+y-P:x+y;}
int dec(int x,int y)  {return x-y< 0?x-y+P:x-y;}
int mul(int x,int y)  {return 1ll*x*y%P;}
int power(int a,int b){
	int ans=1;
	for(;b;b>>=1,a=mul(a,a))  if(b&1)  ans=mul(ans,a);
	return ans;
}
int main(){
	scanf("%d",&n);
	int inv=power(100,P-2);
	for(int i=1;i<=n;++i){
		scanf("%d%d",&a[i],&b[i]);
		a[i]=mul(a[i],inv),b[i]=mul(b[i],inv);
	}
	int A=a[1],B=b[1];
	for(int i=2;i<=n;++i){
		int Inv=power(dec(1,mul(B,b[i])),P-2);
		B=add(b[i],mul(mul(mul(a[i],a[i]),B),Inv));
		A=mul(mul(A,a[i]),Inv);
	}
	printf("%d\n",A);
	return 0;
}