[seqfan] Re: complementing A121558

[Ray Chandler]

> Richard Mathar schrieb:
> > Values are missing after the 3313 = A121558(16).
> > Can someone submit a correction?
> > http://http://research.att.com/~njas/sequences/A121558 .
> >
> > Richard

> It even seems that the terms after 3313 are plain wrong.
> Or did I miss something in the following pari computation:
> 
> P=prod(k=1, 1000, 1+x^k, 1+x*O(x^1000));
> for(n=1,1000,t=n+polcoeff(P,n);if(isprime(t),print1(t,",")))
> 
> 2,3,5,17,23,31,73,97,127,167,191,797,1021,1153,1301,3313,9851,
> 40099,206939,
> 345953,444893,570181,2032409,10327297,15757649,40982701,695455
> 27,281138239,657667789,
> 1333640927,32173820179,48626519377,53868949807,495811829089,41
> 328689179429,
> 232282884905431,406202130845981,2608327191290453,15475586776823633,
> 27574512158560853,36704653322668189,69395901032399983,91840127
> 431010423,
> 105585005382256889,130050574409983361,476517655182416039,87073
> 4084093744939,
> 6162297935619708167,8475550903782346751,23192540574062407799,8
> 4421875547541660201,
> 604490895403729930283,3819342674540204978827,
> 
> Apparently, 281641, 451331, 966527, ... do not occur.
> 
> 
> Hagen

I confirm your terms with the following Mathematica code:
Select[Table[n+PartitionsQ[n],{n,1000}],PrimeQ]

None of the terms following 3313 in the original sequence belong.
Ray