CDOJ 1330 柱爷与远古法阵【高斯消元,卡精度】

柱爷与远古法阵

Time Limit: 125/125MS (Java/Others)     Memory Limit: 240000/240000KB (Java/Others)

Submit Status

众所周知,柱爷的数学非常好,尤其擅长概率论!

某日柱爷在喵哈哈村散步,无意间踏入了远古法阵!

法阵很奇怪,是一个长度为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

2016 UESTC Training for Dynamic Programming

题目链接:http://acm.uestc.edu.cn/#/problem/show/1330

分析:第一次学高斯消元,看着卿学姐的视频长大的QAQ,然后,,,,,卿学姐的代码输出结果,,,,一直为0.0000000,我奋力找了5个小时的bug,最后竟然错在输出的精度问题,要取double,我也不知道为啥QAQ!

下面给出AC代码:

  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 }
时间: 2024-10-28 08:41:04

CDOJ 1330 柱爷与远古法阵【高斯消元,卡精度】的相关文章

UVALive - 6185 Find the Outlier暴力填表+高斯消元+卡eps

https://cn.vjudge.net/problem/UVALive-6185 我真的是服了orz eps 1e5,1e6过不了 开1e2 1e1都能过 题意:给你一个d阶多项式f的f(0),f(1)...f(d+1),f(d+2) 有一个是错误的,问第几个是错的 题解:题目多给了两个方程(约束). 想了一下如果只给一个,是找不出来的. 给两个的话,可以这么考虑: 先取出一个方程X,再取剩下的n个高斯消元一下,将解得的系数带入最后一个方程,if成立,说明X是错的,else再取另一个(说明错

【bzoj2460】[BeiJing2011]元素 贪心+高斯消元求线性基

题目描述 相传,在远古时期,位于西方大陆的 Magic Land 上,人们已经掌握了用魔法矿石炼制法杖的技术.那时人们就认识到,一个法杖的法力取决于使用的矿石.一般地,矿石越多则法力越强,但物极必反:有时,人们为了获取更强的法力而使用了很多矿石,却在炼制过程中发现魔法矿石全部消失了,从而无法炼制出法杖,这个现象被称为“魔法抵消” .特别地,如果在炼制过程中使用超过一块同一种矿石,那么一定会发生“魔法抵消”. 后来,随着人们认知水平的提高,这个现象得到了很好的解释.经过了大量的实验后,著名法师 D

【弱校胡策】2016.4.14 (bzoj2164)最短路+状压DP+矩阵乘法+高斯消元+树链剖分+线段树+背包DP

cyyz&qhyz&lwyz&gryz弱校胡策 命题人:cyyz ws_fqk T3暴力写挫了 50+10+0滚粗辣! 奇妙的约会(appointment.cpp/c/pas) [问题描述] DQS和sxb在网上结识后成为了非常好的朋友,并且都有着惊人 的OI水平.在NOI2333的比赛中,两人均拿到了金牌,并保送进入 HU/PKU.于是两人决定在这喜大普奔的时刻进行面基. NOI2333参赛选手众多,所以安排了n个考点,DQS在1号考点, 而sxb在n号考点.由于是举办全国性赛事

poj_1222_高斯消元

第一次学习使用高斯消元,将灯板化为线性方程组,进行求解. /*######################################################################### # File Name: poj_1222.cpp # Author: CaoLei # Created Time: 2015/7/20 15:48:04 ###################################################################

HDU 4870 Rating(高斯消元)

HDU 4870 Rating 题目链接 题意:一个人注册两个账号,初始rating都是0,他每次拿低分的那个号去打比赛,赢了加50分,输了扣100分,胜率为p,他会打到直到一个号有1000分为止,问比赛场次的期望 思路:f(i, j)表示i >= j,第一个号i分,第二个号j分时候,达到目标的期望,那么可以列出转移为f(i, j) = p f(i', j') + (1 - p) f(i'' + j'') + 1 f(i', j')对应的是赢了加分的状态,f(i'', j'')对应输的扣分的状态

【BZOJ 4171】 4171: Rhl的游戏 (高斯消元)

4171: Rhl的游戏 Time Limit: 20 Sec  Memory Limit: 128 MBSubmit: 74  Solved: 33[Submit][Status][Discuss] Description RHL最近迷上一个小游戏:Flip it.游戏的规则很简单,在一个N*M的格子上,有一些格子是黑色,有一些是白色 .每选择一个格子按一次,格子以及周围边相邻的格子都会翻转颜色(边相邻指至少与该格子有一条公共边的格子 ),黑变白,白变黑.RHL希望把所有格子都变成白色的.不幸

POJ 1830 开关问题 高斯消元,自由变量个数

http://poj.org/problem?id=1830 如果开关s1操作一次,则会有s1(记住自己也会变).和s1连接的开关都会做一次操作. 那么设矩阵a[i][j]表示按下了开关j,开关i会被操作一次,记得a[i][i] = 1是必须的,因为开关i操作一次,本身肯定会变化一次. 所以有n个开关,就有n条方程, 每个开关的操作次数总和是:a[i][1] + a[i][2] + ... + a[i][n] 那么sum % 2就代表它的状态,需要和(en[i] - be[i] + 2) % 2

BZOJ 3105: [cqoi2013]新Nim游戏 [高斯消元XOR 线性基]

以后我也要用传送门! 题意:一些数,选择一个权值最大的异或和不为0的集合 终于有点明白线性基是什么了...等会再整理 求一个权值最大的线性无关子集 线性无关子集满足拟阵的性质,贪心选择权值最大的,用高斯消元判断是否和已选择的线性相关 每一位记录pivot[i]为i用到的行 枚举要加入的数字的每一个二进制为1的位,如果有pivot[i]那么就异或一下(消元),否则pivot[i]=这个数并退出 如果最后异或成0了就说明线性相关... #include <iostream> #include &l

[bzoj1013][JSOI2008]球形空间产生器sphere-题解[高斯消元]

Description 有一个球形空间产生器能够在n维空间中产生一个坚硬的球体.现在,你被困在了这个n维球体中,你只知道球面上n+1个点的坐标,你需要以最快的速度确定这个n维球体的球心坐标,以便于摧毁这个球形空间产生器. Input 第一行是一个整数n(1<=N=10).接下来的n+1行,每行有n个实数,表示球面上一点的n维坐标.每一个实数精确到小数点后6位,且其绝对值都不超过20000. Output 有且只有一行,依次给出球心的n维坐标(n个实数),两个实数之间用一个空格隔开.每个实数精确到