Piの数字を計算する


15

これはやや異なるタスクです。1024桁目の16進数から始まる1024桁のπを計算します。

正式:プログラムは1分以内に完了し、次の出力が生成されます。

25d479d8f6e8def7e3fe501ab6794c3b976ce0bd04c006bac1a94fb6409f60c45e5c9ec2196a246368fb6faf3e6c53b51339b2eb3b52ec6f6dfc511f9b30952ccc814544af5ebd09bee3d004de334afd660f2807192e4bb3c0cba85745c8740fd20b5f39b9d3fbdb5579c0bd1a60320ad6a100c6402c7279679f25fefb1fa3cc8ea5e9f8db3222f83c7516dffd616b152f501ec8ad0552ab323db5fafd23876053317b483e00df829e5c57bbca6f8ca01a87562edf1769dbd542a8f6287effc3ac6732c68c4f5573695b27b0bbca58c8e1ffa35db8f011a010fa3d98fd2183b84afcb56c2dd1d35b9a53e479b6f84565d28e49bc4bfb9790e1ddf2daa4cb7e3362fb1341cee4c6e8ef20cada36774c01d07e9efe2bf11fb495dbda4dae909198eaad8e716b93d5a0d08ed1d0afc725e08e3c5b2f8e7594b78ff6e2fbf2122b648888b812900df01c4fad5ea0688fc31cd1cff191b3a8c1ad2f2f2218be0e1777ea752dfe8b021fa1e5a0cc0fb56f74e818acf3d6ce89e299b4a84fe0fd13e0b77cc43b81d2ada8d9165fa2668095770593cc7314211a1477e6ad206577b5fa86c75442f5fb9d35cfebcdaf0c7b3e89a0d6411bd3ae1e7e4900250e2d2071b35e226800bb57b8e0af2464369bf009b91e5563911d59dfa6aa78c14389d95a537f207d5ba202e5b9c5832603766295cfa911c81a4a94a14a94a14a94a94a4a73a14a94a14a94a14a4a94a14a94a14a7a4a7a94a14a4a7a7a9a4a7a9a4a7a4a7a4a7a4a7a4a7a7a4a7a4a7a4a7a4a7a4a7a4a7a4a7a4a7a4a7a4a3a4a7a4a3a4a3a4a7a7a7a8a8a8a7a7a8a7a8a7b8a9a8a7b8a9b8a9b8e7e4e2501e7e4900250e2d2071b35e226e

最短の長さのプログラムが勝ちます。実行時にすべての数字を計算する必要があります。πを計算するアルゴリズムを実装する必要はありません。言語がすでにその機能を提供している場合は、それを使用できます。


ああ、簡単にできるはずです。(+1、まだ)数秒以内に実行できるに違いない。
Mateen Ulhaq

@muntoo:そして?ソリューションはどこにありますか?
FUZxxl

私はそれをするのを忘れました。:)ところで、速度!=コードゴルフ。
Mateen Ulhaq

@muntoo:わかった。しかし、私はまた、5日間はそのような簡単な作業に適した時期だと思います。
-FUZxxl

回答:


13

セージ、29文字

数字は実行時に計算されるため、これは技術的に不正行為ではありません。とは言っても、それはまだ地獄のように安価です。

hex(floor(pi*2^8192))[1025:]

1
間違いなくごまかしません。
FUZxxl

11
うーん、フロア円。
ブレッドボックス

13

シェルユーティリティ:48

curl -sL ow.ly/5u3hc|grep -Eom 1 '[a-f0-9]{1024}'

  • すべての出力は、実行時に「計算」されます。(ソリューションを投稿するOPに感謝)
  • 1分以内に実行されます。(インターネット接続速度に依存する場合があります)

通常、私はその種の解決策に賛成票を投じます。なぜなら、それらはルールの一般的な乱用であり、もはや面白くないからです。しかし、提供された参照ソリューションを取得して書き込むのが非常に卑劣であるという理由だけで、すべての出力は実行時に「計算」されます。(ソリューションを投稿するOPに感謝)、私はあなたに
賛成票

