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.

A263319 a(n) = pi(n^2)*phi(n)/2, where pi(x) denotes the number of primes not exceeding x, and phi(.) is Euler's totient function given by A000010.

Original entry on oeis.org

0, 1, 4, 6, 18, 11, 45, 36, 66, 50, 150, 68, 234, 132, 192, 216, 488, 198, 648, 312, 510, 460, 1089, 420, 1140, 732, 1161, 822, 2044, 616, 2430, 1376, 1810, 1528, 2400, 1260, 3942, 2052, 2880, 2008, 5260, 1644, 5943, 2950, 3672, 3509, 7567, 2736, 7497, 3670
Offset: 1

Views

Author

Zhi-Wei Sun, Oct 14 2015

Keywords

Comments

Conjecture: (i) All the terms of this sequence are pairwise distinct.
(ii) All the numbers phi(n)*pi(n*(n-1)) (n = 1,2,3,...) are pairwise distinct.
(iii) All the numbers phi(n^2)*pi(n^2) = n*phi(n)*pi(n^2) (n = 1,2,3,...) are pairwise distinct.
We have checked this conjecture via Mathematica. For example, we have verified that a(n) (n = 1..4*10^5) are indeed pairwise distinct.
See also A263325 for a similar conjecture.

Examples

			a(1) = 0 since pi(1^2)*phi(1)/2 = 0*1/2 = 0.
a(2) = 1 since pi(2^2)*phi(2)/2 = 2*1/2 = 1.
a(3) = 4 since pi(3^2)*phi(3)/2 = 4*2/2 = 4.
		

Crossrefs

Programs

  • Magma
    [#PrimesUpTo(n^2)*EulerPhi(n)/2: n in [1..80]]; // Vincenzo Librandi, Oct 15 2015
  • Mathematica
    a[n_]:=a[n]=PrimePi[n^2]*EulerPhi[n]/2
    Do[Print[n," ",a[n]],{n,1,50}]
  • PARI
    a(n) = primepi(n^2)*eulerphi(n)/2; \\ Michel Marcus, Oct 15 2015
    

A263326 Denominator of the rational number Sum_{d|n}1/(d+1).

Original entry on oeis.org

2, 6, 4, 30, 3, 84, 8, 90, 20, 11, 12, 5460, 7, 40, 48, 1530, 9, 7980, 20, 1155, 88, 276, 24, 81900, 78, 189, 35, 1160, 15, 38192, 32, 16830, 51, 315, 72, 3838380, 19, 780, 280, 142065, 21, 132440, 44, 828, 5520, 376, 48, 9746100, 200, 14586
Offset: 1

Views

Author

Zhi-Wei Sun, Oct 14 2015

Keywords

Comments

Conjecture: For any positive integers k and s, all the numbers Sum_{d|n}1/(d+k)^s (n = 1,2,3,...) have pairwise distinct fractional parts, and none of them is an integer.
This implies that a(n) > 1 for all n > 0.
See also A001157 for a similar conjecture involving Sum_{d|n}1/d^s.
I have verified that Sum_{d|n}1/(d+1) (n = 1..2*10^5) indeed have pairwise distinct fractional parts and none of them is an integer. For each k = 2,3,4,5,6 I have verified that Sum_{d|n}1/(d+k) (n = 1..10^5) have pairwise distinct fractional parts and none of them is integral. - Zhi-Wei Sun, Oct 20 2015.

Examples

			a(1) = 2 since sum_{d|1}1/(d+1) = 1/2.
a(2) = 6 since sum_{d|2}1/(d+1) = 1/2 + 1/3 = 5/6.
a(3) = 4 since sum_{d|3}1/(d+1) = 1/2 + 1/4 = 3/4.
		

Crossrefs

Programs

  • Maple
    f:= n -> denom(add(1/(d+1),d=numtheory:-divisors(n))):
    map(f, [$1..100]); # Robert Israel, Oct 20 2015
  • Mathematica
    Dv[n_]:=Dv[n]=Divisors[n]
    a[n_]:=a[n]=Denominator[Sum[1/(Part[Dv[n],i]+1),{i,1,Length[Dv[n]]}]]
    Do[Print[n," ",a[n]],{n,1,50}]
  • PARI
    a(n) = denominator(sumdiv(n, d, 1/(d+1))); \\ Michel Marcus, Oct 15 2015
Showing 1-2 of 2 results.