Обчисліть точку Ферма трикутника


13

Це дещо схоже на Центри трикутника , але з іншою точкою. Точка Ферма - точка Р у трикутнику ABC така, що значення AP + BP + CP мінімізоване. Є два випадки:

Якщо є кут більше 120 градусів, ця вершина є точкою фермату. В іншому випадку намалюйте рівносторонні трикутники на кожній із сторін ABC. З’єднайте далеку вершину кожного рівностороннього трикутника з протилежною вершиною трикутника ABC. Виконання цього для кожного з трьох рівносторонніх трикутників призводить до єдиної загальної точки перетину для всіх трьох прямих, яка є точкою Ферма.

Він повинен працювати протягом 5 секунд на розумній машині.

Введення : Набір з 3 балів, не обов'язково цілих чисел. Це може сприйматися як вкладений масив, рядок, список кортежів тощо (що відповідає вашій мові).

Вихід : Координати точки Ферма, знову ж таки, проте ваша мова найкраще обробляє точки. Неточності з плаваючою точкою не будуть зараховані проти вас.

Випробування :

[[1, 1], [2, 2], [1, 2]] --> [1.2113248654051871, 1.788675134594813]
[[-1, -1], [-2, -1], [0, 0]] --> [-1, -1]
[[-1, -1], [1, -1], [0, 1]] --> [0, -0.42264973081037427]
[[0, 0], [0.5, 0.8660254037844386], [-5, 0]] --> [0, 0]
[[0, 0], [0, -5], [-0.8660254037844386, 0.5]] --> [0, 0]

Це код гольфу, тому найкоротший виграш коду!


1
Чи добре спробувати всі точки з кроком точності з плаваючою точкою та вибрати ту, яка мінімізує загальну відстань?
xnor

1
@xnor Якщо ви можете це зробити протягом 5 секунд.
soktinpk

На скільки значущих цифр повинен бути точний результат? Крім того, чи добре, якщо -0.0виводиться замість деяких 0.0s?
Р. Кап

@R. Кап, я б сказав про 5 чи 6 значущих цифр. Існує не так багато, що помилки округлення повинні бути проблемою. Щодо другого питання, то це здається прекрасним.
soktinpk

Відповіді:


3

Haskell, 346 291 285 байт

infixl 5£
z=zipWith
(?)=z(-)
t[a,b]=[-b,a]
a¤b=sum$z(*)a b
a%b=t a¤b
r a b c=[c%b/a%b,c%a/a%b]
x£y=2*x¤y<= -sqrt(x¤x*y¤y)
f[a,b,c]|a?b£c?b=b|a?c£b?c=c|b?a£c?a=a|[n,m,p,o]<-c?k a b c++a?k b c a=r[m,o][n,p][c%[n,m],a%[p,o]]
k a b c=map(/2)$z(+)a b?map(signum((b?a)%(c?a))*sqrt 3*)(t$b?a)

Той самий код з деякими поясненнями

infixl 5£
z=zipWith

-- operator ? : difference of two vectors
(?)=z(-)            

-- function t : rotate a vector by +90 degrees
t[a,b]=[-b,a]       

-- operator ¤ : scalar product of two vectors ( a¤b = a0 * b0 + a1 * b1 )
a¤b=sum$z(*)a b     

-- operator % : "cross product" of two vectors ( a%b = a0 * b1 - a1 * b0 )
--      this returns actually the z coordinate of the 3d cross vector
--      other coordinates are nul since a and b are in the xy plan
a%b=t a¤b

-- function r : solves the system of two linear equations with two variables x0,x1 :
--      a0*x0 - b0*x1 = c0
--      a1*x0 - b1*x1 = c1
r a b c=[c%b/a%b,c%a/a%b]

-- operator £ : returns true if the angle between two vectors is >= 120 degrees
--      x¤y = ||x|| * ||y|| * cos(xyAngle)
--      so xyAngle>=120° is equivalent to : x¤y / (||x|| * ||y||) <= -0.5
x£y=2*x¤y<= -sqrt(x¤x*y¤y)

-- function k : takes 3 points A B C of a triangle and constructs the point C' 
--              of the equilateral triangle ABC' which is opposite to C:
--              C' = (A+B)/2 - ((+/-) sqrt(3)/2 * t(AB))
--
--      the sign +/- is given by the sign of the cross vector of AB an AC ((b?a)%(c?a))
--      which is >0 if the angle between AB and AC is positive
--      and <0 otherwise.
k a b c=map(/2)$z(+)a b?map(signum((b?a)%(c?a))*sqrt 3*)(t$b?a)

