【codeforces 175D】 Plane of Tanks: Duel

http://codeforces.com/problemset/problem/175/D (题目连接)

题意

　　A,B两人玩坦克大战，坦克有生命值，射击间隔，伤害范围，未命中的几率。问A赢的几率是多少。

Solution

　　直接作并很差作，注意到精度要求只有$10^{-4}$，也就是说当射击次数达到必定上限以后，知足了精度要求咱们就能够不作了。咱们考虑按照时间dp。怎么设计状态呢，若是$f[i][j]$表示A还剩$i$点生命值，B还剩$j$点生命值，这样复杂度过高了，咱们考虑将A和B剩余生命分开dp。

　　$f[i][j]$表示A被射击$i$次后，还剩下$j$点生命值的几率，特殊的，$f[i][0]$表示在射击次数$<=i$时A死亡的几率。B的同理。转移很显然：$$f[i][max(0,j-x)]=f[i-1][j]*(1-p)/(r-l+1)+f[i-1][j]*p$$

　　咱们预处理出$D$次射击后A,B的dp数组。考虑怎么计算答案。若是B在第$K$轮挂掉了，那么要求A在第$K$轮射击以前没有死，那么这个时候A被射击了多少次呢，咱们能够算出来：$$T=((K-1)*dt_A+dt_B-1)/dt_B$$

　　那么A在第$K$轮胜出的几率为：$$ans_K=(1-f_A[T][0])*f_B[K][0]$$

　　最后把全部的$ans$累加起来就能够了。

细节

　　mdzz写了半天，各类坑点，本机1s多CF上竟然TLE了是什么鬼嘛，智力-2。

　　时间设到4000差很少够了。

 　　当$f[i][j]==0$的时候不dp，可是！可是！！要这样写：$f[i][j]<=eps$。。一把辛酸泪啊T_T

　　要特判$p_A=100$和$p_B=100$的状况，并且前者必定要写在后者以前T_T。

代码

// codeforces175D
#include<algorithm>
#include<iostream>
#include<cstdlib>
#include<cstring>
#include<cstdio>
#include<cmath>
#include<queue>
#include<ctime>
#define LL long long
#define eps 1e-8
#define Pi acos(-1.0)
#define free(a) freopen(a".in","r",stdin),freopen(a".out","w",stdout);
using namespace std;

double f[2][4010][210],p[2],ans;
int hp[2],dt[2],L[2],R[2],P[2],n,x;

int main() {
	for (int i=0;i<2;i++) scanf("%d%d%d%d%d",&hp[i],&dt[i],&L[i],&R[i],&P[i]);
	if (P[0]==100) {puts("0");return 0;}  //must first
	if (P[1]==100) {puts("1");return 0;}
	f[0][0][hp[0]]=f[1][0][hp[1]]=1;p[0]=P[0]/100.0;p[1]=P[1]/100.0;
	for (x=1;x<=4000;x++) {
		for (int i=0;i<2;i++)
			for (int l=i^1,j=1;j<=hp[i];j++) {
				if (f[i][x-1][j]<=eps) continue;
				for (int k=L[l];k<=R[l];k++) 
					f[i][x][max(0,j-k)]+=f[i][x-1][j]*(1-p[l])/(R[l]-L[l]+1);
				f[i][x][j]+=f[i][x-1][j]*p[l];
			}
		f[0][x][0]+=f[0][x-1][0];
	}
	for (int i=1;i<=x;i++) {
		int T=((i-1)*dt[0]+dt[1]-1)/dt[1];
		if (T>x) break;
		ans+=(1-f[0][T][0])*f[1][i][0];
	}
	printf("%.6lf\n",ans);
	return 0;
}