Các thách thức

Thử thách này rất đơn giản. Cho bốn điểm 3 chiều, tính diện tích bề mặt của tứ diện mà chúng tạo thành. Đây là mã golf , vì vậy mã ngắn nhất sẽ thắng. Các lỗ hổng tiêu chuẩn được áp dụng, với quy định bổ sung rằng bất kỳ chức năng tích hợp nào để thực hiện nhiệm vụ này với bốn điểm đều bị cấm.

Bạn có thể giả sử tất cả bốn điểm sẽ khác biệt và sẽ được cung cấp qua STDIN, 1 điểm trên mỗi dòng. Mỗi điểm sẽ bao gồm ba số nguyên không dấu 16 bit. Định dạng chính xác của từng điểm có thể được sửa đổi nếu nó làm cho mọi thứ dễ dàng hơn, chẳng hạn như ba số nguyên được phân tách bằng dấu cách. Tuy nhiên, có mỗi điểm trên một dòng riêng biệt là bắt buộc. Đầu ra phải thông qua STDOUT, đến ít nhất 2 chữ số thập phân.

Đối với những người bạn không biết, một khối tứ diện là một khối rắn 3 chiều, được hình thành bởi 4 mặt tam giác.

Thí dụ

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

# output
2932496435.95

Vui lòng để lại một lưu ý nếu bạn nhận thấy toán học của tôi sai.

Python, 198 178 161 ký tự

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

Các định dạng đầu vào như được đưa ra trong câu hỏi.

Nó tính toán độ dài của các cạnh liền kề với mỗi khuôn mặt và sau đó sử dụng công thức của Heron .

Matlab / Octave 103

Tôi giả sử các giá trị được lưu trữ trong biến c. Điều này sử dụng thực tế là diện tích của một hình tam giác là một nửa chiều dài của tích chéo của hai vectơ bên của nó.

%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←⎕⎕⎕-⊂⎕

Hoạt động bằng cách tính toán các sản phẩm chéo

Giải thích
Dòng đầu tiên xác định một hàm có hai đối số (hàm ý được đặt tên ⍺và ⍵), hoàn toàn mong đợi chúng là các mảng số có độ dài 3, coi chúng là các vectơ 3d và tính độ lớn bình phương của sản phẩm chéo của chúng.

                        ⊂⍺   # 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`

Dòng thứ hai làm phần còn lại.

                         ⎕⎕⎕-⊂⎕ # 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

Con trăn 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))

Điều này sử dụng:

• Định lý Pythagore (dưới dạng 3D) để tính ra độ dài của mỗi dòng
• Công thức của Heron để tính diện tích của mỗi tam giác

Toán học 168 154

Điều này tìm thấy độ dài của các cạnh của khối tứ diện và sử dụng công thức của Heron để xác định các khu vực của các mặt.

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}])]

Có một tuyến đường trực tiếp hơn chỉ cần 60 ký tự , nhưng nó vi phạm các quy tắc trong khi nó tính diện tích của mỗi mặt với chức năng tích hợp , Area:

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

Hiền nhân - 103

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

Phần đọc đầu vào được điều chỉnh từ câu trả lời của Keith Randall .

Con trăn - 260

Tôi không chắc chắn về nghi thức khi đăng câu trả lời cho câu hỏi của bạn là gì, nhưng cô ấy là giải pháp của tôi, mà tôi đã sử dụng để xác minh ví dụ của mình, đánh gôn:

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

Nó sử dụng thủ tục tương tự như laurencevs.

C, 303

Không bao gồm khoảng trắng không cần thiết. Tuy nhiên, vẫn còn rất nhiều việc chơi golf được thực hiện ở đây (tôi sẽ cố gắng quay lại và thực hiện sau.) Đây là lần đầu tiên tôi tuyên bố một forvòng lặp trong a #define. Trước đây tôi luôn tìm cách giảm thiểu số vòng lặp.

Tôi đã phải thay đổi từ floatđể doublecó được câu trả lời giống như OP cho trường hợp thử nghiệm. Trước đó, nó là một vòng 300.

scanf hoạt động tương tự cho dù bạn tách đầu vào của mình bằng dấu cách hoặc dòng mới, vì vậy bạn có thể định dạng nó thành nhiều hoặc ít dòng tùy thích.

#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);       
}