Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
login
The OEIS is supported by the many generous donors to the OEIS Foundation.

 

Logo
Hints
(Greetings from The On-Line Encyclopedia of Integer Sequences!)
Search: a097218 -id:a097218
Displaying 1-3 of 3 results found. page 1
     Sort: relevance | references | number | modified | created      Format: long | short | data
A076078 a(n) is the number of nonempty sets of distinct positive integers that have a least common multiple of n. +10
49
1, 2, 2, 4, 2, 10, 2, 8, 4, 10, 2, 44, 2, 10, 10, 16, 2, 44, 2, 44, 10, 10, 2, 184, 4, 10, 8, 44, 2, 218, 2, 32, 10, 10, 10, 400, 2, 10, 10, 184, 2, 218, 2, 44, 44, 10, 2, 752, 4, 44, 10, 44, 2, 184, 10, 184, 10, 10, 2, 3748, 2, 10, 44, 64, 10, 218, 2, 44, 10, 218, 2, 3392, 2, 10 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,2
COMMENTS
a(n)=1 iff n=1, a(p^k)=2^k, a(p*q)=10; where p & q are unique primes. a(n) cannot equal an odd number >1. - Robert G. Wilson v
If m has more divisors than n, then a(m) > a(n). - Matthew Vandermast, Aug 22 2004
If n is of the form p^r*q^s where p & q are distinct primes and r & s are nonnegative integers then a(n)=2^(rs)*(2^(r+s+1) -2^r-2^s+1); for example f(1400846643)=f(3^5*7^8)=2^(5*8)*(2^ (5+8+1)-2^5-2^8+1)=17698838672310272. Also if n=p_1^r_1*p_2^r_2*...*p_k^r_k where p_1,p_2,...,p_k are distinct primes and r_1,r_2,...,r_k are natural numbers then 2^(r_1*r_2*...*r_k)||a(n). - Farideh Firoozbakht, Aug 06 2005
None of terms is divisible by Mersenne numbers 3 or 7. For any n, a(n) is congruent to A008836(n) mod 3. Since A008836(n) is always 1 or -1, this implies that A000225(2)=3 never divides a(n). - Matthew Vandermast, Oct 12 2010
There are terms divisible by larger Mersenne numbers. For example, a(2*3*5*7*11*13*19*23^3) is divisible by 31. - Max Alekseyev, Nov 18 2010
LINKS
FORMULA
2^d(n) - 1 = Sum_{m|n} a(m), where d(n) = A000005(n) is the number of divisors of n, so a(n) = Sum_{m|n} mu(n/m)*(2^d(m) - 1).
a(n) = 2*A069626(n), for n > 1. - Ridouane Oudra, Mar 12 2024
EXAMPLE
a(6) = 10. The sets with LCM 6 are {6}, {1,6}, {2,3}, {2,6}, {3,6}, {1,2,3}, {1,2,6}, {1,3,6}, {2,3,6}, {1,2,3,6}.
MAPLE
with(numtheory): seq(add(mobius(n/d)*(2^tau(d)-1), d in divisors(n)), n=1..80); # Ridouane Oudra, Mar 12 2024
MATHEMATICA
f[n_] := Block[{d = Divisors[n]}, Plus @@ (MoebiusMu[n/d](2^DivisorSigma[0, d] - 1))]; Table[ f[n], {n, 75}] (* Robert G. Wilson v *)
PROG
(PARI) a(n) = local(f, l, s, t, q); f = factor(n); l = matsize(f)[1]; s = 0; forvec(v = vector(l, i, [0, 1]), q = sum(i = 1, l, v[i]); t = (-1)^(l - q)*2^prod(i = 1, l, f[i, 2] + v[i]); s += t); s; \\ Definition corrected by David Wasserman, Dec 26 2007
CROSSREFS
Cf. A069626.
KEYWORD
easy,nonn,nice
AUTHOR
Amarnath Murthy, Oct 05 2002
EXTENSIONS
Edited by Dean Hickerson, Oct 08 2002
Definition corrected by David Wasserman, Dec 26 2007
Edited by Charles R Greathouse IV, Aug 02 2010
Edited by Max Alekseyev, Nov 18 2010
STATUS
approved
A214547 Deficient numbers for which the (absolute value of) abundance is not a divisor. +10
1
3, 5, 7, 9, 11, 13, 14, 15, 17, 19, 21, 22, 23, 25, 26, 27, 29, 31, 33, 34, 35, 37, 38, 39, 41, 43, 45, 46, 47, 49, 50, 51, 52, 53, 55, 57, 58, 59, 61, 62, 63, 65, 67, 68, 69, 71, 73, 74, 75, 76, 77, 79, 81, 82, 83, 85, 86, 87, 89, 91, 92, 93, 94, 95, 97, 98, 99, 101, 103, 105, 106 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
This is to A214408 as deficient numbers are to abundant numbers.
Differs from A097218, which does not contain 105, for example.
The deficient numbers which are *not* in the sequence are 2, 4, 8, 10, 16, 32, 44, 64, 128, 136, 152, 184, 256, 512, 752, 884, 1024, 2048, 2144, 2272, 2528, 4096, 8192, 8384, 12224, 16384, 17176, 18632, 18904, 32768, 32896, 33664, ... the union of powers of 2 and the terms of A060326. - M. F. Hasler, Jul 21 2012
LINKS
FORMULA
Terms A005100(n) such that |A033880(A005100(n))| does not divide A005100(n).
EXAMPLE
7 is in the sequence because 7 is deficient, and its abundance is -6, and |-6| = 6 does not divide 7.
MAPLE
filter:= proc(n) local t;
t:= 2*n-numtheory:-sigma(n);
t > 0 and n mod t <> 0
end proc:
select(filter, [$1..200]); # Robert Israel, Nov 13 2019
MATHEMATICA
q[n_] := Module[{def = 2*n - DivisorSigma[1, n]}, def > 0 && !Divisible[n, def]]; Select[Range[120], q] (* Amiram Eldar, Apr 07 2024 *)
PROG
(PARI) is_A214547(n)={sigma(n)<2*n & n%(2*n-sigma(n))} \\ M. F. Hasler, Jul 21 2012
CROSSREFS
KEYWORD
nonn,easy
AUTHOR
Jonathan Vos Post, Jul 20 2012
EXTENSIONS
Given terms double-checked with the PARI script by M. F. Hasler, Jul 21 2012
STATUS
approved
A097217 Odd numbers n such that A076078(n) > n, where A076078(n) equals the number of sets of distinct positive integers with a least common multiple of n. +10
0
105, 135, 165, 195, 225, 315, 405, 495, 525, 567, 585, 675, 693, 765, 819, 825, 855, 945, 975, 1035, 1071, 1125, 1155, 1197, 1215, 1275, 1287, 1305, 1323, 1365, 1395, 1425, 1449, 1485, 1575, 1665, 1683, 1701, 1725, 1755, 1785, 1827, 1845, 1881, 1925 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
Odd members of A097216.
LINKS
CROSSREFS
KEYWORD
nonn
AUTHOR
Matthew Vandermast, Aug 13 2004
STATUS
approved
page 1

Search completed in 0.004 seconds

Lookup | Welcome | Wiki | Register | Music | Plot 2 | Demos | Index | Browse | More | WebCam
Contribute new seq. or comment | Format | Style Sheet | Transforms | Superseeker | Recents
The OEIS Community | Maintained by The OEIS Foundation Inc.

License Agreements, Terms of Use, Privacy Policy. .

Last modified August 18 23:05 EDT 2024. Contains 375284 sequences. (Running on oeis4.)