Meydan okuma

Bir sayı listesi verildiğinde, listenin popülasyon standart sapmasını hesaplayın.

Popülasyon standart sapmasını hesaplamak için aşağıdaki denklemi kullanın:

Giriş

Giriş, herhangi bir formatta (liste, dize vb.) Bir tamsayı listesi olacaktır. Bazı örnekler:

56,54,89,87
67,54,86,67

Sayılar her zaman tamsayı olur.

Girdi STDIN veya işlev bağımsız değişkenlerine yapılacaktır.

Çıktı

Çıktı bir kayan nokta sayısı olmalıdır.

kurallar

Standart sapmayı bulmak için yerleşik işlevleri kullanabilirsiniz.

Cevabınız tam bir program veya bir fonksiyon olabilir.

Örnekler

10035, 436844, 42463, 44774 => 175656.78441352615

45,67,32,98,11,3 => 32.530327730015607

1,1,1,1,1,1 => 0.0

Kazanan

En kısa program veya işlev kazanır.

Liderler Sıralaması

var QUESTION_ID=60901,OVERRIDE_USER=30525;function answersUrl(e){return"http://api.stackexchange.com/2.2/questions/"+QUESTION_ID+"/answers?page="+e+"&pagesize=100&order=desc&sort=creation&site=codegolf&filter="+ANSWER_FILTER}function commentUrl(e,s){return"http://api.stackexchange.com/2.2/answers/"+s.join(";")+"/comments?page="+e+"&pagesize=100&order=desc&sort=creation&site=codegolf&filter="+COMMENT_FILTER}function getAnswers(){jQuery.ajax({url:answersUrl(answer_page++),method:"get",dataType:"jsonp",crossDomain:!0,success:function(e){answers.push.apply(answers,e.items),answers_hash=[],answer_ids=[],e.items.forEach(function(e){e.comments=[];var s=+e.share_link.match(/\d+/);answer_ids.push(s),answers_hash[s]=e}),e.has_more||(more_answers=!1),comment_page=1,getComments()}})}function getComments(){jQuery.ajax({url:commentUrl(comment_page++,answer_ids),method:"get",dataType:"jsonp",crossDomain:!0,success:function(e){e.items.forEach(function(e){e.owner.user_id===OVERRIDE_USER&&answers_hash[e.post_id].comments.push(e)}),e.has_more?getComments():more_answers?getAnswers():process()}})}function getAuthorName(e){return e.owner.display_name}function process(){var e=[];answers.forEach(function(s){var r=s.body;s.comments.forEach(function(e){OVERRIDE_REG.test(e.body)&&(r="<h1>"+e.body.replace(OVERRIDE_REG,"")+"</h1>")});var a=r.match(SCORE_REG);a&&e.push({user:getAuthorName(s),size:+a[2],language:a[1],link:s.share_link})}),e.sort(function(e,s){var r=e.size,a=s.size;return r-a});var s={},r=1,a=null,n=1;e.forEach(function(e){e.size!=a&&(n=r),a=e.size,++r;var t=jQuery("#answer-template").html();t=t.replace("{{PLACE}}",n+".").replace("{{NAME}}",e.user).replace("{{LANGUAGE}}",e.language).replace("{{SIZE}}",e.size).replace("{{LINK}}",e.link),t=jQuery(t),jQuery("#answers").append(t);var o=e.language;/<a/.test(o)&&(o=jQuery(o).text()),s[o]=s[o]||{lang:e.language,user:e.user,size:e.size,link:e.link}});var t=[];for(var o in s)s.hasOwnProperty(o)&&t.push(s[o]);t.sort(function(e,s){return e.lang>s.lang?1:e.lang<s.lang?-1:0});for(var c=0;c<t.length;++c){var i=jQuery("#language-template").html(),o=t[c];i=i.replace("{{LANGUAGE}}",o.lang).replace("{{NAME}}",o.user).replace("{{SIZE}}",o.size).replace("{{LINK}}",o.link),i=jQuery(i),jQuery("#languages").append(i)}}var ANSWER_FILTER="!t)IWYnsLAZle2tQ3KqrVveCRJfxcRLe",COMMENT_FILTER="!)Q2B_A2kjfAiU78X(md6BoYk",answers=[],answers_hash,answer_ids,answer_page=1,more_answers=!0,comment_page;getAnswers();var SCORE_REG=/<h\d>\s*([^\n,]*[^\s,]),.*?(\d+)(?=[^\n\d<>]*(?:<(?:s>[^\n<>]*<\/s>|[^\n<>]+>)[^\n\d<>]*)*<\/h\d>)/,OVERRIDE_REG=/^Override\s*header:\s*/i;

