HDU 4870 Rating (高斯消元)
题目链接 2014 多校1 Problem J
题意 现在有两个账号,初始$rating$都为$0$,现在每次打分比较低的那个,如果进前$200$那么就涨$50$分,否则跌$100$分。
每一次打进前$200$的概率为$p$,且每一次竞赛是相互独立的。求当一个号打到$1000$分时已经打的期望场数。
把$1000$变成$20$,并且分析出所有有效的状态。
有效的状态一共$211$个。
$(0, 0)$
$(1, 0), (1, 1)$
$(2, 0), (2, 1), (2, 2)$
......
$(19, 0), (19, 1), (19, 2), ..., (19, 19)$
再加一个$(20, 19)$,这是目标状态。
设$f(x, y)$为$(x, y)$状态时所需要到达目标的期望场数。
显然$f(19, 20) = 0$
设$(x, y)$赢了之后状态为$(x1, y1)$, 输了之后状态为$(x2, y2)$
那么$f(x, y) = p*f(x1, y1) + (1-p)*f(x2, y2) + 1$
注意状态$(0, 0)$输了之后还是$(0, 0)$
那么建立211个方程组,高斯消元求解即可。
#include <bits/stdc++.h> using namespace std; #define rep(i, a, b) for (int i(a); i <= (b); ++i)
#define dec(i, a, b) for (int i(a); i >= (b); --i) typedef long long LL; const int N = 235;
const double eps = 1e-12; double a[N][N], x[N], p;
int equ, var;
int c[N][N];
int cnt;
int y; void Gauss(){
int row, col, max_r;
row = col = 0;
while (row < equ && col < var){
max_r = row;
rep(i, row + 1, equ - 1) if (fabs(a[i][col]) - fabs(a[max_r][col]) > eps) max_r = i;
if (max_r != row) rep(j, col, var) swap(a[row][j], a[max_r][j]); if (fabs(a[row][col]) < eps){
col++;
continue;
} rep(i, row + 1, equ - 1){
if (fabs(a[i][col]) > eps){
double t = a[i][col] / a[row][col];
a[i][col] = 0.0;
rep(j, col + 1, var) a[i][j] -= a[row][j] * t;
}
}
row++;
col++;
} dec(i, equ - 1, 0){
if (fabs(a[i][i]) < eps) continue;
double tmp = a[i][var];
rep(j, i + 1, var - 1) tmp -= a[i][j] * x[j];
x[i] = tmp / a[i][i];
}
} int main(){ cnt = 0;
for (int i = 0; i <= 22; ++i){
for (int j = 0; j <= 22; ++j){
c[i][j] = -1e9;
}
} rep(i, 0, 19) rep(j, 0, i) c[i][j] = cnt++;
c[20][19] = cnt++; equ = var = cnt; while (~scanf("%lf", &p)){
memset(a, 0, sizeof a);
rep(i, 0, 19){
rep(j, 0, i - 1){
y = c[i][j];
a[y][y] = 1;
a[y][211] = 1;
a[y][c[i][max(0, j - 2)]] -= (1 - p);
a[y][c[i][j + 1]] -= p;
}
y = c[i][i];
a[y][y] = 1;
a[y][211] = 1;
a[y][c[i][max(0, i - 2)]] -= (1 - p);
a[y][c[i + 1][i]] -= p;
} a[210][210] = 1;
Gauss();
printf("%.6lf\n", x[0]);
} return 0;
}
HDU 4870 Rating (高斯消元)的更多相关文章
- hdu 4870 rating(高斯消元求期望)
Rating Time Limit: 10000/5000 MS (Java/Others) Memory Limit: 32768/32768 K (Java/Others)Total Sub ...
- HDU 3949 XOR 高斯消元
题目大意:给定一个数组,求这些数组通过异或能得到的数中的第k小是多少 首先高斯消元求出线性基,然后将k依照二进制拆分就可以 注意当高斯消元结束后若末尾有0则第1小是0 特判一下然后k-- 然后HDU输 ...
- hdu 5755(高斯消元——模线性方程组模板)
PS. 看了大神的题解,发现确实可以用m个未知数的高斯消元做.因为确定了第一行的情况,之后所有行的情况都可以根据第一行推. 这样复杂度直接变成O(m*m*m) 知道了是高斯消元后,其实只要稍加处理,就 ...
- HDU 3949 XOR [高斯消元XOR 线性基]
3949冰上走 题意: 给你 N个数,从中取出若干个进行异或运算 , 求最后所有可以得到的异或结果中的第k小值 N个数高斯消元求出线性基后,设秩为$r$,那么总共可以组成$2^r$中数字(本题不能不选 ...
- Time travel HDU - 4418(高斯消元)
Agent K is one of the greatest agents in a secret organization called Men in Black. Once he needs to ...
- hdu 6465 线性变换高斯消元
http://acm.hdu.edu.cn/showproblem.php?pid=6465 题意 给你三个点,再给你经过线性变换后的三个点,然后q次询问,给你一个点,需要你输出线性变换后的点 题解 ...
- HDU 3364 Lanterns 高斯消元
Lanterns Problem Description Alice has received a beautiful present from Bob. The present contains ...
- HDU 3915 Game (高斯消元)
题意:有n堆石子,每个人只能从某一堆至少拿走一个,不能拿者败.问事先拿走某些堆的石子,使得先手必败. 析:将石子拆成二进制,未知数为1表示保留该堆石子,为0表示事先拿走该堆石子.最后求自由变元的数目, ...
- HDU 4870 Rating(高斯消元 )
HDU 4870 Rating 这是前几天多校的题目,高了好久突然听旁边的大神推出来说是可以用高斯消元,一直喊着赶快敲模板,对于从来没有接触过高斯消元的我来说根本就是一头雾水,无赖之下这几天做DP ...
- 2014多校第一场J题 || HDU 4870 Rating(DP || 高斯消元)
题目链接 题意 :小女孩注册了两个比赛的帐号,初始分值都为0,每做一次比赛如果排名在前两百名,rating涨50,否则降100,告诉你她每次比赛在前两百名的概率p,如果她每次做题都用两个账号中分数低的 ...
随机推荐
- Codeforces Round #460 (Div. 2)-D. Substring
D. Substring time limit per test3 seconds memory limit per test256 megabytes Problem Description You ...
- IOS开发学习笔记016-Foundation框架
Foundation 框架的学习 一.Foundation 常用结构体 1.NSRange(location,length) typedef struct _NSRange { NSUIntege ...
- NGUI 学习总结
NGUI 学习一段时间了,这里总结一下,用于以后查看. 获取组件 在Awake函数里获取组件,然后就可在Start以及其他函数里使用 lbl = GetComponent<UILabel> ...
- MongoDB快速入门学习笔记7 MongoDB的用户管理操作
1.修改启动MongoDB时要求用户验证加参数 --auth 即可.现在我们把MongoDB服务删除,再重新添加服务 mongod --dbpath "D:\work\MongoDB\dat ...
- [oldboy-django][2深入django]老师管理 -- form表单如何生成多选框标签,多选框的默认值显示,以及多选框数据插入到数据库,多选框数据更改到数据库
1 form表单如何生成多选框(包含了多选框可选择内容) - Form设置班级输入框为 select多选 - 多选 class TeacherForm(Form): name = fields.Cha ...
- Python设计模式之一(单例模式)
单例模式就是告诉你,只有一个对象 (1)单例模式不适用的场景 #单例模式就是告诉你,其实只有一个对象 class Person: def __init__(self,name,age): self.n ...
- w3wp.exe占用cpu特别高
w3wp.exe占用cpu特别高,百度了一下在任务管理器标记出PID可以看到进程号. 试了一下,发现一个xxx网站占用cpu特别高,然后就结束了下进程,再重启网站cpu一下子降下来. 很奇怪,还需要具 ...
- layer 体验
做后端的,前端一直不擅长,提示语以前也只会alert,非常难看,后来在别人代码看到lay.msg() https://pan.baidu.com/s/1eRHH59g <!DOCTYPE htm ...
- BZOJ 3990: [SDOI2015]排序(搜索+剪枝)
[SDOI2015]排序 Description 小A有一个1-2^N的排列A[1..2^N],他希望将A数组从小到大排序,小A可以执行的操作有N种,每种操作最多可以执行一次,对于所有的i(1< ...
- 省选算法学习-数据结构-splay
于是乎,在丧心病狂的noip2017结束之后,我们很快就要迎来更加丧心病狂的省选了-_-|| 所以从写完上一篇博客开始到现在我一直深陷数据结构和网络流的漩涡不能自拔 今天终于想起来写博客(只是懒吧.. ...