2021-05-01から1ヶ月間の記事一覧

オリジナル ディオファントス方程式の解法

不定方程式 a^3+b^5=c^7を解いてください。解答*は乗算,^はべき乗を表します。 まずa=p^5, b=p^3, c=2*p^2 とおきます。するとa^3+b^5=p^15+p^15=(2^7)*(p^14)=c^7となります。まとめるとp^14*(p+p-2^7)=0となります。よって2*p=2^7 ゆえにp=2^6となります。…

エラトステネスの篩

DIM a(1000)FOR b=2 TO 1000 LET a(b)=bNEXT BREM------------------------------------FOR i=2 TO 1000 FOR j=2 TO 1000 IF i<>a(j) AND MOD(a(j),i)=0 THEN LET a(j)=0 NEXT JNEXT IREM ------------------------------LET t=1FOR k=2 TO 1000 IF a(k)<>0 …

銅バッハ数(copperbach number)

10 REM 銅バッハ数とはちょうど2通りの異なる方法で2つの素数の和として表わされる整数です。例えば14=7+7=11+3であるから14は銅バッハ数です。20 REM 銅バッハ数はOEISによると全部で9個のようです。30 REM !の記号以下は注釈化です。40 REM 10行目から170…

キュバン素数

number-pedia.com キュバン素数を求めるプログラム。十進BASIC使用。最初の1000個を求めています。 LET s=0LET t=0FOR a=1 TO 10000 LET b=3*a^2+3*a+1 FOR c=2 TO INT(SQR(b)) IF MOD(b,c)=0 THEN EXIT FOR IF MOD(b,c)<>0 THEN LET s=s+1 NEXT C IF s=INT(…

エジプト式分数

REM 有理数モードで実行 REM エジプト式分数を求めるプログラムです。十進BASIC使用。INPUT n,mLET b=n/m LET t=0DIM a(200)FOR i=2 TO 100000 LET c=b-1/i IF c<=0 THEN GOTO 300 LET t=t+1 LET a(t)=1/i IF numer(c)=1 THEN LET t=t+1 IF numer(c)=1 THEN …

ズッカーマン数

ja.wikipedia.org 十進basicを使っています ズッカーマン数をもとめています。 オプションタブより microsoft 互換モードで実行 LET s=1FOR a=1 TO 1000000 LET aa$=STR$(a) LET b=LEN(aa$) LET product=1 FOR r=1 TO b LET product=product*VAL(mid$(aa$,r,…