π(1E+8) Prime-counting function 1秒投稿者：たろさ 投稿日：2017年 4月23日(日)18時30分31秒 |
|
初級者向けParact BASICによるマルチスレッド化プログラム まだ不慣れなので、間違えていたら教えてください。 自然数1億から素数を数えるプログラム 2 !6n-1 6n+1 Parallel activity Ver.1 DECLARE STRUCTURE abc1: 1 OF NUMERIC DECLARE MESSAGE AMACHANN OF abc1 Paract Par1 OPTION ARITHMETIC NATIVE DECLARE NUMERIC xv START Par2 LET t0=TIME LET k=100000000 LET k2=5761455 !エラトステネスの篩 LET Fu=10007 LET Fm=1230 DIM P(Fu) DIM A(Fm) !素数 MAT P=ZER MAT A=ZER LET A(1)=2 LET H1=1 FOR I=3 TO SQR(Fu) STEP 2 IF P(I)=0 THEN FOR J=I*I TO Fu STEP I LET P(J)=1 NEXT J END IF NEXT I FOR I=3 TO Fu STEP 2 IF P(I)=0 THEN LET H1=H1+1 LET A(H1)=I END IF NEXT I LET Q=6 LET k7=k !篩の計算範囲 LET k5=INT(k7/Q) DIM Au(k5) MAT Au = ZER !(6*n-1) FOR n=3 TO Fm LET Pu=A(n) IF MOD(Pu+1,Q)=0 THEN !(6*n-1) LET ru=(Pu+1)/Q FOR i=1 TO k5 IF Pu*i+ru>k5 THEN EXIT FOR LET Au(Pu*i+ru)=1 NEXT i END IF IF MOD(Pu-1,Q)=0 THEN LET ru=(Pu-1)/Q FOR i=1 TO k5 IF Pu*i-ru>k5 THEN EXIT FOR LET Au(Pu*i-ru)=1 NEXT i END IF NEXT n FOR n=1 TO k5 IF 6*n-1>k7 THEN EXIT FOR IF Au(n)=0 THEN LET cz=cz+1 END IF NEXT n receive from AMACHANN TO xv PRINT "6n-1 ";cz PRINT "pi(x)";cz+xv+2 LET TM=TIME-t0 PRINT USING"####.##":TM; PRINT "秒" END PARACT PARACT Par2 OPTION ARITHMETIC NATIVE !2進モード DECLARE NUMERIC zv !6n+1 素数 LET k=100000007 LET k2=5761456 !エラトステネスの篩 LET Fu=10007 LET Fm=1230 DIM P(Fu) DIM A(Fm) !素数 MAT P=ZER MAT A=ZER LET A(1)=2 LET H1=1 FOR I=3 TO SQR(Fu) STEP 2 IF P(I)=0 THEN FOR J=I*I TO Fu STEP I LET P(J)=1 NEXT J END IF NEXT I FOR I=3 TO Fu STEP 2 IF P(I)=0 THEN LET H1=H1+1 LET A(H1)=I END IF NEXT I LET Q=6 LET k7=k !篩の計算範囲 LET k5=INT(k7/Q) DIM Av(k5) MAT Av = ZER !(6*n+1) FOR n=3 TO Fm LET Pu=A(n) IF MOD(Pu+1,Q)=0 THEN !(6*n+1) LET ru=(Pu+1)/Q FOR i=1 TO k5 IF Pu*i-ru>k5 THEN EXIT FOR LET Av(Pu*i-ru)=1 NEXT i END IF IF MOD(Pu-1,Q)=0 THEN LET ru=(Pu-1)/Q FOR i=1 TO k5 IF Pu*i+ru>k5 THEN EXIT FOR LET Av(Pu*i+ru)=1 NEXT i END IF NEXT n FOR n=1 TO k5 IF 6*n+1>k7 THEN EXIT FOR IF Av(n)=0 THEN LET zv=zv+1 END IF NEXT n PRINT "6n+1 ";zv send TO AMACHANN from zv END PARACT --------------------------------------- !6n-1 6n+1 Parallel activity Ver.2 DECLARE STRUCTURE abc1: 1 OF NUMERIC DECLARE STRUCTURE abc2: NUMERIC(1230) DECLARE SHARED sh OF abc2 DECLARE MESSAGE AMACHANN OF abc1 Paract Par1 OPTION ARITHMETIC NATIVE DECLARE NUMERIC cz DECLARE NUMERIC zv DECLARE NUMERIC A(1230) !素数 LET t0=TIME LET k=100000000 LET k2=5761455 !エラトステネスの篩 LET Fu=10007 LET Fm=1230 DIM P(Fu) MAT P=ZER LET I=2 LET H=1 LET A(H)=I FOR I=3 TO SQR(Fu) STEP 2 IF P(I)=0 THEN FOR J=I*I TO Fu STEP I LET P(J)=1 NEXT J END IF NEXT I FOR I=3 TO Fu STEP 2 IF P(I)=0 THEN LET H=H+1 LET A(H)=I END IF NEXT I PUT TO sh FROM A start par2 START par3 RECEIVE FROM AMACHANN TO cz RECEIVE FROM AMACHANN TO zv PRINT "6n-1 ";cz PRINT "6n+1 ";zv PRINT "pi(x)";cz+zv+2 LET TM=TIME-t0 PRINT USING"###.###":TM; PRINT "秒" END PARACT paract par2 OPTION ARITHMETIC NATIVE DECLARE NUMERIC A(1230) !素数 DECLARE NUMERIC cz GET from sh TO A LET Q=6 LET k7=100000000 !篩の計算範囲 LET k5=INT(k7/Q) DIM Au(k5) MAT Au = ZER !(6*n-1) FOR n=3 TO 1229 LET Pu=A(n) IF MOD(Pu+1,Q)=0 THEN !(6*n-1) LET ru=(Pu+1)/Q FOR i=1 TO k5 IF Pu*i+ru>k5 THEN EXIT FOR LET Au(Pu*i+ru)=1 NEXT i END IF IF MOD(Pu-1,Q)=0 THEN LET ru=(Pu-1)/Q FOR i=1 TO k5 IF Pu*i-ru>k5 THEN EXIT FOR LET Au(Pu*i-ru)=1 NEXT i END IF NEXT n FOR n=1 TO k5 IF 6*n-1>k7 THEN EXIT FOR IF Au(n)=0 THEN LET cz=cz+1 END IF NEXT n send TO AMACHANN from cz END PARACT Paract Par3 DECLARE NUMERIC A(1230) !素数 DECLARE NUMERIC zv GET FROM sh TO A LET Q=6 LET k7=100000000 !篩の計算範囲 LET k5=INT(k7/Q) DIM Av(k5) MAT Av = ZER !(6*n+1) FOR n=3 TO 1229 LET Pu=A(n) IF MOD(Pu+1,Q)=0 THEN !(6*n+1) LET ru=(Pu+1)/Q FOR i=1 TO k5 IF Pu*i-ru>k5 THEN EXIT FOR LET Av(Pu*i-ru)=1 NEXT i END IF IF MOD(Pu-1,Q)=0 THEN LET ru=(Pu-1)/Q FOR i=1 TO k5 IF Pu*i+ru>k5 THEN EXIT FOR LET Av(Pu*i+ru)=1 NEXT i END IF NEXT n FOR n=1 TO k5 IF 6*n+1>k7 THEN EXIT FOR IF Av(n)=0 THEN LET zv=zv+1 END IF NEXT n send TO AMACHANN from zv END PARACT 
|
