cp's OEIS Frontend

This is a front-end for the Online Encyclopedia of Integer Sequences, made by Christian Perfect. The idea is to provide OEIS entries in non-ancient HTML, and then to think about how they're presented visually. The source code is on GitHub.

Showing 1-3 of 3 results.

A306097 Terms of A121707 not in A267999.

Original entry on oeis.org

697, 1241, 1247, 1271, 1513, 2057, 2201, 2329, 2501, 2873, 3053, 3131, 3683, 3689, 3961, 4015, 4061, 4141, 4777, 4859, 4991, 5321, 5921, 5963, 6137, 6851, 6953, 7421, 7769, 7781, 7957, 8471, 8711, 8857, 9017, 9211, 9271, 9401, 9641, 9673, 10217, 10277, 10489, 10795, 11033, 11501
Offset: 1

Views

Author

M. F. Hasler, following remarks from Tomasz Ordowski, Oct 03 2018

Keywords

Comments

Numbers n such that gcd(n, 2^n-2) > 1 and gcd(n, b^n-b) = 1 for some b > 2, b < n.
Or: Numbers n such that gcd(n, 2^n-2) > 1 and for every prime factor p of n, p-1 does not divide n-1.
2057 is the first term not in A008367, nor in A287391. - M. F. Hasler, Oct 04 2018

Examples

			The smallest element of this sequence is a(1) = 697 = 17*41.
		

Crossrefs

