怎么有优化这个，让他在n=1e9的时候可以再1s内通过。

double e=(m+1)/2.0;
for(long long i=2;i<=n;i++){
	e=(floor(e)*e+(ceil(e)+m)*(m-ceil(e)+1)/2.0)/m*1.0;
}

矩阵乘法好像不行