any suitable sieve written in C or Python?

[Shen]

It is well known that if q | 2^p - 1, then q = 2px + 1 and q = 8r +(-) 1.

when I do some research on Fp (Fibonacci prime), I guess :
1: if q | Fp, then q = 2px + (-1)^x
2: F(2n+1) must have a prime factor q, q = 2(2n+1)x + (-1)^x
3: ...

At this moment , I try to factor C261 ( F1565 = F5* F313 * C261), which has 261 digital. If guess is right, then C261 must have a prime factor = 2*1565*x + (-1)^x.

The GIMPS factoring code creates a modified sieve of Eratosthenes with each bit representing a potential 2kp+1 factor. The sieve then eliminates any potential factors that are divisible by prime numbers below 40,000 or so. Also, bits representing potential factors of 3 or 5 mod 8 are cleared. This process eliminates roughly 95% of potential factors.

Since I don't know assemble language, I could not understand such code.
I tried to code it by myself, but the code efficiency is very low.

I have no other good choice but ask for help :
Could someone provide a efficient sieve ?

Thanks a lot!

[carpetpool]

Didn't you get the idea from here?

[axn]

Originally Posted by Shen

At this moment , I try to factor C261 ( F1565 = F5* F313 * C261), which has 261 digital. If guess is right, then C261 must have a prime factor = 2*1565*x + (-1)^x.

Trial factoring is the wrong tool for numbers this size. You should use ECM.

[carpetpool]

Originally Posted by axn

Trial factoring is the wrong tool for numbers this size. You should use ECM.

YAFU?

[Shen]

Thanks for your advice.