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-2 of 2 results.

A080365 Composite numbers k whose smallest and largest prime factors are unitary prime factors.

Original entry on oeis.org

6, 10, 14, 15, 21, 22, 26, 30, 33, 34, 35, 38, 39, 42, 46, 51, 55, 57, 58, 62, 65, 66, 69, 70, 74, 77, 78, 82, 85, 86, 87, 90, 91, 93, 94, 95, 102, 105, 106, 110, 111, 114, 115, 118, 119, 122, 123, 126, 129, 130, 133, 134, 138, 141, 142, 143, 145, 146, 154, 155, 158
Offset: 1

Views

Author

Labos Elemer, Feb 21 2003

Keywords

Examples

			k=90 is not a prime; 90 = 2*3*3*5; extremal prime factors are 2 and 5; gcd(2, 90/2) = gcd(5, 90/5) = 1, so 2 and 5 are unitary prime divisors of 90, thus 90 is in the sequence.
		

Crossrefs

Programs

  • GAP
    D:=List(Filtered([2..160],i->not IsPrime(i)),Factors);;
    a:=[];;  for i in D do if Gcd(i[1],Product(i)/i[1])=1 and Gcd(i[Length(i)],Product(i)/i[Length(i)])=1 then Add(a,Product(i)); fi; od; a;  # Muniru A Asiru, Jul 10 2018~
  • Mathematica
    ma[x_] := Part[Reverse[Flatten[FactorInteger[x]]], 2] mi[x_] := Part[Flatten[FactorInteger[x]], 1] k=0; Do[s=mi[n]; s1=ma[n]; If[Equal[GCD[s, n/s], 1]&&Equal[GCD[s1, n/s1], 1]&&!PrimeQ[n], Print[n]], {n, 2, 256}]
  • PARI
    lista(nn) = {forcomposite(n=1, nn, my(f=factor(n)[,1], p = vecmin(f), q = vecmax(f)); if ((gcd(p, n/p) == 1) && (gcd(q, n/q) == 1), print1(n, ", ")););} \\ Michel Marcus, Jul 09 2018
    

A080366 Numbers k whose least and greatest prime divisors are non-unitary.

Original entry on oeis.org

4, 8, 9, 16, 25, 27, 32, 36, 49, 64, 72, 81, 100, 108, 121, 125, 128, 144, 169, 196, 200, 216, 225, 243, 256, 288, 289, 300, 324, 343, 361, 392, 400, 432, 441, 484, 500, 512, 529, 576, 588, 600, 625, 648, 675, 676, 729, 784, 800, 841, 864, 900, 961, 968, 972
Offset: 1

Views

Author

Labos Elemer, Feb 21 2003

Keywords

Examples

			n=300: it is not a prime, 300 = 2*2*3*5*5; extremal prime factors are 2 and 5; gcd(2, 300/2) > 1 and gcd(5, 300/5) > 1, so neither 2 nor 5 is a unitary prime divisor of 300, thus 300 is in this sequence. - _Labos Elemer_, corrected by _Jeppe Stig Nielsen_, Jun 27 2017
		

Crossrefs

Programs

  • Mathematica
    ma[x_] := Part[Reverse[Flatten[FactorInteger[x]]], 2] mi[x_] := Part[Flatten[FactorInteger[x]], 1] k=0; Do[s=mi[n]; s1=ma[n]; If[ !Equal[GCD[s, n/s], 1]&&!Equal[GCD[s1, n/s1], 1]&&!PrimeQ[n], Print[n]], {n, 2, 1000}]
  • PARI
    isA080366(n) = e=factor(n)[,2];e&&e[1]>1&&e[#e]>1 \\ Jeppe Stig Nielsen, Jun 27 2017
Showing 1-2 of 2 results.