Calculer la somme des n premiers nombres premiers

15

Je suis surpris que ce défi ne soit pas déjà là, car il est si évident. (Ou je suis surpris de ne pas l'avoir trouvé et n'importe qui le marquera comme doublon.)

Tâche

Étant donné un entier non négatif , calculez la somme des n premiers nombres premiers et sortez-la.nn

Exemple 1

Pour , les cinq premiers nombres premiers sont:n=5

  • 2
  • 3
  • 5
  • sept
  • 11

La somme de ces nombres est , donc le programme doit sortir 28 .2+3+5+7+11=2828

Exemple # 2

Pour , les "premiers zéro" premiers sont nuls. Et la somme d'aucun nombre est - bien sûr - 0 .n=00

Règles

  • Vous pouvez utiliser des fonctions intégrées, par exemple, pour vérifier si un nombre est premier.
  • Il s'agit de , donc le plus petit nombre d'octets dans chaque langue gagne!
xanoetux
la source
2
OEIS - A7504 (à part: LOL à ceci dans la section des formules, "a (n) = A033286 (n) - A152535 (n).")
Jonathan Allan
@JonathanAllan: Lié, mais pas équivalent. Je pense que c'est une différence importante si vous cochez des nombres premiers dans la plage ou un certain nombre de nombres premiers. Ce que les deux tâches ont en commun est a) de vérifier si un nombre est premier et b) de résumer les nombres - ce qui est commun à de nombreuses tâches de code-golf ici.
xanoetux

Réponses:

15

6502 routine de code machine , 75 octets

A0 01 84 FD 88 84 FE C4 02 F0 32 E6 FD A0 00 A5 FD C9 04 90 1F 85 64 B1 FB 85
65 A9 00 A2 08 06 64 2A C5 65 90 02 E5 65 CA D0 F4 C9 00 F0 DC C8 C4 FE D0 DB
A5 FD A4 FE 91 FB C8 D0 C8 A9 00 18 A8 C4 FE F0 05 71 FB C8 D0 F7 60

Attend un pointeur vers un stockage temporaire dans $fb/ $fcet le nombre de nombres premiers à résumer $2. Renvoie la somme dans A(le registre accu).

Jamais fait quelques vérifications de base dans le code machine 6502, alors le voici enfin;)

Notez que cela commence à donner des résultats incorrects pour les entrées> = 14. Ceci est dû au débordement, le code fonctionne avec la plage de nombres "naturels" de la plate-forme 8 bits qui est 0 - 255pour non signé .

Démontage commenté

; function to sum the first n primes
;
; input:
;   $fb/$fc: pointer to a buffer for temporary storage of primes
;   $2:      number of primes to sum (n)
; output:
;   A:       sum of the first n primes
; clobbers:
;   $fd:     current number under primality test
;   $fe:     number of primes currently found
;   $64:     temporary numerator for modulo check
;   $65:     temporary divisor for modulo check
;   X, Y
 .primesum:
A0 01       LDY #$01            ; init variable for ...
84 FD       STY $FD             ; next prime number to test
88          DEY                 ; init number of found primes
 .mainloop:
84 FE       STY $FE             ; store current number of found primes
C4 02       CPY $02             ; compare with requested number
F0 32       BEQ .sum            ; enough primes -> calculate their sum
 .mainnext:
E6 FD       INC $FD             ; check next prime number
A0 00       LDY #$00            ; start check against first prime number
 .primecheckloop:
A5 FD       LDA $FD             ; load current number to check
C9 04       CMP #$04            ; smaller than 4?
90 1F       BCC .isprime        ; is a prime (shortcut to get list started)
85 64       STA $64             ; store to temp as numerator
B1 FB       LDA ($FB),Y         ; load from prime number table
85 65       STA $65             ; store to temp as divisor
A9 00       LDA #$00            ; init modulo to 0
A2 08       LDX #$08            ; iterate over 8 bits
 .bitloop:
06 64       ASL $64             ; shift left numerator
2A          ROL A               ; shift carry into modulo
C5 65       CMP $65             ; compare with divisor
90 02       BCC .bitnext        ; smaller -> to next bit
E5 65       SBC $65             ; otherwise subtract divisor
 .bitnext:
