Trójkątne pancerniki (problem z geometrią obliczeniową)

18

Jesteś kapitanem pancernika. W tym roku dział inżynierii robi postępy, więc statek, na którym płyniesz, ma kształt prostego trójkąta.

Wychodzisz na pokład i rozkoszujesz się morską bryzą ... choć nie na długo. Wróg strzelił do ciebie! - ale czy strzał trafi?

Wejście

Możesz napisać funkcję lub pełny program dla tego wyzwania.

Twój program przyjmie 11 liczb całkowitych, z których dziesięć jest sparowanych:

  • Pierwsze trzy pary liczb całkowitych (x 1 , y 1 ), (x 2 , y 2 ), (x 3 , y 3 ) określą wierzchołki twojego statku. Utworzony trójkąt będzie miał niezerowe pole.

  • Kolejna para liczb całkowitych (e x , e y ) określa lokalizację działa przeciwnika. Działo wroga nigdy nie będzie leżeć ani w granicach twojego statku. *

  • Para ( x , y ) po tym określa, dokąd skierowany jest wróg. Będzie to różnić się od (e x , e y ).

  • Ostateczna dodatnia liczba całkowita R określa zasięg strzału wroga

* Byłbyś okropnym kapitanem, gdybyś nawet tego nie zauważył!

Wynik

Musisz wydrukować / zwrócić prawdziwą wartość (np. Prawda, 1), jeśli pancernik zostanie trafiony, w przeciwnym razie wartość fałszu (np. Fałsz, 0).

Co to jest hit?

Strzał wroga to odcinek linii prostej o długości R od (e x , e y ) w kierunku (a x , y y ). Jeśli ten segment linii pokrywa się z dowolną częścią wnętrza twojego trójkątnego pancernika, liczy się to jako trafienie. W przeciwnym razie nie jest to hit.

Strzały, które pasą się lub sięgają tylko do granicy trójkąta, nie liczą się jako trafienie.

Przykłady

0 0 0 1 1 0
1 1
0 0
2

test1

Trafienie: wróg strzelił prosto w środek twojego statku!


2 0 0 2 4 4
0 0
1 1
1

test2

Bez trafienia: zasięg wroga jest zbyt krótki, więc jesteś bezpieczny.


0 0 1 2 3 0
-4 0
0 0
8

test3

Bez trafienia: wróg zranił bok twojego statku, więc nie liczy się to jako trafienie. Szczęściarz!


0 0 -1 3 4 -1
-3 -4
3 4
5

test4

Bez trafienia: Strzał wroga zatrzymuje się tuż przed statkiem, więc jesteś bezpieczny. Gdyby działo przeciwnika miało nawet nieco lepszy zasięg, trafiłbyś w cel! Uff!


-2 -3 -3 6 7 -2
-6 2
1 -4
7

test5

Trafienie: Mimo że strzał nie przeszedł na drugą stronę, nadal jest to trafienie.


-3 2 2 -4 7 -3
-3 -4
-3 0
10

test6

Bez trafienia: Dla przypomnienia jest to kolejny bliski traf .


Dodatkowe przypadki testowe

0 0 6 0 6 8
-6 -8
6 8
20

test7

Bez trafienia: to kolejny pas, ale pod kątem.


0 0 -2 -5 5 3
-3 4
0 0
6

test8

Trafienie: Strzał wszedł przez wierzchołek statku.

Punktacja

To jest , więc wygrywa najkrótszy kod w bajtach. Obowiązują standardowe luki .

