前置知识:条件期望。(应该是对的叭?)
记 s = ∑ a i s=\sum a_i s=∑ai。考虑一个颜色的球有 i i i个,所有球变成它颜色的概率,答案是 p i = i s p_i=\frac{i}{s} pi=si。这是因为,可以将所有球看成两两不同,最终变成每个球的颜色的概率是相同的。
可以考虑固定最后剩下的那个颜色。设 f n f_n fn表示这个颜色的球还剩 n n n个,最后所有球都变成这个颜色期望步数。答案就是 ∑ f a i \sum f_{a_i} ∑fai。
我们可以观察到,一次操作是好的,当且仅当选择了两个颜色不同的球,并且向左和向右的概率是相等的(指下标)。因此递推式如下:
f i = ( s 2 ) i ( s − i ) × i s + 1 2 ( f i − 1 + f i + 1 ) f_i=\frac{\binom{s}{2}}{i(s-i)}\times \frac{i}{s}+\frac{1}{2}(f_{i-1}+f_{i+1}) fi=i(s−i)(2s)×si+21(fi−1+fi+1)
注意计算条件期望的时候,每一步转移的代价要乘上对应的能到达目标状态的概率(感性理解)。可以认为 0 0 0和 s s s都是终止状态,我们要乘上的就是到达 s s s的概率。
记 δ i = f i + 1 − f i \delta_i=f_{i+1}-f_i δi=fi+1−fi,将上式变形可以得到:
δ i = δ i − 1 − s − 1 s − i = δ 0 − ∑ j = 1 i s − 1 s − j \delta_i=\delta_{i-1}-\frac{s-1}{s-i}=\delta_0-\sum_{j=1}^i\frac{s-1}{s-j} δi=δi−1−s−is−1=δ0−j=1∑is−js−1
根据 ∑ i = 0 s − 1 δ i = 0 \sum_{i=0}^{s-1}\delta_i=0 ∑i=0s−1δi=0,事实上已经可以把 δ 0 \delta_0 δ0解出来了。解出来 δ 0 = ( s − 1 ) 2 s \delta_0=\frac{(s-1)^2}{s} δ0=s(s−1)2,从前往后递推即可。
复杂度 O ( max ( a i ) ) O(\max(a_i)) O(max(ai))。
upd on 2024/2/6:没有看懂自己写的题解,警钟长鸣。
#include<bits/stdc++.h>
#define fi first
#define se second
#define ll long long
#define pb push_back
#define db double
#define inf 0x3f3f3f3f3f3f3f3f
using namespace std;
const int mod=1e9+7;
int n,m,a[2505];
ll s,f[100005],g[100005],res;
ll fpow(ll x,ll y){
ll z(1);
for(;y;y>>=1){
if(y&1)z=z*x%mod;
x=x*x%mod;
}return z;
}
int main(){
ios::sync_with_stdio(false);
cin.tie(0),cout.tie(0);
cin>>n;for(int i=1;i<=n;i++)cin>>a[i],m=max(m,a[i]),s+=a[i];
g[0]=(s-1)*(s-1)%mod*fpow(s,mod-2)%mod;
for(int i=1;i<m;i++)g[i]=(g[i-1]-(s-1)*fpow(s-i,mod-2)%mod)%mod;
for(int i=1;i<=m;i++)f[i]=(f[i-1]+g[i-1])%mod;
for(int i=1;i<=n;i++)res=(res+f[a[i]])%mod;
cout<<(res+mod)%mod;
}