Batting Practice LightOJ - 1408(概率dp)


题意:有无限个球,进球的概率为p,问你连续不进k1个球或者连续进k2个球需要使用的球的个数的期望

思路:

\(定义f[i]表示已经连续不进i个球,还需要连续不进k1-i个球的期望\)

\(g[i]表示已经连续进了i个,还需要连续进k2-i个球的期望\)

显然\(f[k1] = g[k2] = 0\)

\(任意0<=i<k1有f[i] = (1-p) \cdot f[i+1] + p \cdot g[1] + 1\)

\(任意0<=i<k2有g[i] = (1-p) \cdot f[1] + p \cdot g[i+1] + 1\)

不好解的样子,列矩阵高斯消元一发,WA,把eps从1e-6改到1e-15还是过不了,用long double交,时限卡的紧,TLE

搜了一发题解,原来上面那个式子可以直接求通项,依次回代就可以求出f[0]和g[0]了

#include<bits/stdc++.h>
#define LL long long
#define P pair<int,int>
using namespace std;
const double eps = 1e-6;
double a[60][60];
int gauss(int n,int m){
int col,i,mxr,j,row;
for(row=col=0;row<=n&&col<=m;row++,col++){
mxr = row;
for(i=row+1;i<=n;i++)
if(fabs(a[i][col])>fabs(a[mxr][col]))
mxr = i;
if(mxr != row) swap(a[row],a[mxr]);
if(fabs(a[row][col]) < eps){
row--;
continue;
}
for(i=0;i<=n;i++)///消成上三角矩阵
if(i!=row&&fabs(a[i][col])>eps)
for(j=m;j>=col;j--)
a[i][j]-=a[row][j]/a[row][col]*a[i][col];
}
row--;
for(int i = row;i>=0;i--){///回代成对角矩阵
for(int j = i + 1;j <= row;j++){
a[i][m] -= a[j][m] * a[i][j];
}
a[i][m] /= a[i][i];
}
return row;
}
int main()
{
int T;
int cas = 1;
cin>>T;
while(T--)
{
int k1,k2;
double p;
scanf("%lf%d%d",&p,&k1,&k2);
printf("Case %d: ",cas++);
/*
memset(a,0,sizeof(a));
int col = k1 + k2 + 2;
for(int i = 0;i < k1;i++){
a[i][i] = 1;
a[i][i+1] = p - 1;
a[i][k1+2] = -p;
a[i][col] = 1;
}
a[k1][k1] = 1,a[k1][col] = 0;
for(int i = 0;i < k2;i++){
a[k1+1+i][k1+1+i] = 1;
a[k1+1+i][k1+1+i+1] = -p;
a[k1+1+i][1] = p - 1;
a[k1+1+i][col] = 1;
}
a[k1+k2+1][k1+k2+1] = 1,a[k1+k2+1][col] = 0;
int row = gauss(k1+k2+1,k1+k2+2);
printf("%.6f\n",a[0][col]);
*/
if(p < eps) printf("%.6f\n",1.0 * k1);
else if(1 - p < eps) printf("%.6f\n",1.0 * k2);
else{
double q;
q=1-p;
double a1=1-pow(q,k1-1),b1=a1/(1-q);
double a2=1-pow(p,k2-1),b2=a2/(1-p);
double t1=(a1*b2+b1)/(1-a1*a2),f1=a2*t1+b2;
printf("%.6f\n",p*f1+q*t1+1);
}
}
return 0;
}

