[go: up one dir, main page]

login
Revision History for A329572 (Bold, blue-underlined text is an addition; faded, red-underlined text is a deletion.)

Showing entries 1-10 | older changes
For all n >= 0, exactly 12 sums are prime among a(n+i) + a(n+j), 0 <= i < j < 7; lexicographically earliest such sequence of distinct nonnegative numbers.
(history; published version)
#12 by Michel Marcus at Wed Feb 19 04:07:22 EST 2020
STATUS

reviewed

approved

#11 by M. F. Hasler at Sat Feb 15 08:37:41 EST 2020
STATUS

proposed

reviewed

#10 by Michel Marcus at Fri Feb 14 23:55:39 EST 2020
STATUS

editing

proposed

Discussion
Sat Feb 15
08:37
M. F. Hasler: Yes of course. (In school I was told to spell out numbers <= 12 but I agree that in math it's useful to see the digits, to I tried to provide both versions, digits in NAME & letters in COMMENT, but I'm completely happy with digits everywhere.)
#9 by Michel Marcus at Fri Feb 14 23:55:31 EST 2020
COMMENTS

That is, there are twelve 12 primes, counted with multiplicity, among the 21 pairwise sums of any 7 consecutive terms.

STATUS

proposed

editing

Discussion
Fri Feb 14
23:55
Michel Marcus: ok ?
#8 by M. F. Hasler at Fri Feb 14 09:37:46 EST 2020
STATUS

editing

proposed

#7 by M. F. Hasler at Fri Feb 14 09:37:42 EST 2020
COMMENTS

Conjectured to be a permutation of the nonnegative integers. See A329573 for the "positive" variant: same definition but with offset 1 and positive terms, leading to a quite different sequence.

#6 by M. F. Hasler at Fri Feb 14 09:34:30 EST 2020
NAME

allocated for M. F. Hasler

For all n >= 0, exactly 12 sums are prime among a(n+i) + a(n+j), 0 <= i < j < 7; lexicographically earliest such sequence of distinct nonnegative numbers.

DATA

0, 1, 2, 5, 6, 11, 12, 17, 26, 35, 36, 47, 24, 54, 77, 7, 43, 60, 13, 30, 96, 4, 67, 97, 16, 133, 34, 3, 40, 27, 63, 100, 10, 20, 171, 9, 8, 51, 21, 22, 52, 15, 32, 38, 75, 141, 56, 41, 71, 122, 152, 45, 68, 29, 59, 14, 39, 44, 50, 23, 53, 57, 74, 107, 170, 176, 93, 134, 137, 86, 177, 65, 476, 62, 87, 92, 101

OFFSET

0,3

COMMENTS

That is, there are twelve primes, counted with multiplicity, among the 21 pairwise sums of any 7 consecutive terms.

This is the theoretical maximum: there can't be more than 12 primes in pairwise sums of 7 distinct numbers > 1. See the wiki page for more details.

Conjectured to be a permutation of the nonnegative integers. See A329573 for the "positive" variant: same definition but with offset 1 and positive terms, leading to a quite different sequence.

For a(3) and a(4) resp. a(5) one must forbid the values < 5 resp. < 11 which would be the greedy choices, in order to get a solution for a(7), but from then on, the greedy choice gives the correct solution, at least for several hundred terms.

LINKS

M. F. Hasler, <a href="/wiki/User:M._F._Hasler/Prime_sums_from_neighboring_terms">Prime sums from neighboring terms</a>, OEIS wiki, Nov. 23, 2019

PROG

