[Lydsy2017年4月月赛]抵制克苏恩题解
Posted Hzoi_joker
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了[Lydsy2017年4月月赛]抵制克苏恩题解相关的知识,希望对你有一定的参考价值。
考试的时候以为就是简单的概率期望题,考完后知道是简单的概率期望DP题,完美爆零。
这道题数据范围很小,很容易让人想到状压,不过貌似没什么可压的。那么只能说明这道题复杂度很高了,状态数组f[o][i][j][k]为第o次攻击后出现有i个1滴血的j个2滴血的,k个3滴血的情况的概率。那么转移方程就明了了。
1.f[o][i][j][k]+=f[o-1][i][j][k]/(i+j+k+1)英雄收到攻击
2.f[o][i][j][k]+=f[o-1][i+1][j][k]*(i+1)/(i+j+k+2)有一个1滴血的奴隶主被击杀。前提是i+j+k+1<=7
3.f[o][i][j][k]+=f[o-1][i-1][j+1][k-1]*(j+1)/(i+j+k)有一个2滴血的奴隶主收到攻击,召唤一个奴隶主,前提是i!=0&&k!=0
4.f[[o][i][j][k]+=f[o-1][i][j-1][k]*(k/i+j+k)有一个3滴血的奴隶主收到攻击,召唤一个奴隶主,前提是j!=0。
5由于仆从满7个后不会再召唤,所以当i+j+k=7时需要特判受到攻击却未召唤的情况,即f[o][i][j][k]+=f[o-1][i-1][j+1][k],f[o][i][j][k]+=f[o-1][i][j-1][k+1],前提分别为i!=0,j!=0。
转移方程写出来别挂精度上就行了,但是要注意我们的状态为受到攻击后场上局面的状态,真正结果为0~k-1中所有状态转移方程中1/仆从+1,一开始写的为1~k身败名裂……
1 #include<iostream> 2 #include<cstdlib> 3 #include<cstdio> 4 #include<cstring> 5 #include<queue> 6 #include<algorithm> 7 #include<cmath> 8 using namespace std; 9 int t,n; 10 double f[55][10][10][10]; 11 int main(){ 12 // freopen("defcthun.in","r",stdin); 13 // freopen("defcthun.out","w",stdout); 14 scanf("%d",&t); 15 while(t--) 16 { 17 memset(f,0,sizeof(f)); 18 scanf("%d",&n); 19 int xx,yy,zz; 20 scanf("%d%d%d",&xx,&yy,&zz); 21 f[0][xx][yy][zz]=1.0; 22 for(int o=1;o<n;o++) 23 { 24 for(int i=0;i<=7;i++) 25 { 26 for(int j=0;j<=7;j++) 27 { 28 for(int k=0;k<=7;k++) 29 { 30 if(i+j+k>7)break; 31 if(i+j+k+1<=7) f[o][i][j][k]+=f[o-1][i+1][j][k]*(double(i+1)/double(i+1+j+k+1)); 32 f[o][i][j][k]+=f[o-1][i][j][k]/double(i+j+k+1); 33 if(i&&k) f[o][i][j][k]+=f[o-1][i-1][j+1][k-1]*(double(j+1)/double(i+k+j)); 34 if(j) f[o][i][j][k]+=f[o-1][i][j-1][k]*(double(k)/double(i+j+k)); 35 if(i+j+k==7) 36 { 37 if(i) f[o][i][j][k]+=f[o-1][i-1][j+1][k]*double(j+1)/double(i+j+k+1); 38 if(j) f[o][i][j][k]+=f[o-1][i][j-1][k+1]*double(k+1)/double(i+j+k+1); 39 } 40 } 41 } 42 } 43 } 44 double sum=0.0; 45 for(int o=0;o<n;o++) 46 { 47 for(int i=0;i<=7;i++) 48 { 49 for(int j=0;j<=7;j++) 50 { 51 for(int k=0;k<=7;k++) 52 { 53 if(i+j+k>7)break; sum+=f[o][i][j][k]*1.0/double(j+k+i+1); 54 } 55 } 56 } 57 } 58 printf("%.2lf\n",sum); 59 } 60 //while(1); 61 return 0; 62 }
以上是关于[Lydsy2017年4月月赛]抵制克苏恩题解的主要内容,如果未能解决你的问题,请参考以下文章
4832: [Lydsy2017年4月月赛]抵制克苏恩]解题报告
[Bzoj4832][Lydsy2017年4月月赛]抵制克苏恩 (期望dp)