-- function f : returns the fermat point of a triangle
f[a,b,c]
    |a?b£c?b=b  -- return B if angle ABC >= 120°
    |a?c£b?c=c  -- return C if angle BCA >= 120°
    |b?a£c?a=a  -- return A if angle CAB >= 120°
    |[n,m,p,o]<-c?k a b c++a?k b c a= -- calculate the two segments C'C and A'A
        r[m,o][n,p][c%[n,m],a%[p,o]]  -- return their intersection

Тести:

main = do 
    print $ f [[1, 1], [2, 2], [1, 2]]
    print $ f [[-1, -1], [-2, -1], [0, 0]]
    print $ f [[-1, -1], [1, -1], [0, 1]]
    print $ f [[0, 0], [0.5, 0.8660254037844386], [-5, 0]]
    print $ f [[0, 0], [0, -5], [-0.8660254037844386, 0.5]]

Вихід:

[1.2113248654051871,1.7886751345948126]
[-1.0,-1.0]
[0.0,-0.42264973081037427]
[0.0,0.0]
[0.0,0.0]

Як ви рахуєте байти? £ і ¤ - 2 байти в UTF-8, і я не знаю компілятора Haskell, який приймає ISO-8859-1. (Однак у вас є безліч безкоштовних 1-байтних операторів ASCII, тому це легко виправити.)
Anders Kaseorg

Я рахую це зі своїм редактором, який насправді рахує символів. Я не знав, що це 2 байти, але все одно, як ви сказали, я міг би замінити його іншими операторами на 1 байт. Цей код компілюється з GHC 7.8.3
Damien

GHC компілює код , коли кодується як UTF-8 з , £а ¤також 2 - байтових операторів, але не тоді , коли кодується як ISO-8859-1 з , £і , ¤як 1 байт операторів. Доступні оператори 1 байт в UTF-8 !, #, %, &, ?. Вам слід замінити 2-х байт-операторів або скорегувати кількість байтів.
Андерс Касеорг

2

Пітон, 475 448 440 байт

Будь-яка допомога в гольфі далі цінується.

from math import *
d=lambda x,y:((x[0]-y[0])**2+(x[1]-y[1])**2)**0.5
s=lambda A,B,C:(d(B,C), d(C,A), d(A,B))
j=lambda a,b,c:acos((b*b+c*c-a*a)/(2*b*c))
t=lambda a,b,c:1/cos(j(a,b,c)-pi/6)
b=lambda A,B,C,p,q,r:[(p*A[i]+q*B[i]+r*C[i])/(p+q+r) for i in [0,1]] 
f=lambda A,B,C:A if j(*s(A,B,C)) >= 2*pi/3 else B if j(*s(B,C,A)) >= 2*pi/3 else C if j(*s(C,A,B)) >= 2*pi/3 else b(A,B,C,d(B,C)*t(*s(A,B,C)),d(C,A)*t(*s(B,C,A)),d(A,B)*t(*s(C,A,B)))

Безголівки:

from math import *
#distance between two points
d = lambda x,y: ((x[0]-y[0])**2+(x[1]-y[1])**2)**0.5

#given the points, returns the sides 
s = lambda A,B,C : (d(B,C), d(C,A), d(A,B))

#given the sides, returns the angle
j = lambda a,b,c : acos((b*b+c*c-a*a)/(2*b*c))

#given the sides, returns secant of that angle
t = lambda a,b,c: 1/cos(j(a,b,c)-pi/6)

#given the sides and the Trilinear co-ordinates, returns the Cartesian co-ordinates
b = lambda A,B,C,p,q,r: [(p*A[i]+q*B[i]+r*C[i])/(p+q+r) for i in [0,1]] 

#this one checks if any of the angle is >= 2π/3 returns that point else computes the point
f = lambda A,B,C: A if j(*s(A,B,C)) >= 2*pi/3 else B if j(*s(B,C,A)) >= 2*pi/3 else C if j(*s(C,A,B)) >= 2*pi/3 else b(A,B,C,d(B,C)*t(*s(A,B,C)),d(C,A)*t(*s(B,C,A)),d(A,B)*t(*s(C,A,B)))

Вхід:

