书 题解

book.in/.out

Hazel有n本书,编号1为n到 ,叠成一堆。当她每次抽出一本书的时候,上方的书会因重力而下落,这本被取出的书则会被放置在书堆顶。

每次有pi的概率抽取编号为i的书。她每次抽书所消耗的体力与这本书在这堆中是第几本成正比。具体地,抽取堆顶的书所耗费体力值为1 ,抽取第二本耗费体力值为2 ,以此类推。

现在 想知道,在很久很久以后(可以认为几乎是无穷的),她每次抽书所耗费的体力的期望值是多少。

最终的答案显然可以表示成a/b的形式,请输出a*(b^-1)模1e9+7的值。

【输入格式】

第一行一个整数n

接下来n行,每行两个整数ai,bi,代表抽取第i本书的概率是ai/bi

保证所有书的概率和等于1

【输出格式】

输出一行一个整数,代表期望值

【输入样例1】

2

227494 333333

105839 333333

【输出样例1】

432679642

【输入样例2】

10

159073 999999

1493 142857

3422 333333

4945 37037

2227 111111

196276 999999

190882 999999

142721 999999

34858 999999

101914 999999

【输出样例2】

871435606

【数据规模与约定】

对于30%的数据,1<=n<=10。

对于100%的数据,1<=n<=1000,0<=ai<=bi,bi!=0。

solution:

  一碰到概率就蒙,这个玩意真的是........,用E[x]表示时间X期望花费的体力值,P[x]表示事件X发生的概率。

  Ans=∑pi*E[抽第i本书]

      =∑pi*(1+∑P[j最后一次被抽到的时间比i的晚])

      =∑ pi*(1+∑pj/(pi+pj))

 1 #include<iostream>
 2 #include<cstdio>
 3 #include<cstring>
 4 #include<algorithm>
 5 using namespace std;
 6 #define int long long
 7 #define mod 1000000007
 8 int read() {
 9     int s=0,f=1;
10     char ch=getchar();
11     while(ch>'9'||ch<'0') {
12         if(ch=='-') {
13             f=-1;
14         }
15         ch=getchar();
16     }
17     while(ch>='0'&&ch<='9') {
18         s=(s<<1)+(s<<3)+(ch^48);
19         ch=getchar();
20     }
21     return s*f;
22 }
23 long long p[100005];
24 int n;
25 int mul(int x) {
26     int ans=1,t=mod-2;
27     for(; t; t>>=1,x=(x*x)%mod) {
28         if(t&1) {
29             ans=(ans*x)%mod;
30         }
31     }
32     return ans;
33 }
34 signed main() {
35     n=read();
36     for(int i=1; i<=n; i++) {
37         int x=read(),y=read();
38         p[i]=x*mul(y)%mod;
39     }
40     int ans=0;
41     for(int i=1; i<=n; i++) {
42         int sum=1;
43         for(int j=1; j<=n; j++) {
44             if(i!=j) {
45                 sum=(sum+p[j]*(mul(p[i]+p[j])%mod)%mod)%mod;
46             }
47         }
48         ans=(ans+p[i]*sum%mod)%mod;
49     }
50     printf("%lld",ans%mod);
51     return 0;
52 }

 

posted @ 2017-08-07 19:50  Forever_goodboy  阅读(272)  评论(0编辑  收藏  举报