裏 RjpWiki

Julia ときどき R, Python によるコンピュータプログラム,コンピュータ・サイエンス,統計学

Hello world! を出力するプログラム

2014年11月18日 | ブログラミング

まあ,2 つの言語である文字列を書けということだけど,むしろ同じ言語でどんなバリエーションがあるかのほうが面白いんだろうなあ。

コマンドラインでやる。一行野郎だ。

R でちょっと凝ったやり方?

Rscript -e 'cat(rawToChar(as.raw(c(73, 32, 108, 111, 118, 101, 32, 121, 111, 117, 46))), "\\n")'

I love you.

AWK はあまり,凝ったやり方が思い浮かばなかったので,ずるっこなやり方でごまかす。

gawk 'BEGIN{print "I love you."}'

I love you.

コメント
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする

ネイピア数(e)の近似式

2014年11月18日 | ブログラミング
A ~ H に 1 から 8 までの数字をあてると、式の結果が自然対数の底、ネイピア数(e=2.7182818284590…)に極めて近い値となる。

まあ,虫食い算ではある。しかし,一筋縄ではいかぬ。

n を大きくすると,e ≒ (1+1/n)^n になるということを利用したもの。
B^(10*C+D)E^(10*F+G) が(ほぼ同じ大きさで)べらぼうに大きな数になるように,決める。A=1 で,H はゴミ。

library(e1071)
d = permutations(8)
colnames(d) = LETTERS[1:8]
a = d[,2]^(d[,3]*10+d[,4])
b = d[,5]^(d[,6]*10+d[,7])
c = abs(a-b)
unique(d[c == 0, 2:7])

     B C D E F G
[1,] 2 7 6 4 3 8
[2,] 2 3 6 4 1 8
[3,] 4 1 8 2 3 6
[4,] 4 3 8 2 7 6

2^76 == 4^38 = 75557863725914323419136
2^36 == 4^18 = 68719476736

どちらでもよいけど,大きい方がよいだろう。

(1 + 2^(-76)) ^ (4^38 + 0.5) でよいだろう。

> library(Rmpfr)
> (1+mpfr(2^-76, 10000))^(4^38+0.5)
1 'mpfr' number of precision  10000   bits
[1] 2.71828182845904523536026948327243884562017542260751355373386591261370824199370367512673498355740221866300150192512577013712402194525729088107603625695791765621277477110787864123095526391463127114520685536306474105328349829970690292025415223570462246776367494240117044146359689778284658748377053631108528524369649708990088499705308087739188743458948765814067112847980775789601385668978369082970614962359622188612525084100685222682523211422142283833500680692953117870857304240655781007179275504357626794755609455702336615836088456321464184491100070840364060583518742678339843281161018070243859730035061281270650049058577898918962636823041601318514133991899068659175614049936572365313687224918238106256134904163514252685340647340456890314801271060284733277103452154217953947794212256610001859188824305819005493101420995707315149924255908483971356622254440811665861497806433425245530414577413150274475908707244099490588998007971471915652104821918366705834682993754472524280666636135840177290124580033865188449321263327246243517903169586281573699595084163491236453339564455102365405626778182944392138305434444792827055459439497524425904694000237296631033220159598665136864124246396854272963872820317637333769237193242884202568303364418046273435762454117040706608172930748692315119997997841619429686669820862061110068904028565437826781936159992926487440700174612607470999978707697207082787868267520682657954314633255397701954003355903287968237899866164531903161262850661513441425018870844607137295088827744560605904851144463845968530092277811965079409662380103667771723922279391991904178680345914824907021940796420741369938414632649949272385757279737779566836807257955131811102861602826936065959718009010543032829786047163352214962220611887790481846213267329885769701504016942449439825975274239480548636777535688113115286363107890608944546004286395664234621620727850509634845140884460424113809699482972677973864783637850054751653120718800418467312735020003534712022848421166354423545617401554314100709757290349151738731506203567718721798205167916146075161513837519243811057993341024401841139508335521579606091693535359204092916638388380674274980856483447553569646936493930811895678168317661540052904197578737606560554891406196968343986453614462329059943329291691296011434571917668237269127335807768840376797480758223346640014402206484032202724578013542472530555454711911214217977945257047510183549101102775973480564986345641104547537979739157453386815737096359570223404164680296807328667655874844462658801682855771290628496762097876608805420529910640789738394606194141138693883902905388877002919338757647512180506028142016641093968849053875292720095305141013750955950014732924418826737952161667833410486719424913741594089595708283399701258904526782568491372754701118896742792090433744556255065162659355857680625978388657011433036136201159872406077243795651120889127932226093230037111704708148334727141931962344248175205214809447454142671052905204264772129669310850553225577837867783253200036289105779201198217449032195


コメント
  • X
  • Facebookでシェアする
  • はてなブックマークに追加する
  • LINEでシェアする

PVアクセスランキング にほんブログ村

PVアクセスランキング にほんブログ村