ゴルフバージョン:curl -sL ow.ly/shKGY|grep -Po \\w{99,}(37)。ダッシュで動作します。Bashには追加のバイトが必要です。
デニス14年

6

J、156、140、137 127

d=:3 :'1|+/4 _2 _1 _1*+/(y&(16^-)%1 4 5 6+8*])"0 i.y+9'
,1([:}.'0123456789abcdef'{~[:|.[:<.[:(],~16*1|{.)^:8 d)"0\1024x+8*i.128

BBP式を使用します。

する NOT分の下で実行されます(:Pが、我々はJの答えを持っています)

πの最初の104桁の例(これは高速に実行されます):

,1([:}.'0123456789abcdef'{~[:|.[:<.[:(],~16*1|{.)^:8 d)"0\8*i.13x

243f6a8885a308d313198a2e03707344a4093822299f31d0082efa98ec4e6c89
452821e638d01377be5466cf34e90c6cc0ac29b7

数字を16進数に変換するために#:を使用しないのはなぜですか?
FUZxxl

どういう意味かわかりません。#:16進数を出力しません。
エルベックス

私見では、#:と、16進数を生成するための変形を使用する方が、現在のアプローチよりも簡単です。
FUZxxl

次のようなもの(... 16 #:) Piですか?桁数が足りないので、とにかく生成する必要があると思います。
エルベックス

1
ところで、hfd数値を16進数に変換する動詞があることがわかりました。
FUZxxl

5

JavaScript、536

(読みやすくするための改行とインデントのみ)

var d='0123456789abcdef',p='',o='',l=3e3,c=0,e='length';d=d+d;
function $(n,r){return n[e]<=r?0:d.indexOf(n[r])}
function g(a,b){for(i=0,t='',s=16;i<l;i++,t+=d[~~(s/b)],s=(s%b)*16);
for(;a--;t=_(t,t,1));return t}
function _(a,b,s){for(i=(a[e]>b[e]?a[e]:b[e])-1,r='',c=0;i>=0;r=(s?
  function(k){c=k>15;return d[k]}($(a,i)+$(b,i)+c):
  function(k){c=k<0;return d[k+16]}($(a,i)-$(b,i)-c))+r,i--);return r}
for(i=0;i<l;i++,p+='2');
for(j=1;j<l;p=_(p,(o+='0')+_(_(_(g(2,8*j+1),g(1,8*j+4)),g(0,8*j+5)),g(0,8*j+6)),1),j++);
console.log(p.slice(1024,2048))

Intel i5コアを使用するラップトップのGoogle Chrome 14では、約25秒かかります。他の誰かがこのコードをゴルフできますか?私はうまくゴルフできません.. :(

以下はゴルフではありません。コメントをすべて削除し、ループをゴルフに変更しました。

については言及しないでくださいfor(;s>=b;s-=b);s*=16;。に変更しましたs=(s%b)*16。:P

/**
Calculate PI-3 to 3000 (3e3) digits.
a : a
b : b
c : carry
d : digits
e : length
f : get from d
g : calculate (2^a)/b.
i,j, : for looping
l : length to calculate
p : pi
r,t : return value
*/
var d='0123456789abcdef',p='',o='',l=3e3,c=0,e='length';
d=d+d;//for carring

function $(n,r){return n[e]<=r?0:d.indexOf(n[r])}
/*
Calculate (2^a)/b. Assume that 2^a < b.
*/
function g(a,b){
    for(i=0,t='',s=16;i<l;i++){t+=d[~~(s/b)];for(;s>=b;s-=b);s*=16;}
    for(;a--;t=_(t,t,1));return t}
/*
Calculate a±b. (+ when s=1, - when s=0) When calculating minus, assume that 1>b>a>0.
*/
function _(a,b,s){
    for(i=(a[e]>b[e]?a[e]:b[e])-1,r='',c=0;i>=0;
        r=(s?function(k){c=k>15;return d[k]}($(a,i)+$(b,i)+c):
            function(k){c=k<0;return d[k+16]}($(a,i)-$(b,i)-c))+r,i--);return r;
}
/*
Using BBP formula. Calc when j=0...
4/1 - 2/4 - 1/5 - 1/6 = 3.22222222.... (b16)
*/
for(i=0;i<l;i++,p+='2');
//Calc when j>0
for(j=1;j<l;p=_(p,(o+='0')+_(_(_(g(2,8*j+1),g(1,8*j+4)),g(0,8*j+5)),g(0,8*j+6)),1),j++);
console.log(p.slice(1024,2048));

編集:完全に未使用の機能を削除しました。(なぜそれを維持したのですか?:/)

PS。PIの最初の100桁

243f6a8885a308d313198a2e03707344a4093822299f31d0082efa98ec4e6c89452821e638d01377be5466cf34e90c6cc0ab


@FUZxxl:ゴルフをしていないコードでは不十分でしたか?.. :(
JiminP

確かにそうでした。しかし、私見では、バックティック構文を使用する代わりにコードの書式設定を使用すると、見た目が良くなります。私が書いたように、もしあなたが嫌い​​なら元に戻してください。
FUZxxl

d='0123456789abcdef',l=3e3,p=Array(l+1).join(2),o='',c=0,e='length';d+=d;function _(a,b,s){for(i=(a[e]>b[e]?a[e]:b[e])-1,r='',c=0;i+1;r=d[Z=F(b,i,1)+c,k=F(a,i,1)+(s?Z:16-Z),c=s?k>15:k<16,k]+r,i--);return r}function F(a,b,f){if(f)f=a[e]>b?d.indexOf(a[b]):0;else{for(i=0,f='',s=16;i++<l;f+=d[~~(s/b)],s=(s%b)*16);while(a--)f=_(f,f,1)}return f}for(j=0;++j<l;p=_(p,(o+='0')+_(_(_(F(2,z=8*j+1),F(1,z+3)),F(0,z+4)),F(0,z+5)),1));console.log(p.slice(1024,2048))
ピーターテイラー

いくつかは非常に大きく見えますが、それは本当にすべてのマイクロ最適化です。最大の節約は_,演算子を支持する途中で2つの匿名関数を削除することです。最も厄介なのは、1つの関数を1つの関数にマージすることで$ありg、それらの間で選択するオプションの引数があります。functionそしてreturn両方とも非常に高価であるため、if(f)...elseとのカップルは,1、合理的なトレードオフです。
ピーターテイラー

4

PHP 116 114バイト

<?for(;$g?$d=0|($$g=$g--/2*$d+($$g?:2)%$g*$f)/$g--:4613^printf($i++>257?'%04x':'',$e+$d/$f=4*$g=16384)^$e=$d%$f;);

このソリューションでは、最大2048桁の16進数字、一度に4桁の16進数字すべてを計算し、それらの後半を出力します。実行時間は5秒未満です。計算に使用される式は次のとおりです。

pi = 2 + 1/3*(2 + 2/5*(2 + 3/7*(2 + 4/9*(2 + 5/11*(2 + 6/13*(2 + 7/15*(2 + ... )))))))

精度は、剰余を配列に格納し、2 ^ 14の各除算を増分的に継続して取得されます。

Python 64バイト

x=p=16385
while~-p:x=p/2*x/p+2*2**8192;p-=2
print('%x'%x)[1025:]

上記と同じ方法。約0.2秒で実行されます。

または、ワンライナーとして 73バイトの

print('%x'%reduce(lambda x,p:p/2*x/p+2*2**8192,range(16387,1,-2)))[1025:]

3

PARI / GP-2.4、141

forstep(d=1024,2047,8,y=frac(apply(x->sum(k=0,d+30,16^(d-k)/(8*k+x)),[1,4,5,6])*[4,-2,-1,-1]~);for(i=0,7,y=16*frac(y);printf("%X",floor(y))))

Bailey–Borwein–Plouffe式を使用します(もちろん)。

1分未満で実行されます。


3

Cコード:

long ki,k,e,d=1024;
int  dig,tD=0,c,Co=0,j,js[4]  ={1,4,5,6};
double res=0.0,tres=0.0,gT,ans[4] ={0.0};
while(tD < 1024)
{while(Co<4){ j= js[Co],gT=0.0,ki= 0;
 for(; ki < d+1;ki++){ k = 8*ki+j,e= d-ki,c=1; while(e--) c = 16*c % k; gT+=((double)(c)/(double)k);}
 ans[Co] = (gT - (int)gT),++Co;}
 double gA = 4*ans[0]-2*ans[1]-ans[2]-ans[3];
 gA = (gA<0) ? gA + -1*(int)gA +1 : gA -(int)gA;
 dig=0;while(dig++ < 6 && tD++ < 1024) gA *=16, printf("%X",gA),gA -= (int)gA;
 d+=6,Co = 0;}

ランタイム= Intelクアッドコアで8.06秒


これはコードゴルフです。短い変数名を使用し、空白を避けて、コードをできるだけ圧縮してください。たとえばprintf("%X",(int)gA)、その長いリストの代わりに使用することで、多くの文字を保存できます。
-FUZxxl

1

PARI / GP-40バイト

このバージョン\xは、結果の16進数字を表示するために使用することで「チート」します。

\p8197
x=Pi<<4^6;x-=x\1;x=(x<<4^6)\1
\xx

このバージョンは、通常の方法で16進数に変換するのに87バイトかかります。

\p8197
x=Pi<<4^6;x-=x\1;concat([Vec("0123456789abcdef")[n+1]|n<-digits((x<<4^6)\1,16)])

両方のバージョンは、ほんの一瞬で実行されます。


1

Perl-59

use ntheory"Pi";say substr int(Pi(3000)<<8192)->as_hex,1027

0.1秒未満。


0

シェル68

ツール:bc -l、tr、cut

echo "scale=2468;obase=16;4*a(1)"|bc -l|tr -d '\\\n'|cut -c1027-2051

シェル64、ツール:bc -l、tr、tail、最後の場所の丸めが異なります

echo "scale=2466;obase=16;4*a(1)"|bc -l|tr -d '\\\n'|tail -c1024

PIの計算方法に関する知識は4 * a(1)にあり、その1ではscale = 2466を使用する必要があることが反復的に調査されたため、不正行為と見なされる場合があります。

cutを使用するアイデアを提供してくれたbreadboxに感謝します。


私はそれがいかに不正行為であると考えられるかわかりません。数字は実行時に計算されます。実行すると、最後の桁(Aではなく7)で出力が異なることに注意する必要があります。ところで、いくつかの文字を保存するためにddコマンドを置き換えることができると思いますtail -c1024
ブレッドボックス

はい、私も違いを観察しました(そしてそれに半時間を費やします:))。bcにx桁のスケールを使用するように指示すると、10進モードでその桁に丸められ、その後16進変換が行われます。したがって、もう1桁を取得すると、7ではなく69が生成されます。ただし、質問で丸めスタイルまたは切り捨てが指定されていません。そして、尾のアイデアに感謝します。:)
ユーザー不明

質問は次のように指定します:指定されたものと同じように出力を終了して生成します。
FUZxxl

@FUZxxl:「する必要があります...」ではなく「する必要があります...」-仕様の一部としてではなく、16進変換がうまくいくかどうか、スキップする文字数を確認するのに役立つはずだと思いました、最後の数字1024に準拠するようにします。しかし、16文字を追加して1文字を置き換えました。
ユーザー不明

1
文が「正式に」という単語で始まることを考えると、OPはおそらくその単語のRFCの意味でそれを意味していたことに同意します。また、あなたはまだの使用置き換えることによって、あなたの新しいソリューションを向上させることができますddしてをcut -c1027-2051。(シェルには、テキストストリームを操作するための多くのツールがあります。)
ブレッドボックス
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.