[seqfan] Re: Possible error in PARI program in A259234?

[aai]

Hi Felix,

I hope it's not too late.

a(490) = 4
with the following PARI program

v=vector(0);for(i=1,500,b=2;p=prime(i); while(Mod(b,p^2)^(p-1)!=1 || 
#select(x->(x==b),v)>0,b++); v=concat(v,[b]);print1(v[#v],", "))



Felix Fröhlich schreef op 15-7-2015 om 19:31:
> Dear Sequence fans
>
> I just discovered a possible error that (unfortunately) I didn't catch
> earlier. According to Michel's comment in A243905, A259234(183) should be 2
> (which is the case). So I wondered what the value of A259234(490) is. I was
> surprised to find in the b-file, A259234(490) = 2 as well, so clearly there
> must be something wrong with my PARI program.
>
> What is the error in the PARI code? I am willing to fix the program and
> compute a new b-file once I understand what the error is.
>
> Regards
> Felix
>
> _______________________________________________
>
> Seqfan Mailing list - http://list.seqfan.eu/