题意:给定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;
}