Sp3000
źródło
Tylko dla pewności, że nie możemy: czy możemy założyć, że statek nie ma dna i że między każdą z dwóch stron jest niewielka szczelina, więc jeśli strzał uda się wejść do statku przez jego róg, liczymy go jako chybienie?
John Dvorak
@JanDvorak Jeśli strzał przebije statek, przechodząc przez wierzchołek, będzie to trafienie, ponieważ segment linii zachodzi na wnętrze statku. W czwartym przykładzie, jeśli zasięg byłby większy niż 5, byłby to trafienie.
Sp3000,
Ile wolno nam grać z argumentami? Czy możemy je grupować, zmieniać kolejność lub wymagać, aby były zmiennoprzecinkowe?
FryAmTheEggman
@FryAmTheEggman W razie potrzeby możesz grupować lub zmieniać kolejność argumentów. Możesz używać liczb zmiennoprzecinkowych, ale program musi działać poprawnie dla małych siatek (powiedzmy do 20x20) bez obawy o precyzję.
Sp3000,
Myślę, że w przykładach brakuje jednego ważnego przypadku, który powoduje, że moje zamierzone rozwiązanie zawodzi: na przykład, gdy statek jest penetrowany przez zakręt 0 0 -1 3 4 -1 -3 -4 3 4 6.
nutki

Odpowiedzi:

3

Python 3, 252 bajty

To z pewnością najbardziej zmienne, których kiedykolwiek używałem w golfie kodowym. : ^ P

from math import*;A=atan2
def h(a,b,c,d,e,f,g,h,i,j,R):
 r=R;_=0
 while r>0:Q=A(j-h,i-g);k,l=g+r*cos(Q),h+r*sin(Q);D=A(d-b,c-a);E=A(f-b,e-a);F=A(l-b,k-a);G=A(b-d,a-c);H=A(f-d,e-c);I=A(l-d,k-c);_=_ or(D<F<E or E<F<D)and(G<I<H or H<I<G);r-=.001
 return _

Nieco golfa z komentarzami:

from math import*
# Parameters:
#  (a,b) (c,d) (e,f) - vertices of the triangle
#  (g,h) - location of cannon
#  (i,j) - aim of cannon
#  R - range of cannon
# Variables within function:
#  _ - was this shot a hit?
#  r - distance 0 < r <= R that we're testing
#  Q - angle between cannon source and destination
#  (k,l) - point that we're testing
#  D,E,F - angles between point 1 and 2,3,test
#  G,H,I - angles between point 2 and 1,3,test
def h(a,b,c,d,e,f,g,h,i,j,R):
    r=R;_=0
    while r>0:
        Q=atan2(j-h,i-g)
        k,l=g+r*cos(Q),h+r*sin(Q)
        D=atan2(d-b,c-a)
        E=atan2(f-b,e-a)
        F=atan2(l-b,k-a)
        G=atan2(b-d,a-c)
        H=atan2(f-d,e-c)
        I=atan2(l-d,k-c)
        _=_ or(D<F<E or E<F<D)and(G<I<H or H<I<G)
        r-=.001
    return _

Jak to działa:

  • Oblicz punkt końcowy strzału.
  • Przetestuj wiele punktów wzdłuż linii od punktu końcowego do położenia armaty:
    • Oblicz kąty od wierzchołka 1 do pozostałych dwóch wierzchołków i do punktu testowego;
    • Oblicz kąty od wierzchołka 2 do pozostałych dwóch wierzchołków i do punktu testowego;
    • Jeśli kąt punktu testowego znajduje się pomiędzy pozostałymi dwoma kątami, w obu przypadkach punkt testowy znajduje się w trójkącie i statek został trafiony.

Przykładowe przebiegi:

>>> h(0,0,0,1,1,0,1,1,0,0,2)
True
>>> h(0,0,1,2,3,0,-4,0,0,0,8)
False
>>> h(0,0,-1,3,4,-1,-3,-4,3,4,5)
False
>>> h(-2,-3,-3,6,7,-2,-6,2,1,-4,7)
True
DLosc
źródło
2

Python 2.7, 235 bajtów

from numpy import*
X=cross
h=lambda q,w,a,s,y,x,c,v,d,f,r:max([(X([a-q,s-w],[c+k*(d-c)-q,v+k*(f-v)-w])>0)==(X([y-a,x-s],[c+k*(d-c)-a,v+k*(f-v)-s])>0)==(X([q-y,w-x],[c+k*(d-c)-y,v+k*(f-v)-x])>0)for k in arange(0,r/hypot(d-c,f-v),1e-4)])

