找回密码
 欢迎注册
查看: 16478|回复: 8

[分享] superabundant问题的数学及算法描述

[复制链接]
发表于 2008-1-9 15:52:46 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?欢迎注册

×
精华
问题转自 http://topic.csdn.net/u/20080103 ... html?seed=157425768
The abundance of an integer n is the sum of the divisors of n (including n itself), divided
by n. Integer n is k-abundant if its abundance is at least k.

For example, the sum of the divisors of 6 is 6+3+2+1=12, and 12/6=2, so 6 is 2-abundant.
As another example, the sum of the divisors of 120 is

120+60+40+30+24+20+15+12+10+8+6+5+4+3+2+1=360

so 120 is 3-abundant. It happens that 6 is the smallest 2-abundant number and 120 is the smallest
3-abundant number. They happen to be exactly 2- and 3-abundant, respectively, but it is generally
possible that the smallest k-abundant number has abundance greater than k.

The task is to write a program that finds the smallest k-abundant number for k=1,2,..
How high can you go?

附件中是我给的相关数学分析以及算法描述 abundance.zip (64.43 KB, 下载次数: 6)
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2008-1-9 17:32:18 | 显示全部楼层
一个更新版本,修改了一些错误,
并且删除了使用拉格朗日极值法求极值,而是使用琴生不等式
abundance.zip (70.46 KB, 下载次数: 25)
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2008-1-10 07:52:33 | 显示全部楼层
可以对原贴直接进行编辑,附件也是可以重新上传更新的(必要时提醒网友重新下载即可)。
这个论坛功能够强吧

上传附件时,最好在“描述”栏中写些摘要或说明;
如果附件是图片,则有可能被首页自动调用显示,具体见:关于“首页调用”图片的几点说明
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2008-2-5 15:46:36 | 显示全部楼层
CSDN中的yaos提出了一个相关的问题,就是如何找出一些abundance是整数的数,当然是越多越好。
http://topic.csdn.net/u/20080204 ... 3-90ce2d8f2c8d.html
这个问题我们可以描述为给定一个整数K,请找出尽量多的大整数使得它的abundance正好是K.
关于这个问题,我有一些想法,觉得使用以后用计算机应该可以将搜索范围裁剪得很小,但是没有真正试验过。
游客,如果您要查看本帖隐藏内容请回复
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2008-3-7 11:06:17 | 显示全部楼层
很难很难

至今还无法想到如何做
关键是起始方幂无法确定
一旦错误, 就得不到结果
而且2, 3, 5, 7的幂都是基本关联的
大的2 ^ n不会对应小的3 ^ m

估计试探法对高次幂成功概率不会超过1/10
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
 楼主| 发表于 2008-3-7 12:47:22 | 显示全部楼层
1/10对于计算机来说是个很大的概率了。
比如计算一个素数的原根,通常就用试探法,而理论上试探成功率下限大概在$1/log(p)$,但是算法已经非常有效了。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2009-8-20 13:02:56 | 显示全部楼层
一个更新版本,修改了一些错误,
并且删除了使用拉格朗日极值法求极值,而是使用琴生不等式
15
mathe 发表于 2008-1-9 17:32

下载链接失效了。
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-8-22 21:44:37 | 显示全部楼层
同一个问题,就是下载链接失效了,另外就是想请教下mathe
  一个数的重数公式:

p1-p1^(-a1)        pt - pt^(-at)
-----------      *...* --------------           是如推导得出的
  p1-1                 pt - 1
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
发表于 2013-8-24 08:09:48 | 显示全部楼层
其实这个问题还有个难点在于要在程序里集成椭圆曲线跟二次筛整数分解算法

当然可以用现成的库
毋因群疑而阻独见  毋任己意而废人言
毋私小惠而伤大体  毋借公论以快私情
您需要登录后才可以回帖 登录 | 欢迎注册

本版积分规则

小黑屋|手机版|数学研发网 ( 苏ICP备07505100号 )

GMT+8, 2024-5-7 13:39 , Processed in 0.063435 second(s), 19 queries .

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表