print('{}'.format(f([1, 1], [2, 2], [1, 2])))
print('{}'.format(f([-1, -1], [-2, -1], [0, 0])))
print('{}'.format(f([-1, -1], [1, -1], [0, 1])))
print('{}'.format(f([0, 0], [0.5, 0.8660254037844386], [-5, 0])))
print('{}'.format(f([0, 0], [0, -5], [-0.8660254037844386, 0.5])))

Вихід:

[1.2113248652983113, 1.7886751347016887]
[-1, -1]
[0.0, -0.42264973086764884]
[0, 0]
[0, 0]

2
from math import*- досить поширений гольф. Це також дозволить вам використовувати piзамість жорсткого кодування (однакова довжина для 2*pi/3). Ви також можете залишити багато прогалин , як: d=lambda x,y:(....
FryAmTheEggman

2

Python 3.5, 1019 1016 998 982 969 953 байт:

from math import*
def H(z,a,b):c=complex;T=lambda A,B:abs(c(*A)-c(*B));d=T(z,a);e=T(z,b);f=T(a,b);g=[d,e,f];h=max(g);g.remove(h);i=acos((sum(i*i for i in g)-(h*h))/(2*g[0]*g[-1]));_=[[z,a],[z,b],[a,b]];j,s,t=cos,sin,atan;N=[[b,a]for a,b in zip([b,a,z],[max(i,key=i.get)if i!=''else''for i in[{(g[0]+(h*j(t(l))),g[1]+(h*s(t(l)))):T(k,(g[0]+(h*j(t(l))),g[1]+(h*s(t(l))))),(g[0]-(h*j(t(l))),g[1]-(h*s(t(l)))):T(k,(g[0]-(h*j(t(l))),g[1]-(h*s(t(l)))))}if l else{(g[0]+h,g[1]):T(k,(g[0]+h,g[1])),(g[0]-h,g[1]):T(k,(g[0]-h,g[1]))}if l==0else''for g,h,l,k in zip([((a[0]+b[0])/2,(a[1]+b[1])/2)for a,b in _],[(3**0.5)*(i/2)for i in[d,e,f]],[-1/p if p else''if p==0else 0for p in[((a[1]-b[1])/(a[0]-b[0]))if a[0]-b[0]else''for a,b in _]],[b,a,z])]])if b!=''];I=N[0][0][1];J=N[0][0][0];K=N[1][0][1];G=N[1][0][0];A=(N[0][1][1]-I)/(N[0][1][0]-J);B=I-(A*J);C=(K-N[1][1][1])/(G-N[1][1][0]);D=K-(C*G);X=(D-B)/(A-C);Y=(A*X)+B;return[[X,Y],[[a,b][h==d],z][h==f]][i>2.0943]

Неймовірно довго в порівнянні з іншими відповідями, але ей, принаймні, це працює! Я не міг бути щасливішим з результатом, який я отримав, оскільки це має бути однією з найважчих проблем, яку я коли-небудь робив. Я просто такий щасливий, що це насправді працює! : D Тепер, до більш технічних зауважень:

  • Ця функція приймає кожну впорядковану пару як список або кортеж. Наприклад, H((1,1),(2,2),(1,2))буде працювати, але так буде H([1,1],[2,2],[1,2]).
  • Виводить координати точок у списку цілих чи плаваючих точок залежно від того, чи існує один кут, більший або рівний 120º.
  • Це може вивести -0.0замість 0.0деяких входів. Наприклад, вихід для вводу [-1, -1], [1, -1], [0, 1]є [-0.0, -0.4226497308103744]. Я сподіваюся, що це нормально, хоча якщо це не так, я його зміню, хоча це коштуватиме мені ще декількох байтів. Це нормально, як підтвердили ОП .
  • Повинно бути точним, щонайменше, 13до 14значних цифр.

Я з часом спробую пограти ще більше. Пояснення, можливо, дуже довге, незабаром.

Спробуйте в Інтернеті! (Ідеон)


1

Математика, 39 байт

Sum[Norm[p-{x,y}],{p,#}]~NArgMin~{x,y}&

Побудує рівняння на основі відстані між вершинами і точкою {x,y}. Потім використовує NArgMinфункцію, щоб знайти глобальний мінімум для цього рівняння, який буде точкою Ферма за визначенням.

Приклад


1
39 байт, коли наступна найкоротша відповідь - 285 ...
Белінт
Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.