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

n^2 + n + 41

n^2 - 79n + 1601

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个质数

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

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++)
        {
...