2020年河南省CCPC大学生程序设计竞赛 --- 太阳轰炸

最新推荐文章于 2023-04-19 19:16:49 发布

chstor

最新推荐文章于 2023-04-19 19:16:49 发布

阅读量859

点赞数 1

分类专栏： ------Online Judge-------

本文链接：https://blog.csdn.net/weixin_45821690/article/details/116493665

版权

------Online Judge------- 专栏收录该内容

3 篇文章 0 订阅

订阅专栏

该博客探讨了一种概率计算问题，涉及虚空碎片摧毁的条件。当给定一定的半径和高度时，计算在多次投掷中摧毁虚空碎片所需的最小次数及其概率。使用了二项分布、费马小定理和逆元来解决这个问题，并提供了C++代码实现。文章强调了在计算过程中注意取模和数组大小限制的重要性。

摘要由CSDN通过智能技术生成

题目链接：太阳轰炸

在这里插入图片描述

答案对质数109+7 取模的定义：设 M=109+7，可以证明答案可表示为一个既约分数 p/q，其中 , 均为整数且 q 模 M 不余 0。输出满足 0≤x<M 且x⋅q≡p(mod M) 的整数 x。

题目分析：

$min(1,{(r1+r)\over{r2^2}}),需要砸至少[h/a]次（向上取整）$

1、 $令 t = [h / a] ，当 t > n 时，不可能摧毁，概率为 0$

2、 $当 r 2 < = r 1 + r 时，伤害范围一定与虚空碎片相交，概率为 1$

3、根据二项分布，计算所有可能事件的概率和

注意：

用费马小定理求逆元
预处理阶乘，逆元
乘法、加法要记得取模
数组多了不要开long long，会超出空间限制

#include<cstring>
#include<cstdio>
#include<iostream>
#include<algorithm>
 
using namespace std;
 
const int N = 5e6 + 10;
const int mod = 1e9 + 7;
 
typedef long long LL;
 
LL n,r1,r2,r,a,h;
int fact[N],inv_fact[N];
int fact_p[N],fact_q[N];
 
int qmi(int a,int b){
    int res = 1 % mod;
    while(b){
        if(b & 1) res = 1LL * res * a % mod;
        a = 1LL * a * a % mod;
        b >>= 1;
    }
    return res;
}
 
int main(){
    cin >> n >> r1 >> r2 >> r >> a >> h;
    LL t = (h + a - 1) / a;
    if(t > n){
        puts("0");
        return 0;
    }
    if(r2 <= r1 + r){
        puts("1");
        return 0;
    }
    LL p = (r1 + r) * (r1 + r) % mod;
    LL q = (r2 * r2 - (r1 + r) * (r1 + r)) % mod;
    LL R2 = r2 * r2 % mod;
    fact[0] = 1;
    for(int i = 1; i < N; i ++ ) {
        fact[i] = 1LL * fact[i - 1] * i % mod;
    }
     
    inv_fact[n] = qmi(fact[n],mod - 2);
    for(int i = n - 1; i >= 0; i --){
        inv_fact[i] = 1LL * inv_fact[i + 1] * (i + 1) % mod;
    }
     
    R2 = qmi(R2,mod - 2);
    fact_p[0] = fact_q[0] = 1;
    for(int i = 1; i <= n; i ++ ){
        fact_p[i] = 1LL * fact_p[i - 1] * p % mod * R2 % mod;
        fact_q[i] = 1LL * fact_q[i - 1] * q % mod * R2 % mod;
    }
    LL ans = 0;
    for(int i = t; i <= n; i ++ ){
        ans = (1LL * fact[n] % mod * inv_fact[i] % mod * inv_fact[n - i] % mod * fact_p[i] % mod * fact_q[n - i] % mod + ans) % mod;
    }
    printf("%lld",ans);
    return 0;
}