Factoring all 2n-1 up to 400.
addprimes([ 26986333437777017 , 61676882198695257501367 , 1155685395246619182673033 , 120226360536848498024035943 , 31133636305610209482201109050392404721 , 40122362455616221971122353 , 2647649373910205158468946067671 , 3857194764289141165278097 , 4868122671322098041565641 , 14143189112952632419639 , 34720396273212657799920861294559 , 2927455476800301964116805545194017 , 370906580744492785430299503112990447 , 301311116540899114446723859201 , 503823044204581129045587727 , 15174923558680812616818436353130417 , 814132872808522587940886856743 , 1234904213576000272542841146073 ])
addprimes([ 179058312604392742511009 , 36450568206770608791178096385783 ])
addprimes([ 51441563151591093599 , 365505823711978039310711 , 1953272766780718501831 , 39044358788825633753 , 33157029794959983067039 , 126901141805369975317583 , 6834040335349578249140287 , 25806248225716242845491832244899635927231330561 , 5560125493425335999 ])
for(q=1 , 200 , i=q ; f=factorint(2^i+1) ; print(i , " " , gettime , " " , f) ; d=matsize(f) ; for(j=1 , d[1] , if(f[j , 1]>10^15 , addprimes(f[j , 1]))))
for(q=1 , 400 , i=q ; f=factorint(2^i-1) ; print(i , " " , gettime , " " , f) ; d=matsize(f) ; for(j=1 , d[1] , if(f[j , 1]>10^15 , addprimes(f[j , 1]))))
No comments :
Post a Comment