本文主要是介绍BZOJ 4832 [Lydsy2017年4月月赛]抵制克苏恩 期望dp，希望对大家解决编程问题提供一定的参考价值，需要的开发者们随着小编来一起学习吧！
Description
小Q同学现在沉迷炉石传说不能自拔。他发现一张名为克苏恩的牌很不公平。如果你不玩炉石传说，不必担心，小Q
同学会告诉你所有相关的细节。炉石传说是这样的一个游戏，每个玩家拥有一个 30 点血量的英雄，并且可以用牌
召唤至多 7 个随从帮助玩家攻击对手，其中每个随从也拥有自己的血量和攻击力。小Q同学有很多次游戏失败都是
因为对手使用了克苏恩这张牌，所以他想找到一些方法来抵御克苏恩。他去求助职业炉石传说玩家椎名真白，真白
告诉他使用奴隶主这张牌就可以啦。如果你不明白我上面在说什么，不必担心，小Q同学会告诉你他想让你做什么
。现在小Q同学会给出克苏恩的攻击力是 K ，表示克苏恩会攻击 K 次，每次会从对方场上的英雄和随从中随机选
择一个并对其产生 1 点伤害。现在对方有一名克苏恩，你有一些奴隶主作为随从，每名奴隶主的血量是给定的。
如果克苏恩攻击了你的一名奴隶主，那么这名奴隶主的血量会减少 1 点，当其血量小于等于 0 时会死亡，如果受
到攻击后不死亡，并且你的随从数量没有达到 7 ，这名奴隶主会召唤一个拥有 3 点血量的新奴隶主作为你的随从
；如果克苏恩攻击了你的英雄，你的英雄会记录受到 1 点伤害。你应该注意到了，每当克苏恩进行一次攻击，你
场上的随从可能发生很大的变化。小Q同学为你假设了克苏恩的攻击力，你场上分别有 1 点、 2 点、 3 点血量的
奴隶主数量，你可以计算出你的英雄受到的总伤害的期望值是多少吗？
Input
输入包含多局游戏。
第一行包含一个整数 T (T<100) ，表示游戏的局数。
每局游戏仅占一行，包含四个非负整数 K, A, B 和 C ，表示克苏恩的攻击力是 K ，你有 A 个 1 点血量的奴隶
主， B 个 2 点血量的奴隶主， C 个 3 点血量的奴隶主。
保证 K 是小于 50 的正数， A+B+C 不超过 7 。
Output
对于每局游戏，输出一个数字表示总伤害的期望值，保留两位小数。
Sample Input
1
1 1 1 1
1 1 1 1
Sample Output
0.25
HINT
传送门
炉石传说骨灰级玩家2333
虽说一眼能大致看出期望dp以及状态，但我对期望dp的理解真的不够啊=c=
大概的感觉就是f[i][a][b][c]表示攻击了i点，1血奴隶主a个，2血奴隶主b个，3血奴隶主c个，
这个东西的概率
不怎么会直接用一个f……那么再用一个E[i][a][b][c]表示……的期望，
能够知道，很明显地，一血，二血，三血奴隶主不管有几个效果都是一样的，
那么根据打到哪个单位的概率，
能够把f转移掉了。。具体不说了自己yy吧。。
然后同时把E一起转移了，只要在打到人的时候加上1*f[i][a][b][c]就行了。。
期望dp真的不会啊。。。
#include<bits/stdc++.h>
using namespace std;
int K,A,B,C;
double f[55][8][8][8],E[55][8][8][8];
void solve(){memset(f,0,sizeof(f));memset(E,0,sizeof(E));scanf("%d%d%d%d",&K,&A,&B,&C);f[0][A][B][C]=1.0;for (int i=0;i<K;i++)for (int a=0;a<=7;a++)for (int b=0;b<=7;b++)for (int c=0;c<=7;c++){if (a+b+c>7) break;if (f[i][a][b][c]==0) continue;if (a) f[i+1][a-1][b][c]+=f[i][a][b][c]*a/(double)(1+a+b+c),E[i+1][a-1][b][c]+=E[i][a][b][c]*a/(double)(1+a+b+c);if (a+b+c<=6){if (b) f[i+1][a+1][b-1][c+1]+=f[i][a][b][c]*b/(double)(1+a+b+c),E[i+1][a+1][b-1][c+1]+=E[i][a][b][c]*b/(double)(1+a+b+c);if (c) f[i+1][a][b+1][c]+=f[i][a][b][c]*c/(double)(1+a+b+c),E[i+1][a][b+1][c]+=E[i][a][b][c]*c/(double)(1+a+b+c);}if (a+b+c==7){if (b) f[i+1][a+1][b-1][c]+=f[i][a][b][c]*b/(double)(1+a+b+c),E[i+1][a+1][b-1][c]+=E[i][a][b][c]*b/(double)(1+a+b+c);if (c) f[i+1][a][b+1][c-1]+=f[i][a][b][c]*c/(double)(1+a+b+c),E[i+1][a][b+1][c-1]+=E[i][a][b][c]*c/(double)(1+a+b+c);}f[i+1][a][b][c]+=f[i][a][b][c]/(double)(1+a+b+c);E[i+1][a][b][c]+=(E[i][a][b][c]+f[i][a][b][c])/(double)(1+a+b+c);}double ans=0.0;for (int a=0;a<=7;a++)for (int b=0;b<=7;b++)for (int c=0;c<=7;c++)if (a+b+c>7) break;else ans+=E[K][a][b][c];printf("%.2lf\n",ans);
}
int main(){int cas;scanf("%d",&cas);while (cas--) solve();return 0;
}
这篇关于BZOJ 4832 [Lydsy2017年4月月赛]抵制克苏恩 期望dp的文章就介绍到这儿，希望我们推荐的文章对编程师们有所帮助！
