[BJOI2019] 光线

最新推荐文章于 2023-04-20 12:39:00 发布

weixin_30550081

最新推荐文章于 2023-04-20 12:39:00 发布

阅读量50

点赞数

原文链接：http://www.cnblogs.com/PaperCloud/p/11202344.html

版权

传送门

Description

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

有\(n\)层玻璃，已知每层玻璃的发射率和折射率，求从最上面一层射入\(1\)单位的光，有多少单位的光将会穿透这些玻璃

Solution

推式子。

Code

#include<bits/stdc++.h>
#define reg register
#define ll long long
#define db double
using namespace std;
inline int read()
{
    int x=0,f=1;char ch=getchar();
    while(ch<'0'||ch>'9'){if(ch=='-')f=-1;ch=getchar();}
    while(ch>='0'&&ch<='9'){x=(x<<3)+(x<<1)+ch-'0';ch=getchar();}
    return x*f;
}
const int MN=5e5+5,P=1e9+7;
int n,a[MN],b[MN],c[MN],f[MN];
int Mul(int x,int y){return (1ll*x*y%P);}
int Add(int x,int y){return (x+y)%P;}
int inv(int x)
{
    int m=P-2,r=1;
    for(;m;m>>=1,x=Mul(x,x))
        if(m&1) r=Mul(r,x);
    return r;
}
int main()
{
    n=read();
    register int i;int inv100=inv(100);
    for(i=1;i<=n;++i) a[i]=Mul(read(),inv100),b[i]=Mul(read(),inv100);
    f[1]=a[1];c[1]=b[1];
    for(i=2;i<=n;++i)
    {
        c[i]=Add(b[i],Mul(Mul(c[i-1],Mul(a[i],a[i])),inv(Add(1,P-Mul(c[i-1],b[i])))));
        f[i]=Mul(f[i-1],Mul(a[i],inv(Add(1,P-Mul(b[i],c[i-1])))));
    }
    printf("%d\n",f[n]);
    return 0;
}

Blog来自PaperCloud，未经允许，请勿转载，TKS！

转载于:https://www.cnblogs.com/PaperCloud/p/11202344.html