- 最后登录
- 2017-10-10
- 在线时间
- 88 小时
- 阅读权限
- 20
- 注册时间
- 2008-3-19
- 积分
- 421
- 帖子
- 233
- 精华
- 2
- UID
- 25681
- 性别
- 保密

- 积分
- 421
- 帖子
- 233
- 精华
- 2
- UID
- 25681
- 性别
- 保密
|
再来提供一点线索,以及一个比较理想的近似公式,呵呵!
1/N*(S-1)/(m-1) (m=S~N之和)
在N比较大时,以Pmin=(S-1)/N*ln(N/(S-1))为下界;以Pmax=(S-1)/N*ln(N/(S-2))为上界。
再给出一个近似计算公式,刚才验证了一下,可以保证到N=500000,准确率100%,可以按照这个公式的思路继续往下构造,保证公式的100%准确。
S=floor(N/exp(1)+1.31605935892)-
floor(1/(abs(N-97)*abs(N-24586)+1))+
floor(1/(abs(N-443900)+1));
如果认为0的0次方为1(有些计算软件如此规定,不过通常认为无意义,或其值不定),公式还可以写为:
S=floor(N/exp(1)+1.31605935892)-0^(abs(N-97)*abs(N-24586))+0^abs(N-443900);
上述公式是编程写法,写成数学形式为:
S=[ N/e+1.31605935892]-[1/(|N-97|*|N-24586|+1)]+ [1/(|N-443900|+1)];
或S=[ N/e+1.31605935892]-0^(|N-97|*|N-24586|)+ 0^|N-443900|; |
|