UVA 10886 Standard Deviation

https://vjudge.net/problem/UVA-10886

计算标准差

碰到这种题将式子展开

#include<cmath>
#include<cstdio>
//#include<iostream>
using namespace std;
unsigned long long seed;
long double gen()
{
    static const long double Z = ( long double )1.0 / (1LL<<);
    seed >>= ;
    seed &= ( 1ULL <<  ) - ;
    seed *= seed;
    return seed * Z;
}
int main()
{
    int T,n;
    double xi,xi2,y,x;
    scanf("%d",&T);
    for(int t=;t<=T;t++)
    {
        xi=xi2=;
        scanf("%d%llu",&n,&seed);
        //cin>>n>>seed;
        for(int i=;i<=n;i++)
        {
            y=gen();
            xi+=y;
            xi2+=y*y;
        }
        x=xi/n;
        printf("Case #%d: %.5lf\n",t,sqrt((xi2-*x*xi+n*x*x)/n));
    }
}