A286163 -id:A286163 - cp's OEIS frontend

A286152 Compound filter: a(n) = T(A051953(n), A046523(n)), where T(n,k) is sequence A000027 used as a pairing function.

0, 2, 2, 12, 2, 40, 2, 59, 18, 61, 2, 179, 2, 86, 73, 261, 2, 265, 2, 265, 100, 148, 2, 757, 33, 185, 129, 367, 2, 1297, 2, 1097, 166, 271, 131, 1735, 2, 320, 205, 1105, 2, 1741, 2, 619, 517, 430, 2, 3113, 52, 850, 295, 769, 2, 1747, 205, 1517, 346, 625, 2, 5297, 2, 698, 730, 4497, 248, 2821, 2, 1117, 460, 2821, 2, 7069, 2, 941, 1070, 1315, 248, 3457, 2, 4513

Offset: 1

Antti Karttunen, May 04 2017

nonn

Antti Karttunen, Table of n, a(n) for n = 1..10000
MathWorld, Pairing Function

Cf. A000027, A046523, A051953, A285729, A286142, A286143, A286144, A286149, A286154, A286160, A286161, A286162, A286163, A286164.

Table[(2 + (#1 + #2)^2 - #1 - 3 #2)/2 & @@ {n - EulerPhi@ n, Times @@ MapIndexed[Prime[First@ #2]^#1 &, Sort[FactorInteger[n][[All, -1]], Greater]] - Boole[n == 1]}, {n, 80}] (* Michael De Vlieger, May 04 2017 *)

A051953(n) = (n - eulerphi(n));
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); };  \\ This function from Charles R Greathouse IV, Aug 17 2011
A286152(n) = (2 + ((A051953(n)+A046523(n))^2) - A051953(n) - 3*A046523(n))/2;
for(n=1, 10000, write("b286152.txt", n, " ", A286152(n)));

from sympy import factorint, totient
def T(n, m): return ((n + m)**2 - n - 3*m + 2)/2
def P(n):
    f = factorint(n)
    return sorted([f[i] for i in f])
def a046523(n):
    x=1
    while True:
        if P(n) == P(x): return x
        else: x+=1
def a(n): return T(n - totient(n), a046523(n)) # Indranil Ghosh, May 05 2017

(define (A286152 n) (* (/ 1 2) (+ (expt (+ (A051953 n) (A046523 n)) 2) (- (A051953 n)) (- (* 3 (A046523 n))) 2)))

a(n) = (1/2)*(2 + ((A051953(n)+A046523(n))^2) - A051953(n) - 3*A046523(n)).

A286154 Compound filter: a(n) = T(A055396(n), A000010(n)), where T(n,k) is sequence A000027 used as a pairing function.

Original entry on oeis.org

0, 1, 5, 2, 18, 2, 40, 7, 23, 7, 96, 7, 142, 16, 38, 29, 238, 16, 308, 29, 80, 46, 444, 29, 234, 67, 173, 67, 676, 29, 791, 121, 212, 121, 328, 67, 1093, 154, 302, 121, 1339, 67, 1499, 191, 302, 232, 1785, 121, 994, 191, 530, 277, 2227, 154, 864, 277, 668, 379, 2718, 121, 2944, 436, 668, 497, 1228, 191, 3505, 497, 992, 277, 3936, 277, 4207, 631, 822, 631

Offset: 1

Antti Karttunen, May 04 2017

nonn

Antti Karttunen, Table of n, a(n) for n = 1..10000
MathWorld, Pairing Function

Cf. A000010, A000027, A055396, A286142, A286143, A286144, A286152, A286160, A286161, A286162, A286163, A286164.

Table[(2 + (#1 + #2)^2 - #1 - 3 #2)/2 & @@ {If[n == 1, 0, PrimePi[ FactorInteger[n][[1, 1]] ]], EulerPhi@ n}, {n, 76}] (* Michael De Vlieger, May 04 2017 *)

A000010(n) = eulerphi(n);
A055396(n) = if(n==1, 0, primepi(factor(n)[1, 1])); \\ This function from Charles R Greathouse IV, Apr 23 2015
A286154(n) = (2 + ((A055396(n)+A000010(n))^2) - A055396(n) - 3*A000010(n))/2;
for(n=1, 10000, write("b286154.txt", n, " ", A286154(n)));

from sympy import primepi, isprime, primefactors, totient
def a049084(n): return primepi(n)*(1*isprime(n))
def a055396(n): return 0 if n==1 else a049084(min(primefactors(n)))
def T(n, m): return ((n + m)**2 - n - 3*m + 2)/2
def a(n): return T(a055396(n), totient(n)) # Indranil Ghosh, May 05 2017

(define (A286154 n) (* (/ 1 2) (+ (expt (+ (A055396 n) (A000010 n)) 2) (- (A055396 n)) (- (* 3 (A000010 n))) 2)))

a(n) = (1/2)*(2 + ((A055396(n)+A000010(n))^2) - A055396(n) - 3*A000010(n)).

A319705 Filter sequence which for primes p records a distinct value for each distinct multiset formed from the lengths of 1-runs in its binary representation [A286622(p)], and for all other numbers assigns a unique number.

Original entry on oeis.org

1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 11, 13, 14, 15, 5, 16, 11, 17, 18, 19, 20, 21, 22, 23, 24, 25, 20, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 33, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 38, 48, 49, 50, 51, 52, 53, 54, 42, 55, 56, 57, 58, 59, 11, 60, 61, 62, 20, 63, 33, 64, 65, 66, 67, 68, 42, 69, 70, 71, 38, 72, 73, 74, 75, 76, 38, 77, 78, 79, 80, 81, 82, 83, 11, 84, 85, 86, 38

Offset: 1

Antti Karttunen, Sep 26 2018

nonn

Restricted growth sequence transform of function f defined as f(n) = A278222(n) when n is a prime, otherwise -n.
After its initial term 3, Fermat primes (A019434) gives the positions of 5 in this sequence, while the Mersenne primes (A000668) are each assigned to their own singleton equivalence class.
For all i, j:
  a(i) = a(j) => A305900(i) = A305900(j),
  a(i) = a(j) => A286622(i) = A286622(j),
  a(i) = a(j) => A305795(i) = A305795(j).

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

Cf. A000668, A019434, A278222, A286163, A286622, A305795, A305900.

Cf. also A319704, A319706.

up_to = 100000;
rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om,invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om,invec[i],i); outvec[i] = u; u++ )); outvec; };
A005940(n) = { my(p=2, t=1); n--; until(!n\=2, if((n%2), (t*=p), p=nextprime(p+1))); t };
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); };  \\ From A046523
A278222(n) = A046523(A005940(1+n));
A319705aux(n) = if(isprime(n),A278222(n),-n);
v319705 = rgs_transform(vector(up_to,n,A319705aux(n)));
A319705(n) = v319705[n];