Oblicza iloczyn krzyżowy AB x APmiędzy narożnikami A, B i punktem P. Jeśli wszystkie trzy mają ten sam znak, wówczas punkt znajduje się w trójkącie.

Nie golfowany:

from numpy import *
def i(q,w,a,s,y,x,e,r): # helper-function, checks whether ER is inside the triangle QW-AS-YX
  t=cross([a-q,s-w],[e-q,r-w])>0
  g=cross([y-a,x-s],[e-a,r-s])>0
  b=cross([q-y,w-x],[e-y,r-x])>0
  return t==g==b

def h(q,w,a,s,y,x,c,v,d,f,r):
  R=arange(0,r/hypot(d-c,f-v),1e-3)
  return max([i(q,w,a,s,y,x,c+k*(d-c),v+k*(f-v)) for k in R])

Testy:

In : h(0,0,0,1,1,0,1,1,0,0,2)
Out: True

In : h(-3,2,2,-4,7,-3,-3,-4,-3,0,10)
Out: False

In : h(0,0,1,2,3,0,-4,0,0,0,8)
Out: True
     Grazes may count as hits...
In : h(1,2,0,0,3,0,-4,0,0,0,8)
Out: False
     ...or not, depending on the order of edges
DenDenDo
źródło
1

C, 247 bajtów

Zdecydowanie jeszcze nie dość golfa.

#include<math.h>
int z(float*k){float s=1e-3,t=s,p=k[8]-k[6],q=k[9]-k[7],r=k[10]/hypot(p,q);int w=0;for(;t<1;t+=s){float x=k[6]-k[0]+p*r*t,y=k[7]-k[1]+q*r*t,b=k[2]*k[5]-k[3]*k[4],d=(x*k[5]-y*k[4])/b,e=(x*k[3]-y*k[2])/b;w|=d>0&e<0&d-e<1;}return w;}

Obecnie wykorzystuje to podejście podobne do rozwiązania DLosc, tzn. Iteruje wszystkie możliwe współrzędne na segmencie linii, aby ustalić, czy przecina się on z trójkątem. (Więc nie powiedzie się, jeśli zakres jest większy niż 1000). Jednak używa wzoru z http://mathworld.wolfram.com/TriangleInterior.html, aby ustalić, czy punkt znajduje się wewnątrz trójkąta. Pozwala to uniknąć szeregu funkcji trygonometrycznych.


Przykładowa kontrola powinna zostać wydrukowana 1 0 0 0 1 0.

#include <stdio.h>
int main() {
    {
        float arr[] = {0,0,0,1,1,0,1,1,0,0,2};
        printf("%d\n", z(arr));
    }

    {
        float arr[] = {2,0,0,2,4,4,0,0,1,1,1};
        printf("%d\n", z(arr));
    }

    {
        float arr[] = {0,0,1,2,3,0,-4,0,0,0,8};
        printf("%d\n", z(arr));
    }

    {
        float arr[] = {0,0,-1,3,4,-1,-3,-4,3,4,5};
        printf("%d\n", z(arr));
    }

    {
        float arr[] = {-2,-3,-3,6,7,-2,-6,2,1,-4,7};
        printf("%d\n", z(arr));
    }

    {
        float arr[] = {-3,2,2,-4,7,-3,-3,-4,-3,0,10};
        printf("%d\n", z(arr));
    }
}
kennytm
źródło
1

JavaScript (ES6) 320 448 522 627

(Czy można jeszcze grać w golfa?)

Kroki:

  1. Znajdź rzeczywisty trafiony cel (punkt na odległość r na linii od wroga do celu)
  2. Trafienie: jeśli odcinek od wroga do celu przecina jedną ze stron statku, ale nie w punktach końcowych
  3. Traf też: jeśli cel znajduje się wewnątrz statku - nawet jeśli strzał trafił w wierzchołek - przypadek testowy 8

