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!)
A273201 Integers which are prime power-like but not prime powers. 3
21, 33, 39, 65, 85, 95, 115, 133, 145, 155, 161, 185, 203, 205, 215, 217, 235, 259, 261, 265, 279, 287, 295, 301, 305, 329, 335, 341, 355, 365, 371, 395, 407, 413, 415, 427, 445, 451, 469, 473, 481, 485, 497 (list; graph; refs; listen; history; text; internal format)
OFFSET
1,1
COMMENTS
For an integer n>0 and not the unity we define DTD(n) to be the difference table of the divisors of n. We say that DTD(n) is positive if all entries in the table are positive and we call DTD(n) monotone if all rows and all columns of the table are nondecreasing (reading from left to right and from top to bottom).
We define an integer n to be prime power-like if and only if DTD(n) is positive and monotone. All prime powers (in the sense of A246655 (but not in the sense of A000961)) are prime power-like integers. Sequence A273200 provides the prime power-like integers. This sequence (A273201) lists the integers which are prime power-like but not prime powers.
Thus we have the inclusions A000040 < A246655 < A273200 and the union A273200 = A273201 U A246655. Integers which have a positive but not monotone DTD are listed in A273199. Integers with a positive DTD are listed in A273130.
LINKS
EXAMPLE
95 is in this sequence because the DTD of 95 has positive entries and nondecreasing rows and columns:
[ 1 5 19 95]
[ 4 14 76]
[10 62]
[52]
MATHEMATICA
pplikeQ[n_] := Module[{T, DTD, DTD2}, If[n == 1 || PrimePowerQ[n], Return[False]]; T = Divisors[n]; DTD = Table[Differences[T, k], {k, 0, Length[T]-1}]; If[AnyTrue[Flatten[DTD], NonPositive], Return[False]]; DTD2 = Transpose[PadRight[#, Length[T], Infinity]& /@ DTD]; AllTrue[DTD, OrderedQ] && AllTrue[DTD2, OrderedQ]];
Select[Range[500], pplikeQ] (* Jean-François Alcover, Jun 28 2019 *)
PROG
(Sage) # uses[is_prime_power_like from A273200]
# Compare script in A273200.
def is_A273201(n):
return not is_prime_power(n) and is_prime_power_like(n)
print(list(filter(is_A273201, range(1, 500))))
CROSSREFS
Sequence in context: A368949 A298855 A354813 * A260730 A119973 A279229
KEYWORD
nonn
AUTHOR
Peter Luschny, May 17 2016
STATUS
approved

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 April 19 12:14 EDT 2024. Contains 371792 sequences. (Running on oeis4.)