VK Cup 2016 - Round 2 D. Little Artem and Random Variable(已知两个撒子掷出的点数较大为1,2,...,n的概率,较小为1,2,...,n概率)
2016-05-09 10:04
573 查看
D. Little Artem and Random Variable
time limit per test
2 seconds
memory limit per test
256 megabytes
input
standard input
output
standard output
Little Artyom decided to study probability theory. He found a book with a lot of nice exercises and now wants you to help him with one of them.
Consider two dices. When thrown each dice shows some integer from 1 to n inclusive.
For each dice the probability of each outcome is given (of course, their sum is 1), and different dices may have different probability distributions.
We throw both dices simultaneously and then calculate values max(a, b) and min(a, b),
where a is equal to the outcome of the first dice, while b is
equal to the outcome of the second dice. You don't know the probability distributions for particular values on each dice, but you know the probability distributions for max(a, b) and min(a, b).
That is, for each x from 1 to n you
know the probability thatmax(a, b) would be equal to x and
the probability that min(a, b) would be equal to x.
Find any valid probability distribution for values on the dices. It's guaranteed that the input data is consistent, that is, at least one solution exists.
Input
First line contains the integer n (1 ≤ n ≤ 100 000) —
the number of different values for both dices.
Second line contains an array consisting of n real values with up to 8 digits after the decimal point — probability distribution
formax(a, b), the i-th
of these values equals to the probability that max(a, b) = i. It's guaranteed that the sum of these values for one dice is 1.
The third line contains the description of the distribution min(a, b) in the same format.
Output
Output two descriptions of the probability distribution for a on the first line and for b on
the second line.
The answer will be considered correct if each value of max(a, b) and min(a, b)
probability distribution values does not differ by more than10 - 6 from
ones given in input. Also, probabilities should be non-negative and their sums should differ from 1 by no more than 10 - 6.
Examples
input
output
input
output
题意:
有两个不同的撒子,每个撒子都可以掷出1-n,而且掷出每个点的概率都是未知的,
现在已知两个撒子掷出的点数较大为1,2,...,n的概率,较小为1,2,...,n的概率
求每个撒子掷出1,2,...,n的概率
思路:
设P(a<=k)为x1,P(b<=k)为x2
P(max(a,b)<=k)=P(a<=k)*P(b<=k)=x1*x2=第一行的前缀和
P(min(a,b)>=k+1)=(1-P(a<=k))*(1-P(b<=k))=1+x1*x2-x1-x2=第二行的后缀和
根据韦达定理,可以求出关于x1,x2的二元方程组的b,c,a默认为1,解方程便可以了。
#include<bits/stdc++.h>
using namespace std;
const int maxn=100100;
double ans1[maxn],ans2[maxn];
double pre[maxn],suf[maxn];
void solve(double b,double c,double &x1,double &x2){
double x=max(b*b-4*c,0.0);
x1=(-b-sqrt(x))*0.5;
x2=(-b+sqrt(x))*0.5;
}
int main(){
int n;
double a;
scanf("%d",&n);
for(int i=1;i<=n;i++){
scanf("%lf",&a);
pre[i]=pre[i-1]+a;
}
for(int i=1;i<=n;i++)
scanf("%lf",&suf[i]);
for(int i=n;i>=1;i--)
suf[i]+=suf[i+1];
double x1,x2;
for(int i=1;i<=n;i++){
solve(-pre[i]-1+suf[i+1],pre[i],x1,x2);
ans1[i]=x1,ans2[i]=x2;
}
for(int i=1;i<=n;i++)
printf("%.10f ",ans1[i]-ans1[i-1]);
printf("\n");
for(int i=1;i<=n;i++)
printf("%.10f ",ans2[i]-ans2[i-1]);
printf("\n");
}
time limit per test
2 seconds
memory limit per test
256 megabytes
input
standard input
output
standard output
Little Artyom decided to study probability theory. He found a book with a lot of nice exercises and now wants you to help him with one of them.
Consider two dices. When thrown each dice shows some integer from 1 to n inclusive.
For each dice the probability of each outcome is given (of course, their sum is 1), and different dices may have different probability distributions.
We throw both dices simultaneously and then calculate values max(a, b) and min(a, b),
where a is equal to the outcome of the first dice, while b is
equal to the outcome of the second dice. You don't know the probability distributions for particular values on each dice, but you know the probability distributions for max(a, b) and min(a, b).
That is, for each x from 1 to n you
know the probability thatmax(a, b) would be equal to x and
the probability that min(a, b) would be equal to x.
Find any valid probability distribution for values on the dices. It's guaranteed that the input data is consistent, that is, at least one solution exists.
Input
First line contains the integer n (1 ≤ n ≤ 100 000) —
the number of different values for both dices.
Second line contains an array consisting of n real values with up to 8 digits after the decimal point — probability distribution
formax(a, b), the i-th
of these values equals to the probability that max(a, b) = i. It's guaranteed that the sum of these values for one dice is 1.
The third line contains the description of the distribution min(a, b) in the same format.
Output
Output two descriptions of the probability distribution for a on the first line and for b on
the second line.
The answer will be considered correct if each value of max(a, b) and min(a, b)
probability distribution values does not differ by more than10 - 6 from
ones given in input. Also, probabilities should be non-negative and their sums should differ from 1 by no more than 10 - 6.
Examples
input
2 0.25 0.75 0.75 0.25
output
0.5 0.5 0.5 0.5
input
3 0.125 0.25 0.625 0.625 0.25 0.125
output
0.25 0.25 0.5 0.5 0.25 0.25
题意:
有两个不同的撒子,每个撒子都可以掷出1-n,而且掷出每个点的概率都是未知的,
现在已知两个撒子掷出的点数较大为1,2,...,n的概率,较小为1,2,...,n的概率
求每个撒子掷出1,2,...,n的概率
思路:
设P(a<=k)为x1,P(b<=k)为x2
P(max(a,b)<=k)=P(a<=k)*P(b<=k)=x1*x2=第一行的前缀和
P(min(a,b)>=k+1)=(1-P(a<=k))*(1-P(b<=k))=1+x1*x2-x1-x2=第二行的后缀和
根据韦达定理,可以求出关于x1,x2的二元方程组的b,c,a默认为1,解方程便可以了。
#include<bits/stdc++.h>
using namespace std;
const int maxn=100100;
double ans1[maxn],ans2[maxn];
double pre[maxn],suf[maxn];
void solve(double b,double c,double &x1,double &x2){
double x=max(b*b-4*c,0.0);
x1=(-b-sqrt(x))*0.5;
x2=(-b+sqrt(x))*0.5;
}
int main(){
int n;
double a;
scanf("%d",&n);
for(int i=1;i<=n;i++){
scanf("%lf",&a);
pre[i]=pre[i-1]+a;
}
for(int i=1;i<=n;i++)
scanf("%lf",&suf[i]);
for(int i=n;i>=1;i--)
suf[i]+=suf[i+1];
double x1,x2;
for(int i=1;i<=n;i++){
solve(-pre[i]-1+suf[i+1],pre[i],x1,x2);
ans1[i]=x1,ans2[i]=x2;
}
for(int i=1;i<=n;i++)
printf("%.10f ",ans1[i]-ans1[i-1]);
printf("\n");
for(int i=1;i<=n;i++)
printf("%.10f ",ans2[i]-ans2[i-1]);
printf("\n");
}
相关文章推荐
- sftp 限制用户登陆指定目录(家目录)
- 已超过传入消息(65536)的最大消息大小配额
- jQuery AJAX实现调用页面后台方法
- 【回文子串的个数】HDU1544Palindromes【字符缩距/奇偶判断】
- Dream_Spark-----Spark 定制版:005~贯通Spark Streaming流计算框架的运行源码
- Bootstrap进阶七:LESS语法详解
- MBProgressHUD 提示自动换行
- TabBar高度修改以及UITabBarItem坐标移动
- MyBatis 入门(五)--typeHandlers
- 用collectionview实现瀑布流-转(后面附demo,供参考)
- (整理)SQLServer_DBA 工具
- PHP 使用 password_hash() 给密码加密
- 界面设计,事件处理--GUI编程的核心--
- 启动tomcat时 错误: 代理抛出异常 : java.rmi.server.ExportException: Port already in use: 1099的解决办法
- Linux内核中的PID散列表实例
- iOS多线程开发――NSThread浅析
- Win32 常用数据类型
- quartz储存方式之JDBC JobStoreTX
- 判断一个单链表是否有环及环的连接点
- 百度地图API