BZOJ4832抵制克苏恩

我们是在从后往前推 ,即我们是在用当前局推上一局

i:表示还有i次没打，a:表示上一局血量为1的还有多少个，b:2，c:3

f[i]中a,b,c考虑这一局和上一局比a,b,c的变化

double k=1/(1+a+b+c) //我们要转移状态选择每一个人物的概率（因为有一个英雄，所以+1）
int tot=a+b+c;
f[i+1][a][b][c]+=(f[i][a][b][c]+1)*k;//攻击英雄
if(a) f[i+1][a][b][c]+=f[i][a-1][b][c]*a*k;//攻击血量为1的随从，死亡
if(b)
{
    if(tot<7) f[i+1][a][b][c]+=f[i][a+1][b-1][c+1]*b*k;//可以增加奴隶主 
    else f[i+1][a][b][c]+=f[i][a+1][b-1][c]*b*k;
} 
if(c)
{
    if(tot<7) f[i+1][a][b][c]+=f[i][a][b+1][c]*c*k;
    else f[i+1][a][b][c]+=f[i][a][b+1][c-1]*c*k;
}

这就是期望dp的核心

嗯，也没什么好说的了，代码：

#include<cstdio>
#include<iostream>
using namespace std;
double f[60][10][10][10];
int t;
int read()
{
    int num=0,f=1;
    char ch=getchar();
    while(!isdigit(ch)) {if(ch=='-') f=-1; ch=getchar();}
    while(isdigit(ch)) num=(num<<3)+(num<<1)+(ch^48),ch=getchar();
    return num*f;
}
void work()
{
    for(int i=0;i<50;i++)
        for(int a=0;a<=7;a++)
            for(int b=0;b<=7-a;b++)
                for(int c=0;c<=7-a-b;c++)
                {
                    double k=1.0/(1+a+b+c);
                    int tot=a+b+c;
                    f[i+1][a][b][c]+=(f[i][a][b][c]+1)*k;
                    if(a) f[i+1][a][b][c]+=f[i][a-1][b][c]*a*k;
                    if(b)
                    {
                        if(tot<7) f[i+1][a][b][c]+=f[i][a+1][b-1][c+1]*b*k; 
                        else f[i+1][a][b][c]+=f[i][a+1][b-1][c]*b*k;
                    } 
                    if(c)
                    {
                        if(tot<7) f[i+1][a][b][c]+=f[i][a][b+1][c]*c*k;
                        else f[i+1][a][b][c]+=f[i][a][b+1][c-1]*c*k;
                    }
                }
}
int main()
{
    work();
    t=read();
    while(t--)
    {
        int k=read(),a=read(),b=read(),c=read();
        printf("%.2lf\n",f[k][a][b][c]);
    }
    return 0;
}