博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
LOJ#3093. 「BJOI2019」光线(递推+概率期望)
阅读量:5995 次
发布时间:2019-06-20

本文共 2139 字,大约阅读时间需要 7 分钟。

题面

题解

\(a_i\)\(b_i\)都变成小数的形式,记\(f_i\)表示\(1\)单位的光打到第\(i\)个玻璃上,能从第\(n\)个玻璃下面出来的光有多少,记\(g_i\)表示能从第\(i\)块玻璃反射出来的光有多少,,递推式的话,我们枚举一下这束光在\(i\)\(i+1\)块玻璃之间反射了几次就可以了

\[ \begin{aligned} f_i &=a_i\left(f_{i+1}+g_{i+1}\times b_i\times f_{i+1}+g_{i+1}\times b_i\times g_{i+1}\times b_i\times f_{i+1}+...\right)\\ &=a_if_{i+1}\sum\limits_{k=0}^\infty (b_i\times g_{i+1})^k\\ &=a_if_{i+1}{1\over 1-b_i\times g_{i+1}}\\ g_i &=b_i+a_ig_{i+1}a_i+a_ig_{i+1}b_ig_{i+1}a_i+...\\ &=b_i+a_i^2g_i\sum\limits_{k=0}^\infty (b_i\times g_{i+1})^k\\ &=b_i+a_i^2g_{i+1}{1\over 1-b_ig_{i+1}} \end{aligned} \]

然而问题来了,如果\(b_i=0\)的特殊情况该怎么办?

我们发现在这种情况下,手玩出来的和代入柿子计算的值似乎是一样的?

所以直接代入柿子并没有问题的说……

递推即可

//minamoto#include
#define R register#define inline __inline__ __attribute__((always_inline))#define fp(i,a,b) for(R int i=(a),I=(b)+1;i
I;--i)#define go(u) for(int i=head[u],v=e[i].v;i;i=e[i].nx,v=e[i].v)template
inline bool cmax(T&a,const T&b){return a
inline bool cmin(T&a,const T&b){return a>b?a=b,1:0;}using namespace std;char buf[1<<21],*p1=buf,*p2=buf;inline char getc(){return p1==p2&&(p2=(p1=buf)+fread(buf,1,1<<21,stdin),p1==p2)?EOF:*p1++;}int read(){ R int res,f=1;R char ch; while((ch=getc())>'9'||ch<'0')(ch=='-')&&(f=-1); for(res=ch-'0';(ch=getc())>='0'&&ch<='9';res=res*10+ch-'0'); return res*f;}const int N=5e5+5,P=1e9+7,inv=570000004;inline int add(R int x,R int y){return x+y>=P?x+y-P:x+y;}inline int dec(R int x,R int y){return x-y<0?x-y+P:x-y;}inline int mul(R int x,R int y){return 1ll*x*y-1ll*x*y/P*P;}int ksm(R int x,R int y){ R int res=1; for(;y;y>>=1,x=mul(x,x))(y&1)?res=mul(res,x):0; return res;}int a[N],b[N],f[N],g[N],res,n,tmp;int main(){// freopen("testdata.in","r",stdin); n=read(); fp(i,1,n)a[i]=read(),b[i]=read(),a[i]=mul(a[i],inv),b[i]=mul(b[i],inv); f[n]=a[n],g[n]=b[n]; fd(i,n-1,1){ tmp=ksm(P+1-mul(g[i+1],b[i]),P-2); f[i]=1ll*a[i]*f[i+1]%P*tmp%P; g[i]=add(b[i],1ll*a[i]*a[i]%P*g[i+1]%P*tmp%P); } printf("%d\n",f[1]); return 0;}

转载于:https://www.cnblogs.com/bztMinamoto/p/10754755.html

你可能感兴趣的文章
微软BI 之SSIS 系列 - 导出数据到 Excel 2013 的实现
查看>>
php 学习路线 赵兴壮2014年4月28 日 加油
查看>>
高精确度且线程分离的定时器——多媒体定时器
查看>>
Linux命令工具基础04 磁盘管理
查看>>
设计模式---建造者模式Builder(创建型)
查看>>
SVG
查看>>
maven web配置发布路径 cargo自动部署项目到tomcat
查看>>
linxu select 返回值
查看>>
代码中特殊的注释技术——TODO、FIXME和XXX的用处
查看>>
Android开发(20)--RadioGroup的使用
查看>>
iphone开发之获取网卡的MAC地址和IP地址
查看>>
【网站国际化必备】Asp.Net MVC 集成Paypal(贝宝)快速结账 支付接口 ,附源码demo...
查看>>
java中不常见的keyword:strictfp,transient
查看>>
INDEX--创建索引和删除索引时的SCH_M锁
查看>>
linux C(hello world)
查看>>
微信平台BAE
查看>>
Java程序编译和运行的过程
查看>>
数学图形之牟合方盖
查看>>
Objective-C-类(static)方法、实例方法、overwrite(覆写)、属性(property)复习...
查看>>
PHP多次调用Mysql存储过程报错解决办法
查看>>