平方数问题一则

无心人 · 发表于 2009-1-31 08:46:20

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

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

×

假设完全平方数$a$的十进制表示的各位数字和是$b$ 如果$b$也是完全平方数，我们称$a$是二重平方数如果$b$也是二重平方数，则我们称$a$是三重平方数依次类推现在请求出$10^24$内的所有二重平方数

无心人 · 发表于 2009-1-31 14:26:11

看来别人么兴趣做啊

无心人 · 发表于 2009-1-31 14:33:11

代码 let square = [s| a <- [1..16], let s = a * a] let isSquare = listArray (1, 256) [s | a <- 1..256], let s = (/= []) \$ findIndices (==a) square] let digitsSum n = sum \$ zipWith (*) [0..9] \$ map (flip (-) 1) \$ map length \$ group \$ sort \$ (show n) ++ "0123456789" let test n m = [(i, s, ss) | i <- [n..m], let s = i * i, let ss = digitsSum s, isSquare!(toInteger ss)] 得到10^12内的结果并保存 writeFile "Result.Txt" \$ show \$ test 1 1000000

无心人 · 发表于 2009-1-31 14:37:56

呵呵 10^12内耗时几秒共160264组解解的数量太多了

result24.part1.rar (439.45 KB, 下载次数: 1)

result24.part2.rar (345.27 KB, 下载次数: 1)

kofeffect · 发表于 2009-1-31 22:02:36

使用Haskell算的啊耗时多少？确实太多了你只求$10^12$内的而已

Clear[a, b, s, i, j, k, t];
k = 1;
For[i = 1, i < 10^12 + 1, i++,
a = i^2;
(*二重平方数*)
len = IntegerLength[a];
s = 0;
b = a;
For[j = len, j > 0, j--,
s += IntegerPart[a/(10^(j - 1))];
a = Mod[a, (10^(j - 1))];
];
t = IntegerPart[Sqrt[s]];
If[t^2 != s, Continue[],
PutAppend[NO : k | b | s, "2SquareNum.txt"]; k += 1;];
] // Timing
Print["No.", k]

复制代码

无心人 · 发表于 2009-1-31 22:05:00

多谢提醒，当时有点糊涂了呵呵，年纪大了，人一过10岁就有点迷糊了呵呵

kofeffect · 发表于 2009-1-31 22:06:19

不是吧，我刚回复，你就看到贴了，难道你也一直刷新帖子啊

无心人 · 发表于 2009-1-31 22:06:36

这个东西有优化算法在那个平方数字和帖子上的代码可以修改来做大概不会很长时间就能得到10^24内的

无心人 · 发表于 2009-1-31 22:07:05

Haskell代码是用最笨的算法写的

kofeffect · 发表于 2009-1-31 22:09:58

呵呵我现在还看不懂Haskell 只会mathematica穷举

账号		自动登录	找回密码
密码			欢迎注册

[讨论] 平方数问题一则

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