RosettaCodeData/Task/Entropy/PARI-GP/entropy.parigp

2 lines
101 B
Text
Raw Permalink Normal View History

2023-07-01 11:58:00 -04:00
entropy(s)=s=Vec(s);my(v=vecsort(s,,8));-sum(i=1,#v,(x->x*log(x))(sum(j=1,#s,v[i]==s[j])/#s))/log(2)