题目链接:Problem - 453A - Codeforces
简化题意: m个面的骰子,第i个面点数为i(1<=i<=m),掷到每个面的概率均为1/m。求独立掷n次的最大点数的期望。(1<=m,n<=1e5)
分析:期望就是所有可能的情况乘以其概率再求和,设P(X=i)表示独立掷n次的最大点数为i时的概率。
那么我们就是要求,我们先来看一下P(X=i)怎么求,当P(X=i)不好直接求的时候我们可以考虑求P(X<=i),然后由P(X=i)=P(X<=i)-P(X<=i-1)求出P(X=i)。很显然有P(X<=i)=(i/m)^n,那么对应的就有P(X=i)=( i / m ) ^ n - ( ( i - 1 ) / m ) ^ n
最后我们只需要for循环遍历一边就可以求出答案。
下面是代码:
#include<iostream>
#include<algorithm>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<vector>
#include<queue>
using namespace std;
const int N=1e5+10;
double qpow(double a,long long b)
{
if(fabs(a-0.0)<1e-6) return 0;
double ans=1.0;
while(b)
{
if(b&1) ans*=a;
a*=a;
b>>=1;
}
return ans;
}
int main()
{
int m,n;
cin>>m>>n;
double ans=0;
for(int i=1;i<=m;i++)
ans+=i*(qpow(i*1.0/m,n)-qpow((i-1)*1.0/m,n));
printf("%.10lf",ans);
return 0;
}