997651

#!/usr/bin/env python

def generate_primes(n):
    isprimes=[True]*n
    for i in range(2,n):
        if isprimes[i]:
            for j in range(2*i,n,i):
                isprimes[j] = False
    primes=[i for i in range(3,n,2) if isprimes[i]]
    primes.insert(0,2)
    return primes


def euler50(num):
    p=generate_primes(num)
    max_n=1
    while sum(p[:max_n])         max_n+=1
   
    for i in range(max_n,2,-1):
        for j in range(len(p)-i):
            tmp=sum(p[j:j+i])
            if tmp>num:break
            if tmp in p:
                return tmp
        
if __name__ == "__main__":
    print euler50(1000000)

% Elapsed time is 0.083633 seconds.
% ans =
%       997651
function result = euler50()
tic;
result = 0;
maxLen = 0;

ap = primes(1e6);

sumlist = cumsum(ap);
sumlist(sumlist>=1e6) = [];

maxLen = 21;
for i=1:length(sumlist)
    for j=i+maxLen:length(sumlist)
        temp = sumlist(j)-sumlist(i);
        if isprime(temp) && (j-i)>maxLen
            maxLen = j-i;
            result = temp;
        end
    end
end
toc;
end