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

24小时热门版块排行榜    

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

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 个回答

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的回帖

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的回帖
查看全部 16 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂288 +6 ioodiiij 2026-04-10 8/400 2026-04-10 21:07 by zhouxiaoyu
[考研] 085410-273求调剂 +4 X1999 2026-04-10 4/200 2026-04-10 15:29 by hemengdong
[考研] 0702物理学学硕299求调剂 +6 祁柒连 2026-04-06 6/300 2026-04-10 11:10 by Roomoo
[考研] 机械专硕270求调剂，接受跨专业 +10 老师看看我吧aba 2026-04-09 11/550 2026-04-10 10:41 by cal0306
[考研] 一志愿双非085400电子信息344 求调剂，对材料和化学方向也感兴趣 +8 无情的小羊 2026-04-09 9/450 2026-04-10 09:30 by 松花缸1201
[考研] 电子信息270求调剂 +11 terminal469 2026-04-07 11/550 2026-04-09 18:28 by hy861222
[论文投稿] 求助文献原文 10+3 18500821399 2026-04-08 3/150 2026-04-09 16:56 by 北京莱茵润色
[考研] 085400电子信息类（川大控制工程）求调剂可跨专业 求老师联系 +3 626776879 2026-04-08 3/150 2026-04-09 16:05 by 猪会飞
[考研] 322求调剂，08工科 +3 今天是个小号 2026-04-08 3/150 2026-04-09 15:53 by wp06
[考研] 070300化学学硕311分求调剂 +18 梁富贵险中求 2026-04-04 20/1000 2026-04-09 11:18 by 哒哒哒呱呱呱
[考博] 申博 +8 IQwQl 2026-04-04 8/400 2026-04-08 09:43 by 0608104024
[考研] 计算机408｜在校多次国家级竞赛获奖｜申请调剂 +4 东山大白鹅 2026-04-05 4/200 2026-04-08 00:18 by chongya
[考研] 388求调剂 +6 四川王涛 2026-04-07 8/400 2026-04-08 00:17 by JourneyLucky
[考研] 调剂求助（生物与医药） +6 @6952 2026-04-06 6/300 2026-04-07 23:52 by lys0704
[考研] 机械调剂 +3 zzzbcb 2026-04-07 3/150 2026-04-07 22:19 by hemengdong
[考研] 生物与医药求调剂 +7 heguanhua 2026-04-05 8/400 2026-04-06 18:41 by macy2011
[考研] 326求调剂 +3 顾若浮生 2026-04-05 3/150 2026-04-05 18:32 by 蓝云思雨
[考研] 一志愿北交大材料工程总分358求调剂 +6 cs0106 2026-04-05 6/300 2026-04-05 16:34 by imissbao
[考研] 085701求调剂 +7 龚禹铭 2026-04-04 8/400 2026-04-04 13:49 by 小小树2024
[考研] 求生物学专业调剂-332分 +5 云朵遛弯指南 2026-04-04 5/250 2026-04-04 10:05 by rzh123456
信息提示
关闭
请填处理意见
关闭