チャレンジ

この課題は非常に簡単です。4つの3次元ポイントが与えられたら、それらが形成する四面体の表面積を計算します。これはcode-golfなので、最短のコードが優先されます。標準的な抜け穴が適用されますが、4つのポイントが与えられた場合にこのタスクを実行する組み込み関数は禁止されるという追加の規定があります。

4つのポイントすべてが明確であり、STDINを介して1行に1ポイント与えられると仮定できます。各ポイントは、3つの16ビット符号なし整数で構成されます。3つのスペースで区切られた整数など、作業を簡単にする場合は、各ポイントの正確な形式を変更できます。ただし、各ポイントを別々の行に置くことは必須です。出力は、STDOUTを介して、少なくとも小数点以下2桁にする必要があります。

知らない人にとっては、四面体は4つの三角形の面で構成される3次元の立体です。

例

# input (format is up to you, see clarification above)
[23822, 47484, 57901]
[3305, 23847, 42159]
[19804, 11366, 14013]
[52278, 28626, 52757]

# output
2932496435.95

私の数学が間違っていることに気づいたら、メモを残してください。

Python、198 178 161文字

V=eval('input(),'*4)
A=0
for i in range(4):F=V[:i]+V[i+1:];a,b,c=map(lambda e:sum((a-b)**2for a,b in zip(*e)),zip(F,F[1:]+F));A+=(4*a*b-(a+b-c)**2)**.5
print A/4

入力形式は質問のとおりです。

各面に隣接するエッジの長さを計算し、Heronの式を使用します。

Matlab / Octave 103

値は変数に格納されると仮定しますc。これは、三角形の面積がその2つの側面ベクトルの外積の半分の長さであるという事実を使用しています。

%input
[23822, 47484, 57901;
3305, 23847, 42159;
19804, 11366, 14013;
52278, 28626, 52757]



%actual code
c=input('');
a=0;
for i=1:4;
    d=c;d(i,:)=[];
    d=d(1:2,:)-[1 1]'*d(3,:);
    a=a+norm(cross(d(1,:),d(2,:)))/2;
end
a

APL、59

f←{+.×⍨⊃1 2-.⌽(⊂⍵)×1 2⌽¨⊂⍺}
.5×.5+.*⍨(f/2-/x),2f/4⍴x←⎕⎕⎕-⊂⎕

クロス積を計算して動作します

説明
1行目は、2つの引数（暗黙の名前付き⍺および⍵）を取り、暗黙的に長さ3の数値配列であると想定し、3dベクトルとして扱い、それらの外積の2乗の大きさを計算する関数を定義します。

                        ⊂⍺   # Wrap the argument in a scalar
                   1 2⌽¨     # Create an array of 2 arrays, by rotating `⊂⍺` by 1 and 2 places
             (⊂⍵)×           # Coordinate-wise multiply each of them with the other argument
        1 2-.⌽               # This is a shorthand for:
        1 2  ⌽               #   Rotate the first array item by 1 and the second by 2
           -.                #   Then subtract the second from the first, coordinate-wise
       ⊃                     # Unwrap the resulting scalar to get the (sorta) cross product
   +.×                       # Calculate the dot product of that...
      ⍨                      # ...with itself
f←{+.×⍨⊃1 2-.⌽(⊂⍵)×1 2⌽¨⊂⍺} # Assign function to `f`

2行目は残りを行います。

                         ⎕⎕⎕-⊂⎕ # Take 4 array inputs, create an array of arrays by subtracting one of them from the other 3
                       x←        # Assign that to x
                     4⍴          # Duplicate the first item and append to the end
                  2f/            # Apply f to each consecutive pair
            2-/x                 # Apply subtraction to consecutive pairs in x
          f/                     # Apply f to the 2 resulting arrays
         (f/2-/x),2f/4⍴x←⎕⎕⎕-⊂⎕ # Concatenate to an array of 4 squared cross products
   .5+.*⍨                        # Again a shorthand for:
   .5  *⍨                        #   Take square root of each element (by raising to 0.5)
     +.                          #   And sum the results
.5×                              # Finally, divide by 2 to get the answer

Python 3、308 298 292 279 258 254

from itertools import*
def a(t,u,v):w=(t+u+v)/2;return(w*(w-t)*(w-u)*(w-v))**.5
z,x,c,v,b,n=((lambda i,j:(sum((i[x]-j[x])**2for x in[0,1,2]))**.5)(x[0],x[1])for*x,in combinations([eval(input())for i in">"*4],2))
print(a(z,x,v)+a(z,c,b)+a(b,v,n)+a(x,c,n))

これは以下を使用します：

• 各行の長さを計算するためのピタゴラスの定理（3D）
• 各三角形の面積を計算するヘロンの式

Mathematica 168 154

これは、四面体のエッジの長さを見つけて使用します ヘロンの式をて面の面積を決定します。

t = Subsets; p = Table[Input[], {4}];
f@{a_, b_, c_} := Module[{s = (a + b + c)/2}, N[Sqrt[s (s - #) (s - #2) (s -#3)] &[a, b, c], 25]]
  Tr[f /@ (EuclideanDistance @@@ t[#, {2}] & /@ t[p, {3}])]

60文字のみを必要とするより直接的なルートがありますが、組み込み関数で各面の面積を計算する限り、ルールに違反しますArea：

p = Table[Input[], {4}];
N[Tr[Area /@ Polygon /@ Subsets[p, {3}]], 25]

セージ– 103

print sum((x*x*y*y-x*y*x*y)^.5for x,y in map(differences,Combinations(eval('vector(input()),'*4),3)))/2

入力読み取り部分は、Keith Randallの回答に基づいています。

パイソン-260

私はあなた自身の質問への回答を投稿する際のエチケットが何であるかわかりませんが、彼女は私の解決策であり、私の例を検証するために使用しました。

import copy,math
P=[input()for i in"1234"]
def e(a, b):return math.sqrt(sum([(b[i]-a[i])**2 for i in range(3)]))
o=0
for j in range(4):p=copy.copy(P);p.pop(j);a,b,c=[e(p[i],p[(i+1)%3])for i in range(3)];s=(a+b+c)/2;A=math.sqrt(s*(s-a)*(s-b)*(s-c));o+=A
print o

laurencevsと同じ手順を使用します。

C、303

不要な空白を除外します。ただし、ここで行うべきゴルフはまだたくさんあります（後で戻ってやってみようと思います。）でforループを宣言したのは初めてです。#define。でです。私はいつもループの数を最小化する方法を見つけました。

テストケースのOPと同じ答えを得るためにからに変更floatするdouble必要がありました。それ以前は、ラウンド300でした。

scanf 入力をスペースまたは改行で区切っても同じように機能するので、好きなだけ行をフォーマットすることができます。

#define F ;for(i=0;i<12;i++)
#define D(k) (q[i]-q[(i+k)%12])
double q[12],r[12],s[4],p,n;

main(i){
  F scanf("%lf",&q[i])
  F r[i/3*3]+=D(3)*D(3),r[i/3*3+1]+=D(6)*D(6)
  F r[i]=sqrt(r[i])
  F i%3||(s[i/3]=r[(i+3)%12]/2),s[i/3]+=r[i]/2
  F i%3||(p=s[i/3]-r[(i+3)%12]),p*=s[i/3]-r[i],n+=(i%3>1)*sqrt(p)   
  ;printf("%lf",n);       
}