51nod1836-战忽局的手段【期望dp,矩阵乘法】

1|0正题

题目连接:http://www.51nod.com/Challenge/Problem.html#problemId=1836

---

1|1题目大意

$n$个点$m$次随机选择一个点标记（可以重复），求最后被标记点的期望个数。

$1\leq n,m\leq 10^{18}$

---

1|2解题思路

额开始拿方案数推了半天后面发现要斯特林数就放弃了，然后换了种方法发现很简单？

设$i$轮之后被标记点的期望个数是$f_i$，那么有

$$f_i=f_{i-1}+\frac{n-f_{i-1}}{n}$$

$$f_i=f_{i-1}\frac{n-1}{n}+1$$

然后矩阵乘法就好了。

有一说一我第一次用期望值来算概率（（（

时间复杂度$O(T\log n)$

---

1|3code

#include<cstdio>
#include<cstring>
#include<algorithm>
using namespace std;
const int S=2;
struct Matrix{
	__float128 a[S][S];
}f,ans,c;
long long T,n,m;
Matrix operator*(const Matrix &a,const Matrix &b){
	c.a[0][0]=c.a[0][1]=c.a[1][0]=c.a[1][1]=0;
	for(int i=0;i<S;i++)
		for(int j=0;j<S;j++)
			for(int k=0;k<S;k++)
				c.a[i][j]+=a.a[i][k]*b.a[k][j];
	return c;
}
int main()
{
	scanf("%lld",&T);
	while(T--){
		scanf("%lld%lld",&n,&m);
		f.a[1][1]=(__float128)(n-1)/n;
		f.a[0][1]=f.a[0][0]=1;f.a[1][0]=0;
		ans.a[0][0]=1;ans.a[0][1]=0;
		while(m){
			if(m&1)ans=ans*f;
			f=f*f;m>>=1;
		}
		printf("%.12lf\n",(double)ans.a[0][1]);
	}
	return 0;
}

__EOF__

本文作者：QuantAsk
本文链接：https://www.cnblogs.com/QuantAsk/p/15087088.html
关于博主：退役OIer，GD划水选手
版权声明：本博客所有文章除特别声明外，均采用 BY-NC-SA 许可协议。转载请注明出处！
声援博主：如果您觉得文章对您有帮助，可以点击文章右下角【推荐】一下。您的鼓励是博主的最大动力！