“蔚来杯“2022牛客暑期多校训练营4-N-Particle Arts

"蔚来杯"2022牛客暑期多校训练营4-N-Particle Arts

原题题面：https://ac.nowcoder.com/acm/contest/33189/N

题目大意

空间中有$n(2\le n\le 10^5)$个点，每个点具有热量$a_i(0\le a_i\le 2^{15})$焦耳，每两个点可能随机碰撞，如果两个带有能量$a$焦耳与$b$焦耳的粒子碰撞，将形成$a|b$焦耳与$a$&$b$焦耳
两个能量的粒子。
已证这$n$个粒子的能量方差终会收敛，求最终收敛的方差。
求方差公式如下：
$${\sigma }^2=\frac{1}{n}\sum _{i=1}^n(x_i-\mu {)}^2$$
而
$$\mu =\frac{1}{n}\sum _{i=1}^n{x}_i$$
最终结果以$a/b$的最简分数形式输出(当$a=0$时，$b=1$)。

解题思路

发现两粒子碰撞后，二进制位为$1$的趋向于能量为$a|b$的粒子，二进制位为$0$的趋向于能量为$a\&b$的粒子，经过有限次碰撞，$1$不断进入最大能量的粒子，$0$不断进入携带较小能量的粒子
如图：
$\begin{array}{ccccc}1& 0& 1& 1& \\ 0& 0& 1& 1& \\ 0& 1& 0& 0& \\ 1& 0& 1& 1& \\ 0& 0& 0& 1& \end{array}$
变为：
$\begin{array}{ccccc}0& 0& 0& 0& \\ 0& 0& 0& 1& \\ 0& 0& 1& 1& \\ 1& 0& 1& 1& \\ 1& 1& 1& 1& \end{array}$
显然，只有任意两个粒子二进制都具有包含关系，碰撞时才不会使二进制位交换，此时能量值的方差收敛。

代码实现

代码只需构造出如上关系即可。

#include<bits/stdc++.h>
#define ll long long
using namespace std;
const int N=1e5+5;
ll n,a[N],b[100],m;
ll gcd(ll a,ll b){
    if(a%b==0)return b;
    return gcd(b,a%b);
}
int main(){
    cin>>n;
    for(int i=1;i<=n;i++){
        cin>>a[i];
        ll t=0;
        while(a[i]){
            t++;
            b[t]+=a[i]&1;
            m=max(m,t);
            a[i]>>=1;
        }
    }
    for(int i=1;i<=m;i++){
        while(b[i]){
            a[b[i]]|=1<<i-1;
            b[i]--;
        }
    }
    ll a1=0,a2=0;
    for(int i=1;i<=n;i++){
        a1+=a[i]*a[i],a2+=a[i];
    }
    ll ans=a1*n-a2*a2;
    cout<<ans/gcd(ans,n*n)<<'/'<<n*n/gcd(ans,n*n);
}

另

末尾所用公式为：

$\begin{array}{rl}设一个数组的& 方差为D(x)，平均数为E(x)。\\ 则D(x)& =E(x^2)-E^2(x)\\ 证明：D(x)& =\frac{1}{n}\sum _{i=1}^n(x_i-E(x){)}^2\\ 展开可得D(x)& =\frac{1}{n}\sum _{i=1}^n(x_i^2-2x_{i}E(x)+E^2(x))\\ & =\frac{1}{n}\times n\times E^2(x)+\frac{1}{n}\times \sum _{i=1}^n(x_i^2-2x_{i}E(x))\\ & =E^2(x)+\frac{1}{n}\times \sum _{i=1}^n(x_i^2-2x_{i}E(x))\\ & =E^2(x)+E(x^2)-\frac{1}{n}\times \sum _{i=1}^{n}2x_{i}E(x)\\ & =E^2(x)+E(x^2)-2\times E^2(x)\\ & =E(x^2)-E^2(x)\\ & 证毕。\end{array}$