Programs

  • PARI
    is(n,p)={for(i=1, #p=factor(n)[,1], (n-1)%(p[i]-1)||return); gcd(n, lift(Mod(2,n)^n-2))>1}

Formula

A008367 Composite but smallest prime factor >= 17.

Original entry on oeis.org

289, 323, 361, 391, 437, 493, 527, 529, 551, 589, 629, 667, 697, 703, 713, 731, 779, 799, 817, 841, 851, 893, 899, 901, 943, 961, 989, 1003, 1007, 1037, 1073, 1081, 1121, 1139, 1147, 1159, 1189, 1207, 1219, 1241, 1247, 1271, 1273, 1333, 1343, 1349, 1357, 1363, 1369, 1387, 1403, 1411, 1457
Offset: 1

Views

Author

Keywords

Comments

Composite numbers k such that k^720 mod 30030 = 1. - Gary Detlefs, May 02 2012
The asymptotic density of this sequence is 192/1001. - Amiram Eldar, Mar 22 2021

Crossrefs

Intersection of A002808 and A008366.
Cf. A287391.

Programs

  • GAP
    Filtered([17..1500],n->PowerMod(n,720,30030)=1 and not IsPrime(n)); # Muniru A Asiru, Nov 24 2018
  • Maple
    for i from 1 to 2000 do if gcd(i,30030) = 1 and not isprime(i) then print(i); fi; od;
  • Mathematica
    Select[ Range[ 1500 ], (GCD[ #1, 30030 ]==1&&!PrimeQ[ #1 ])& ]
    Select[Range[2000], ! PrimeQ[#] && FactorInteger[#][[1, 1]] >= 17 &] (* T. D. Noe, Mar 16 2013 *)
  • PARI
    is(n)={gcd(n,30030)==1 && !ispseudoprime(n)} \\ M. F. Hasler, Oct 04 2018
    

Formula

For 1 <= n < 107, a(n) = A287391(n+2); then a(107) = 2329, a(108) = 2363 are not in A287391, but again a(n) = A287391(n) for 108 < n < 120. - M. F. Hasler, Oct 04 2018

A287917 Number of distinct primorials A002110(k) > A285784(n) such that the primorials are coprime to A285784(n).

Original entry on oeis.org

1, 1, 2, 1, 1, 1, 1, 2, 1, 2, 3, 1, 2, 1, 3, 1, 2, 2, 4, 1, 3, 1, 3, 1, 2, 4, 1, 2, 3, 4, 2, 1, 3, 1, 2, 3, 5, 4, 1, 3, 5, 2, 1, 4, 1, 6, 4, 2, 3, 1, 2, 5, 1, 4, 3, 2, 6, 1, 3, 5, 2, 4, 2, 5, 1, 6, 3, 1, 6, 1, 2, 3, 4, 5, 7, 3, 1, 4, 2, 1, 6, 1, 3, 2, 7, 5, 4
Offset: 2

Views

Author

Jamie Morken and Michael De Vlieger, Jun 08 2017

Keywords

Comments

Terms of A285784 that have a(n) = 1 appear in A287390.
Terms of A285784 that have a(n) > 1 appear in A287391.
From Michael De Vlieger, Jun 09 2017: (Start)
Let primorial p_n# = A002110(n) and let m be a nonzero positive number called a totative such that gcd(t, p_n#) = 1. This sequence concerns nonprime m. A285784 is the sequence that lists unique nonprime totatives m of primorials p_n#.
For A285784(1), a(n) = infinity, since 1 is the empty product and a totative of (i.e., coprime to) all numbers. Hence the offset of a(n) is 2 and for this reason hereinafter we only consider composite totatives m.
Consider the composite totative m in A285784. For a given composite term in A285784, there is a least primorial p_a# to which m is coprime. Such m < p_a# are products of prime totatives q > p_a, the gpf of p_a#. Therefore m "appears" when there are prime totatives q < sqrt(p_a#). The smallest a for which we have this condition is a = 4, as q = 11 is less than sqrt(210). For the same reason the first composite term of A285784 is 11^2 = 121.
For n >= 2, m is coprime to a finite range of primorials p_a# .. p_b#. If m is coprime to p_b#, then it must be coprime to all primorials p_a# .. p_b# by the definition of primorial. m is no longer coprime to p_(b+1)# since at least one of its prime divisors p_(b+1) also divides p_(b+1)#. This sequence gives the range b - a + 1.
To generate data that includes all the terms of A285784 less than a limit x, we can write a while statement that operates so long as there is at least 1 totative m < x of p_n#. Since primorial p_n# is the product of the smallest n primes, fewer numbers less than x are coprime to p_n# as n increases, until exhaustion. Thus we can produce a list of unique m < x (i.e., terms of A285784 less than x) for relatively large primorials p_n#. Then we can count the instances of terms of A285784 for a list of lists of totatives m < x for primorials p_1# .. p_n# and obtain certainty about the number of instances of terms of A285784.
First position of values of a(n): {2, 4, 12, 20, 38, 47, 76, 96, 111, 139, 228, 241, 300, 339, 363, 434, 482, 566, 689, 752, 790, 862, 902, 973, 1264, 1361, 1506, 1562, 1816, ...}
Terms of A285784 that set records in a(n): {121, 169, 361, 529, 841, 961, 1369, 1681, 1849, 2209, 3481, 3721, 4489, 5041, 5329, 6241, 6889, 7921, 9409, 10201, 10609, 11449, 11881, 12769, ...}
(End)

Examples

			The sequence starts:
   n  A285784(n)  a(n)
   2      121     1
   3      143     1
   4      169     2
   5      187     1
   6      209     1
   7      221     1
   8      247     1
   9      289     2
  10      299     1
  11      323     2
  12      361     3
  13      377     1
  14      391     2
  15      403     1
  16      437     3
  17      481     1
  18      493     2   ...
		

Crossrefs

Programs

  • Mathematica
    Block[{nn = 1600, k = 1, P = 2, a}, a = Most@ Reap[While[Nand[k > 3, Length@ Sow@ Rest@ Select[Range[Min[P, nn]], And[! PrimeQ@ #, CoprimeQ[#, P]] &] == 0], k++; P *= Prime@ k]][[-1, 1]]; Function[b, Map[Count[b, #] &, Union@ b]]@ Flatten@ a] (* Michael De Vlieger, Jun 09 2017 *)

Extensions

Edited by Michael De Vlieger, Jun 09 2017
Showing 1-3 of 3 results.