A286149 Compound filter: a(n) = T(A046523(n), A109395(n)), where T(n,k) is sequence A000027 used as a pairing function.

Original entry on oeis.org

1, 5, 8, 14, 17, 34, 30, 44, 19, 51, 68, 103, 93, 72, 196, 152, 155, 103, 192, 132, 72, 126, 278, 349, 32, 159, 53, 165, 437, 976, 498, 560, 709, 237, 786, 739, 705, 282, 159, 402, 863, 660, 948, 243, 337, 384, 1130, 1273, 49, 132, 1546, 288, 1433, 349, 126, 459, 282, 567, 1772, 2761, 1893, 636, 165, 2144, 2421, 1921, 2280, 390, 2707, 2046, 2558, 2773, 2703

Offset: 1

Antti Karttunen, May 04 2017

nonn

Antti Karttunen, Table of n, a(n) for n = 1..10000
MathWorld, Pairing Function

Cf. A000027, A046523, A109395, A285729, A286142, A286143, A286144, A286152, A286154, A286160, A286161, A286162, A286163, A286164.

Table[(2 + (#1 + #2)^2 - #1 - 3 #2)/2 & @@ {Times @@ MapIndexed[ Prime[First@ #2]^#1 &, Sort[FactorInteger[n][[All, -1]], Greater]] - Boole[n == 1], Denominator[EulerPhi[n]/n]}, {n, 73}] (* Michael De Vlieger, May 04 2017 *)

A109395(n) = n/gcd(n, eulerphi(n));
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); };  \\ This function from Charles R Greathouse IV, Aug 17 2011
A286149(n) = (1/2)*(2 + ((A046523(n)+A109395(n))^2) - A046523(n) - 3*A109395(n));
for(n=1, 10000, write("b286149.txt", n, " ", A286149(n)));

from sympy import factorint, totient, gcd
def T(n, m): return ((n + m)**2 - n - 3*m + 2)/2
def P(n):
    f = factorint(n)
    return sorted([f[i] for i in f])
def a046523(n):
    x=1
    while True:
        if P(n) == P(x): return x
        else: x+=1
def a(n): return T(a046523(n), n/gcd(n, totient(n))) # Indranil Ghosh, May 05 2017

(define (A286149 n) (* (/ 1 2) (+ (expt (+ (A046523 n) (A109395 n)) 2) (- (A046523 n)) (- (* 3 (A109395 n))) 2)))

a(n) = (1/2)*(2 + ((A046523(n)+A109395(n))^2) - A046523(n) - 3*A109395(n)).

A291762 Restricted growth sequence transform of ((-1)^A000120(n))*A046523(n); filter combining the parity of binary weight with the prime signature of n.

Original entry on oeis.org

1, 2, 3, 4, 3, 5, 2, 6, 7, 5, 2, 8, 2, 9, 5, 10, 3, 8, 2, 8, 9, 9, 3, 11, 4, 9, 12, 13, 3, 14, 2, 15, 5, 5, 9, 16, 2, 9, 5, 11, 2, 17, 3, 13, 8, 5, 2, 18, 4, 13, 5, 13, 3, 11, 9, 19, 5, 5, 2, 20, 2, 9, 8, 21, 5, 14, 2, 8, 9, 17, 3, 22, 2, 9, 8, 13, 5, 14, 2, 18, 10, 9, 3, 23, 5, 5, 9, 19, 3, 20, 9, 8, 9, 9, 5, 24, 2, 13, 8, 25, 3, 14, 2, 19, 14, 5, 2, 22, 2, 17

Offset: 1

Antti Karttunen, Sep 11 2017

nonn

Equally, restricted growth sequence transform of sequence b defined as b(1) = 1; b(n) = A046523(n) + A010060(n) for n > 1, which starts as 1, 3, 2, 5, 2, 6, 3, 9, 4, 6, 3, 12, 3, 7, 6, 17, 2, 12, 3, 12, 7, 7, ...

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

Cf. A000120, A010060, A046523.
Cf. A101296, A286163, A291761 (related or similar filtering sequences).
Cf. A027697 (positions of 2's), A027699 (of 3's), A130593 (of 5's and 7's), A230095 (of 9's).
Cf. also A231431, A235001.

rgs_transform(invec) = { my(om = Map(), outvec = vector(length(invec)), u=1); for(i=1, length(invec), if(mapisdefined(om,invec[i]), my(pp = mapget(om, invec[i])); outvec[i] = outvec[pp] , mapput(om,invec[i],i); outvec[i] = u; u++ )); outvec; };
write_to_bfile(start_offset,vec,bfilename) = { for(n=1, length(vec), write(bfilename, (n+start_offset)-1, " ", vec[n])); }
A046523(n) = { my(f=vecsort(factor(n)[, 2], , 4), p); prod(i=1, #f, (p=nextprime(p+1))^f[i]); };  \\ This function from Charles R Greathouse IV, Aug 17 2011
write_to_bfile(1,rgs_transform(vector(65537,n,((-1)^hammingweight(n))*A046523(n))),"b291762_upto65537.txt");
\\ Or alternatively:
A010060(n) = (hammingweight(n)%2);
f(n) = if(1==n,n,A046523(n)+A010060(n));
write_to_bfile(1,rgs_transform(vector(16385,n,f(n))),"b291762.txt");

cp's OEIS Frontend

A286152 Compound filter: a(n) = T(A051953(n), A046523(n)), where T(n,k) is sequence A000027 used as a pairing function.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Python

Scheme

Formula

A286154 Compound filter: a(n) = T(A055396(n), A000010(n)), where T(n,k) is sequence A000027 used as a pairing function.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Python

Scheme

Formula

A319705 Filter sequence which for primes p records a distinct value for each distinct multiset formed from the lengths of 1-runs in its binary representation [A286622(p)], and for all other numbers assigns a unique number.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI

A286149 Compound filter: a(n) = T(A046523(n), A109395(n)), where T(n,k) is sequence A000027 used as a pairing function.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Python

Scheme

Formula

A291762 Restricted growth sequence transform of ((-1)^A000120(n))*A046523(n); filter combining the parity of binary weight with the prime signature of n.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

PARI