Ref:
Przecięcie segmentu
Punkt wewnątrz trójkąta
Punkt w odcinku o określonej odległości

Przetestuj w przeglądarce Firefox

C=(i,j,k,l,m,n,g,h,a,b,r,
  d=a-g,e=b-h,f=r/Math.sqrt(d*d+e*e),
  p=g+f*d,q=h+f*e,
  z=j*(m-k)+i*(l-n)+k*n-l*m,
  s=(j*m-i*n+(n-j)*p+(i-m)*q)/z,
  t=(i*l-j*k+(j-l)*p+(k-i)*q)/z,
  S=(i,j,k,l,
     a=k-i,b=l-j,c=p-g,d=q-h,e=i-g,f=j-h,
     s=a*f-b*e,t=c*f-d*e,m=a*d-c*b)=>
     m&&((s/=m)>0&s<1&(t/=m)>0&t<1)
)=>s>0&t>0&s+t<1|S(i,j,k,l)|S(i,j,m,n)|S(m,n,k,l)

// Test
MyOutput.innerHTML = ['Test1', C(0,0, 0,1, 1,0, 1,1, 0,0, 2),
'<br>Test2', C(2,0, 0,2, 4,4, 0,0, 1,1, 1),
'<br>Test3', C(0,0, 1,2, 3,0, -4,0, 0,0, 8),
'<br>Test4', C(0,0, -1,3, 4,-1, -3,-4, 3,4, 5),
'<br>Test5', C(-2,-3, -3,6, 7,-2, -6,2, 1,-4, 7),
'<br>Test6', C(-3,2, 2,-4, 7,-3, -3,-4, -3,0 ,10),
'<br>Test7', C(0,0, 6,0, 6,8, -6,-8, 6,8, 20),
'<br>Test8', C(0,0,-2,-5, 5,3, -3,4, 0,0, 6)];
<div id="MyOutput"></div>

Bez golfa

function check(p0x, p0y, p1x, p1y, p2x, p2y, ex, ey, ax, xy, r)
{
  var sec = function(p0x, p0y, p1x, p1y, p2x, p2y, p3x, p3y)
  {
      var s10x = p1x - p0x, s10y = p1y - p0y, 
          s32x = p3x - p2x, s32y = p3y - p2y,
          s02x = p0x - p2x, s02y = p0y - p2y,
          s = s10x * s02y - s10y * s02x, t = s32x * s02y - s32y * s02x,
          d = s10x * s32y - s32x * s10y;
      return d && (s/=d) > 0 && s<1 && (t/=d) > 0 && t < 1 && [p0x + (t * s10x), p0y + (t * s10y)];
  }
  var pr = function(p0x, p0y, p1x, p1y, r)
  {
      var dx = (p1x-p0x), dy = (p1y-p0y), f = r/Math.sqrt(dx*dx+dy*dy);
      return [p0x + f*dx, p0y+f*dy];
  }
  var inside = function(p0x, p0y, p1x, p1y, p2x, p2y, px, py)
  {
      var area2 = (-p1y*p2x + p0y*(-p1x + p2x) + p0x*(p1y - p2y) + p1x*p2y),
          s = (p0y*p2x - p0x*p2y + (p2y - p0y)*px + (p0x - p2x)*py)/area2,
          t = (p0x*p1y - p0y*p1x + (p0y - p1y)*px + (p1x - p0x)*py)/area2;
      return s > 0 && t > 0 && s+t < 1;
  }
  var tx, xy;
  [tx, ty] = pr(ex, ey, ax, ay, r);

  return inside(p0x, p0y, p1x, p1y, p2x, p2y, tx,ty)
  || sec(p0x, p0y, p1x, p1y, ex, ey, tx, ty)
  || sec(p0x, p0y, p2x, p2y, ex, ey, tx, ty)
  || sec(p2x, p2y, p1x, p1y, ex, ey, tx, ty);
}  
edc65
źródło