(PARI) {A329572(n, show=0, o=0, N=12, M=6, D=[3, 5, 4, 6, 5, 11], p=[], u=o, U)=for(n=o+1, n, show>0&& print1(o", "); show<0&& listput(L, o); U+=1<<(o-u); U>>=-u+u+=valuation(U+1, 2); p=concat(if(#p>=M, p[^1], p), o); D&& D[1]==n&& [o=D[2], D=D[3..-1]]&& next; my(c=N-sum(i=2, #p, sum(j=1, i-1, isprime(p[i]+p[j])))); for(k=u, oo, bittest(U, k-u)|| min(c-#[0|p<-p, isprime(p+k)], #p>=M)|| [o=k, break])); show&&print([u]); o} \\ optional args: show=1: print a(o..n-1), show=-1: append them on global list L, in both cases print [least unused number] at the end. See the wiki page for more.

CROSSREFS

Cf. A055273 (analog starting with a(1) = 1), A055265 & A128280 (1 prime using 2 terms), A055266 & A253074 (0 primes using 2 terms), A329405 - A329416, A329425, A329333, A329449 - A329456, A329563 - A329581.

KEYWORD

allocated

nonn,changed

AUTHOR

M. F. Hasler, Feb 09 2020

STATUS

approved

editing

#5 by M. F. Hasler at Fri Feb 14 09:22:33 EST 2020
NAME

For every n >= 0, exactly 10 sums are prime among a(n+i) + a(n+j), 0 <= i < j < 7; lexicographically earliest such sequence of distinct nonnegative numbers.

allocated for M. F. Hasler

DATA

0, 1, 2, 3, 4, 5, 8, 9, 10, 14, 33, 15, 20, 27, 26, 11, 32, 16, 41, 21, 57, 116, 22, 51, 38, 23, 50, 63, 86, 6, 17, 24, 77, 65, 18, 13, 114, 25, 36, 28, 35, 43, 12, 31, 61, 66, 40, 19, 47, 42, 90, 241, 7, 52, 37, 34, 45, 30, 55, 49, 394, 58, 73, 39, 48, 64, 109, 115

OFFSET

0,3

COMMENTS

That is, there are 10 primes, counted with multiplicity, among the 21 pairwise sums of any 7 consecutive terms.

Is this a permutation of the nonnegative integers?

If so, then the restriction to [1..oo) is a permutation of the positive integers, but maybe not the lexicographically earliest one with this property.

This is the first example of a sequence of this type for which the greedy choice of a(n) is frequently incorrect beyond the initial terms, see Examples.

LINKS

M. F. Hasler, <a href="/wiki/User:M._F._Hasler/Prime_sums_from_neighboring_terms">Prime sums from neighboring terms</a>, OEIS Wiki, Nov. 23, 2019, updated Feb. 2020.

EXAMPLE

At the beginning of the sequence, we must avoid the choice of 6 or 7 for a(6): these choices would be possible w.r.t. the constraints, but then make it impossible to find a successor term.

The same happens again for a(37) and a(58), where the value 19 resp. 46 must be avoided.

PROG

(PARI) {A329572(n, show=0, o=0, N=10, M=6, X=[[6, 6], [6, 7], [37, 19], [58, 46]], p=[], u=o, U)=for(n=o+1, n, show>0&& print1(o", "); show<0&& listput(L, o); U+=1<<(o-u); U>>=-u+u+=valuation(U+1, 2); p=concat(if(#p>=M, p[^1], p), o); my(c=N-sum(i=2, #p, sum(j=1, i-1, isprime(p[i]+p[j])))); for(k=u, oo, bittest(U, k-u)|| min(c-#[0|x<-p, isprime(x+k)], #p>=M)|| setsearch(X, [n, k])|| [o=k, break])); show&&print([u]); o} \\ optional args: show=1: print a(o..n-1), show=-1: append them on global list L, in both cases print [least unused number] at the end. Parameters N, M, o, ... allow to get other variants, see the wiki page for more.

CROSSREFS

Cf. A055265, A128280 (1 prime from 2 terms), A329333 (1 prime from 3 terms), A329405, ..., A329416 (N primes from M terms >= 1), A329425, A329449, ..., A329581 (N primes from M terms >= 0).

KEYWORD

nonn,changed

allocated

AUTHOR

M. F. Hasler, Feb 10 2020

STATUS

editing

approved

#4 by M. F. Hasler at Mon Feb 10 21:46:52 EST 2020
NAME

allocated for M. F. Hasler

For every n >= 0, exactly 10 sums are prime among a(n+i) + a(n+j), 0 <= i < j < 7; lexicographically earliest such sequence of distinct nonnegative numbers.

DATA

0, 1, 2, 3, 4, 5, 8, 9, 10, 14, 33, 15, 20, 27, 26, 11, 32, 16, 41, 21, 57, 116, 22, 51, 38, 23, 50, 63, 86, 6, 17, 24, 77, 65, 18, 13, 114, 25, 36, 28, 35, 43, 12, 31, 61, 66, 40, 19, 47, 42, 90, 241, 7, 52, 37, 34, 45, 30, 55, 49, 394, 58, 73, 39, 48, 64, 109, 115

OFFSET

0,3

COMMENTS

That is, there are 10 primes, counted with multiplicity, among the 21 pairwise sums of any 7 consecutive terms.

Is this a permutation of the nonnegative integers?

If so, then the restriction to [1..oo) is a permutation of the positive integers, but maybe not the lexicographically earliest one with this property.

This is the first example of a sequence of this type for which the greedy choice of a(n) is frequently incorrect beyond the initial terms, see Examples.

LINKS

M. F. Hasler, <a href="/wiki/User:M._F._Hasler/Prime_sums_from_neighboring_terms">Prime sums from neighboring terms</a>, OEIS Wiki, Nov. 23, 2019, updated Feb. 2020.

EXAMPLE

At the beginning of the sequence, we must avoid the choice of 6 or 7 for a(6): these choices would be possible w.r.t. the constraints, but then make it impossible to find a successor term.

The same happens again for a(37) and a(58), where the value 19 resp. 46 must be avoided.

PROG

(PARI) {A329572(n, show=0, o=0, N=10, M=6, X=[[6, 6], [6, 7], [37, 19], [58, 46]], p=[], u=o, U)=for(n=o+1, n, show>0&& print1(o", "); show<0&& listput(L, o); U+=1<<(o-u); U>>=-u+u+=valuation(U+1, 2); p=concat(if(#p>=M, p[^1], p), o); my(c=N-sum(i=2, #p, sum(j=1, i-1, isprime(p[i]+p[j])))); for(k=u, oo, bittest(U, k-u)|| min(c-#[0|x<-p, isprime(x+k)], #p>=M)|| setsearch(X, [n, k])|| [o=k, break])); show&&print([u]); o} \\ optional args: show=1: print a(o..n-1), show=-1: append them on global list L, in both cases print [least unused number] at the end. Parameters N, M, o, ... allow to get other variants, see the wiki page for more.

CROSSREFS

Cf. A055265, A128280 (1 prime from 2 terms), A329333 (1 prime from 3 terms), A329405, ..., A329416 (N primes from M terms >= 1), A329425, A329449, ..., A329581 (N primes from M terms >= 0).

KEYWORD

allocated

nonn

AUTHOR

M. F. Hasler, Feb 10 2020

STATUS

approved

editing

Discussion
Fri Feb 14
09:22
M. F. Hasler: (moved this draft to A329574 to use this number A329572 for M=7,N=12.)
#3 by M. F. Hasler at Sat Jan 04 00:11:57 EST 2020
NAME

a(n) = Product_{prime p} p^floor(log_p n).

allocated for M. F. Hasler

DATA

1, 2, 6, 12, 60, 60, 420, 840, 2520, 2520, 27720, 27720, 360360, 360360, 360360, 720720, 12252240, 12252240, 232792560, 232792560, 232792560, 232792560, 5354228880, 5354228880, 26771144400, 26771144400, 80313433200, 80313433200, 2329089562800, 2329089562800, 72201776446800, 144403552893600

OFFSET

1,2

COMMENTS

Related to the inequality (54) in Ramanujan's paper about highly composite numbers A002182, also used in A199337: a(A329570(m))^2 is a (not minimal) bound above which all highly composite numbers are divisible by m, according to the right part of that inequality.

Like the highly composite numbers A002182, all terms in this sequence are a product of primorials.

LINKS

S. Ramanujan, <a href="https://doi.org/10.1112/plms/s2_14.1.347">Highly composite numbers</a>, Proceedings of the London Mathematical Society ser. 2, vol. XIV, no.. 1 (1915): 347-409. (DOI: 10.1112/plms/s2_14.1.347, a variant of better quality with an additional footnote is available at http://ramanujan.sirinudi.org/Volumes/published/ram15.html)

PROG

(PARI) apply( {A329572(n)=vecprod([p^logint(n, p)|p<-primes([2, n])])}, [1..33])

CROSSREFS
KEYWORD

nonn,changed

allocated

AUTHOR

M. F. Hasler, Jan 03 2020

STATUS

editing

approved