CodeForces 453 A.Little Pony and Expected Maximum（快速幂）

Description

掷一个$m$面的骰子（设$m$个面上的点数分别是$1,2,...,m$,且掷到每一个面的概率相同均为$\frac{1}{m}$）$n$次，问$n$次中最大点数的期望值

Input

两个整数$m$和$n$ $(1\leq m,n\leq 10^5)$

Output

输出问$n$次中最大点数的期望值，结果与精确值的绝对误差和相对误差均需小于$10^{-4}$

Sample Input

6 1

Sample Output

3.500000000000

Solution

点数最大值不小于$i$的概率为$1-(\frac{i-1}{m})^n,i=1,...,m$

点数最大值为$i$的概率=点数最大值不小于$i$的概率-点数最大值不小于$i+1$的概率=$1-(\frac{i-1}{m})^n-(1-(\frac{i}{m})^n)=(\frac{i}{m})^n-(\frac{i-1}{m})^n$

点数最大值的期望$ans=\sum\limits_{i=1}^mi\cdot((\frac{i}{m})^n-(\frac{i-1}{m})^n)=m-\sum\limits_{i=1}^{m-1}(\frac{i}{m})^n$

用快速幂加速求$(\frac{i}{m})^n$的过程，时间复杂度$O(mlogn)$

Code

#include<cstdio>
#include<iostream>
#include<cstring>
#include<algorithm>
#include<cmath>
#include<vector>
#include<queue>
#include<map>
#include<set>
#include<ctime>
using namespace std;
typedef long long ll;
#define INF 0x3f3f3f3f
double Pow(double a,int n)
{
    double ans=1;
    while(n)
    {
        if(n&1)ans=ans*a;
        a=a*a;
        n>>=1;
    }
    return ans;
}
int main()
{
    int n,m;
    while(~scanf("%d%d",&m,&n))
    {
        double ans=m;
        for(int i=1;i<m;i++)ans-=Pow(1.0*i/m,n);
        printf("%.12f\n",ans);
    }
    return 0;
}