Znajdź Emirps!

20

Emirp jest non-palindromiczna prime, który po odwróceniu, jest również podstawowym.

Lista 10 podstawowych emirps znajduje się w OEIS . Pierwsze sześć to:

13, 17, 31, 37, 71, 73

Jednak ze względu na zasadę odwrócenia emirps są różne w każdej bazie. Na przykład pierwsze sześć binarnych emirps to:

Bin  | 1011, 1101, 10111, 11101, 101001, 100101
Dec  | (11 , 13  , 23   , 29   , 37    , 41   ) 

... a w systemie szesnastkowym są to:

Hex |  17, 1F, 35, 3B, 3D, 53
Dec | (23, 31, 53, 59, 61, 83)

Ciekawostka: nie ma emirps w jednoskładnikowa jak każda liczba jest palindrom.


Wyzwanie

Twoim zadaniem jest utworzenie funkcji (lub pełnego programu), która pobierze dwa parametry, n i b , i wygeneruje listę pierwszych n emirps w bazie b .

Zasady / szczegóły:

  • n ib są dodatnimi liczbami całkowitymi większymi niż0 .
  • Możesz założyć 2b16 : to znaczy, że podstawa będzie między wartością binarną a szesnastkową.
  • Powinieneś być w stanie obliczyć dla wartości od n do 100 . 100
  • Wygenerowana lista może znajdować się w bazie b lub standardowej podstawie liczb całkowitych twojego języka, o ile podasz to w odpowiedzi.
  • Wbudowane testy emirp są niedozwolone (wbudowane testy pierwotności są w porządku)
  • Nie można na stałe zakodować emirps ani odczytać z plików zewnętrznych.
  • Standardowe luki są jak zawsze zakazane.
  • To jest , więc wygrywa najkrótsza odpowiedź (w bajtach).

Przypadki testowe

Dla każdego przypadku testowego zawarłem listę w bazie bi jej 10 ekwiwalentów.

B = 2, N = 10

BIN: [1011, 1101, 10111, 11101, 100101, 101001, 101011, 101111, 110101, 111101]
DEC: [11, 13, 23, 29, 37, 41, 43, 47, 53, 61] 


B = 3, N = 5

BASE3: [12, 21, 102, 201, 1011]
DEC:   [5, 7, 11, 19, 31]


B = 12, N = 7

BASE12: [15, 51, 57, 5B, 75, B5, 107]
DEC: [17, 61, 67, 71, 89, 137, 151]


B = 16, N = 4

HEX: [17, 1F, 35, 3B]
DEC: [23, 31, 53, 59] 

Możesz dalej testować swój program na moim (nie golfowym) przykładzie Pythona na repl.it

FlipTack
źródło

Odpowiedzi:

6

Galaretka , 16 bajtów

bµU,ḅ⁹QÆPḄ=3
⁸ç#

TryItOnline!

W jaki sposób?

bµU,ḅ⁹QÆPḄ=3 - Link 1, in-sequence test: n, b
b            - convert n to base b - a list
 µ           - monadic chain separation
  U          - reverse the list
   ,         - pair with the list
     ⁹       - link's right argument, b
    ḅ        - convert each of the two lists from base b
      Q      - get unique values (if palindromic a list of only one item)
       ÆP    - test if prime(s) - 1 if prime, 0 if not
         Ḅ   - convert to binary
          =3 - equal to 3? (i.e. [reverse is prime, forward is prime]=[1,1])

⁸ç# - Main link: b, N
  # - count up from b *see note, and find the first N matches (n=b, n=b+1, ...) for:
 ç  - last link (1) as a dyad with left argument n and right argument
⁸   - left argument, b

* Uwaga bw bazie bjest [1,0], która po odwróceniu jest [0,1]tym 1, co jest, a nie liczbą pierwszą; cokolwiek mniejszego niż bjedna cyfra w bazie, ba zatem palindromiczna.

Jonathan Allan
źródło
Gratulujemy wygranej!
FlipTack
8

05AB1E , 17 bajtów

Wykorzystuje kodowanie CP-1252 .

Kolejność wprowadzania to n, b
Wyjście jest w bazie-10.

µN²BÂD²öpŠÊNpPD–½

Wypróbuj online!

Wyjaśnienie

                    # implicit input a,b
µ                   # loop until counter is a
 N²B                # convert current iteration number to base b
    ÂD              # create 2 reversed copies
      ²ö            # convert one reversed copy to base 10
        p           # check for primality
         ŠÊ         # compare the normal and reversed number in base b for inequality
           Np       # check current iteration number for primality
             P      # product of all
              D     # duplicate
               –    # if 1, print current iteration number
                ½   # if 1, increase counter
Emigna
źródło
4

Mathematica, 70 bajtów