CA          DEX                 ; next bit
D0 F4       BNE .bitloop
C9 00       CMP #$00            ; compare modulo with 0
F0 DC       BEQ .mainnext       ; equal? -> no prime number
C8          INY                 ; next index in prime number table
C4 FE       CPY $FE             ; checked against all prime numbers?
D0 DB       BNE .primecheckloop ; no -> check next
 .isprime:
A5 FD       LDA $FD             ; prime found
A4 FE       LDY $FE             ; then store in table
91 FB       STA ($FB),Y
C8          INY                 ; increment number of primes found
D0 C8       BNE .mainloop       ; and repeat whole process
 .sum:
A9 00       LDA #$00            ; initialize sum to 0
18          CLC
A8          TAY                 ; start adding table from position 0
 .sumloop:
C4 FE       CPY $FE             ; whole table added?
F0 05       BEQ .done           ; yes -> we're done
71 FB       ADC ($FB),Y         ; add current entry
C8          INY                 ; increment index
D0 F7       BNE .sumloop        ; and repeat
 .done:
60          RTS

Exemple de programme assembleur C64 utilisant la routine:

Démo en ligne

Code dans la syntaxe ca65 :

.import primesum   ; link with routine above

.segment "BHDR" ; BASIC header
                .word   $0801           ; load address
                .word   $080b           ; pointer next BASIC line
                .word   2018            ; line number
                .byte   $9e             ; BASIC token "SYS"
                .byte   "2061",$0,$0,$0 ; 2061 ($080d) and terminating 0 bytes

.bss
linebuf:        .res    4               ; maximum length of a valid unsigned
                                        ; 8-bit number input
convbuf:        .res    3               ; 3 BCD digits for unsigned 8-bit
                                        ; number conversion
primebuf:       .res    $100            ; buffer for primesum function

.data
prompt:         .byte   "> ", $0
errmsg:         .byte   "Error parsing number, try again.", $d, $0

.code
                lda     #$17            ; set upper/lower mode
                sta     $d018

input:
                lda     #<prompt        ; display prompt
                ldy     #>prompt
                jsr     $ab1e

                lda     #<linebuf       ; read string into buffer
                ldy     #>linebuf
                ldx     #4
                jsr     readline

                lda     linebuf         ; empty line?
                beq     input           ; try again

                lda     #<linebuf       ; convert input to int8
                ldy     #>linebuf
                jsr     touint8
                bcc     numok           ; successful -> start processing
                lda     #<errmsg        ; else show error message and repeat
                ldy     #>errmsg
                jsr     $ab1e
                bcs     input

numok:          
                sta     $2
                lda     #<primebuf
                sta     $fb
                lda     #>primebuf
                sta     $fc
                jsr     primesum        ; call function to sum primes
                tax                     ; and ...
                lda     #$0             ; 
                jmp     $bdcd           ; .. print result

; read a line of input from keyboard, terminate it with 0
; expects pointer to input buffer in A/Y, buffer length in X
.proc readline
                dex
                stx     $fb
                sta     $fc
                sty     $fd
                ldy     #$0
                sty     $cc             ; enable cursor blinking
                sty     $fe             ; temporary for loop variable
getkey:         jsr     $f142           ; get character from keyboard
                beq     getkey
                sta     $2              ; save to temporary
                and     #$7f
                cmp     #$20            ; check for control character
                bcs     checkout        ; no -> check buffer size
                cmp     #$d             ; was it enter/return?
                beq     prepout         ; -> normal flow
                cmp     #$14            ; was it backspace/delete?
                bne     getkey          ; if not, get next char
                lda     $fe             ; check current index
                beq     getkey          ; zero -> backspace not possible
                bne     prepout         ; skip checking buffer size for bs
checkout:       lda     $fe             ; buffer index
                cmp     $fb             ; check against buffer size
                beq     getkey          ; if it would overflow, loop again
prepout:        sei                     ; no interrupts
                ldy     $d3             ; get current screen column
                lda     ($d1),y         ; and clear 
                and     #$7f            ;   cursor in
                sta     ($d1),y         ;   current row
output:         lda     $2              ; load character
                jsr     $e716           ;   and output
                ldx     $cf             ; check cursor phase
                beq     store           ; invisible -> to store
                ldy     $d3             ; get current screen column
                lda     ($d1),y         ; and show
                ora     #$80            ;   cursor in
                sta     ($d1),y         ;   current row
                lda     $2              ; load character
