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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » 其它 » Euler 工程 第廿七题：系数的积
51/1返回列表
查看: 1459  |  回复: 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的回帖

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的回帖
查看全部 16 个回答

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的回帖
查看全部 16 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料与化工调剂 10+11 下一站上岸@ 2026-04-10 33/1650 2026-04-11 00:27 by success5733
[考研] 284求调剂 +19 梵@@ 2026-04-06 21/1050 2026-04-10 21:12 by zhouxiaoyu
[考研] 284求调剂 +9 让我上岸吧阿西 2026-04-09 11/550 2026-04-10 19:18 by 靖jing
[考研] 22408调剂求助 +5 毂12 2026-04-09 7/350 2026-04-10 16:32 by 高维春
[考研] 一志愿京区985，085401电子信息，本科电子信息 +3 阳光开朗的男孩 2026-04-10 3/150 2026-04-10 16:29 by sophia_93
[考研] 085402通信工程调剂，有4项学科竞赛国奖（电赛国二），硕士研究生调剂自荐信。 +4 m永o不v言o弃m 2026-04-09 4/200 2026-04-10 11:30 by asy1wn
[考研] 一志愿中南大学物理学，英一66，求调剂 +4 长烟旖旎 2026-04-08 5/250 2026-04-10 10:31 by 颖果儿
[考研] 材料专硕283求调剂 +18 试试看呗 2026-04-04 19/950 2026-04-10 10:05 by may_新宇
[考研] 材料调剂 +10 18815505510 2026-04-09 11/550 2026-04-09 17:07 by 544594351
[论文投稿] 求助文献原文 10+3 18500821399 2026-04-08 3/150 2026-04-09 16:56 by 北京莱茵润色
[考研] 353求调剂 +8 晴空万里air 2026-04-07 8/400 2026-04-09 00:18 by GouQ
[考研] 机械专硕273请求调剂 +6 庚申壬申 2026-04-07 6/300 2026-04-08 22:41 by bljnqdcc
[考研] 一志愿吉大化学327求调剂 +12 王王白石 2026-04-06 13/650 2026-04-08 16:05 by luoyongfeng
[考研] 265求调剂 +19 小木虫085600 2026-04-06 21/1050 2026-04-08 10:38 by 逆水乘风
[考研] 316求调剂 +4 15318418673 2026-04-07 4/200 2026-04-07 22:12 by hemengdong
[考研] 软工学硕299求调剂 +6 useryy 2026-04-07 6/300 2026-04-07 09:50 by vgtyfty
[考研] 285求调剂 +15 哦呦呼o 2026-04-04 17/850 2026-04-06 23:02 by chenzhimin
[考研] 考研调剂 +3 Wwwwwww哇 2026-04-06 3/150 2026-04-06 20:55 by lbsjt
[考研] 285求调剂 +5 mapmath 2026-04-06 6/300 2026-04-06 17:18 by 蓝云思雨
[考研] 本科211，专业085404，293分请求调剂 +5 莲菜就是藕吧 2026-04-04 5/250 2026-04-04 14:08 by 这是一个无聊的�
信息提示
关闭
请填处理意见
关闭