A179645 -id:A179645 - cp's OEIS frontend

A182944 Square array A(i,j), i >= 1, j >= 1, of prime powers prime(i)^j, by descending antidiagonals.

2, 4, 3, 8, 9, 5, 16, 27, 25, 7, 32, 81, 125, 49, 11, 64, 243, 625, 343, 121, 13, 128, 729, 3125, 2401, 1331, 169, 17, 256, 2187, 15625, 16807, 14641, 2197, 289, 19, 512, 6561, 78125, 117649, 161051, 28561, 4913, 361, 23

Offset: 1

Clark Kimberling, Dec 14 2010

We alternatively refer to this sequence as a triangle T(.,.), with T(n,k) = A(k,n-k+1) = prime(k)^(n-k+1).
The monotonic ordering of this sequence, prefixed by 1, is A000961.
The joint-rank array of this sequence is A182869.
Main diagonal gives A062457. - Omar E. Pol, Sep 11 2018

			Square array A(i,j) begins:
  i \ j: 1      2      3      4      5  ...
  ---\-------------------------------------
  1:     2,     4,     8,    16,    32, ...
  2:     3,     9,    27,    81,   243, ...
  3:     5,    25,   125,   625,  3125, ...
  4:     7,    49,   343,  2401, 16807, ...
  ...
The triangle T(n,k) begins:
  n\k:  1     2     3     4     5     6  ...
  1:    2
  2:    4     3
  3:    8     9     5
  4:   16    27    25     7
  5:   32    81   125    49    11
  6:   64   243   625   343   121    13
  ...

Michael De Vlieger, Table of n, a(n) for n = 1..11325 (rows 1..150, flattened)
Michael De Vlieger, Diagram showing row n of triangle in a semicircle as noted, with a color function associated with the magnitude of T(n,k) compared to 2^n in light blue, where prime(n) is the smallest and the prime power indicated in red the largest in the row.

Cf. A000961, A006939 (row products of triangle), A062457, A182945, A332979 (row maxima of triangle).
Columns: A000040 (1), A001248 (2), A030078 (3), A030514 (4), A050997 (5), A030516 (6), A092759 (7), A179645 (8), A179665 (9), A030629 (10).
A319075 extends the array with 0th powers.
Subtable of A242378, A284457, A329332.

TableForm[Table[Prime[n]^j,{n,1,14},{j,1,8}]]

From Peter Munn, Dec 29 2019: (Start)
A(i,j) = A182945(j,i) = A319075(j,i).
A(i,j) = A242378(i-1,2^j) = A329332(2^(i-1),j).
A(i,i) = A062457(i).
(End)

Clarified in respect of alternate reading as a triangle by Peter Munn, Aug 28 2022

A189990 Numbers with prime factorization p^2*q^6.

Original entry on oeis.org

576, 1600, 2916, 3136, 7744, 10816, 18225, 18496, 23104, 33856, 35721, 53824, 61504, 62500, 87616, 88209, 107584, 118336, 123201, 140625, 141376, 179776, 210681, 222784, 238144, 263169, 287296, 322624, 341056, 385641, 399424, 440896, 470596

Offset: 1

Vladimir Joseph Stephan Orlovsky, May 03 2011

nonn

T. D. Noe, Table of n, a(n) for n = 1..1000
Will Nicholes, List of Prime Signatures
Index to sequences related to prime signature

Subsequence of A137484.
Cf. A179645, A179664.
Cf. A065036, A085548, A085966, A085968.

f[n_]:=Sort[Last/@FactorInteger[n]]=={2,6}; Select[Range[800000],f]