body{text-align:left!important}#answer-list,#language-list{padding:10px;width:290px;float:left}table thead{font-weight:700}table td{padding:5px}

<script src="https://ajax.googleapis.com/ajax/libs/jquery/2.1.1/jquery.min.js"></script> <link rel="stylesheet" type="text/css" href="//cdn.sstatic.net/codegolf/all.css?v=83c949450c8b"> <div id="answer-list"> <h2>Leaderboard</h2> <table class="answer-list"> <thead> <tr><td></td><td>Author</td><td>Language</td><td>Size</td></tr></thead> <tbody id="answers"> </tbody> </table> </div><div id="language-list"> <h2>Winners by Language</h2> <table class="language-list"> <thead> <tr><td>Language</td><td>User</td><td>Score</td></tr></thead> <tbody id="languages"> </tbody> </table> </div><table style="display: none"> <tbody id="answer-template"> <tr><td>{{PLACE}}</td><td>{{NAME}}</td><td>{{LANGUAGE}}</td><td>{{SIZE}}</td><td><a href="{{LINK}}">Link</a></td></tr></tbody> </table> <table style="display: none"> <tbody id="language-template"> <tr><td>{{LANGUAGE}}</td><td>{{NAME}}</td><td>{{SIZE}}</td><td><a href="{{LINK}}">Link</a></td></tr></tbody> </table>

Klip , 3

.sk

.sstandart sapmadır, kgirdiyi formda ayrıştırır {1,2,3}.

Mathematica, 24 22 bayt

Güzel, Mathematica'da StandardDevi...popülasyon standart sapmasını değil, örnek standart sapmayı hesaplayan yerleşik bir oh vardır.

Ama ya kullanırsak Variance... oh ... aynı anlaşma.

Ancak başka bir ilgili yerleşik var:

CentralMoment[#,2]^.5&

Yaşasın. :)

Bu aynı zamanda 22 bayt için de geçerlidir:

Mean[(#-Mean@#)^2]^.5&

Ve bu 27 için:

N@RootMeanSquare[#-Mean@#]&

Oktav, 14 bayt

g=@(a)std(a,1)

İdeone üzerinde deneyin .

kdb + , 3 bayt

dev

APL derviatlarından birinin buna yerleşik olarak sahip olması gerekiyordu .

Test sürüşü

q)dev 56, 54, 89, 87
16.53028
q)f:dev
q)f 10035, 436844, 42463, 44774
175656.8
q)f 45,67,32,98,11,3
32.53033

Dyalog APL, 24 23 21 20 19 17 bayt

*∘.5∘M×⍨-M×M←+/÷≢

Bu, aşağıdaki işleve eşdeğer olan adsız, monadik bir işlev dizisi tanımlar.

{.5*⍨M(×⍨⍵)-M⍵×(M←{(+/⍵)÷≢⍵})⍵}

TryAPL'de çevrimiçi deneyin .

Nasıl çalışır

Kod birkaç trenden oluşur.

M←+/÷≢

Bu , doğru argüman için ( tüm öğelerin toplamı) ve (uzunluk) Myürüten , daha sonra sonuçlara (bölme) uygulayarak girişin aritmetik ortalamasını döndüren monadik bir 3-treni (çatal) tanımlar .+/≢÷

M×M

Bu, Mdoğru argüman için geçerli olan, bunu ikinci kez tekrarlayan ×ve sonuçlara (ürün) μ ² döndüren başka bir çatal .

×⍨-(M×M)

Bu, daha önce açıklandığı gibi aritmetik ortalamanın karesini hesaplayan, ×⍨doğru argümana (kendisiyle ürün) ve son -olarak sonuçlara (fark) uygulayan başka bir çataldır.

Giriş (x ₁ ,…, x _N ) için bu işlev geri döner (x ₁ - μ ² ,…, x _N - μ ² ) .

*∘.5∘M

Bu oluşturulan işlev, Mdoğru argümanı için geçerlidir *∘.5. İkincisi kullanımları doğru argüman tımar haritası girişini uygulamak aiçin a*0.5(kare kökü a).

(*∘.5∘M)(×⍨-(M×M))

Son olarak, önce doğru işlevi, sonra solu sonuca uygulayan ve standart sapmayı aşağıdaki gibi hesaplayan bu monadik 2-tren (üstüne) var.

