login
a(n) = Sum_{1<=k<=n, k is not a divisor of n and k is not coprime to n} k.
5

%I #41 Jan 06 2024 23:35:38

%S 0,0,0,0,0,4,0,6,6,18,0,27,0,40,37,42,0,79,0,89,74,108,0,145,45,154,

%T 96,183,0,274,0,210,184,270,163,360,0,340,257,411,0,556,0,467,418,504,

%U 0,669,140,683,439,657,0,880,369,805,548,810,0,1183,0,928,779,930,502

%N a(n) = Sum_{1<=k<=n, k is not a divisor of n and k is not coprime to n} k.

%C This function can be used to prove no p^k is perfect or multi-perfect.

%H Michael De Vlieger, <a href="/A066760/b066760.txt">Table of n, a(n) for n = 1..10000</a>

%F a(n) = 1 + n*(n+1)/2 - sigma(n) - n*phi(n)/2.

%F a(n) = 0 if and only if n = 1, 4 or a prime. - _Robert G. Wilson v_, Jul 31 2004

%F a(n) = 1 + A067392(n) - A000203(n). - _Amiram Eldar_, Dec 03 2023

%F a(n) = Sum_{k=1..n} k * (1 - floor(1/gcd(n,k))) * (ceiling(n/k) - floor(n/k)). - _Wesley Ivan Hurt_, Jan 06 2024

%e There are three integers that satisfy this definition for n = 12, namely 8, 9 and 10. These sum to 27, hence a(12) = 27.

%p f:= n -> 1 + n*(n+1)/2 - numtheory:-sigma(n) - n*numtheory:-phi(n)/2;

%p 0, seq(f(n),n=2..100); # _Robert Israel_, Nov 02 2014

%t Table[n(n + 1)/2 + 1 - EulerPhi[n] * n/2 - DivisorSigma[1, n], {n, 2, 65}] (* _Robert G. Wilson v_, Jul 31 2004 *)

%t Table[Sum[k * Boole[Not[Divisible[n, k]]] * Boole[GCD[n, k] > 1], {k, n - 1}], {n, 65}] (* _Alonso del Arte_, Nov 02 2014 *)

%o (PARI) a(n)=n*(n + 1 - eulerphi(n))\2 + 1 - sigma(n) \\ _Charles R Greathouse IV_, Nov 02 2014

%o (Magma) [0] cat [1 + n*(n+1)/2 - SumOfDivisors(n) - n*EulerPhi(n)/2: n in [2..70]]; // _Vincenzo Librandi_, Nov 03 2014

%Y Cf. A000010, A000203, A000217, A023896, A024816, A045763, A067392.

%K nonn,easy

%O 1,6

%A _Jon Perry_, Jan 17 2002

%E Offset corrected to 1 by _Michael De Vlieger_, Jul 05 2014