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.

A308643 Odd squarefree composite numbers k, divisible by the sum of their prime factors, sopfr (A001414).

Original entry on oeis.org

105, 231, 627, 805, 897, 1581, 2967, 3055, 4543, 5487, 6461, 6745, 7881, 9717, 10707, 14231, 15015, 16377, 21091, 26331, 29607, 33495, 33901, 33915, 35905, 37411, 38843, 40587, 42211, 45885, 49335, 50505, 51051, 53295, 55581, 60297
Offset: 1

Views

Author

David James Sycamore, Jun 13 2019

Keywords

Comments

Every term has an odd number of prime divisors (A001221(k) is odd), since if not, sopfr(k) would be even, and so not divide k, which is odd.
Some Carmichael numbers appear in this sequence, the first of which is 3240392401.
From Robert Israel, Jul 05 2019: (Start)
Includes p*q*r if p and q are distinct odd primes and r=(p-1)*(q-1)-1 is prime. Dickson's conjecture implies that there are infinitely many such terms for each odd prime p. Thus for p=3, q is in A063908 (except 3), for p=5, q is in A156300 (except 2), and for p=7, q is in A153135 (except 2). (End)

Examples

			105=3*5*7; sum of prime factors = 15 and 105 = 7*15, so 105 is a term.
		

Crossrefs

Programs

  • Magma
    [k:k in [2*d+1: d in [1..35000]]|IsSquarefree(k) and not IsPrime(k) and k mod &+PrimeDivisors(k) eq 0]; // Marius A. Burtea, Jun 19 2019
  • Maple
    with(NumberTheory);
    N := 500;
    for n from 2 to N do
    S := PrimeFactors(n);
    X := add(S);
    if IsSquareFree(n) and not mod(n, 2) = 0 and not isprime(n) and mod(n, X) = 0 then print(n);
    end if:
    end do:
  • Mathematica
    aQ[n_] := Module[{f = FactorInteger[n]}, p=f[[;;,1]]; e=f[[;;,2]]; Length[e] > 1 && Max[e]==1 && Divisible[n, Plus@@(p^e)]]; Select[Range[1, 61000, 2], aQ] (* Amiram Eldar, Jul 04 2019 *)

A289556 Primes p such that both 5*p - 4 and 4*p - 5 are prime.

Original entry on oeis.org

3, 7, 13, 43, 67, 109, 127, 151, 163, 211, 277, 307, 373, 457, 463, 601, 613, 673, 727, 853, 919, 967, 1021, 1117, 1171, 1231, 1399, 1471, 1483, 1747, 1789, 1933, 2029, 2251, 2311, 2389, 2503, 2521, 2557, 2659, 2851, 2857, 3019, 3067, 3121, 3229, 3583, 3613, 3637, 3691, 3697
Offset: 1

Views

Author

David James Sycamore, Aug 02 2017

Keywords

Comments

The terms of this sequence belong to two disjoint subsequences, namely those for which |A(5*p) - A(4*p)| = 9; (3,7,13,43,67,127,163,211,277,307,457,...), and those for which 5*A(4*p) - 3*A(5*p) = 3, (109,151,373,673,919,...), where A = A288814.
Note: A288814(n) = A056240(n) for all composite n.

Examples

			P=7: 5*7 - 4 = 31, 4*7 - 5 = 23, both prime so 7 is in this sequence, and belongs to the subsequence of terms satisfying A(4*p) - A(3*p) = 9.
P=109: 5*109 - 4 = 541, 4*109 - 5 = 431, both prime so 109 is in this sequence, and belongs to the subsequence of terms satisfying 5*A(4*p) - 3*A(5*p) = 3.
		

Crossrefs

Intersection of A136051 and A156300. - Michel Marcus, Aug 04 2017

Programs

  • Mathematica
    Select[Prime@ Range@ 516, Times @@ Boole@ Map[PrimeQ, {5 # - 4, 4 # - 5}] > 0 &] (* Michael De Vlieger, Aug 02 2017 *)

Extensions

More terms from Altug Alkan, Aug 02 2017
Showing 1-2 of 2 results.