R, 41 40 39 36 30 28 bayt

kod

Sayesinde beher , Alex A. ve MickyT çok bayt için.

cat(sd(c(v=scan(),mean(v))))   

eski kodlar

v=scan();n=length(v);sd(v)/(n/(n-1))**0.5
m=scan();cat(sqrt(sum(mean((m-mean(m))^2))))
m=scan();cat(mean((m-mean(m))^2)^.5) 

Bu, popülasyon standart sapmasını sağlamalıdır.

Pyth, 20 19 17 13 bayt

@.O^R2-R.OQQ2

4 bayt golf için @FryAmTheEggman teşekkürler!

Çevrimiçi deneyin.

Nasıl çalışır

        .OQ    Compute the arithmetic mean of the input (Q).
      -R   Q   Subtract the arithmetic mean of all elements of Q.
   ^R2         Square each resulting difference.
 .O            Compute the arithmetic mean of the squared differences.
@           2  Apply square root.

CJam, 24 22 21 bayt

q~_,_@_:+d@/f-:mh\mq/

@Aditsu'ya 1 baytlık golf için teşekkürler!

CJam yorumlayıcısında çevrimiçi deneyin .

Nasıl çalışır

q~                    e# Read all input and evaluate it.
  _,                  e# Copy the array and push its length.
    _@                e# Copy the length and rotate the array on top.
      _:+d            e# Copy the array and compute its sum. Cast to Double.
          @/          e# Rotate the length on top and divide the sum by it.
            f-        e# Subtract the result (μ) from the array's elements.
              :mh     e# Reduce by hypotenuse.
                      e# a b mh -> sqrt(a^2 + b^2)
                      e# sqrt(a^2 + b^2) c mh -> sqrt(sqrt(a^2 + b^2)^2 + c^2)
                      e#                           = sqrt(a^2 + b^2 + c^2)
                      e# ⋮
                 \mq/ e# Divide the result by the square root of the length.

APL, 24 bayt

{.5*⍨+/(2*⍨⍵-+/⍵÷≢⍵)÷≢⍵}

Dennis'in Dyalog APL çözümünden biraz farklı bir yaklaşım . Bu, herhangi bir APL uygulamasıyla çalışmalıdır.

Bu, vektörü ( x - µ ) hesaplayan adsız bir monadik işlev oluşturur ² olarak 2*⍨⍵-+/⍵÷≢⍵ile bölme bu N ( ÷≢⍵), kullanılarak, bu vektörden toplamını alır +/ve daha sonra kare kökünü (alır .5*⍨).

Çevrimiçi deneyin

Julia, 26 19 bayt

x->std([x;mean(x)])

Bu, bir diziyi kabul eden ve bir kayan nokta döndüren adsız bir işlev oluşturur.

Kurtulmamış, sanırım:

function f(x::Array{Int,1})
    # Return the sample standard deviation (denominator N-1) of
    # the input with the mean of the input appended to the end.
    # This corrects the denominator to N without affecting the
    # mean.
    std([x; mean(x)])
end

TI-BASIC, 7 bayt

