UVA10529 Dumb Bones -（数学期望，DP）

额，看见很多大佬写了DP，然后这里提供一个记忆化搜索的方法

最开始的想法

对于一段区间[i，j]，往里面插入第$j-i+2$块骨牌，最后复杂度是$O(n^3)$，转移方程：

$f[i][j]=min(f[i][k-1]\times \frac{1-pr}{1-pl-pr}+f[k+1][j]\times \frac{1-pl}{1-pl-pr}+\frac{1}{1-pl-pr})(i<k<j)$

如何优化

得到第一个$O(n^3)$的做法之后，我们发现其实i到j并不重要，只需考虑i到j的骨牌数量就好了，复杂度为$O(n^2)$，转移方程：

$f[i]=min(ans,f[j]\times \frac{1-pr}{1-pl-pr}+f[i-j-1]\times \frac{1-pl}{1-pl-pr}+\frac{1}{1-pl-pr})(0\le j\le i)$

再优化

利用决策单调性，优化成$O(n)$，由于这题n只有1000，这里就不过多赘述（其实是我也不会呀QAQ）

放代码

//DP代码
#include<bits/stdc++.h>
using namespace std;
int n;
double pl,pr,f[1005];
int main(){
	do{
		memset(f,0,sizeof(f));
	    scanf("%d",&n);
	    if(n==0) break;
		scanf("%lf%lf",&pl,&pr);
		f[1]=1/(1-pl-pr);
		for(int i=2;i<=n;i++){
			f[i]=(1-pl)/(1-pl-pr)*f[i-1]+f[1];
		    for(int j=1;j<i;j++)
		        f[i]=min(f[i],(1-pr)/(1-pl-pr)*f[j]+(1-pl)/(1-pl-pr)*f[i-j-1]+f[1]);
		}
		printf("%.2f\n",f[n]);
	}while(1);
	return 0;
}

//记忆化搜索代码
#include<bits/stdc++.h>
using namespace std;
int n;
double pl,pr,f[1005];
double calc(int x){
	if(x==1) return 1/(1-pl-pr);
	if(x==0) return 0;
	if(f[x]) return f[x];
	double ans=(1-pl)/(1-pl-pr)*calc(x-1)+1/(1-pl-pr);
	for(int i=1;i<x;i++)
	    ans=min(ans,(1-pr)/(1-pl-pr)*calc(i)+(1-pl)/(1-pl-pr)*calc(x-i-1)+1/(1-pl-pr));
	f[x]=ans;
	return ans;
}
int main(){
	do{
		memset(f,0,sizeof(f));
	    scanf("%d",&n);
	    if(n==0) break;
		scanf("%lf%lf",&pl,&pr);
		printf("%.2f\n",calc(n));
	}while(1);
	return 0;
}