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\n",ans);

    }

    return 0;

}