Batting Practice LightOJ - 1408的更多相关文章

  1. lightoj 1408 Batting Practice (概率问题,求期望,推公式)

    题意:一个人若连续进k1个球或连续不进k2个球,游戏结束,给出这个人不进球的概率p(注意:是不进球!!!),求到游戏结束时这个投球个数的期望. 不进球概率为p,进概率 q=1-p.设 f[i] 表示连 ...

  2. lightoj 1408 Batting Practice

    题意:一个人若连续进k1个球或连续不进k2个球,游戏结束,给出这个人进球的概率p,求到游戏结束时这个投球个数的期望. 进球概率为p,不进概率 q=1-p 设 f[i] 表示连续 i 次不进距离连续k2 ...

  3. LightOj_1408 Batting Practice

    题目链接 题意: 击球训练中, 你击中一个球的概率为p,连续击中k1个球, 或者连续击空k2个球, 则训练结束. 求结束训练所击球次数的期望. 思路: 设f[x]为连续击中x个球, 距离结束训练所需要 ...

  4. lightoj 1408 概率dp

    https://blog.csdn.net/moon_sky1999/article/details/98097470 博主在此,牛逼神犇 #include<bits/stdc++.h> ...

  5. LightOJ 1341 唯一分解定理

    Aladdin and the Flying Carpet Time Limit:3000MS     Memory Limit:32768KB     64bit IO Format:%lld &a ...

  6. lightoj 1370 欧拉函数

    A - Bi-shoe and Phi-shoe Time Limit:2000MS     Memory Limit:32768KB     64bit IO Format:%lld & % ...

  7. lightoj 1074 spfa判断负环

     Extended Traffic Time Limit:2000MS     Memory Limit:32768KB     64bit IO Format:%lld & %llu Sub ...

  8. lightoj.1048.Conquering Keokradong(二分 + 贪心)

    Conquering Keokradong Time Limit:1000MS     Memory Limit:32768KB     64bit IO Format:%lld & %llu ...

  9. LightOJ 1234 Harmonic Number

    D - Harmonic Number Time Limit:3000MS     Memory Limit:32768KB     64bit IO Format:%lld & %llu S ...

随机推荐

  1. python基础数据类型之字典的操作

    一. 字典的简单介绍字典(dict)是python中唯一的一个映射类型.他是以{ }括起来的键值对组成. 在dict中key是唯一的. 在保存的时候, 根据key来计算出一个内存地址. 然后将key- ...

  2. 给网站添加icon图标

    只需制成ico结尾的图片即可

  3. 以CentOS为操作系统的vps或服务器安装lnmp运行环境的方法

    安装步骤: 1.使用putty或类似的SSH工具登陆VPS或服务器: 登陆后运行:screen -S lnmp 如果提示screen: command not found 命令不存在可以执行:yum ...

  4. windows下的node.js和npm的安装步骤详解

    一.使用之前,我们先来掌握3个东西是用来干什么的. npm: Nodejs下的包管理器. webpack: 它主要的用途是通过CommonJS的语法把所有浏览器端需要发布的静态资源做相应的准备,比如资 ...

  5. PHP 二维数组按某一个键值排序

    一.前言 在某个项目中,需要读取某个文件夹下的所有文件,在本地的 Windows 环境下时,读取出来的二维数组的文件名称和在 Windows 文件夹的文件排序一致, 但是项目上线后,环境为 Linux ...

  6. 【异常】The server time zone value 'Öйú±ê׼ʱ¼ä' is unrecognized or represents more than one time zone.

    异常错误:The server time zone value 'Öйú±ê׼ʱ¼ä' is unrecognized or represents more than one time zone ...

  7. POJ:2674-Linear world(名字交换碰撞)

    Linear world Time Limit: 3000MS Memory Limit: 65536K Total Submissions: 4514 Accepted: 1025 Descript ...

  8. Diycode开源项目 如何解决InputMethodManager造成的内存泄漏问题

    1.内存泄漏的状况及原因 1.1.利用LeakCanary查看内存泄漏的状况 1.2.内存泄漏怎么产生的呢? InputMethodManager.mServicedView持有一个最后聚焦View的 ...

  9. 10,before_request 和 after_request

    Flask我们已经学习很多基础知识了,现在有一个问题 我们现在有一个 Flask 程序其中有3个路由和视图函数,如下: from flask import Flask app = Flask(__na ...

  10. 常用doc 命令

    开始-->运行 regedit 进入注册表补充些: 1. gpedit.msc-----组策略 2. sndrec32-------录音机 3. Nslookup-------IP地址侦测器 4 ...