store:          cli                     ; enable interrupts
                cmp     #$14            ; was it backspace/delete?
                beq     backspace       ; to backspace handling code
                cmp     #$d             ; was it enter/return?
                beq     done            ; then we're done.
                ldy     $fe             ; load buffer index
                sta     ($fc),y         ; store character in buffer
                iny                     ; advance buffer index
                sty     $fe
                bne     getkey          ; not zero -> ok
done:           lda     #$0             ; terminate string in buffer with zero
                ldy     $fe             ; get buffer index
                sta     ($fc),y         ; store terminator in buffer
                sei                     ; no interrupts
                ldy     $d3             ; get current screen column
                lda     ($d1),y         ; and clear 
                and     #$7f            ;   cursor in
                sta     ($d1),y         ;   current row
                inc     $cc             ; disable cursor blinking
                cli                     ; enable interrupts
                rts                     ; return
backspace:      dec     $fe             ; decrement buffer index
                bcs     getkey          ; and get next key
.endproc

; parse / convert uint8 number using a BCD representation and double-dabble
.proc touint8
                sta     $fb
                sty     $fc
                ldy     #$0
                sty     convbuf
                sty     convbuf+1
                sty     convbuf+2
scanloop:       lda     ($fb),y
                beq     copy
                iny
                cmp     #$20
                beq     scanloop
                cmp     #$30
                bcc     error
                cmp     #$3a
                bcs     error
                bcc     scanloop
error:          sec
                rts
copy:           dey
                bmi     error
                ldx     #$2
copyloop:       lda     ($fb),y
                cmp     #$30
                bcc     copynext
                cmp     #$3a
                bcs     copynext
                sec
                sbc     #$30
                sta     convbuf,x
                dex
copynext:       dey
                bpl     copyloop
                lda     #$0
                sta     $fb
                ldx     #$8
loop:           lsr     convbuf
                lda     convbuf+1
                bcc     skipbit1
                ora     #$10
skipbit1:       lsr     a
                sta     convbuf+1
                lda     convbuf+2
                bcc     skipbit2
                ora     #$10
skipbit2:       lsr     a
                sta     convbuf+2
                ror     $fb
                dex
                beq     done
                lda     convbuf
                cmp     #$8
                bmi     nosub1
                sbc     #$3
                sta     convbuf
nosub1:         lda     convbuf+1
                cmp     #$8
                bmi     nosub2
                sbc     #$3
                sta     convbuf+1
nosub2:         lda     convbuf+2
                cmp     #$8
                bmi     loop
                sbc     #$3
                sta     convbuf+2
                bcs     loop
done:           lda     $fb
                clc
                rts
.endproc
Felix Palmen
la source
4
J'apprécie tellement plus que le flux constant de langues de golf (je peux ou non porter un t-shirt MOS 6502 aujourd'hui).
Matt Lacey
1
@MattLacey merci :) Je suis juste trop paresseux pour apprendre toutes ces langues ... et faire des puzzles dans le code 6502 semble un peu "naturel" car économiser de l'espace est en fait une pratique de programmation standard sur cette puce :)
Felix Palmen
J'ai besoin d'acheter un T-shirt MOS 6502.
Titus
8

Python 2 , 49 octets

f=lambda n,t=1,p=1:n and p%t*t+f(n-p%t,t+1,p*t*t)

