2698: 太阳轰炸


击中碎片的概率 p=(R1+r)2R22p = \frac {(R1+r)^2}{R2^2}p=R22(R1+r)2​
未击中的概率为 q=R22−(R1+r)2R22q = \frac {R2^2-(R1+r)^2} {R2^2}q=R22R22−(R1+r)2​
摧毁碎片的概率

x=C(n,k)∗pk∗qn−k+C(n,k+1)∗pk+1∗qn−k−1+…+C(n,n)∗pn∗q0x = C(n,k) * p^k * q^{n-k} + C(n,k+1) * p^{k+1} * q^{n-k-1} + … + C(n,n)*p^n*q^0x=C(n,k)∗pk∗qn−k+C(n,k+1)∗pk+1∗qn−k−1+…+C(n,n)∗pn∗q0

#include<cstdio>
const int N = 5000010, p = 1e9+7;
typedef long long LL;
LL fac[N],inv[N];
int ksm(int a,int b)
{int res=1;for(;b;b>>=1){if(b&1) res=(LL)res*a%p;a=(LL)a*a%p;}return res;
}
int C(int a,int b)
{return fac[a]*inv[b]%p*inv[a-b]%p;
}
int main()
{LL n,R1,R2,r,a,h,k;scanf("%lld%lld%lld%lld%lld%lld",&n,&R1,&R2,&r,&a,&h);if(h%a==0) k=h/a;else k=h/a+1;//上取整if(k>n)//k最大到n,>n时不可能,概率为0{puts("0");return 0; }if(R2<=R1+r)//概率最大为1{puts("1");return 0;}fac[0]=1;for(int i=1;i<=n+1;i++) fac[i]=fac[i-1]*i%p;//预处理出阶乘、阶乘的逆元inv[n+1]=ksm(fac[n+1],p-2);for(int i=n;i>=0;i--) inv[i]=inv[i+1]*(i+1)%p;int rn=ksm(R2*R2%p,p-2);//r2的逆元int pp=(R1+r)*(R1+r)%p*rn%p;//即概率pint qq=(R2*R2%p-(R1+r)*(R1+r)%p+p)%p*rn%p;//即概率qint pk=ksm(pp,k);//p^kint qnk=ksm(qq,n-k);//q^(n-k)LL res=0;for(int i=k;i<=n;i++){res=res+(LL)C(n,i)*pk%p*qnk%p;res%=p;pk=(LL)pk*pp%p;qnk=(LL)qnk*ksm(qq,p-2)%p;}printf("%lld\n",res);return 0;
}

ZZULIOJ 2698: 太阳轰炸相关推荐

  1. 2020河南省第二届CCPC真题解析(附榜单)

    这里写目录标题 A:班委竞选 题目: 分析: AC代码: B:广告投放 题目: 分析: AC代码: C:我得重新集结部队 题目: 分析: AC代码: D:园艺大师 题目: 分析: AC代码: E:发通 ...

  2. 阴霾“强势霸屏”长三角 天际上演“失踪的太阳”

    图为民众大雾中出行.孟德龙 摄 中新社南京1月14日电 (杨颜慈)在能见度不足200米的强浓雾中,江苏14日开始上演新一周的"迷雾重重".开年以来,雾霾.雨雪轮番"轰炸& ...

  3. NASA望远镜拍摄太阳表面高能X射线爆发瞬间

    2003年11月4日拍摄的照片,展示了巨型太阳黑子486(右下).太阳黑子486产生一次猛烈的耀斑.此次大爆发发生后不久,致电离辐射袭击地球大气层,导致北美洲出现大面积无线电中断. 新浪科技讯 北京时 ...

  4. 洛谷 P1142 轰炸

    洛谷 P1142 轰炸 题目描述 "我该怎么办?"飞行员klux向你求助. 事实上,klux面对的是一个很简单的问题,但是他实在太菜了. klux要想轰炸某个区域内的一些地方,它们 ...

  5. 太阳直射点纬度计算公式_高中地理——每日讲1题(二十四节气、正午太阳高度、日出时间)...

    知识点 (1)二十四节气 地球绕太阳公转,公转一周接近360°,将它分为24份,即为24节气.因此从立春开始,地球每绕太阳公转15°,即度过一个节气. 地球公转一周的周期为365日5时48分46秒,分 ...

  6. cesium 渲染分析(以太阳为例),实现卫星任务规划中地球赤道平面绘制

    分析太阳的绘制,对绘制赤道和其它天体有很大的启发. 太阳的定义在文件sun.js 中 function Sun() {/*** Determines if the sun will be shown. ...

  7. LSTM模型(基于Keras框架)预测特定城市或者区域的太阳光照量实战

    LSTM模型(基于Keras框架)预测特定城市或者区域的太阳光照量实战 LSTM在解决序列预测的问题时非常强大,因为它们能够存储之前的信息. LSTM是一种时间递归神经网络,它出现的原因是为了解决RN ...

  8. OpenGL(glut)模拟太阳系:太阳,地球,月球,火星

    GitHub项目地址OpenGL(glut)模拟太阳系:太阳,地球,月球,火星 核心代码如下: // orgin.cpp : 此文件包含 "main" 函数.程序执行将在此处开始并 ...

  9. 记《浪潮之巅》-第一版-12.短暂的春秋--与机会失之交臂的公司之一太阳公司(sun microsystems)...

    1.从斯坦福大学孵化出的高科技公司首推SUN. 它的辉煌: 最早进入中国市场并直接与中国政府开展技术合作的计算机公司 全球拥有5万雇员 市值超过2千亿美元,而且远超今天市值排名第二的google的17 ...

最新文章

  1. Intellij IDEA 高效使用教程
  2. gear s3刷android wear,3星gear|3星gear s3自动选择手表刷机办法图文详细教程以及风险介绍...
  3. 很久没写东西了,留个言。
  4. spring-boot-starter-actuator(健康监控)配置和使用
  5. go 使用sarama写入kafka数据时间戳问题
  6. cad怎么设置线的粗细_CAD软件中怎么设置CAD线宽?
  7. mysql查看版本号_十分钟了解MySQL事务机制
  8. Go程序:利用命令行参数做四则运算
  9. 原子操作、互斥锁、读写锁
  10. java poi 页码,apache poi页码
  11. 剑指offer面试题[20]-顺时针打印矩阵
  12. [Network Architecture]ResNext论文笔记(转)
  13. Nvivo的使用学习
  14. html基础教程全集零基础入门到精通汇总整理(附详细的学习路线图)
  15. 小米盒子3增强版 android,小米盒子3增强版拆机评测 怎么样?好不好?
  16. 当前车牌识别相机技术发展现状
  17. Windows中redis设置密码
  18. 自然语言处理--加载使用facebook的预训练 fastText 模型wiki-news-300d-1M.vec
  19. 电脑连接打印机出现的问题(打印机显示未指定)解决方法
  20. excel函数去重_Excel 2010中去除重复项的几种常用技巧

热门文章

  1. typedef给类型起别名
  2. 图数据结构学习:无向图
  3. 5月12日沪深机构资金增仓前50名
  4. 十进制数转换成N进制数
  5. Webpack 基础配置 Vue+Ts
  6. java实现潜艇大战游戏_java潜艇大战游戏源码项目
  7. ubuntu20.0.4+RTX3090配置深度学习环境的踩坑
  8. 三星Exyos4412
  9. 数据库表字段命名规范
  10. js replace函数