[HDU4652] Dice

扔一个有m个面的骰子，每个面有一个数字，这些数字互不相同，求连续扔到n个数相同的期望步数或连续扔到n个数不同的期望步数

相同的：

　　E[i]表示已经有i连续个相同，到达目标的期望步数

　　E[i]=E[i+1]/m+(1-1/m)E[1].................... 1

　　E[i+1]=E[i+2]/m+(1-1/m)E[1]................ 2

　　1-2得 E[i]-E[i+1]=(E[i+1]+E[i+2])/m

　　设s[i]=E[i]-E[i+1]，

　　则s[0]=E[0]-E[1]=1，s[i+1]=m*s[i]，ans=sigma(s[i])(0<=i<n)　　

不同的：

　　E[i]表示已经有i连续个不同，到达目标的期望步数

　　E[i]=simga(E[j])/m+((m-i)/m)E[i+1],(1<=j<=i).................. 1

　　E[i+1]=sigma(E[j])/m+((m-i-1)/m)E[i+2],(1<=j<=i+1)....... 2

　　1-2得 E[i]-E[i+1]=(m-i-1)/m*(E[i+1]-E[i+2])

　　设s[i]=E[i]-E[i+1]，

　　则s[0]=E[0]-E[1]=1，s[i+1]=m/(m-i-1)s[i]，ans=sigma(s[i])(0<=i<n)

#include<bits/stdc++.h>
using namespace std;
#define maxn 1000005
int T,op,n;
double m,s[maxn];
void solve0(){
    s[0]=1;
    double ans=s[0];
    for(int i=0;i<n-1;i++)
        s[i+1]=s[i]*m,ans+=s[i+1];
    printf("%.9lf\n",ans);
}
void solve1(){
    s[0]=1;
    double ans=s[0];
    for(int i=0;i<n-1;i++)
        s[i+1]=m/(m-i-1)*s[i],ans+=s[i+1];
    printf("%.9lf\n",ans);
}
int main(){
    scanf("%d",&T);
    while(T--){
        scanf("%d%lf%d",&op,&m,&n);
        if(op==0)solve0();
        else solve1();
    }
    return 0;
}