Utilise le théorème de Wilson , (tel qu'introduit sur le site par xnor, je crois ici )

Essayez-le en ligne!

La fonction fest récursive, avec une entrée initiale de net une queue lorsqu'elle natteint zéro, ce qui donne ce zéro (en raison de la logique and); nest décrémenté chaque fois tqu'un nombre de test qui s'incrémente à chaque appel à fest premier. Le premier test est alors de savoir si pour laquelle nous gardons une trace d'un carré de la factorielle en.(n1)!  1(modn)p

Jonathan Allan
la source
J'adaptais l' une des fonctions d'aide courantes de Lynn et j'ai réalisé exactement la même chose.
M. Xcoder
... ah donc le théorème a été introduit sur le site par xnor. Bon poste de référence, merci!
Jonathan Allan
6

05AB1E , 3 octets

ÅpO

Essayez-le en ligne.

Explication:

Åp     # List of the first N primes (N being the implicit input)
       #  i.e. 5 → [2,3,5,7,11]
  O    # Sum of that list
       #  i.e. [2,3,5,7,11] → 28
Kevin Cruijssen
la source
6

Java 8, 89 octets

n->{int r=0,i=2,t,x;for(;n>0;r+=t>1?t+0*n--:0)for(t=i++,x=2;x<t;t=t%x++<1?0:t);return r;}

Essayez-le en ligne.

Explication:

n->{               // Method with integer as both parameter and return-type
  int r=0,         //  Result-sum, starting at 0
      i=2,         //  Prime-integer, starting at 2
      t,x;         //  Temp integers
  for(;n>0         //  Loop as long as `n` is still larger than 0
      ;            //    After every iteration:
       r+=t>1?     //     If `t` is larger than 1 (which means `t` is a prime):
           t       //      Increase `r` by `t`
           +0*n--  //      And decrease `n` by 1
          :        //     Else:
           0)      //      Both `r` and `n` remain the same
    for(t=i++,     //   Set `t` to the current `i`, and increase `i` by 1 afterwards
        x=2;       //   Set `x` to 2
        x<t;       //   Loop as long as `x` is still smaller than `t`
      t=t%x++<1?   //    If `t` is divisible by `x`:
         0         //     Set `t` to 0
        :          //    Else:
         t);       //     `t` remains the same
                   //   If `t` is still the same after this loop, it means it's a prime
  return r;}       //  Return the result-sum
Kevin Cruijssen
la source
5

Perl 6 , 31 octets

{sum grep(&is-prime,2..*)[^$_]}

Essayez-le en ligne!

L' is-primeintégré est malheureusement long.

Jo King
la source
5

Brachylog , 8 7 octets

~lṗᵐ≠≜+

Essayez-le en ligne!

1 octet enregistré grâce à @sundar.

Explication

~l        Create a list of length input
  ṗᵐ      Each element of the list must be prime
    ≠     All elements must be distinct
     ≜    Find values that match those constraints
      +   Sum
Fatalize
la source
~lṗᵐ≠≜+semble fonctionner, pour 7 octets (Aussi, je suis curieux de savoir pourquoi il donne la sortie 2 * entrée + 1 si exécuté sans l'étiquetage.)
sundar - Réinstallez Monica
2
@sundar J'ai vérifié en utilisant le débogueur et j'ai trouvé pourquoi: il ne choisit pas de valeurs pour les nombres premiers, mais il sait toujours que chacun doit être [2,+inf)évidemment. Par conséquent, il sait que la somme de 5 nombres premiers (si l'entrée est 5) doit être au moins 10, et il sait en partie que parce que les éléments doivent être différents, ils ne peuvent pas tous être 2 donc c'est au moins 11. L'implémentation TL; DR de la labellisation implicite n'est pas assez forte.
Fatalize
C'est très intéressant. J'aime que la raison ne soit pas une bizarrerie de syntaxe ou un accident aléatoire d'implémentation, mais quelque chose qui a du sens en fonction des contraintes. Merci d'avoir jeté un coup d'œil!
sundar
2

Rétine , 41 octets

K`_
"$+"{`$
$%"_
)/¶(__+)\1+$/+`$
_
^_

_

Essayez-le en ligne! Je voulais continuer à ajouter 1 jusqu'à ce que j'aie trouvé des nnombres premiers, mais je ne pouvais pas savoir comment faire cela dans Retina, j'ai donc eu recours à une boucle imbriquée. Explication:

K`_

Commencez par 1.

"$+"{`

nTemps de boucle .

$
$%"_

Faites une copie de la valeur précédente et incrémentez-la.

)/¶(__+)\1+$/+`$
_

Continuez à l'incrémenter pendant qu'il est composite. (Le )ferme la boucle extérieure.)

^_

Supprimez l'original 1.

_

Additionnez et convertissez en décimal.

Neil
la source
2

MATL , 4 octets

:Yqs

Essayez-le en ligne!

Explication:

       % Implicit input: 5
:      % Range: [1, 2, 3, 4, 5]
 Yq    % The n'th primes: [2, 3, 5, 7, 11]
   s   % Sum: 28
Stewie Griffin
la source
2

PHP, 66 octets

en utilisant à nouveau ma propre fonction principale ...

for(;$k<$argn;$i-1||$s+=$n+!++$k)for($i=++$n;--$i&&$n%$i;);echo$s;

Exécuter en tant que pipe avec -nrou l' essayer en ligne .

panne

for(;$k<$argn;      # while counter < argument
    $i-1||              # 3. if divisor is 1 (e.g. $n is prime)
        $s+=$n              # add $n to sum
        +!++$k              # and increment counter
    )
    for($i=++$n;        # 1. increment $n
        --$i&&$n%$i;);  # 2. find largest divisor of $n smaller than $n:
echo$s;             # print sum
Titus
la source
même longueur, une variable en moins:for(;$argn;$i-1||$s+=$n+!$argn--)for($i=++$n;--$i&&$n%$i;);echo$s;
Titus
2

Haskell , 48 octets

sum.(`take`[p|p<-[2..],all((>0).mod p)[2..p-1]])

Essayez-le en ligne!

\p-> all((>0).mod p)[2..p-1]True0,12

ბიმო
la source
2

C (gcc) , 70 octets

f(n,i,j,s){s=0;for(i=2;n;i++)for(j=2;j/i?s+=i,n--,0:i%j++;);return s;}

Essayez-le en ligne!

Curtis Bechtel
la source
Suggérer à la n=splace dereturn s
plafondcat
2

C, C ++, D: 147 142 octets

int p(int a){if(a<4)return 1;for(int i=2;i<a;++i)if(!(a%i))return 0;return 1;}int f(int n){int c=0,v=1;while(n)if(p(++v)){c+=v;--n;}return c;}

Optimisation 5 octets pour C et C ++:

-2 octets grâce à Zacharý

#define R return
int p(int a){if(a<4)R 1;for(int i=2;i<a;++i)if(!(a%i))R 0;R 1;}int f(int n){int c=0,v=1;while(n)if(p(++v))c+=v,--n;R c;}

pteste si un nombre est un nombre premier, frésume lan premiers nombres

Code utilisé pour tester:

C / C ++:

for (int i = 0; i < 10; ++i)
    printf("%d => %d\n", i, f(i));

D Réponse optimisée de Zacharý , 133 131 octets

D a un système de gabarits golfy

T p(T)(T a){if(a<4)return 1;for(T i=2;i<a;)if(!(a%i++))return 0;return 1;}T f(T)(T n){T c,v=1;while(n)if(p(++v))c+=v,--n;return c;}
HatsuPointerKun
la source
1
T p(T)(T a){if(a<4)return 1;for(T i=2;i<a;)if(!(a%i++))return 0;return 1;}T f(T)(T n){T c,v=1;while(n)if(p(++v)){c+=v;--n;}return c;}. En outre, le C / C ++ / D peut être int p(int a){if(a<4)return 1;for(int i=2;i<a;++i)if(!(a%i))return 0;return 1;}int f(int n){int c=0,v=1;while(n)if(p(++v)){c+=v;--n;}return c;}(identique à l'optimisation C / C ++, en ajustant simplement l'algorithme)
Zacharý
Peut-être pour toutes les réponses, vous pourriez utiliser la virgule pour faire {c+=v;--n;}être c+=v,--n;?
Zacharý
En voici une autre pour D (et peut-être aussi pour C / C ++, si elle revient à ints):T p(T)(T a){T r=1,i=2;for(;i<a;)r=a%i++?r:0;return r;}T f(T)(T n){T c,v=1;while(n)if(p(++v))c+=v,--n;return c;}
Zacharý
Suggérer à la a>3&i<aplace de i<aet supprimerif(a<4)...
plafondcat
2

Japt -x , 11 octets

;@_j}a°X}hA

