应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » Euler 工程 第廿七题：系数的积
161/2返回列表12下一页
查看: 1423  |  回复: 15
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

holmescn

金虫 (正式写手)

[交流] Euler 工程 第廿七题：系数的积 已有4人参与

Euler大牛给出了一个很牛的二次公式:
引用回帖:
n^2 + n + 41

这个很牛的公式, 当n 从0取到39的时候,能给出40个质数. 可是当n=40的时候,就失灵了.

使用计算机, 我们又得到一个更牛的公式
引用回帖:
n^2 - 79n + 1601

这哥们,当n从0取到79的时候,能给出80个质数.

如果我们定义这样的一个二次公式: n^2 + an + b
a 和 b 的绝对值都小于1000, 当这个公式能产生最多的质数的时候, 给出a和b的积.

致歉：
开始译的时候，我理解错了，结果给出错误的表述，让大家产生了误解，在这里说声对不起了。

[ Last edited by holmescn on 2011-6-16 at 19:54 ]
回复此楼

» 猜你喜欢
1楼 2011-06-16 09:35:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sudo

木虫 (正式写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
余泽成(金币+3): 鼓励交流！ 2011-06-18 16:07:37
“产生最多的质数”这个说法有点模糊呢，看例子，是不是指n从0开始取，然后递增1，直到式子n^2 + an + b不再为质数为止，这个过程中n的个数呢？

然后那个80个质数的例子是暗示一个上限吗：

80^2 + 1000*80 + 1000 = 87400 （使用的质数表中，最大的质数小于这个数~）
赞一下(2人) 回复此楼
2楼2011-06-16 10:43:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huycwork

金虫 (著名写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
余泽成(金币+3): 鼓励交流！ 2011-06-18 16:07:47
引用回帖:
Originally posted by sudo at 2011-06-16 10:43:17:
“产生最多的质数”这个说法有点模糊呢，看例子，是不是指n从0开始取，然后递增1，直到式子n^2 + an + b不再为质数为止，这个过程中n的个数呢？

然后那个80个质数的例子是暗示一个上限吗：

80^2 + 1000*8 ...

应该是没有什么暗示的吧。要找的是从[0~x)自然数区间映射到素数空间的一个函数映射f(n)=n(n+a)+b，要求0~x这个区间最长。

a取正数的时候n+a肯定不能超过b，x的取值就是0~(b-a)，a取负数的时候似乎只能达到|a|，函数形状是对称的，能到达|a|纯属巧合，真正的产生素数的部分是0~|a/2|这个部分，x所在的区间应该是0~|a|。不过再往下也不是没可能，最可靠的估计还是0~b。

a的搜索区间是-1000~1000，b的搜索区间则是0~1000内的素数，算法看起来需要O(n*n/Inn)的复杂度，多项式时间可解的搜索问题吧。

[ Last edited by huycwork on 2011-6-16 at 12:26 ]
赞一下(2人) 回复此楼
漩涡的中心有一块空地，空空的。
3楼2011-06-16 11:50:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军
★ ★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
余泽成(金币+5): 鼓励交流！ 2011-06-18 16:08:02
这个只能暴力解吧,
CODE:
#include
#include
#include

bool isPrime(int n)
{
        int i;
        bool flag = true;
        for(i=2;i         {
                if(n%i==0)
                {
                        flag = false;
                        break;
                }
        }
        return flag;
}

int main(int args,char* argv[])
{
        int a=0,b=0,i,j;
        int maxlen=0, curlen,n;

        for(i=-999;i<1000;i++)
        {
                for(j=-999;j<1000;j++)
                {
                        if(!isPrime(j)) // n^2+a*n+b, b must be a prime while n==0
                                continue;

                        curlen = 1; // n==0
                        for(n=1;n<79;n++)
                        {
                                if(!isPrime(n*n+i*n+j))
                                        break;
                                curlen += 1;
                        }

                        if(curlen>maxlen)
                        {
                                maxlen = curlen;
                                a = i;
                                b = j;
                        }
                }
        }

        printf("While %d*%d=%d, (n^2+(%d)*n+%d) produces %d primes.\n",a,b,a*b,a,b,maxlen);

        return 0;
}

结果
CODE:
% While -61*971=-59231, (n^2+(-61)*n+971) produces 72 primes.
% Elapsed time is 1.578 seconds.

[ Last edited by libralibra on 2011-6-16 at 16:46 ]
赞一下(2人) 回复此楼
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
4楼2011-06-16 16:23:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)
★ ★ ★
余泽成(金币+3): 鼓励交流！ 2011-06-18 16:08:15
我还以为表达清楚了呢. 其实就是n从0取到m,这m+1个数都是质数.然后看a*b等于多少.

Matlab版的穷举法:
CODE:
tic
maxn = 0;
maxp = [0 0];
for a = -1000:1000
    for b = -1000:1000
        n = 0;
        while (n^2 + a*n + b) > 0 && isprime(n^2 + a*n +b)
            n = n + 1;
        end
        if n > maxn
            maxn = n;
            maxp = [a b];
            fprintf('maxn=%d\n', maxn);
        end
    end
end
fprintf('a=%d,b=%d, a*b=%d\n', maxp(1), maxp(2), maxp(1)*maxp(2));
toc

结果:
引用回帖:
a = -61,  b = 971,  a*b = -59231
用时 150 秒. 共有72个质数

[ Last edited by holmescn on 2011-6-16 at 16:37 ]
赞一下(1人) 回复此楼
5楼2011-06-16 16:36:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huycwork

金虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
Originally posted by sudo at 2011-06-16 10:43:17:
“产生最多的质数”这个说法有点模糊呢，看例子，是不是指n从0开始取，然后递增1，直到式子n^2 + an + b不再为质数为止，这个过程中n的个数呢？

然后那个80个质数的例子是暗示一个上限吗：

80^2 + 1000*8 ...

很高的洞察力哈~佩服佩服~
阁下莫非是根据1601超过1000判断这个上限的？
赞一下(1人) 回复此楼
漩涡的中心有一块空地，空空的。
6楼2011-06-16 16:49:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)
★ ★
余泽成(金币+2): 辛苦了！ 2011-06-18 16:08:40
修改以后, 假设b是质数,这样时间变成原来的三分之一了
CODE:
tic
bprime = primes(1000);
maxn = 0;
maxp = [0 0];
for a = -1000:1000
    for i = 1:length(bprime)
        n = 0;
        b = bprime(i);
        while (n^2 + a*n + b) > 0 && isprime(n^2 + a*n +b)
            n = n + 1;
        end
        if n > maxn
            maxn = n;
            maxp = [a b];
        end
    end
end
fprintf('a=%d,b=%d, a*b=%d\n', maxp(1), maxp(2), maxp(1)*maxp(2));
toc
赞一下(1人) 回复此楼
7楼2011-06-16 16:50:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

holmescn

金虫 (正式写手)
引用回帖:
Originally posted by libralibra at 2011-06-16 16:23:25:
这个只能暴力解吧,

[code] #include <stdio.h>
#include <stdlib.h>
#include <math.h>

bool isPrime(int n)
{
        int i;
        bool flag = true;
        for(i=2;i<sqrt(abs(n));i++)
        {
...

如果先给b生成一个质数表呢, 我想还能再快, 可能不到1秒了.
赞一下 回复此楼
8楼2011-06-16 16:55:23
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fatpig8832

铁杆木虫 (著名写手)
★
小木虫(金币+0.5):给个红包，谢谢回帖
汗，题目写得太令人费解了...

n^2 + n + 41 能给出前40个质数

我还以为是指从2开始的40个质数呢...后来一看不是，原来是从41开始的40个质数...后来一看又不是，原来，是指 n^2 + n + 41 这条公式给出的前40个解都是质数...
赞一下(1人) 回复此楼
9楼2011-06-16 17:35:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

libralibra

至尊木虫 (著名写手)

骠骑将军
★
小木虫(金币+0.5):给个红包，谢谢回帖
引用回帖:
Originally posted by fatpig8832 at 2011-06-16 17:35:39:
汗，题目写得太令人费解了...

n^2 + n + 41 能给出前40个质数

我还以为是指从2开始的40个质数呢...后来一看不是，原来是从41开始的40个质数...后来一看又不是，原来，是指 n^2 + n + 41 这条公式给出的前40 ...

哈哈,这个题目英文说明很清楚的,

Euler published the remarkable quadratic formula:

n² + n + 41

It turns out that the formula will produce 40 primes for the consecutive values n = 0 to 39. However, when n = 40, 402 + 40 + 41 = 40(40 + 1) + 41 is divisible by 41, and certainly when n = 41, 41² + 41 + 41 is clearly divisible by 41.

Using computers, the incredible formula  n² − 79n + 1601 was discovered, which produces 80 primes for the consecutive values n = 0 to 79. The product of the coefficients, −79 and 1601, is −126479.

Considering quadratics of the form:

    n² + an + b, where |a| < 1000 and |b| < 1000

    where |n| is the modulus/absolute value of n
    e.g. |11| = 11 and |−4| = 4

Find the product of the coefficients, a and b, for the quadratic expression that produces the maximum number of primes for consecutive values of n, starting with n = 0.
赞一下(1人) 回复此楼
matlab/VB/python/c++/Java写程序请发QQ邮件:790404545@qq.com
10楼2011-06-16 17:41:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 holmescn 的主题更新
161/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 317求调剂 +4 申子申申 2026-03-19 8/400 2026-03-20 11:20 by 申子申申
[考研] 材料080500调剂求收留 +6 一颗meteor 2026-03-13 6/300 2026-03-20 10:41 by EBSD
[考研] 一志愿武汉理工材料工程专硕调剂 +7 Doleres 2026-03-19 7/350 2026-03-20 10:39 by EBSD
[考研] 324分 085600材料化工求调剂 +3 llllkkkhh 2026-03-18 3/150 2026-03-19 14:22 by houyaoxu
[考研] 求调剂，一志愿:南京航空航天大学大学 ，080500材料科学与工程学硕，总分289分 +3 @taotao 2026-03-19 3/150 2026-03-19 14:07 by peike
[考研] 一志愿天大材料与化工（085600）总分338 +5 蔡大美女 2026-03-13 5/250 2026-03-19 10:44 by 是小刘呀～
[考研] 0703化学调剂，求各位老师收留 +10 秋有木北 2026-03-14 10/500 2026-03-19 05:52 by anny19840123
[考研] 304求调剂 +6 司空. 2026-03-18 6/300 2026-03-18 23:03 by 星空星月
[考研] 311求调剂 +4 冬十三 2026-03-18 4/200 2026-03-18 21:47 by 尽舜尧1
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 08工科 320总分 求调剂 +5 梨花珞晚风 2026-03-17 5/250 2026-03-18 14:49 by haxia
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 0703化学调剂 +3 妮妮ninicgb 2026-03-17 3/150 2026-03-18 10:29 by macy2011
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考研] 东南大学364求调剂 +5 JasonYuiui 2026-03-15 5/250 2026-03-16 21:28 by 木瓜膏
[考研] 318求调剂 +3 Yanyali 2026-03-15 3/150 2026-03-16 16:41 by houyaoxu
[考研] 0703 物理化学调剂 +3 我可以上岸的对� 2026-03-13 5/250 2026-03-16 10:50 by 我可以上岸的对�
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 求老师收留调剂 +4 jiang姜66 2026-03-14 5/250 2026-03-15 20:11 by Winj1e
信息提示
关闭
请填处理意见
关闭