UESTC 1330 柱爷与远古法阵【高斯消元】
题目链接【http://acm.uestc.edu.cn/#/problem/show/1330】
题意:有一个长度为L(L <= 300)的长廊,有一人站在最左边,他要到最右边去,他每次可以走1 ~ 6 步,每一次走的步数是随机发生的。并且某些位置有传送门(a,b)表示a位置有个传送门,可以直接到达b,然后求出到达最右边的期望是多少?
题解:首先,我们先列出公式,设dp[i]表示从 i 点到达最右边的期望,则1、dp[L] =0 。2、如果a有传送门b那么dp[a] = dp[b] 。否则 dp[x] = (dp[x+1]+dp[x+2] + dp[x+y])/6+1;( x + y <= L && y <= 6)。化简后就得到,6 * dp[x] - dp[x+1] - dp[x+2] - ... - dp[x+y] == 6。那么我们一共可以列出L个方程,一共有L个未知量,只要用高斯消元解出来就可以了。
注意:1、每个点最多只有一个传送门,2、这道题存在无解的情况,3、这道题精度卡的很高。
#include<bits/stdc++.h>
using namespace std;
typedef long double LB;
const LB eps = 1e-14;
const int maxn = 350 + 15;
LB a[maxn][maxn], x[maxn];
int Guess(int equ, int var)
{
int i, j, k, col, max_r;
for(k = 1, col = 1; k <= equ && col <= var; k++, col++)
{
max_r = k;
for(i = k + 1; i <= equ; i++)
if(fabs(a[i][col] ) > fabs(a[max_r][col])) max_r = i;
if(fabs(a[max_r][col]) < eps ) return 0;//无解
if(k != max_r)
{
for(j = col; j <= var; j++) swap(a[k][j], a[max_r][j]);
swap(x[k], x[max_r]);
}
x[k] /= a[k][col];
for(j = col + 1; j <= var; j++) a[k][j] /= a[k][col];
a[k][col] = 1;
for(i = 1; i <= equ; i++)
if(i != k)
{
x[i] -= x[k] * a[i][col];
for(j = col + 1; j <= var; j++) a[i][j] -= a[k][j] * a[i][col];
a[i][col] = 0;
}
}
return 1;
}
int N, M;
int f[maxn];//表示某点有没有传送门
int main ()
{
scanf("%d %d", &N, &M);
memset(a, 0, sizeof(a));
for(int i = 1; i <= N; i++) f[i] = 0;
for(int i = 1; i <= M; i++)
{
int u, v;
scanf("%d %d", &u, &v);
f[u] = v;
}
for(int i = 1; i < N; i++)
{
a[i][i] = 6.0;
if(f[i])
{
a[i][f[i]] = -6.0;
continue;
}
x[i] = 6.0;
for(int j = 1; j <= 6; j++)
{
if(i + j <= N) a[i][i + j] -= 1.0;
else a[i][i] -= 1.0;
}
}
a[N][N] = 1.0, x[N] = 0;
if(!Guess(N, N)) printf("-1\n");
else printf("%.12Lf\n", x[1]);
return 0;
}
UESTC 1330 柱爷与远古法阵【高斯消元】的更多相关文章
- CDOJ 1330 柱爷与远古法阵(高斯消元)
CDOJ 1330 柱爷与远古法阵(高斯消元) 柱爷与远古法阵 Time Limit: 125/125MS (Java/Others) Memory Limit: 240000/240000K ...
- CDOJ 1330 柱爷与远古法阵【高斯消元,卡精度】
柱爷与远古法阵 Time Limit: 125/125MS (Java/Others) Memory Limit: 240000/240000KB (Java/Others) Submit S ...
- ACM学习历程—UESTC 1219 Ba Gua Zhen(dfs && 独立回路 && xor高斯消元)
题目链接:http://acm.uestc.edu.cn/#/problem/show/1219 题目大意是给了一张图,然后要求一个点通过路径回到这个点,使得xor和最大. 这是CCPC南阳站的一道题 ...
- hdu3949 XOR xor高斯消元
XOR Time Limit: 2000/1000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Submiss ...
- 【BZOJ-3143】游走 高斯消元 + 概率期望
3143: [Hnoi2013]游走 Time Limit: 10 Sec Memory Limit: 128 MBSubmit: 2264 Solved: 987[Submit][Status] ...
- 【BZOJ-3270】博物馆 高斯消元 + 概率期望
3270: 博物馆 Time Limit: 30 Sec Memory Limit: 128 MBSubmit: 292 Solved: 158[Submit][Status][Discuss] ...
- *POJ 1222 高斯消元
EXTENDED LIGHTS OUT Time Limit: 1000MS Memory Limit: 10000K Total Submissions: 9612 Accepted: 62 ...
- [bzoj1013][JSOI2008][球形空间产生器sphere] (高斯消元)
Description 有一个球形空间产生器能够在n维空间中产生一个坚硬的球体.现在,你被困在了这个n维球体中,你只知道球 面上n+1个点的坐标,你需要以最快的速度确定这个n维球体的球心坐标,以便于摧 ...
- hihoCoder 1196 高斯消元·二
Description 一个黑白网格,点一次会改变这个以及与其连通的其他方格的颜色,求最少点击次数使得所有全部变成黑色. Sol 高斯消元解异或方程组. 先建立一个方程组. \(x_i\) 表示这个点 ...
随机推荐
- JavaScript设置粘贴板
设置复制 document.body.oncopy = function(){ alert('不许复制'); return false; }; 设置粘贴 document.getElementById ...
- 巧用margin/padding的百分比值实现高度自适应
原文:https://segmentfault.com/a/1190000004231995 一个基础却又容易混淆的css知识点 本文依赖于一个基础却又容易混淆的css知识点:当margin/padd ...
- 使用showplan.sql分析sql Performance
在HelloDBA网站找到一个分析sql性能的工具-showplan,记录一下 showplan.sql下载路径:http://www.HelloDBA.com/Download/showplan.z ...
- Ubuntu命令设置ip网关dns
本文系转载,介绍Ubuntu如何设置IP和网络来连接网络 如果是在虚拟机中使用Ubuntu,那么设置之前请先参照我的上一遍文章虚拟机Net方式设置连接外网中的网络设置部分,先设置好主机的网络,然后配置 ...
- MySQL Warning: Using a password on the command line interface can be insecure.解决办法
转自 http://www.cnblogs.com/sunss/p/6256706.html 被一个小朋友问到,直接公布答案: If your MySQL client/server version ...
- 简单实现JS上传图片预览功能
HTML代码 <div class="upload"> <input type="button" class="btn" ...
- 338.Counting Bits---位运算---《剑指offer》32
题目链接:https://leetcode.com/problems/counting-bits/description/ 题目大意:求解从0到num的所有数的二进制表示中所有1的个数. 法一:暴力解 ...
- JavaScript 跳转 页面
* window.location.href , self.location, window.location 出现问题不能跳转 Chome 不能本页跳转, IE 有时可以
- PHP的instanceof关键字
PHP5的另一个新成员是instdnceof关键字.使用这个关键字可以确定一个对象是类的实例.类的子类,还是实现了某个特定接口,并进行相应的操作.在某些情况下,我们希望确定某个类是否特定的类型,或者是 ...
- 自家人不认识自家人——考你一道有趣的Javascript小题目
今天的内容很简单,给大家分享一个有趣的Javascript小题目. 题目很简单,就是填空: var a = ______; var b = a; alert(a==b); // alert " ...