Essayez-le en ligne!

Plusieurs octets enregistrés grâce à une nouvelle fonction de langue.

Explication:

;@      }hA    :Get the first n numbers in the sequence:
     a         : Get the smallest number
      °X       : Which is greater than the previous result
  _j}          : And is prime
               :Implicitly output the sum
Kamil Drakari
la source
1

JavaScript (ES6), 55 octets

f=(k,n=2)=>k&&(g=d=>n%--d?g(d):d<2&&k--&&n)(n)+f(k,n+1)

Essayez-le en ligne!

Arnauld
la source
1

APL (Dyalog Unicode) , 7 + 9 = 16 octets

+/pco∘⍳

Essayez-le en ligne!

9 octets supplémentaires pour importer le pcoDfn (et tous les autres):⎕CY'dfns'

Comment:

+/pco∘⍳
        Generate the range from 1 to the argument
        Compose
  pco    P-colon (p:); without a left argument, it generates the first <right_arg> primes.
+/       Sum
J. Sallé
la source
N'avez-vous pas besoin d'ajouter encore un autre octet? import X(nouvelle ligne) X.something()en python est compté avec la nouvelle ligne.
Zacharý
1

Rubis, 22 + 7 = 29 octets

Courir avec ruby -rprime(+7)

->n{Prime.take(n).sum}
Piccolo
la source
1

