hdu5810Balls and Boxes

链接:http://acm.hdu.edu.cn/showproblem.php?pid=5810

题意:给定n个不同的球放入m个不同的盒子中,设第i个盒子个数为xi个平均数为X,求方差V=sigma(xi-X)^2/m的期望。

分析:碰到概率期望的题只能懵逼,刚好补一补。首先我们知道E[v]=E[sigma(xi-X)^2/m]=E[(xi-X)^2]=E[xi^2-2*xi*X+X^2],然后因为期望的独立性我们可以分开算那么E[xi^2-2*xi*X+X^2]=E[xi^2]-2*E[xi]*X+E[X^2],显然后面两部分-2*E[xi]*X+E[X*2]=-X^2,那么就只需要求E[xi^2]了,第i个盒子的个数的平方的期望,又因为每个球之间也是相互独立的,那么我们设yj表示第j个球是否在第i个盒子中,那么E[xi^2]=E[(sigma(yj))^2]=E[sigma(yj^2)]+E[sigma(yj)*(yk)]=n/m+n*(n-1)/m^2,代回去得E[v]=n*(m-1)/m^2。

代码:

#include<map>
#include<set>
#include<cmath>
#include<queue>
#include<bitset>
#include<math.h>
#include<vector>
#include<string>
#include<stdio.h>
#include<cstring>
#include<iostream>
#include<algorithm>
#pragma comment(linker, "/STACK:102400000,102400000")
using namespace std;
const int N=1000010;
const int mod=100000000;
const int MOD1=1000000007;
const int MOD2=1000000009;
const double EPS=0.00000001;
typedef long long ll;
const ll MOD=1000000007;
const int INF=2000000010;
const ll MAX=1ll<<55;
const double pi=acos(-1.0);
const double eps=1e-5;
typedef double db;
typedef unsigned long long ull;
ll gcd(ll a,ll b) {
    return b==0 ? a:gcd(b,a%b);
}
int main()
{
    ll n,m,g,fz,fm;
    while (scanf("%I64d%I64d", &n, &m)&&(n||m)) {
        fz=n*(m-1);fm=m*m;
        g=gcd(fz,fm);
        fz/=g;fm/=g;
        printf("%I64d/%I64d\n", fz, fm);
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值