期望
题目大意: 有一个 m m m面的骰子,一共投 n n n次。设投 n n n次的权值为这 n n n次出现的最大值,求权值的期望。
首先答案是 ∑ i ∗ p ( i ) \sum i*p(i) ∑i∗p(i)的式子,主要问题就是算 p ( i ) p(i) p(i)。而 p ( i ) = ( i m ) n − ( i − 1 m ) n p(i)=(\frac im)^n-(\frac {i-1}m)^n p(i)=(mi)n−(mi−1)n。前一项表示这 n n n次中投的点数都不超过 i i i的概率,后一个为 i − 1 i-1 i−1次的概率。两个相减就表示投 n n n次至少出现一次 i i i且都不大于 i i i的概率,快速幂算一算就好了。
代码:
#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
typedef double DB;
int n,m; DB ans;
DB ksm(DB a,int b){ DB ret=1; for (;b;a=a*a,b>>=1) if (b&1) ret=ret*a; return ret; }
int main(){
scanf("%d%d",&m,&n);
for (int i=1;i<=m;i++)
ans+=(ksm(1.0*i/m,n)-ksm(1.0*(i-1)/m,n))*i;
return printf("%.10lf\n",ans),0;
}