JAEL , 5 octets

#&kȦ

Explication (générée automatiquement):

./jael --explain '#&kȦ'
ORIGINAL CODE:  #&kȦ

EXPANDING EXPLANATION:
Ȧ => .a!

EXPANDED CODE:  #&k.a!,

#     ,                 repeat (p1) times:
 &                              push number of iterations of this loop
  k                             push nth prime
   .                            push the value under the tape head
    a                           push p1 + p2
     !                          write p1 to the tape head
       ␄                print machine state
Eduardo Hoefel
la source
0

Python 2 , 63 59 56 51 octets

f=lambda n:n and prime(n)+f(n-1)
from sympy import*

Essayez-le en ligne!


Enregistré:

  • -5 octets, merci à Jonathan Allan

Sans libs:

Python 2 , 83 octets

n,s=input(),0
x=2
while n:
 if all(x%i for i in range(2,x)):n-=1;s+=x
 x+=1
print s

Essayez-le en ligne!

TFeld
la source
f=lambda n:n and prime(n)+f(n-1)sauve cinq (il pourrait être golfable plus loin aussi)
Jonathan Allan
0

CJam , 21 octets

0{{T1+:Tmp!}gT}li*]:+


Explanation:
0{{T1+:Tmp!}gT}li*]:+ Original code

 {            }li*    Repeat n times
  {        }          Block
   T                  Get variable T | stack: T
    1+                Add one | Stack: T+1 
      :T              Store in variable T | Stack: T+1
        mp!           Is T not prime?     | Stack: !(isprime(T))
            g         Do while condition at top of stack is true, pop condition
             T        Push T onto the stack | Stack: Primes so far
0                 ]   Make everything on stack into an array, starting with 0 (to not throw error if n = 0)| Stack: array with 0 and all primes up to n
                   :+ Add everything in array

Essayez-le en ligne!

lolade
la source
0

F #, 111 octets

let f n=Seq.initInfinite id|>Seq.filter(fun p->p>1&&Seq.exists(fun x->p%x=0){2..p-1}|>not)|>Seq.take n|>Seq.sum

Essayez-le en ligne!

Seq.initInfinitecrée une séquence infiniment longue avec une fonction de générateur qui prend, comme paramètre, l'index de l'article. Dans ce cas, la fonction de générateur est juste la fonction d'identitéid .

Seq.filter sélectionne uniquement les nombres créés par la séquence infinie qui sont premiers.

Seq.take prend le premier n éléments de cette séquence.

Et enfin, les Seq.sumrésume.

Ciaran_McCarthy
la source
0

cQuents , 3 octets

;pz

Essayez-le en ligne!

Explication

;    sum of first n terms for input n
 pz  each term is the next prime after the previous term
Stephen
la source
Notez que la version actuelle utilise Zau lieu dez
Stephen
0

MY , 4 octets

⎕ṀΣ↵

Essayez-le en ligne!

Regrettant toujours aucune entrée / sortie implicite dans ce langage poubelle, cela aurait été deux octets sinon.

  • = entrée
  • = 1er ... nième prime inclus
  • Σ = somme
  • = sortie
Zacharý
la source
0

APL (NARS), 27 caractères, 54 octets

{⍵=0:0⋄+/{⍵=1:2⋄¯2π⍵-1}¨⍳⍵}

{¯2π⍵} ici renverrait le n premier différent de 2. Donc {⍵ = 1: 2⋄¯2π⍵-1} renverrait le n premier 2 en nombre ...

RosLuP
la source