A296302 -id:A296302 - cp's OEIS frontend

A299152 Denominators of the positive solution to 2^(n-1) = Sum_{d|n} a(d) * a(n/d).

1, 1, 1, 2, 1, 1, 1, 2, 1, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 8, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 16, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 4, 1, 1, 1, 1, 1, 1, 1, 1, 1

Offset: 1

Gus Wiseman, Feb 03 2018

			Sequence begins: 1, 1, 2, 7/2, 8, 14, 32, 121/2, 126, 248, 512, 1003, 2048, 4064, 8176, 130539/8, 32768.

Antti Karttunen, Table of n, a(n) for n = 1..65537

Cf. A000005, A000010, A000740, A018804, A023900, A029935, A034691, A046643, A059966, A228369, A257098, A296302, A298971, A299119, A299149, A299151.

nn=50;
sys=Table[2^(n-1)==Sum[a[d]*a[n/d],{d,Divisors[n]}],{n,nn}];
Denominator[Array[a,nn]/.Solve[sys,Array[a,nn]][[2]]]

up_to = 65537;
prepareA299151perA299152(up_to) = { my(vmemo = vector(up_to)); for(n=1,up_to, vmemo[n] = if(1==n,n,(2^(n-1)-sumdiv(n,d,if((d>1)&&(dA299152 = prepareA299151perA299152(up_to);
A299151perA299152(n) = v299151perA299152[n];
\\ Or without memoization as:
A299151perA299152(n) = if(1==n,n,(2^(n-1)-sumdiv(n,d,if((d>1)&&(dA299151perA299152(d)*A299151perA299152(n/d),0)))/2);
A299152(n) = denominator(A299151perA299152(n)); \\ Antti Karttunen, Jul 29 2018

More terms from Antti Karttunen, Jul 29 2018

A298971 Number of compositions of n that are proper powers of Lyndon words.

Original entry on oeis.org

0, 1, 1, 2, 1, 4, 1, 5, 3, 8, 1, 16, 1, 20, 9, 35, 1, 69, 1, 110, 21, 188, 1, 381, 7, 632, 59, 1184, 1, 2300, 1, 4115, 189, 7712, 25, 14939, 1, 27596, 633, 52517, 1, 101050, 1, 190748, 2247, 364724, 1, 703331, 19, 1342283, 7713, 2581430, 1, 4985609, 193

Offset: 1

Gus Wiseman, Jan 30 2018

nonn

a(n) is the number of compositions of n that are not Lyndon words but are of the form p * p * ... * p where * is concatenation and p is a Lyndon word.

			The a(12) = 16 compositions: 111111111111, 1111211112, 11131113, 112112112, 11221122, 114114, 12121212, 123123, 131313, 132132, 1515, 222222, 2424, 3333, 444, 66.

Cf. A000005, A000031, A000740, A000961, A001045, A008965, A019536, A034691, A051953, A052823, A059966, A060223, A178472, A185700, A296302, A296373.

Table[Sum[DivisorSum[d,MoebiusMu[d/#]*(2^#-1)&]/d,{d,Most@Divisors[n]}],{n,100}]

a(n) = sumdiv(n, d, (2^d-1)*(eulerphi(n/d)-moebius(n/d))/n); \\ Michel Marcus, Jan 31 2018

a(n) = Sum_{d|n} (2^d-1)*(phi(n/d)-mu(n/d))/n.

a(n) = A008965(n) - A059966(n).

A299119 Positive solution to 2^(n-1) = (1/n) * Sum_{d|n} a(d) * a(n/d).

Original entry on oeis.org

1, 2, 6, 14, 40, 84, 224, 484, 1134, 2480, 5632, 12036, 26624, 56896, 122640, 261078, 557056, 1176876, 2490368, 5237360, 11008704, 23057408, 48234496, 100635144, 209714400, 436154368, 905962860, 1878931264, 3892314112, 8052800160, 16642998272, 34359209436

Offset: 1

Gus Wiseman, Feb 03 2018

nonn

For prime p, a(p) = 2^(p-2)*p. - Jon E. Schoenfield, Feb 03 2018

Alois P. Heinz, Table of n, a(n) for n = 1..1000

Cf. A000005, A000010, A000740, A001787, A018804, A023900, A029935, A034691, A046643, A059966, A228369, A257098, A296302, A298971, A299149, A299151.

with(numtheory):
a:= proc(n) option remember; `if`(n=1, 1, n*2^(n-2)-
       add(a(d)*a(n/d), d=divisors(n) minus {1, n})/2)
    end:
seq(a(n), n=1..35);  # Alois P. Heinz, Mar 07 2018

nn=50;
sys=Table[2^(n-1)*n==Sum[a[d]*a[n/d],{d,Divisors[n]}],{n,nn}];
Array[a,nn]/.Solve[sys,Array[a,nn]][[2]]

A319054 Maximum product of an aperiodic integer partition of n.

Original entry on oeis.org

1, 2, 3, 4, 6, 8, 12, 18, 24, 36, 54, 72, 108, 162, 216, 324, 486, 648, 972, 1458, 1944, 2916, 4374, 5832, 8748, 13122, 17496, 26244, 39366, 52488, 78732, 118098, 157464, 236196, 354294, 472392, 708588, 1062882, 1417176, 2125764, 3188646, 4251528, 6377292

Offset: 1

Gus Wiseman, Sep 09 2018

nonn

An integer partition is aperiodic if its multiplicities are relatively prime.

			Among the aperiodic partitions of 9, those with maximum product are (432) and (3222), so a(9) = 24. If periodic partitions were allowed, we would have (333) with product 27.

Cf. A000740, A000792, A000793, A000837, A007916, A100953, A281116, A289508, A289509, A296302.

Table[Max[Times@@@Select[IntegerPartitions[n],GCD@@Length/@Split[#]==1&]],{n,30}]

cp's OEIS Frontend

A299152 Denominators of the positive solution to 2^(n-1) = Sum_{d|n} a(d) * a(n/d).

Views

Author

Keywords

Examples

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A298971 Number of compositions of n that are proper powers of Lyndon words.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica

PARI

Formula

A299119 Positive solution to 2^(n-1) = (1/n) * Sum_{d|n} a(d) * a(n/d).

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Maple

Mathematica

A319054 Maximum product of an aperiodic integer partition of n.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

Mathematica