[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月月赛]抵制克苏恩题解的主要内容,如果未能解决你的问题,请参考以下文章

[Lydsy2017年4月月赛]抵制克苏恩

4832: [Lydsy2017年4月月赛]抵制克苏恩]解题报告

[Bzoj4832][Lydsy2017年4月月赛]抵制克苏恩 (期望dp)

BZOJ 4832 4832: [Lydsy2017年4月月赛]抵制克苏恩 (期望DP)

Lydsy2017年4月月赛 抵制克苏恩

BZOJ4832[Lydsy2017年4月月赛]抵制克苏恩 概率与期望