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.

A210452 Number of integers k

Original entry on oeis.org

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

Views

Author

Zhi-Wei Sun, Jan 20 2013

Keywords

Comments

Conjecture: a(n)>0 for all n>4.
This implies the twin prime conjecture since k*p is not practical for any prime p>sigma(k)+1.
Zhi-Wei Sun also made the following conjectures:
(1) For each integer n>197, there is a practical number k
(2) For every n=9,10,... there is a practical number k
(3) For any integer n>26863, the interval [1,n] contains five consecutive integers m-2, m-1, m, m+1, m+2 with m-1 and m+1 both prime, and m-2, m, m+2, m*n all practical.

Examples

			a(11)=1 since 5 and 7 are twin primes, and 6 and 6*11 are both practical.
		

Programs

  • Mathematica
    f[n_]:=f[n]=FactorInteger[n]
    Pow[n_, i_]:=Pow[n, i]=Part[Part[f[n], i], 1]^(Part[Part[f[n], i], 2])
    Con[n_]:=Con[n]=Sum[If[Part[Part[f[n], s+1], 1]<=DivisorSigma[1, Product[Pow[n, i], {i, 1, s}]]+1, 0, 1], {s, 1, Length[f[n]]-1}]
    pr[n_]:=pr[n]=n>0&&(n<3||Mod[n, 2]+Con[n]==0)
    a[n_]:=a[n]=Sum[If[PrimeQ[k-1]==True&&PrimeQ[k+1]==True&&pr[k]==True&&pr[k*n]==True,1,0],{k,1,n-1}]
    Do[Print[n," ",a[n]],{n,1,100}]