[go: up one dir, main page]

login

Year-end appeal: Please make a donation to the OEIS Foundation to support ongoing development and maintenance of the OEIS. We are now in our 61st year, we have over 378,000 sequences, and we’ve reached 11,000 citations (which often say “discovered thanks to the OEIS”).

Search: a107448 -id:a107448
     Sort: relevance | references | number | modified | created      Format: long | short | data
Least p for which min { x >= 0 | p + (2n+1)*x + x^2 is composite } reaches the (local) maximum given in A273770.
+10
4
41, 43, 47, 53, 61, 71, 83, 97, 113, 131, 151, 173, 197, 223, 251, 281, 313, 347, 383, 421, 461, 503, 547, 593, 641, 691, 743, 797, 73303, 73361, 73421, 73483, 3443897, 3071069, 3071137, 15949847, 76553693, 365462323, 365462399, 2204597, 9721, 1842719, 246407633, 246407719, 246407807, 246407897, 246407989
OFFSET
0,1
COMMENTS
All terms are prime, since this is necessary and sufficient to get a prime for x = 0.
The values given in A273770 are the number of consecutive primes obtained for x = 0, 1, 2, ....
Sequence A273595 is the subsequence of terms for which 2n+1 is prime.
For even coefficients of the linear term, the answer would always be q=2, the only choice that yields a prime for x=0 and also for x=1 if (coefficient of the linear term)+3 is prime.
The initial term a(n=0) = 41 corresponds to Euler's famous prime-generating polynomial 41+x+x^2. Some subsequent terms are equal to the primes this polynomial takes for x=1,2,3,.... This stems from the fact that adding 2 to the coefficient of the linear term is equivalent to shifting the x-variable by 1. Since here we require x >= 0, we find a reduced subset of the previous sequence of primes, missing the first one, starting with q equal to the second one. (It is known that there is no better prime-generating polynomial of this form than Euler's, see the MathWorld page and A014556. "Better" means a larger p producing p-1 primes in a row. However, the prime k-tuple conjecture suggests that there should be arbitrarily long runs of primes of this form (for much larger p), i.e., longer than 41, but certainly much less than the respective p. Therefore we speak of local maxima.)
LINKS
Eric Weisstein's World of Mathematics, Prime-Generating Polynomial
PROG
(PARI) A273756(n, p=2*n+1, L=10^(5+n\10), m=0, Q)={forprime(q=1, L, for(x=1, oo, ispseudoprime(q+p*x+x^2)&& next; x>m&& [Q=q, m=x]; break)); Q}
CROSSREFS
Cf. also A002837 (n such that n^2-n+41 is prime), A007634 (n such that n^2+n+41 is composite), A005846 (primes of form n^2+n+41), A097823, A144051, A187057 ... A187060, A190800, A191456 ff.
The first line of data coincides with that of A202018, A107448, A155884 (and also A140755, A142719, except for some initial terms), which are all related.
KEYWORD
nonn
AUTHOR
M. F. Hasler, May 26 2016
EXTENSIONS
Edited, following a remark by Don Reble, by M. F. Hasler, Jan 23 2018
a(27) corrected and more terms from Don Reble, Feb 15 2018
STATUS
approved
Irregular triangle T(n, k) = 10 - ( (b(n) + k^2 + k + 1) mod 10 ), where b(n) = A056486(n-1) - (1/2)*[n=1], for n >= 1 and 1 <= k <= b(n) - 1, read by rows.
+10
2
5, 3, 9, 3, 7, 3, 7, 9, 9, 7, 3, 7, 9, 1, 7, 1, 3, 3, 1, 7, 1, 3, 3, 1, 7, 1, 3, 3, 7, 3, 7, 9, 9, 7, 3, 7, 9, 9, 7, 3, 7, 9, 9, 7, 3, 7, 9, 9, 7, 3, 7, 9, 9, 7, 3, 7, 9, 9, 7, 3, 7, 9, 9, 7, 3, 7, 9, 3, 9, 3, 5, 5, 3, 9, 3, 5, 5, 3, 9, 3, 5, 5, 3, 9, 3, 5, 5, 3, 9, 3, 5, 5, 3, 9, 3, 5, 5, 3, 9, 3, 5, 5, 3, 9, 3
OFFSET
1,1
LINKS
FORMULA
T(n, k) = 10 - (b(n) + k^2 + k + 1) mod 10, where b(n) = A056486(n-1) - (1/2)*[n=1], for n >= 1 and 1 <= k <= b(n) - 1. - G. C. Greubel, Mar 24 2024
EXAMPLE
The irregular triangle begins as:
5;
3, 9, 3;
7, 3, 7, 9, 9, 7, 3, 7, 9;
1, 7, 1, 3, 3, 1, 7, 1, 3, 3, 1, 7, 1, 3, 3;
MATHEMATICA
b[n_]:= 2^(n-3)*(9-(-1)^n) -Boole[n==1]/2;
T[n_, k_]:= 10 -Mod[k^2+k+1+b[n], 10];
Table[T[n, k], {n, 8}, {k, b[n]-1}]//Flatten (* G. C. Greubel, Mar 24 2024 *)
PROG
(Magma)
b:= func< n | n eq 1 select 2 else 2^(n-3)*(9-(-1)^n) >;
A107448:= func< n, k | 10 - ((b(n) +k^2 +k +1) mod 10) >;
[5, 3, 9, 3] cat [A107448(n, k): k in [1..b(n)-1], n in [3..8]]; // G. C. Greubel, Mar 24 2024
(SageMath)
def b(n): return 2^(n-3)*(9-(-1)^n) - int(n==1)/2
def A107449(n, k): return 10 - ((b(n) + k^2+k+1)%10);
flatten([[A107449(n, k) for k in range(1, b(n))] for n in range(1, 8)]) # G. C. Greubel, Mar 24 2024
CROSSREFS
KEYWORD
nonn,less
AUTHOR
Roger L. Bagula, May 26 2005
EXTENSIONS
Edited by G. C. Greubel, Mar 24 2024
STATUS
approved

Search completed in 0.008 seconds