6502.org

Code: Select all

PROC INMOS.ENTRY.POINT (CHAN OF ANY from.link, to.link)
  -- DO NOT CHANGE PROCEDURE NAME, IT IS IMPORTANT.

  #USE "\toolset\flibs"      -- File server library
  #USE "\toolset\realio"     -- IO-conversions

  PROC is.it.prim(VAL INT testnum, BOOL result)
    VAL step IS 2:
    VAL init IS 3:
    INT i:
    SEQ
      i := init
      WHILE (( (i TIMES i) < testnum) AND (( testnum REM i) <> 0))
        i := i PLUS step
      result := (testnum < (i TIMES i))
  :

  VAL Close.Option IS 0:
  VAL f.ok IS 0:

  INT result, screen.id, key.id:
  SEQ
    open.output.stream(from.link, to.link, 0, screen.id, result)
    open.input.stream(from.link, to.link, 0, key.id, result)
    IF
      result = f.ok
        [20]BYTE count.str:

        INT len, loprim, hiprim, diff, lauf, time1, time2:
        TIMER clock:
        BOOL error,prim:
        VAL step IS 2:
        VAL limit IS ((MOSTPOS INT) >> 1):
        SEQ
          SEQ a = 0 FOR 20
            count.str[a] := ' '   -- init count.string

          write.block(from.link, to.link, screen.id,
                      "minimale Primzahldifferenz ?  :", len, result)

          result := f.ok - 1
          WHILE result <> f.ok
            read.block(from.link, to.link, key.id, 20, count.str, len, result)
          STRINGTOINT(error, diff, [count.str FROM 0 FOR len])

          clock ? time1
          loprim := 3
          hiprim := 7

          WHILE (diff < 300)
            SEQ           
              WHILE ((hiprim MINUS loprim) < diff)
                SEQ
                  lauf := hiprim MINUS loprim
                  WHILE (lauf > 0)
                    SEQ
                      is.it.prim(loprim PLUS lauf,prim)
                      IF
                        prim
                          SEQ
                            loprim := lauf PLUS loprim
                            hiprim := loprim
                            lauf := 0
                        TRUE
                          lauf := lauf MINUS step
                  IF
                    hiprim = loprim
                      SEQ
                        hiprim := (loprim PLUS diff) MINUS step
                    TRUE
                      SEQ
                        lauf := step
                        hiprim := hiprim PLUS step
                        WHILE (lauf > 0)
                          SEQ
                            is.it.prim(hiprim,prim)
                            IF
                              prim
                                SEQ
                                  lauf := 0
                              TRUE
                                hiprim := hiprim PLUS step



              clock ? time2

              SEQ a = 0 FOR 20
                count.str[a] := ' '   -- init count.string
              INTTOSTRING(len, count.str, hiprim)
              write.block(from.link, to.link, screen.id,
                          count.str, len, result)
              SEQ a = 0 FOR 20
                count.str[a] := ' '   -- init count.string
              INTTOSTRING(len, count.str, loprim)
              write.block(from.link, to.link, screen.id,
                          count.str, len, result)
              SEQ a = 0 FOR 20
                count.str[a] := ' '   -- init count.string
              INTTOSTRING(len, count.str, (hiprim-loprim))
              write.block(from.link, to.link, screen.id,
                          [count.str FROM 0 FOR len], len, result)
              time2 := time2 - time1
              SEQ a = 0 FOR 20
                count.str[a] := ' '   -- init count.string
              REAL32TOSTRING(len, count.str,
                            (( REAL32 ROUND time2 ) / 15625.0 (REAL32)), 5, 5)
    --          INTTOSTRING(len, count.str, time2)
              write.block(from.link, to.link, screen.id,
                          [count.str FROM 0 FOR len], len, result)
              write.block(from.link, to.link, screen.id,
                          " Sek. *N", len, result)

              diff := (hiprim-loprim)+2
    
          terminate.filer(from.link, to.link, result)
:

Code: Select all

include "inc/cmdsys.plh"

var a, b, i,loprim,hiprim

def getnum
    var str, num
    num = 0
    str = gets(':'|$80) + 1
    while ^str >= '0' and ^str <= '9'
        num = num * 10 + (^str - '0')
        str++
    loop
    return num
end

def prim(x)
    var i
    i = 3
    while i*i < x and x % i; i = i + 2; loop
    return x < i*i
end

puts("Geben Sie den oberen Grenzwert ein")
a = getnum
puts("Geben Sie die minimale Differenz ein")
b = getnum
loprim, hiprim, i = 1, 1, 3
while i <= a and hiprim-loprim < b
    if prim(i)
        loprim = hiprim
        hiprim = i
    fin
    i = i + 2
loop
if hiprim-loprim < b
    puts("keine Loesung gefunden !")
else
    puti(hiprim); putc(' '); puti(loprim); putc(' '); puti(hiprim-loprim)
fin
putln
done

Code: Select all

A: 5.72s
B: 9.30s
C: 118.36s

Code: Select all

(truncate 5 2) 
2
1

Code: Select all

(multiple-value-bind (div mod)
    (truncate 5 2)
  (format t "Div: ~a, Mod: ~a~%" div mod))

Code: Select all

A <1ms
B <1ms
C 2ms
D 4ms
E 7ms
F 183ms
G (gap 222) - 122164747, 122164969, 222, 627.332 seconds

Code: Select all

#include <stdio.h>
#include <stdint.h>
#include <sys/time.h>

void findGap(const uint8_t target)
{
    struct timeval begin,end;
    uint32_t c=3, z=3, l=2, ll=4;

    gettimeofday(&begin, NULL);

    for(;;) {
        if(c > ll) {
            ll += l+l+1;
            l++;
        }
        for(uint32_t d=3; d < l; d += 2)
            if(!(c%d))
                goto notPrime;
        if(c-z >= target)
            break;
        z = c;

    notPrime:
        c += 2;
    }

    gettimeofday(&end, NULL);
    printf("%u\t%u\t%u\t%.3f\n", z, c, c-z, (end.tv_sec - begin.tv_sec) + 0.000001 * (end.tv_usec - begin.tv_usec));
}

int main(void)
{
    findGap(20);
    findGap(30);
    findGap(35);
    findGap(50);
    findGap(70);
    findGap(100);
    findGap(222);

    return 0;
}

Code: Select all

        A       B       C       D       E       F
int     2.48    4.17    53.61   139.26  262.22  7653.10
long*   13.24   22.61   312.63  826.16  1570.51 46712.30