HDU 3076 ssworld VS DDD 概率dp,无穷级数,oj错误题目 难度:2

http://acm.hdu.edu.cn/showproblem.php?pid=3076

不可思议的题目,总之血量越少胜率越高,所以读取时把两人的血量交换一下

明显每一轮的胜率和负率都是固定的,所以设psc为胜率,pls为负率,peq为平率,

则在每一局中的胜率负率平率可以确定,

而在有结果的一个阶段中的胜率和负率则各是一个无穷级数

psc(new)=1*psc+peq*psc+peq*peq*psc.......=lim(n->正无穷)(1-peq^n)*psc/(1-peq)=psc/(1-peq)

pls(new)=1*pls+peq*pls+peq*peq*pls.......=lim(n->正无穷)(1-peq^n)*pls/(1-peq)=pls/(1-peq)

设p[i][j]为A的hp为i,B的hp为j的概率

很明显转移公式是:

当i>0,j>0时

p[i-1][j]+=p[i][j]*pls(new)

p[i][j-1]+=p[i][j]*psc(new)

#include <cstdio>
#include <cstring>
#include <cmath>
#include <algorithm>
#include <queue>
using namespace std;
const int maxn=2e3+3;
double p[maxn];
double pp[2][6];
double psc,pls,peq;
int a,b;
double ans;

void calc(){
    psc=pls=peq=0;
    ans=0;
    memset(p,0,sizeof p);
    for(int i=0;i<6;i++){
        for(int j=0;j<6;j++){
            if(i==j)peq+=pp[0][i]*pp[1][j];
            else if(i<j)pls+=pp[0][i]*pp[1][j];
            else psc+=pp[0][i]*pp[1][j];
        }
    }
    if(peq!=1){
        psc/=(1-peq);
        pls/=(1-peq);
    }
    p[b]=1;
    for(int i=a;i>0;i--){
        for(int j=b;j>0;j--){
            p[j-1]+=p[j]*psc;
        }
        ans+=p[1]*psc;
        for(int j=b;j>0;j--){
            p[j]*=pls;
        }
    }
}

int main(){
    while(scanf("%d%d",&b,&a)==2){
        for(int i=0;i<6;i++)scanf("%lf",pp[0]+i);
        for(int i=0;i<6;i++)scanf("%lf",pp[1]+i);
        calc();
        printf("%.6f
",ans);
    }

    return 0;
}