CDOJ 1330 柱爷与远古法阵(高斯消元)
CDOJ 1330 柱爷与远古法阵(高斯消元)
柱爷与远古法阵
Time Limit: 125/125MS (Java/Others) Memory Limit: 240000/240000KB (Java/Others)
众所周知,柱爷的数学非常好,尤其擅长概率论!
某日柱爷在喵哈哈村散步,无意间踏入了远古法阵!
法阵很奇怪,是一个长度为NN的走廊,初始时柱爷在最左边,现在柱爷要到最右边去!
柱爷的行动方式如下:
每个回合柱爷会投一次骰子,根据骰子上的点数每个回合柱爷会投一次骰子,根据骰子上的点数X,柱爷会相应的往右边移动,柱爷会相应的往右边移动X步.步.
骰子的数值是骰子的数值是1到到6,取到每面的概率相同,取到每面的概率相同
在某些位置可能有传送门,一旦柱爷在该回合结束后在这个位置上,会被强制传送到传送门的另外一边在某些位置可能有传送门,一旦柱爷在该回合结束后在这个位置上,会被强制传送到传送门的另外一边
传送门是单向的,同时每个位置不会有超过1个传送门,同时不会存在a→b,b→c这种情况传送门是单向的,同时每个位置不会有超过1个传送门,同时不会存在a→b,b→c这种情况
在任意时刻柱爷都必须保证在法阵内,也就说如果在这一回合结束后柱爷的位置在法阵外,那么这回合柱爷将什么都不做在任意时刻柱爷都必须保证在法阵内,也就说如果在这一回合结束后柱爷的位置在法阵外,那么这回合柱爷将什么都不做
那么请问柱爷到达最右边的期望回合数是多少呢?或者是永远都无法到达?
Input
第一行两个整数NN,MM,分别表示法阵的长度和传送门的数量
接下来MM行,每行两个整数uu,vv,表示从uu到vv有一扇传送门
数据保证:
1≤N≤3001≤N≤300
0≤M≤[N−22]0≤M≤[N−22]
1<u<N,1≤v≤N,u≠v1<u<N,1≤v≤N,u≠v
Output
输出仅一行,表示期望的回合数,如果永远不能到达,输出−1−1.
答案误差在10−610−6以内将被忽略
Sample input and output
Sample Input | Sample Output |
---|---|
100 0 |
33.0476190476 |
100 2 2 3 99 100 |
29.8571428571 |
Hint
你可能需要一些概率论 & 线性代数的知识才能解决本题!
Source
1 #include <bits/stdc++.h> 2 using namespace std; 3 const int maxn=305; 4 const long double eps=1e-14; 5 long double a[maxn][maxn];//构造的高斯消元的矩阵,代表第i个方程式的第j个系数是多少 ,精度要求很高 6 int n,m,f[maxn],x,y; 7 inline int read()//读入优化 8 { 9 int x=0,f=1; 10 char ch=getchar(); 11 while(ch<'0'||ch>'9') 12 { 13 if(ch=='-') 14 f=-1; 15 ch=getchar(); 16 } 17 while(ch>='0'&&ch<='9') 18 { 19 x=x*10+ch-'0'; 20 ch=getchar(); 21 } 22 return x*f; 23 } 24 inline void write(int x)//输出优化 25 { 26 if(x<0) 27 { 28 putchar('-'); 29 x=-x; 30 } 31 if(x>9) 32 write(x/10); 33 putchar(x%10+'0'); 34 } 35 int main() 36 { 37 n=read(); 38 m=read(); 39 for(int i=1;i<=n;i++) 40 f[i]=i; 41 for(int i=1;i<=m;i++)//如果有传送的话,到哪里 42 f[read()]=read(); 43 //建立增广矩阵的过程 44 for(int i=1;i<n;i++) 45 { 46 a[i][i]=6;//第一个方程 47 if(f[i]!=i) 48 a[i][f[i]]=-6;//如果有传送门 系数直接抵消 x-y=0 相当于 x=y 49 else 50 { 51 a[i][n+1]=6;//方程右边的常数 52 for(int j=1;j<=6;j++) 53 { 54 if(i+j<=n) 55 a[i][i+j]-=1.0; 56 else 57 a[i][i]-=1.0;//另外一个方程 58 } 59 } 60 } 61 a[n][n]=1.0;//最后的方程 62 a[n][n+1]=0; 63 //高斯消元的过程 64 for(int i=1;i<=n;i++) 65 { 66 int p=i; 67 for(int j=i+1;j<=n;j++) 68 { 69 if(fabs(a[j][i])>eps)//向下查找第j个系数不为0的方程 70 p=j; 71 } 72 if(fabs(a[p][i])>eps) 73 { 74 for(int j=i;j<=n+1;j++) 75 swap(a[i][j],a[p][j]);//把方程移上来 76 for(int j=i+1;j<=n;j++)//向下消元 同时除去其他的系数 77 { 78 if(fabs(a[j][i])>eps) 79 { 80 long double k=a[j][i]/a[i][i];//消元 81 for(int t=i;t<=n+1;t++) 82 a[j][t]-=a[i][t]*k;//系数相减 83 } 84 } 85 } 86 } 87 //回代过程 88 for(int i=n;i>=1;i--) 89 { 90 for(int j=i+1;j<=n;j++) 91 { 92 if(fabs(a[i][j])>eps) 93 a[i][n+1]-=a[i][j]*a[j][n+1];//用已知的解求未知解 94 } 95 if(abs(a[i][i])<=eps&&abs(a[i][n+1])>eps)//如果出现矛盾 96 { 97 printf("-1\n"); 98 return 0; 99 } 100 a[i][n+1]/=a[i][i];//求出当前的解 101 } 102 printf("%.12lf\n",(double)a[1][n+1]);//a[i][n+1]就是第i个未知数的解 103 return 0; 104 }
![](https://gss1.bdstatic.com/9vo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D220/sign=6dd87712f11f3a295ec8d2cca924bce3/f9dcd100baa1cd119307a4ecb812c8fcc2ce2d4d.jpg)
![](https://gss0.bdstatic.com/94o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D220/sign=740fcbc8359b033b2888fbd825cf3620/37d12f2eb9389b50efe7369e8435e5dde6116ee6.jpg)
![](https://gss2.bdstatic.com/-fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D220/sign=337804ed91ef76c6d4d2fc29ad14fdf6/241f95cad1c8a7861fb29cc26609c93d71cf5066.jpg)
![](https://gss2.bdstatic.com/9fo3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D220/sign=1bed4ec85882b2b7a39f3ec601afcb0a/5d6034a85edf8db171baaaa80823dd54574e7462.jpg)
![](https://gss0.bdstatic.com/94o3dSag_xI4khGkpoWK1HF6hhy/baike/s%3D220/sign=1fa208041b4c510faac4e518505b2528/7dd98d1001e93901b9d95bf27aec54e737d19627.jpg)
参考:CDOJ 1330 柱爷与远古法阵【高斯消元,卡精度】 - Angel_Kitty - 博客园
https://www.cnblogs.com/ECJTUACM-873284962/p/7016987.html
转载于:https://www.cnblogs.com/Renyi-Fan/p/8137260.html
CDOJ 1330 柱爷与远古法阵(高斯消元)相关推荐
- 【bzoj2460】[BeiJing2011]元素 贪心+高斯消元求线性基
题目描述 相传,在远古时期,位于西方大陆的 Magic Land 上,人们已经掌握了用魔法矿石炼制法杖的技术.那时人们就认识到,一个法杖的法力取决于使用的矿石. 一般地,矿石越多则法力越强,但物极必反 ...
- Rocksdb Ribbon Filter : 结合 XOR-filter 以及 高斯消元算法 实现的 高效filter
文章目录 前言 XOR-filter 实现原理 xor filter 的构造原理 xor filter 构造总结 XOR-filter 和 ADD-filter对比 XOR-filter 在计算上的优 ...
- poj 1681 Painter#39;s Problem(高斯消元)
http://poj.org/problem? id=1681 求最少经过的步数使得输入的矩阵全变为y. 思路:高斯消元求出自由变元.然后枚举自由变元,求出最优值. 注意依据自由变元求其它解及求最优值 ...
- AC自动机 + 概率dp + 高斯消元 --- HDU 5955 or 2016年沈阳icpc H [AC自动机 + 概率dp + 高斯消元]详解
题目链接 题目大意: 就是有NNN个人,每个人都会猜一个长度为LLL的只包含{1,2,3,4,5,6}\{1,2,3,4,5,6\}{1,2,3,4,5,6}的序列,现在裁判开始投掷骰子,并且把每次的 ...
- ICPC 2005 hangzhou Generator (UVA1358)KMP + 期望DP / 高斯消元
整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 Generator Weblink https://www.luogu.com.cn/problem/ ...
- 2020 ACM / ICPC 济南 A Matrix Equation (高斯消元、乘法原理)
整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 题目链接 给你定义两种 010101 矩阵上的运算: Xi,j×Yi,j=(∑k=1NXi,kYk,j ...
- luogu P4035 [JSOI2008]球形空间产生器(高斯消元 / 模拟退火)
整理的算法模板合集: ACM模板 点我看算法全家桶系列!!! 实际上是一个全新的精炼模板整合计划 数据范围只开到了10,而且是经典的力学结构,所以我们可以用模拟退火,可以做一下 nnn 维的正交分解h ...
- BZOJ 2707: [SDOI2012]走迷宫 [高斯消元 scc缩点]
2707: [SDOI2012]走迷宫 题意:求s走到t期望步数,\(n \le 10^4\),保证\(|SCC| \le 100\) 求scc缩点,每个scc高斯消元,scc之间直接DP 注意每次清 ...
- [BZOJ 3143][Hnoi2013]游走(高斯消元+期望)
Description 一个无向连通图,顶点从1编号到N,边从1编号到M. 小Z在该图上进行随机游走,初始时小Z在1号顶点,每一步小Z以相等的概率随机选 择当前顶点的某条边,沿着这条边走到下一个顶点, ...
- [Sdoi2017]硬币游戏 [高斯消元 KMP]
[Sdoi2017]硬币游戏 题意:硬币序列,H T等概率出现,\(n \le 300\)个人猜了一个长为$ m \le 300$的字符串,出现即获胜游戏结束.求每个人获胜概率 考场用了[1444: ...
最新文章
- fibonacci的几种实现及尾递归
- 物料主数据(MM03)跳转函数
- linux 添加本地源,linux 添加本地yum源
- 网易云音乐的算法有什么特点_当算法设计音乐厅时会发生什么?
- python中控脚本_Python3.7实现中控考勤机自动连接
- 教你如何在google上查阅一个方位的经纬度
- Centos VNC 远程桌面程序的安装使用
- 使用官网支持工具制作富士施乐打印机一键安装的程序
- 机器学习——LASSO算法
- 动画过渡transition
- android开发 高德地图自定义地图应用【最新】
- Android开发之传感器(加速度传感器、方向传感器)
- 《互联网的那些事之时代》第二回:春江潮水连海平,海上明月共潮生
- ipad怎样和计算机连接网络,ipad怎样连接电脑itunes
- Jenkins与DevOps持续交付详解
- 医疗大数据平台的标准化通迅协议构建和架构
- Python和C语言哪个更容易学,感觉学了C语言有点难,只学过C语言的大学生很迷茫?
- 批量升级320k百度音乐(java版本)
- LED驱动电源EMI整改方案
- Encyclopaedia Britannica Ultimate 2014电子版下载|大不列颠百科全书
热门文章
- OPPO R9s成为情侣最喜欢礼物,今年情人节就送它了!
- 2015.4.7-C#入门基础(一)
- PowerTip of the Day-Removing Empty Things
- 超级简单:ASP.NET输出缓存
- ADSL拨号网络环境下实时视频广播的完美实现
- oracle alter database,alter database ... create datafile的原理及用途
- mysql数据库集群版_MySQL数据库集群实战
- mysql怎么判断多行数据日期是否连续_MySQL学习笔记(一)
- spark rdd详解二(transformation与action操作)
- 计算机指针知识,指针_计算机基础知识142页.ppt