Séquence de nombres denses

12

OEIS: A167171

Un nombre dense est un nombre qui a exactement autant de diviseurs premiers que les diviseurs non premiers (y compris 1 et lui-même comme diviseurs). De manière équivalente, c'est soit un nombre premier, soit un produit de deux nombres premiers distincts. Les 100 premiers nombres denses sont:

2, 3, 5, 6, 7, 10, 11, 13, 14, 15, 17, 19, 21, 22, 23, 26, 29, 31, 33, 34, 35, 37, 38, 39, 41, 43, 46, 47, 51, 53, 55, 57, 58, 59, 61, 62, 65, 67, 69, 71, 73, 74, 77, 79, 82, 83, 85, 86, 87, 89, 91, 93, 94, 95, 97, 101, 103, 106, 107, 109, 111, 113, 115, 118, 119, 122, 123, 127, 129, 131, 133, 134, 137, 139, 141, 142, 143, 145, 146, 149, 151, 155, 157, 158, 159, 161, 163, 166, 167, 173, 177, 178, 179, 181, 183, 185, 187, 191, 193, 194

Étant donné un entier non négatif n, sortie dense(n). npeut être indexé 0 ou indexé 1.

Implémentation de référence (Sage)

import itertools

def dense_numbers():
    n = 1
    while True:
        prime_divisors = [x for x in divisors(n) if x.is_prime()]
        non_prime_divisors = [x for x in divisors(n) if not x.is_prime()]
        if len(prime_divisors) == len(non_prime_divisors):
            yield n
        n += 1

N = 20

print itertools.islice(dense_numbers(), N, N+1).next()

Essayez-le en ligne

Mego
la source
Tant de séquences de nombres premiers ... je ne savais pas tant existé
Beta Decay
2
@ βετѧΛєҫαγ Il existe également des nombres premiers appelés Sexy Primes (͡ ° ͜ʖ ͡ °).
Adnan
@Adnan Oh myy; D
Beta Decay
À quoi sert la valeur maximale n?
R. Kap
@ R.Kap Aussi haut que la langue de votre choix peut aller.
Mego

Réponses:

3

Gelée , 9 octets

ÆE²Sḍ2µ#Ṫ

Lit à partir de STDIN et utilise l'indexation basée sur 1. Essayez-le en ligne!

Comment ça fonctionne

ÆE²Sḍ2µ#Ṫ  Main link. No arguments. Implicit argument: 0

      µ#   Read an integer n from STDIN and execute the chain to the left for
           k = 0, 1, 2, ... until n of them return a truthy value.
           Return the array of matches.
ÆE           Compute the exponents of k's prime factorization.
  ²          Square each exponent.
   S         Compute the sum of all squares.
    ḍ2       Test if 2 is divisible by the result (true iff  the sum is 1 or 2).
        Ṫ  Tail; extract the last (n-th) matching value of k.
Dennis
la source
2

En fait , 12 octets

Tous les crédits à Dennis pour son algorithme .

`w♂N;*2%Y`╓N

Essayez-le en ligne!

`w♂N;*2%Y`╓N

`        `     define a function
 w             prime factorization in exponent form:
               18 = (2^1)*(3^2) becomes [[2,1],[3,2]]
  ♂N           get the last element (exponent) of each sublist
    ;*         dot-product with self; equivalent to squaring
               each item and then taking the sum
      2%Y      test divisibility by 2
          ╓    first (input) solutions to the above function
           N   get the last element.
Leaky Nun
la source
1

05AB1E, 12 11 octets

1 indexé

µ             # while counter != input
 NÑ           # get divisors of current number
   p          # check if prime
    D         # duplicate
     O        # sum one copy
      s_O     # invert and sum the other copy
         Q½   # if equal increase counter

Essayez-le en ligne

Emigna
la source
1

Brachylog , 17 octets

:1yt.
1<.=$p#dl<3

Essayez-le en ligne!

Prédicat 0 (prédicat principal)

:1yt.
:1y     Find the first (input) solutions of predicate 1
   t    Last element
    .   Unify with output

Prédicat 1 (prédicat auxiliaire)

1<.=$p#dl<3
1<.            1 < output
  .=           assign a value to output
  . $p#d       output's prime factorization contains no duplicate
        l      and the length
         <3    is less than three
Leaky Nun
la source
0

R, 93 octets

dense=function(n){a=b=0;for(i in which(!n%%1:n))if(which(!i%%2:i)+1==i)a=a+1 else b=b+1;a==b}

Il a tendance à lancer un avertissement. Ce n'est pas vraiment un problème. Autoriser l'avertissement me fait gagner 5 octets.

Non golfé

dense=function(n){
     a=b=0                                #Initializing
     factors = which(!n%%1:n)             #Finds all factors
     for(i in factors)                    #Loops through factors
         prime = which(!i%%2:i)+1==i      #Tests if current factor is prime. If it is -- the first term in this vector will be TRUE. Otherwise, it will be false.
           if (prime) a=a+1 else b=b+1    #If first term is true, add 1 to a. Else add one to b. 
      return(a==b)                        #Test equality of a and b.
}
user5957401
la source
Vous ne pouvez pas utiliser l' +=opérateur pour enregistrer 2 octets?
R. Kap
Malheureusement, R n'a pas d'opérateurs d'incrémentation utiles comme +=ou a++. Parfois, il peut y avoir des moyens plus courts (profitant principalement de la structure de la boucle), mais je n'en connais pas ici.
user5957401
0

Python, 79 octets

f=lambda n,k=2:n<1or-~f(n-(sum((k%i<1)+2*(k%i**2<1)for i in range(2,k))<3),k+1)

Utilise l'indexation basée sur 1. Testez-le sur Ideone .

Dennis
la source
0

PHP, 118 octets

for($i=1;!$o=$s[$argn];$r[2]?:$t+=2*$$i=1,$t?:$s[]=$i)for($t=0,$r=[],$n=++$i;$n;$n--)$i%$n?:$t+=${$r[]=$n}?:-1;echo$o;

Essayez-le en ligne!

Jörg Hülsermann
la source
0

Axiome, 102 octets

f(n:PI):PI==(i:=1;repeat(i:=i+1;a:=divisors(i);2*#[x for x in a|prime?(x)]=#a=>(n=1=>break;n:=n-1));i)

ungolf et résultat

-- 1 base Indexed: return the n_th number i that has 2*#divisorsPrimeOf(i)=#divisors(i)
ff(n:PI):PI==
     i:=1
     repeat
        i:=i+1
        a:=divisors(i)
        2*#[x for x in a|prime?(x)]=#a=>(n=1=>break;n:=n-1)
     i

(3) -> [f(i)  for i in 1..23]
   (3)  [2,3,5,6,7,10,11,13,14,15,17,19,21,22,23,26,29,31,33,34,35,37,38]
                                               Type: List PositiveInteger
RosLuP
la source