链接:https://ac.nowcoder.com/acm/contest/275/A
题目描述
你在一栋楼房下面,楼房一共有n层,第i层每秒有p
i的概率会扔下一个东西并砸到你
求第一秒内你被砸到的概率
求第一秒内你被砸到的概率
输入描述:
第一行一个整数ni
之后有n行,第i+1行有两个整数a
,bi
,表示
![](https://i-blog.csdnimg.cn/blog_migrate/3049abc2a5fada4d1765005c93d2a1c0.png)
输出描述:
设答案为
![](https://i-blog.csdnimg.cn/blog_migrate/0ee16c0d8cd722cfbba4a2b26a755879.png)
,你只需要找到一个最小的非负整数T,使得
![](https://i-blog.csdnimg.cn/blog_migrate/c6dbde14810dc79b6e8d139c8e72d127.png)
输出这个T就行了
示例1
备注:
数据范围5
0 ≤ n ≤ 10
i
1 ≤ a
≤ bi
≤ 105
不愧是签到题,思路一看就知道,但是非得加一个乘法逆元,于是签不了了。
百度乘法逆元代码:
inline long long extend_gcd(long long a,long long b,long long &x,long long &y)
{
if(a==0&&b==0)
return -1ll;
if(b==0)
{
x=1ll;
y=0ll;
return a;
}
long long d=extend_gcd(b,a%b,y,x);
y-=a/b*x;
return d;
}
inline long long mod_reverse(long long a,long long n)
{
long long x,y,d=extend_gcd(a,n,x,y);
if(d==1)
{
if(x%n<=0)return x%n+n;
else return x%n;
}
else
return -1ll;
}
于是摘抄就过了;
#include<iostream>
long long mod=1000000007;
using namespace std;
inline long long extend_gcd(long long a,long long b,long long &x,long long &y)
{
if(a==0&&b==0)
return -1ll;
if(b==0)
{
x=1ll;
y=0ll;
return a;
}
long long d=extend_gcd(b,a%b,y,x);
y-=a/b*x;
return d;
}
inline long long mod_reverse(long long a,long long n)
{
long long x,y,d=extend_gcd(a,n,x,y);
if(d==1)
{
if(x%n<=0)return x%n+n;
else return x%n;
}
else
return -1ll;
}
long long GCD(long long a,long long b)
{
long long c=1;
while(c)
{
c=a%b;
a=b;
b=c;
}
return a;
}
int main()
{
int n;
cin>>n;
long long a,b,c,sum_a=1,sum_b=1;
for(int i=0;i<n;i++)
{
cin>>a>>b;
// c=GCD(a,b);
// a/=c;
// b/=c; //要不要没多大影响
sum_a*=b-a;
sum_b*=b;
sum_b%=mod;
sum_a%=mod; //避免超内存(浮点错误)
}
c=GCD(sum_a,sum_b);
sum_a/=c;
sum_b/=c;
sum_a=(sum_b-sum_a+mod)%mod;//由于取余过所以sum_b不一定大于sum_a
c=mod_reverse(sum_b,mod);
cout<<c*sum_a%mod;
}