CodeForces 717 D.Dexterina’s Lab(博弈论+dp+矩阵快速幂)
来源:互联网 发布:绝地求生a卡优化补丁 编辑:程序博客网 时间:2024/04/30 07:01
Description
n堆石子,每堆石子的数量是i的概率都是dp[i],两人轮流取石子,每次可以选任一堆取走任意数量的石子,但不能不取,谁没得取谁输,问在双方采取最佳策略的前提下先手必胜的概率
Input
第一行两整数n和x分别表示石子堆数和每堆数量上限,之后x+1个实数p[i]表示一堆石子有i个的概率
(1<=n<=1e9,1<=x<=100,0<=p[i]<=1,p[0]+p[1]+…+p[x]=1)
Output
输出先手必胜的概率
Sample Input
2 2
0.500000 0.250000 0.250000
Sample Output
0.62500000
Solution
dp[i][j]表示前i堆石子数量异或和为j的概率,那么答案为1-dp[n][0],在前i-1堆石子数量异或和确定的情况下根据第i堆石子的数量得到转移方程dp[i][k^j]=sum(dp[i-1][k]*p[j]),0<=k<=127,0<=j<=x,用矩阵快速幂加速该转移即可
Code
#include<cstdio>#include<iostream>#include<cstring>#include<algorithm>#include<cmath>#include<vector>#include<queue>#include<map>#include<set>#include<ctime>using namespace std;#define maxn 133struct Mat{ double mat[maxn][maxn];//矩阵 int row,col;//矩阵行列数 };Mat Mul(Mat a,Mat b)//矩阵乘法 { Mat ans; ans.row=a.row; ans.col=b.col; memset(ans.mat,0,sizeof(ans.mat)); for(int i=0;i<ans.row;i++) for(int k=0;k<a.col;k++) for(int j=0;j<ans.col;j++) ans.mat[i][j]+=a.mat[i][k]*b.mat[k][j]; return ans;}Mat Pow(Mat a,int k)//矩阵快速幂 { Mat ans; ans.row=a.row; ans.col=a.col; for(int i=0;i<a.row;i++) for(int j=0;j<a.col;j++) ans.mat[i][j]=(i==j); while(k) { if(k&1)ans=Mul(ans,a); a=Mul(a,a); k>>=1; } return ans;}Mat A;int n,x;double P[maxn];int main(){ while(~scanf("%d%d",&n,&x)) { for(int i=0;i<=x;i++)scanf("%lf",&P[i]); A.row=A.col=128; for(int i=0;i<128;i++) for(int j=0;j<=x;j++) A.mat[i^j][i]=P[j]; A=Pow(A,n); printf("%.10f\n",1.0-A.mat[0][0]); } return 0;}
0 0
- CodeForces 717 D.Dexterina’s Lab(博弈论+dp+矩阵快速幂)
- [矩阵快速幂加速DP] Codeforces 717D Bubble Cup 9 - Finals D. Dexterina’s Lab
- [矩阵快速幂加速DP]Codeforces 717D Bubble Cup 9 -Finals D. Dexterina’s Lab
- codeforces 717 D. Dexterina’s Lab(矩阵快速幂加速dp,好题)
- 【矩阵乘法优化DP】Codeforces 717D Dexterina’s Lab
- 【Nim游戏+DP+矩阵乘法】Codeforces717D[Dexterina’s Lab]题解
- CodeForces 696D AC自动机+DP+矩阵快速幂
- CodeForces 551 D.GukiZ and Binary Operations(dp+矩阵快速幂)
- Codeforces 696D Legen...(AC自动机+矩阵快速幂)
- codeforces 696D Legen...(AC自动机+矩阵快速幂)
- HDU5863 cjj's string game (dp+矩阵快速幂)
- HDU5863cjj's string game(DP+矩阵快速幂)
- CodeForces 678D Iterated Linear Function 矩阵快速幂
- codeforces 118D Caesar's Legions(DP进阶)
- Codeforces 118D Caesar's Legions (dp)
- [CodeForces 118D] Caesar's Legions (计数DP)
- codeforces 118-D. Caesar's Legions(dp计数)
- 矩阵快速幂DP Darth Vader and Tree : CodeForces
- css学习交流(一)
- 2017年京东实习编程题 ----Java
- 把程序和大自然进行类比
- postgresql copy的使用一例
- 删除排序链表中的重复元素
- CodeForces 717 D.Dexterina’s Lab(博弈论+dp+矩阵快速幂)
- matlab 2015 ubuntu 14 15 16 崩溃 创建java对象找不到类
- 1000 A+B Problem
- 【web性能优化】使用SVG中的Symbol元素制作Icon,解决图片分辨率适配问题
- 定时清理Oracle数据库中不活动的session的脚本,解决ora-12516问题
- Unity3D学习笔记(3)——Unity脚本编译顺序
- HTML(1)——body和title
- C++中的基本控制结构
- StarUML类图创建教程