list(lim)=my(v=List(),t);forprime(p=2, (lim\4)^(1/6), t=p^6;forprime(q=2, sqrt(lim\t), if(p==q, next);listput(v,t*q^2))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 24 2011

from math import isqrt
from sympy import primepi, integer_nthroot, primerange
def A189990(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x-sum(primepi(isqrt(x//p**6)) for p in primerange(integer_nthroot(x,6)[0]+1))+primepi(integer_nthroot(x,8)[0])
    return bisection(f,n,n) # Chai Wah Wu, Feb 21 2025

Sum_{n>=1} 1/a(n) = P(2)*P(6) - P(8) = A085548 * A085966 - A085968 = 0.003658..., where P is the prime zeta function. - Amiram Eldar, Jul 06 2020

a(n) = A065036(n)^2. - Chai Wah Wu, Mar 27 2025

A115975 Numbers of the form p^k, where p is a prime and k is a Fibonacci number.

Original entry on oeis.org

1, 2, 3, 4, 5, 7, 8, 9, 11, 13, 17, 19, 23, 25, 27, 29, 31, 32, 37, 41, 43, 47, 49, 53, 59, 61, 67, 71, 73, 79, 83, 89, 97, 101, 103, 107, 109, 113, 121, 125, 127, 131, 137, 139, 149, 151, 157, 163, 167, 169, 173, 179, 181, 191, 193, 197, 199, 211, 223, 227, 229, 233

Offset: 1

Giovanni Teofilatto, Mar 15 2006; corrected Apr 23 2006

Amiram Eldar, Table of n, a(n) for n = 1..10000

Subsequence of A000961 (powers of primes).
Cf. A117245 (partial sums).
Subsequences: A000040, A001248, A030078, A050997, A138031, A179645.

With[{nn=60},Take[Join[{1},Union[First[#]^Last[#]&/@Union[Flatten[ Outer[List,Prime[Range[nn]],Fibonacci[Range[nn/6]]],1]]]],70]] (* Harvey P. Dale, Jun 05 2012 *)
fib[lim_] := Module[{s = {}, f = 1, k = 2}, While[f <= lim, AppendTo[s, f]; k++; f = Fibonacci[k]]; s]; seq[max_] := Module[{s = {1}, p = 2, e = 1, f = {}}, While[e > 0, e = Floor[Log[p, max]]; If[f == {}, f = fib[e], f = Select[f, # <= e &]]; s = Join[s, p^f]; p = NextPrime[p]]; Sort[s]]; seq[250] (* Amiram Eldar, Aug 09 2024 *)

{m=240;v=Set([]);forprime(p=2,m,i=0;while((s=p^fibonacci(i))

Edited and corrected by Klaus Brockhaus, Apr 25 2006

A369938 Numbers whose maximal exponent in their prime factorization is a power of 2.

Original entry on oeis.org

2, 3, 4, 5, 6, 7, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 25, 26, 28, 29, 30, 31, 33, 34, 35, 36, 37, 38, 39, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 55, 57, 58, 59, 60, 61, 62, 63, 65, 66, 67, 68, 69, 70, 71, 73, 74, 75, 76, 77

Offset: 1

Amiram Eldar, Feb 06 2024

First differs from its subsequence A138302 \ {1} at n = 378: a(378) = 432 = 2^4 * 3^3 is not a term of A138302.
First differs from A096432, A220218 \ {1}, A335275 \ {1} and A337052 \ {1} at n = 56, and from A270428 \ {1} at n = 113.
Numbers k such that A051903(k) is a power of 2.
The asymptotic density of this sequence is 1/zeta(3) + Sum_{k>=2} (1/zeta(2^k+1) - 1/zeta(2^k)) = 0.87442038669659566330... .

Amiram Eldar, Table of n, a(n) for n = 1..10000

Cf. A000079, A002117, A051903.
Subsequences: A001146, A001248, A005117, A030514, A030635, A050376, A062503, A113849, A138302 \ {1}, A179645.
Similar sequences: A368714, A369937, A369939.
Cf. A096432, A220218, A270428, A335275, A337052.

pow2Q[n_] := n == 2^IntegerExponent[n, 2];
Select[Range[2, 100], pow2Q[Max[FactorInteger[#][[;; , 2]]]] &]
Select[Range[2,80],IntegerQ[Log2[Max[FactorInteger[#][[;;,2]]]]]&] (* Harvey P. Dale, Nov 06 2024 *)

ispow2(n) = n >> valuation(n, 2) == 1;
is(n) = n > 1 && ispow2(vecmax(factor(n)[, 2]));

A179689 Numbers with prime signature {7,2}, i.e., of form p^7*q^2 with p and q distinct primes.

Original entry on oeis.org

1152, 3200, 6272, 8748, 15488, 21632, 36992, 46208, 54675, 67712, 107163, 107648, 123008, 175232, 215168, 236672, 264627, 282752, 312500, 359552, 369603, 445568, 476288, 574592, 632043, 645248, 682112, 703125, 789507, 798848, 881792, 1013888

Offset: 1

Vladimir Joseph Stephan Orlovsky, Jul 24 2010

nonn

T. D. Noe, Table of n, a(n) for n = 1..1000
OEIS Wiki, Numbers with same prime signature.
Will Nicholes, Prime Signatures

Cf. A006881, A007304, A065036, A085986, A085987, A092759, A178739, A179642, A179643, A179644, A179645, A179646, A179664, A179665, A179666, A179667, A179668, A179669, A179670, A179671, A179672, A179688.

Cf. A085548, A085967, A085969.

a:= proc(n) option remember; local k;
      for k from 1+ `if` (n=1, 1, a(n-1))
        while sort (map (x-> x[2], ifactors(k)[2]), `>`)<>[7, 2]
      do od; k
    end:
seq (a(n), n=1..32);  # Alois P. Heinz, Jan 23 2011

f[n_]:=Sort[Last/@FactorInteger[n]]=={2,7}; Select[Range[10^6], f]

list(lim)=my(v=List(),t);forprime(p=2, (lim\4)^(1/7), t=p^7;forprime(q=2, sqrt(lim\t), if(p==q, next);listput(v,t*q^2))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 20 2011

from math import isqrt
from sympy import primepi, integer_nthroot, primerange
def A179689(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x-sum(primepi(isqrt(x//p**7)) for p in primerange(integer_nthroot(x,7)[0]+1))+primepi(integer_nthroot(x,9)[0])
    return bisection(f,n,n) # Chai Wah Wu, Feb 21 2025

Sum_{n>=1} 1/a(n) = P(2)*P(7) - P(9) = A085548 * A085967 - A085969 = 0.001741..., where P is the prime zeta function. - Amiram Eldar, Jul 06 2020

Title edited by Daniel Forgues, Jan 22 2011

A179696 Numbers with prime signature {7,1,1}, i.e., of form p^7qr with p, q and r distinct primes.

Original entry on oeis.org

1920, 2688, 4224, 4480, 4992, 6528, 7040, 7296, 8320, 8832, 9856, 10880, 11136, 11648, 11904, 12160, 14208, 14720, 15232, 15744, 16512, 17024, 18048, 18304, 18560, 19840, 20352, 20608, 21870, 22656, 23424, 23680, 23936, 25728, 25984, 26240, 26752, 27264

Offset: 1

Vladimir Joseph Stephan Orlovsky, Jul 24 2010

nonn

T. D. Noe, Table of n, a(n) for n = 1..1000
Will Nicholes List of Prime Signatures
OEIS Wiki, Numbers with same prime signature.

Cf. A006881, A007304, A065036, A085986, A085987, A092759, A178739, A179642, A179643, A179644, A179645, A179646, A179664, A179665, A179666, A179667, A179668, A179669, A179670, A179671, A179672, A179688, A179689, A179690, A179691, A179692, A179693, A179694, A179695.

a:= proc(n) option remember; local k;
      for k from 1+ `if` (n=1, 1, a(n-1))
        while sort (map (x-> x[2], ifactors(k)[2]), `>`)<>[7, 1, 1]
      do od; k
    end:
seq (a(n), n=1..40); # Alois P. Heinz, Jan 23 2011

f[n_]:=Sort[Last/@FactorInteger[n]]=={1,1,7}; Select[Range[30000], f]

list(lim)=my(v=List(),t1,t2);forprime(p=2, (lim\6)^(1/7), t1=p^7;forprime(q=2, lim\t1, if(p==q, next);t2=t1*q;forprime(r=q+1, lim\t2, if(p==r,next);listput(v,t2*r)))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 20 2011

from math import isqrt
from sympy import primerange, primepi, integer_nthroot
def A179696(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x+sum((t:=primepi(s:=isqrt(y:=x//r**7)))+(t*(t-1)>>1)-sum(primepi(y//k) for k in primerange(1, s+1)) for r in primerange(integer_nthroot(x,7)[0]+1))+sum(primepi(x//p**8) for p in primerange(integer_nthroot(x,8)[0]+1))-primepi(integer_nthroot(x,9)[0])
    return bisection(f,n,n) # Chai Wah Wu, Mar 27 2025

Title edited by Daniel Forgues, Jan 22 2011

A319075 Square array T(n,k) read by antidiagonal upwards in which row n lists the n-th powers of primes, hence column k lists the powers of the k-th prime, n >= 0, k >= 1.

Original entry on oeis.org

1, 2, 1, 4, 3, 1, 8, 9, 5, 1, 16, 27, 25, 7, 1, 32, 81, 125, 49, 11, 1, 64, 243, 625, 343, 121, 13, 1, 128, 729, 3125, 2401, 1331, 169, 17, 1, 256, 2187, 15625, 16807, 14641, 2197, 289, 19, 1, 512, 6561, 78125, 117649, 161051, 28561, 4913, 361, 23, 1, 1024, 19683, 390625, 823543, 1771561, 371293

Offset: 0

Omar E. Pol, Sep 09 2018

If n = p - 1 where p is prime, then row n lists the numbers with p divisors.

The partial sums of column k give the column k of A319076.

			The corner of the square array is as follows:
         A000079 A000244 A000351  A000420    A001020    A001022     A001026
A000012        1,      1,      1,       1,         1,         1,          1, ...
A000040        2,      3,      5,       7,        11,        13,         17, ...
A001248        4,      9,     25,      49,       121,       169,        289, ...
A030078        8,     27,    125,     343,      1331,      2197,       4913, ...
A030514       16,     81,    625,    2401,     14641,     28561,      83521, ...
A050997       32,    243,   3125,   16807,    161051,    371293,    1419857, ...
A030516       64,    729,  15625,  117649,   1771561,   4826809,   24137569, ...
A092759      128,   2187,  78125,  823543,  19487171,  62748517,  410338673, ...
A179645      256,   6561, 390625, 5764801, 214358881, 815730721, 6975757441, ...
...

Rows 0-13: A000012, A000040, A001248, A030078, A030514, A050997, A030516, A092759, A179645, A179665, A030629, A079395, A030631, A138031.
Other rows n: A030635 (n=16), A030637 (n=18), A137486 (n=22), A137492 (n=28), A139571 (n=30), A139572 (n=36), A139573 (n=40), A139574 (n=42), A139575 (n=46), A173533 (n=52), A183062 (n=58), A183085 (n=60), A261700 (n=100).
Columns 1-15: A000079, A000244, A000351, A000420, A001020, A001022, A001026, A001029, A009967, A009973, A009975, A009981, A009985, A009987, A009991.
Main diagonal gives A093360.
Second diagonal gives A062457.
Third diagonal gives A197987.
Removing the 1's we have A182944/ A182945.
Cf. A006093, A319074, A319076.

PARI
```
T(n, k) = prime(k)^n;
```

T(n,k) = A000040(k)^n, n >= 0, k >= 1.

A179694 Numbers of the form p^6*q^3 where p and q are distinct primes.

Original entry on oeis.org

1728, 5832, 8000, 21952, 85184, 91125, 125000, 140608, 250047, 314432, 421875, 438976, 778688, 941192, 970299, 1560896, 1601613, 1906624, 3176523, 3241792, 3581577, 4410944, 5000211, 5088448, 5359375, 6644672

Offset: 1

Vladimir Joseph Stephan Orlovsky, Jul 24 2010

nonn

T. D. Noe, Table of n, a(n) for n = 1..1000
Will Nicholes, Prime Signatures

Cf. A006881, A007304, A065036, A085986, A085987, A092759, A178739, A179642, A179643, A179644, A179645, A179646, A179664, A179665, A179666, A179667, A179668, A179669, A179670, A179671, A179672, A179688, A179689, A179690, A179691, A179692, A179693.

Cf. A085541, A085966, A085969.

f[n_]:=Sort[Last/@FactorInteger[n]]=={3,6}; Select[Range[10^6], f]

list(lim)=my(v=List(),t);forprime(p=2, (lim\8)^(1/6), t=p^6;forprime(q=2, (lim\t)^(1/3), if(p==q, next);listput(v,t*q^3))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 24 2011

from sympy import primepi, integer_nthroot, primerange
def A179694(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x-sum(primepi(integer_nthroot(x//p**6,3)[0]) for p in primerange(integer_nthroot(x,6)[0]+1))+primepi(integer_nthroot(x,9)[0])
    return bisection(f,n,n) # Chai Wah Wu, Feb 21 2025

Sum_{n>=1} 1/a(n) = P(3)*P(6) - P(9) = A085541 * A085966 - A085969 = 0.000978..., where P is the prime zeta function. - Amiram Eldar, Jul 06 2020

a(n) = A054753(n)^3. - R. J. Mathar, May 05 2023

A179702 Numbers of the form p^4*q^5 where p and q are two distinct primes.

Original entry on oeis.org

2592, 3888, 20000, 50000, 76832, 151875, 253125, 268912, 468512, 583443, 913952, 1361367, 2576816, 2672672, 3557763, 4170272, 5940688, 6940323, 7503125, 8954912, 10504375, 13045131, 20295603, 22632992, 22717712, 29552672, 30074733

Offset: 1

Vladimir Joseph Stephan Orlovsky, Jul 24 2010

nonn

Subsequence of A046312 and of A137493. - R. J. Mathar, Jul 27 2010

T. D. Noe, Table of n, a(n) for n = 1..1000
Will Nicholes, List of Prime Signatures

Cf. A006881, A007304, A065036, A085986, A085987, A092759, A178739, A179642, A179643, A179644, A179645, A179646, A179664, A179665, A179666, A179667, A179668, A179669, A179670, A179671, A179672, A179688, A179689, A179690, A179691, A179692, A179693, A179694, A179695, A179696, A179698, A179699, A179700.

Cf. A085964, A085965, A085969.

fQ[n_] := Sort[Last /@ FactorInteger @n] == {4, 5}; Select[ Range@ 31668000, fQ] (* fixed by Robert G. Wilson v, Aug 26 2010 *)
lst = {}; Do[ If[p != q, AppendTo[lst, Prime@p^4*Prime@q^5]], {p, 12}, {q, 10}]; Take[ Sort@ Flatten@ lst, 27] (* Robert G. Wilson v, Aug 26 2010 *)
Take[Union[First[#]^4 Last[#]^5&/@Flatten[Permutations/@Subsets[ Prime[ Range[30]],{2}],1]],30] (* Harvey P. Dale, Jan 01 2012 *)

list(lim)=my(v=List(),t);forprime(p=2, (lim\16)^(1/5), t=p^5;forprime(q=2, (lim\t)^(1/4), if(p==q, next);listput(v,t*q^4))); vecsort(Vec(v)) \\ Charles R Greathouse IV, Jul 20 2011

from sympy import primepi, integer_nthroot, primerange
def A179702(n):
    def bisection(f,kmin=0,kmax=1):
        while f(kmax) > kmax: kmax <<= 1
        kmin = kmax >> 1
        while kmax-kmin > 1:
            kmid = kmax+kmin>>1
            if f(kmid) <= kmid:
                kmax = kmid
            else:
                kmin = kmid
        return kmax
    def f(x): return n+x-sum(primepi(integer_nthroot(x//p**5,4)[0]) for p in primerange(integer_nthroot(x,5)[0]+1))+primepi(integer_nthroot(x,9)[0])
    return bisection(f,n,n) # Chai Wah Wu, Mar 27 2025

Sum_{n>=1} 1/a(n) = P(4)*P(5) - P(9) = A085964 * A085965 - A085969 = 0.000748..., where P is the prime zeta function. - Amiram Eldar, Jul 06 2020

Edited and extended by Ray Chandler and R. J. Mathar, Jul 26 2010

A377654 Numbers m^2 for which the center part (containing the diagonal) of its symmetric representation of sigma, SRS(m^2), has width 1 and area m.

Original entry on oeis.org

1, 9, 25, 49, 81, 121, 169, 289, 361, 441, 529, 625, 729, 841, 961, 1089, 1369, 1521, 1681, 1849, 2209, 2401, 2601, 2809, 3025, 3249, 3481, 3721, 4225, 4489, 4761, 5041, 5329, 6241, 6561, 6889, 7225, 7569, 7921, 8649, 9025, 9409, 10201, 10609, 11449, 11881, 12321, 12769, 13225, 14161, 14641, 15129, 15625

Offset: 1

Hartmut F. W. Hoft, Nov 03 2024

nonn

Since for numbers m^2 in the sequence the width at the diagonal of SRS(m^2) is 1, the area m of its center part is odd so that this sequence is a proper subsequence of A016754 and since SRS(m^2) has an odd number of parts it is a proper subsequence of A319529. The smallest odd square not in this sequence is 225 = 15^2.  SRS(225) is {113, 177, 113}, its center part has maximum width 2, its width at the diagonal is 1.
The k+1 parts of SRS(p^(2k)), p an odd prime and k >= 0, through the diagonal including the center part have areas (p^(2k-i) + p^i)/2 for 0 <= i <= k. They form a strictly decreasing sequence. Since p^(2k) has 2k+1 divisors and SRS(p^(2k)) has 2k+1 parts, all of width 1 (A357581), the even powers of odd primes form a proper subsequence of A244579. For the subsequence of squares of odd primes p, SRS(p^2) consists of the 3 parts { (p^2 + 1)/2, p, (p^2 + 1)/2 } see A001248, A247687 and A357581.
The areas of the parts of SRS(m^2) need not be in descending order through the diagonal as a(112) = 275^2 = 75625 with SRS(75625) = (37813, 7565, 3443, 1525, 715, 738, 275, 738, 715, 1525, 3443, 7565, 37813) demonstrates.
An equivalent description of the sequence is: The center part of SRS(m^2) has width 1, m is odd, and A249223(m^2, m-1) = 0.
Conjectures (true for all a(n) <= 10^8):
(1) The central part of SRS(a(n)) is the minimum of all parts of SRS(a(n)), 1 <= n.
(2) The terms in this sequence are the squares of the terms in A244579.

			The center part of SRS(a(3)) = SRS(25) has area 5, all 3 parts have width 1, and 25 with 3 divisors also belongs to A244579.
The center part of SRS(a(7)) = SRS(169) has area 13, all 3 parts have width 1, and 169 with 3 divisors also belongs to A244579.
The center part of SRS(a(10)) = SRS(441) has area 21 and width 1, but the maximum width of SRS(441) is 2. Number 441 has 9 divisors and SRS(441) has 7 parts while 21 has 4 divisors and SRS(21) has 4 parts so that 21 is in A244579 while 441 is not.

Cf. A001248, A030514, A030516, A030629, A030631, A179645 (even prime power sequences).

Cf. A003056, A016754, A237591, A237593, A244579, A247687, A249223, A319529, A341969, A357581.

(* t237591 and partsSRS compute rows in A237270 and A237591, respectively *)
(* t249223 and widthPattern are also defined in A376829 *)
row[n_] := Floor[(Sqrt[8 n+1]-1)/2]
t237591[n_] := Map[Ceiling[(n+1)/#-(#+1)/2]-Ceiling[(n+1)/(#+1)-(#+2)/2]&, Range[row[n]]]
partsSRS[n_] := Module[{widths=t249223[n], legs=t237591[n], parts, srs}, parts=widths legs; srs=Map[Apply[Plus, #]&, Select[SplitBy[Join[parts, Reverse[parts]], #!=0&], First[#]!=0&]]; srs[[Ceiling[Length[srs]/2]]]-=Last[widths]; srs]
t249223[n_] := FoldList[#1+(-1)^(#2+1)KroneckerDelta[Mod[n-#2 (#2+1)/2, #2]]&, 1, Range[2, row[n]]]
widthPattern[n_] := Map[First, Split[Join[t249223[n], Reverse[t249223[n]]]]]
centerQ[n_] := Module[{pS=partsSRS[n]}, Sqrt[n]==pS[[(Length[pS]+1)/2]]]/;OddQ[n]
widthQ[n_] := Module[{wP=SplitBy[widthPattern[n], #!=0&]}, wP[[(Length[wP]+1)/2]]]=={1}/;OddQ[n]
a377654[m_, n_] := Select[Map[#^2&, Range[m, n, 2]], centerQ[#]&&widthQ[#]&]/;OddQ[m]
a377654[1, 125]

cp's OEIS Frontend

A182944 Square array A(i,j), i >= 1, j >= 1, of prime powers prime(i)^j, by descending antidiagonals.

Views

Author

Keywords

Comments

Examples

Links

Crossrefs

Programs

Mathematica

Formula

Extensions

A189990 Numbers with prime factorization p^2*q^6.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Python

Formula

A115975 Numbers of the form p^k, where p is a prime and k is a Fibonacci number.

Views

Author

Keywords

Links

Crossrefs

Programs

Mathematica

PARI

Extensions

A369938 Numbers whose maximal exponent in their prime factorization is a power of 2.

Views

Author

Keywords

Comments

Links

Crossrefs

Programs

Mathematica

PARI

A179689 Numbers with prime signature {7,2}, i.e., of form p^7*q^2 with p and q distinct primes.

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Formula

Extensions

A179696 Numbers with prime signature {7,1,1}, i.e., of form p^7*q*r with p, q and r distinct primes.

Views

Author

Keywords

Links

Crossrefs

Programs

Maple

Mathematica

PARI

Python

Extensions

A319075 Square array T(n,k) read by antidiagonal upwards in which row n lists the n-th powers of primes, hence column k lists the powers of the k-th prime, n >= 0, k >= 1.

Views

Author

Keywords

Comments

Examples

Crossrefs

Programs

PARI

Formula

A179694 Numbers of the form p^6*q^3 where p and q are distinct primes.

A179696 Numbers with prime signature {7,1,1}, i.e., of form p^7qr with p, q and r distinct primes.