stdDev(augment(Ans,{mean(Ans

Buradan örnek standart sapmadan popülasyon standart sapmasını elde etmek için algoritmayı ödünç aldım .

Ben olmadan bulabildiğim en kısa çözüm augment(9 bayt:

stdDev(Ans√(1-1/dim(Ans

Haskell, 61 bayt

d n=1/sum(n>>[1])
f a=sqrt$d a*sum(map((^2).(-)(d a*sum a))a)

sum(n>>[1])Haskell'in katı tip sistemini kandırmak için belki de benim özel uzunluk fonksiyonum hariç .

Python 3.4+, 30 bayt

from statistics import*;pstdev

Yerleşik işlevi içe aktarır pstdev, ör.

>>> pstdev([56,54,89,87])
16.53027525481654

JavaScript (ES6), 73 bayt

a=>Math.sqrt(a.reduce((b,c)=>b+(d=c-eval(a.join`+`)/(l=a.length))*d,0)/l)

Jöle , rakipsiz

11 byte Bu cevap yarışmayı zorlaştırır, çünkü meydan okumayı destekleyen bir dil kullanır.

S÷L
Dz_²ÇN½

Bu, Jelly'e APL cevabımın doğrudan çevirisidir .Çevrimiçi deneyin!

Nasıl çalışır

S÷L        Helper link. Argument: z (vector)

S          Compute the sum of z.
  L        Compute the length of z.
 ÷         Divide the former by the latter.
           This computes the mean of z.

Dz_²ÇN½    Main link. Argument: z (vector)

Ç          Apply the previous link, i.e., compute the mean of z.
 ²         Square the mean.
   ²       Square all number in z.
  _        Subtract each squared number from the squared mean.
    Ç      Take the mean of the resulting vector.
     N     Multiply it by -1.
      ½    Take the square root of the result.

J, 18 bayt

[:%:@M*:-M*M=:+/%#

Bu benim APL cevabımın doğrudan çevirisidir J.

Çevrimiçi deneyin!

Simpleks s.0.5 , 43 bayt

Sadece Cuz. Gerçekten bu bir bayt daha golf gerekiyor.

t[@u@RvR]lR1RD@wA@T@{j@@SR2ERpR}u@vR@TR1UEo   
t[      ]                                     ~~ Applies inner function to entire strip (left-to-right)
  @                                           ~~ Copies current value to register
   u                                          ~~ Goes up a strip level
    @                                         ~~ Dumps the register on the current byte
     R                                        ~~ Proceeds right (s1)
      v                                       ~~ Goes back down
       R                                      ~~ Proceeds right (s0)
                                              ~~ Go right until an empty byte is found
         lR1RD                                ~~ Push length, 1, and divide.
              @                               ~~ Store result in register (1/N)
               wA                             ~~ Applies A (add) to each byte, (right-to-left)
                 @T@                          ~~ Puts 1/N down, multiplies it, and copies it to the register
                    {          }              ~~ Repeats until a zero-byte is met
                     j@@                      ~~ inserts a new byte and places register on it
                        SR                    ~~ Subtract it from the current byte and moves right
                          2E                  ~~ Squares result
                            RpR               ~~ Moves to the recently-created cell, deletes it, and continues
                                u@v           ~~ takes 1/N again into register
                                   R@T        ~~ multiplies it by the new sum
                                      R1UE    ~~ takes the square root of previous
                                          o   ~~ output as number

Prolog (SWI), 119 bayt

Kod:

q(U,X,A):-A is(X-U)^2.
p(L):-sumlist(L,S),length(L,I),U is S/I,maplist(q(U),L,A),sumlist(A,B),C is sqrt(B/I),write(C).

Açıklama:

q(U,X,A):-A is(X-U)^2.   % calc squared difference of X and U
p(L):-sumlist(L,S),      % sum input list
      length(L,I),       % length of input list
      U is S/I,          % set U to the mean value of input list
      maplist(q(U),L,A), % set A to the list of squared differences of input and mean
      sumlist(A,B),      % sum squared differences list
      C is sqrt(B/I),    % divide sum of squares by length of list
      write(C).          % print answer

Misal:

p([10035, 436844, 42463, 44774]).
175656.78441352615

Burada çevrimiçi deneyin

Perl5, 39 38

 Komut dosyası için 16
için 22 Manahtar
için +1 Edüğmesi
39 =

perl -MStatistics::Lite=:all -E"say stddevp@ARGV" .1 .2 300

Çilekte Test Edildi 5.20.2.

Oh, ama sonra cevaplarımızın programlar yerine fonksiyonlar olabileceğini söyledin. Bu durumda,

{use Statistics::Lite":all";stddevp@_}

sadece 38. Çilek 5.20.2'de test edilmiştir.

print sub{use Statistics::Lite":all";stddevp@_}->( .1, .2, 300)

Python, 57 bayt

lambda l:(sum((x-sum(l)/len(l))**2for x in l)/len(l))**.5

Girişi liste olarak alır

Teşekkürler @xnor

PowerShell, 122

:\>type stddev.ps1
$y=0;$z=$args -split",";$a=($z|?{$_});$c=$a.Count;$a|%{$y+=$_};$b=$y/$c;$a|%{$x+
=(($_-$b)*($_-$b))/$c};[math]::pow($x,0.5)

açıklama

<#
$y=0                            init
$z=$args -split","              split delim ,
$a=($z|? {$_})                  remove empty items
$c=$a.Count                     count items
$a|%{$y+=$_}                    sum
$b=$y/$c                        average
$a|%{$x+=(($_-$b)*($_-$b))/$c}  sum of squares/count
[math]::pow($x,0.5)             result
#>

sonuç

:\>powershell -nologo -f stddev.ps1 45,67,32,98,11,3
32.5303277300156

:\>powershell -nologo -f stddev.ps1 45,  67,32,98,11,3
32.5303277300156

:\>powershell -nologo -f stddev.ps1 45,  67,32, 98 ,11,3
32.5303277300156

:\>powershell -nologo -f stddev.ps1 10035, 436844, 42463, 44774
175656.784413526

:\>powershell -nologo -f stddev.ps1 1,1,1,1,1,1
0

Fortran, 138 bayt

Fortran'daki denklemin basit bir uygulaması:

double precision function std(x)
integer,dimension(:),intent(in) :: x
std = norm2(dble(x-sum(x)/size(x)))/sqrt(dble(size(x)))
end function

SmileBASIC, 105 bayt (işlev olarak)

Bir fonksiyon olmasına izin verildiğini fark ettim. Hata! Cevabımı önemli ölçüde azaltır. Bu, bir Sdizi alan ve popülasyon standart sapmasını döndüren bir işlevi tanımlar . Bir açıklama için diğerini okuyun, ancak ayrıştırma bölümünü atlayın. Bir daha yapmak istemiyorum.

DEF S(L)N=LEN(L)FOR I=0TO N-1U=U+L[I]NEXT
U=1/N*U FOR I=0TO N-1T=T+POW(L[I]-U,2)NEXT RETURN SQR(1/N*T)END

Bir program olarak, 212 bayt

Ne yazık ki, giriş listesini bir dize olarak almak ve kendim ayrıştırmak zorunda. Bu, yanıta 100 bayt ekler, bu nedenle virgülle ayrılmış bir listeden başka bir giriş biçimine izin verilirse, bunu duymaktan memnuniyet duyarım. Ayrıca VAL, buggy olduğu için virgülden önce bir boşluğa sahip olmanın veya dizenin sonunun programı kırdığını unutmayın. Virgülden sonra veya dizenin başlangıcında sorun yoktur.

DIM L[0]LINPUT L$@L I=INSTR(O,L$,",")IF I>-1THEN PUSH L,VAL(MID$(L$,O,I-O))O=I+1GOTO@L ELSE PUSH L,VAL(MID$(L$,O,LEN(L$)-O))
N=LEN(L)FOR I=0TO N-1U=U+L[I]NEXT
U=1/N*U FOR I=0TO N-1T=T+POW(L[I]-U,2)NEXT?SQR(1/N*T)

Açık ve açık:

DIM L[0]  'define our array
LINPUT L$ 'grab string from input

'parse list
'could've used something cleaner, like a REPEAT, but this was shorter
@L
I=INSTR(O,L$,",")                 'find next comma
IF I>-1 THEN                      'we have a comma
 PUSH L,VAL(MID$(L$,O,I-O))       'get substring of number, parse & store
 O=I+1                            'set next search location
 GOTO @L                          'go again
ELSE                              'we don't have a comma
 PUSH L,VAL(MID$(L$,O,LEN(L$)-O)) 'eat rest of string, parse & store
ENDIF                             'end

N=LEN(L) 'how many numbers we have

'find U
'sum all of the numbers, mult by 1/N
FOR I=0 TO N-1
 U=U+L[I]
NEXT
U=1/N*U

'calculate our popstdev
'sum(pow(x-u,2))
FOR I=0 TO N-1
 T=T+POW(L[I]-U,2)
NEXT
PRINT SQR(1/N*T) 'sqrt(1/n*sum)

Aksiyom, 137 bayt

m(a:List Float):Complex Float==(#a=0=>%i;reduce(+,a)/#a)
s(a:List Float):Complex Float==(#a=0=>%i;n:=m(a);sqrt(m([(x-n)^2 for x in a])))

M () işlevi girişteki listenin ortalamasını döndürür. Hatadaki her iki işlev de% i sanal sabit sqrt (-1) değerini döndürür. Test kodu ve sonuçları. [ama sonuç iyi ise, karmaşık bir sayının gerçek kısmıdır]

(6) -> s([45,67,32,98,11,3])
   (6)  32.5303277300 15604966

(7) -> s([10035,436844,42463,44774])
   (7)  175656.7844135261 4035

(8) -> s([1,1,1,1,1,1])
   (8)  0.0

Python 3 , 49 bayt

lambda l,N:(sum((i-sum(l)/N)**2for i in l)/N)**.5

Çevrimiçi deneyin!

İçinde Alır lsayının bir liste, ve N, tam sayılar mevcut sayısını.

Pyt , 13 bayt

←Đ↔Ł↔е-²Ʃ⇹/√

Standart sapma formülünü uygular