Cases[Prime@Range@437,p_/;(r=p~IntegerReverse~#2)!=p&&PrimeQ@r]~Take~#&

Działa dla 0 <= n <= 100i 2 <= b <= 16. Z listy Prime@Range@437pierwszych 437liczb pierwszych, znaleźć Cases pgdzie IntegerReverse rod pw bazie #2nie jest równa p, a także jest pierwsza, a następnie podjąć pierwszą #takie p.

Oto 95-bajtowe rozwiązanie, które działa dla dowolnego n>=0i b>=2:

(For[i=1;a={},Length@a<#,If[(r=IntegerReverse[p=Prime@i,#2])!=p&&PrimeQ@r,a~AppendTo~p],i++];a)&
ngenisis
źródło
+1 IntegerReverse. Oczywiście! Ładny.
DavidC
79 bajtów na rozwiązanie arbitralne-nb; 77 bajtów, jeśli Reaping jest dozwolone w stopce:For[i=j=0,j<#,If[(r=IntegerReverse[p=Prime@++i,#2])!=p&&PrimeQ@r,j++;Sow@p]]&
Roman
3

Perl, 262 bajtów

($b,$n)=@ARGV;$,=',';sub c{my$z;for($_=pop;$_;$z=(0..9,a..z)[$_%$b].$z,$_=($_-$_%$b)/$b){};$z}sub d{my$z;for(;c(++$z)ne@_[0];){}$z}for($p=2;@a<$n;$p++){$r=qr/^1?$|^(11+?)\1+$/;(c($p)eq reverse c$p)||((1x$p)=~$r)||(1x d($x=reverse c($p)))=~$r?1:push@a,c($p);}say@a

Czytelny:

($b,$n)=@ARGV;
$,=',';
sub c{
    my$z;
    for($_=pop;$_;$z=(0..9,a..z)[$_%$b].$z,$_=($_-$_%$b)/$b){};
    $z
}
sub d{
    my$z;
    for(;c(++$z)ne@_[0];){}
    $z
}
for($p=2;@a<$n;$p++){
    $r=qr/^1?$|^(11+?)\1+$/;
    (c($p)eq reverse c$p)||((1x$p)=~$r)||(1x d($x=reverse c($p)))=~$r?1:push@a,c($p)
}
say@a

ckonwertuje daną liczbę na bazę $bi dkonwertuje daną liczbę z bazy z $bpowrotem na dziesiętną, znajdując pierwszą liczbę, która zwraca wspomnianą $bliczbę bazową po przekazaniu do c. Pętla for sprawdza następnie, czy jest to palindrom i czy obie liczby są pierwsze, używając wyrażenia regularnego złożonego.

Gabriel Benamy
źródło
3

Matematyka 112 bajtów

Cases[Table[Prime@n~IntegerDigits~#2,{n,500}],x_/;x!=(z=Reverse@x)&&PrimeQ[z~(f=FromDigits)~#2]:>x~f~#2]~Take~#&

Przykład

Znajdź pierwsze 10 Emipsów w heksie; zwróć je w postaci dziesiętnej.

Cases[Table[Prime@n~IntegerDigits~#2, {n, 500}], 
x_ /; x != (z = Reverse@x) && PrimeQ[z~(f = FromDigits)~#2] :> x~f~#2]~Take~# &[10, 16]


{23, 31, 53, 59, 61, 83, 89, 113, 149, 179}

Nie golfił

Take[Cases[                                             (* take #1 cases; #1 is the first input argument *)
   Table[IntegerDigits[Prime[n], #2], {n, 500}],        (* from a list of the first 500 primes, each displayed as a list of digits in base #2 [second argument] *) 
   x_ /;                                                (* x, a list of digits, such that *)
   x != (z = Reverse[x]) && PrimeQ[FromDigits[z, #2]]   (* the reverse of the digits is not the same as the list of digits; and the reverse list, when composed, also constitutes a prime *)
   :> FromDigits[x, #2]],                               (* and return the prime *)
   #1] &                                                (* [this is where #1 goes, stating how many cases to Take] *)
DavidC
źródło
2

Perl 6 , 91 bajtów

->\n,\b{(grep {.is-prime&&{$_ ne.flip &&.parse-base(b).is-prime}(.base(b).flip)},1..*)[^n]}

Zwraca listę emirps w bazie 10.

Sean
źródło
81 bajtów
Jo King
2

Python 3 , 232 214 191 188 bajtów

p=lambda n:all(n%i for i in range(2,n))
def f(b,n):
 s=lambda n:(''if n<b else s(n//b))+f'{n%b:X}';l=[];i=3
 while n:i+=1;c=s(i);d=c[::-1];a=(c!=d)*p(i)*p(int(d,b));l+=[c]*a;n-=a
 return l

Wypróbuj online!

movatica
źródło
1
200 bajtów
Herman L
Dobry chwyt! Sprowadzono to do 191 bajtów
movatica
Niezły, @JoKing!
movatica
2

C 293 286 261 bajtów

Poprawiony przez @ceilingcat , 261 bajtów:

v,t,i,j,c,g,s[9],r[9],b;main(n,a)int**a;{for(b=n=atoi(a[1]);g^atoi(a[2]);t|v|!wcscmp(s,r)||printf("%u ",n,++g)){i=j=0;for(c=++n;s[i]=c;c/=b)s[i++]=c%b+1;for(;r[j]=i;)r[j++]=s[--i];p(n);for(t=v;r[i];)c+=~-r[i]*pow(b,i++);p(c);}}p(n){for(j=1,v=0;++j<n;n%j||v++);}

Wypróbuj online!

(Ta osoba jest jak ciągłe śledzenie mnie wokół PPCG i poprawianie moich rzeczy w komentarzach, a jak tylko odpowiem, aby mu podziękować, po prostu usuwa komentarz i znika lol. Welp, jeszcze raz dziękuję!)


Poprawiony przez @movatica , 286 bajtów:

u,v,t,i,j,c,n,g;main(int x,char**a){char s[9],r[9],b=n=atoi(a[1]);x=atoi(a[2]);for(;g^x;){i=j=0;for(c=++n;c;c/=b)s[i++]=c%b+1;s[i]=c=0;for(;i;r[j++]=s[--i]);r[j]=0;p(n);t=v;for(;r[i];)c+=(r[i]-1)*pow(b,i++);p(c);t|v|!strcmp(s,r)?:printf("%u ",n,++g);}}p(n){for(u=1,v=0;++u<n;n%u?:v++);}

Wypróbuj online!


Moja oryginalna odpowiedź, 293 bajty:

u,v,t,i,j,c,n,g;main(int x,char**a){char s[9],r[9],b=n=atoi(a[1]);x=atoi(a[2]);for(++n;g^x;++n){i=j=0;for(c=n;c;c/=b)s[i++]=c%b+1;s[i]=c=0;for(;i;r[j++]=s[--i]);r[j]=0;p(n);t=v;for(--i;r[++i];)c+=(r[i]-1)*pow(b,i);p(c);t|v|!strcmp(s,r)?:printf("%u ",n,++g);}}p(n){for(u=1,v=0;++u<n;n%u?:v++);}

Kompiluj gcc emirp.c -o emirp -lmi uruchamiaj z ./emirp <b> <n>. Wyświetla emirps rozdzielone spacjami w bazie-10.

OverclockedSanic
źródło
@FlipTack Masz rację. Będę musiał to naprawić jutro.
OverclockedSanic
@FlipTack Naprawiono i przetestowano, aby upewnić się, że pomyślnie przejdzie testy. Czy to jest dobre?
OverclockedSanic
Jasne, że tak! Witamy w golfa code.
FlipTack
1
Dobra robota! Przesunąłem niektórych operatorów przyrostowych, aby obniżyć cię do 286
movatica
1
@movatica Awesome! Dodałem twoje ulepszenia do mojej odpowiedzi. Dzięki!
OverclockedSanic
1

JavaScript (ES6), 149 148 141 140 bajtów

Zwraca rozdzieloną spacjami listę emirps w bazie b. (Może być o 2 bajty krótszy, zwracając listę dziesiętną.)

f=(b,n,i=2)=>n?((p=(n,k=n)=>--k<2?k:n%k&&p(n,k))(i)&p(k=parseInt([...j=i.toString(b)].reverse().join``,b))&&k-i&&n--?j+' ':'')+f(b,n,i+1):''

Przypadki testowe

Arnauld
źródło
1

Python 2 , 133 bajty

p=lambda n:all(n%i for i in range(2,n))
b,n=input()
i=b
while n:
 j=i=i+1;r=0
 while j:r=r*b+j%b;j/=b
 if(i-r)*p(i)*p(r):print i;n-=1

Wypróbuj online!

Wyświetla każdą liczbę w nowej linii, w bazie 10

ujemna siódemka
źródło
0

APL (NARS), 87 znaków, 174 bajty

r←a f w;i
i←1⋄r←⍬
→2×⍳∼{∼0π⍵:0⋄k≡v←⌽k←{(a⍴⍨⌊1+a⍟⍵)⊤⍵}⍵:0⋄0πa⊥v:1⋄0}i+←1⋄r←r,i⋄→2×⍳w>≢r

Wynik będzie w bazie 10. Test i wyniki:

  3 f 1
5 
  2 f 10
11 13 23 29 37 41 43 47 53 61 
  3 f 5
5 7 11 19 31 
  12 f 7
17 61 67 71 89 137 151 
  16 f 4
23 31 53 59 

{(⍺⍴⍨⌊1+⍺⍟⍵)⊤⍵}zrobi konwersja w bazie , tablicy całkowitej wyniku; 0π⍵zwróci true [1], jeśli jest liczbą pierwszą, zwróci 0.

RosLuP
źródło