Площадь треугольника


16

Еще один легкий вызов для вас.

Твое задание

Напишите программу или функцию, которая принимает входные данные, которая содержит 3 пары x- и y-координат и вычисляет площадь треугольника, образованного внутри них. Для тех, кто не помнит, как рассчитать его, вы можете найти его здесь .

Пример:

1,2,4,2,3,7       # input as x1,y1,x2,y2,x3,y3
7.5               # output

Посмотреть это на Wolfram Alpha

Некоторые соображения:

  • Входные данные будут шесть целых 10 положительных чисел.
  • Вы можете предположить, что ввод в любом разумном формате .
  • Точки всегда образуют правильный треугольник.
  • Вы можете предположить, что входные данные уже сохранены в такой переменной, как t.
  • Самый короткий код в байтах побеждает!

Изменить: чтобы избежать путаницы, я упростил, как следует обрабатывать ввод, не ставя под угрозу любой из текущих кодов.

Помните, что ваша программа / функция должна выводить правильную область, поэтому она не может дать отрицательное число в качестве вывода


1
Re: ваше редактирование. Означает ли это, что я могу иметь фактический массив пар (например, [[1, 2], [4, 2], [3, 7]]) в T?
Деннис

4
Я все еще в замешательстве. В посте все еще говорится «3 пары» и «шесть ... целых чисел». Обратите внимание, что удаление любого из них приведет к аннулированию некоторых ответов.
xnor

1
Мне не нравится видеть изменение вопроса после публикации и ответа. Но на этот раз я могу сохранить еще 2 байта, так что все в порядке
edc65

1
Если мы можем принять их как три пары, можем ли мы принять их как многомерный массив? То есть [1 2;4 2;3 7](используя синтаксис Юлии)?
Глен О

2
@YiminRong Площадь треугольника не может быть отрицательной по определению. Неважно, в каком порядке находятся очки.
Rainbolt

Ответы:


16

CJam, 18 16 байтов

T(f.-~(+.*:-z.5*

Попробуйте онлайн в интерпретаторе CJam .

идея

Как уже упоминалось в Википедии , площадь треугольника [[0 0] [x y] [z w]]может быть рассчитана как |det([[x y] [z w]])| / 2 = |xw-yz| / 2.

Для общего треугольника [[a b] [c d] [e f]] мы можем перевести его первую вершину в начало координат, получив таким образом треугольник [[0 0] [c-a d-b] [e-a f-b]], площадь которого можно рассчитать по приведенной выше формуле.

Код

T                  e# Push T.
                   e# [[a b] [c d] [e f]]
   (               e# Shift out the first pair.
                   e# [[c d] [e f]] [a b]
    f.-            e# For [c d] and [e f], perform vectorized
                   e# subtraction with [a b].
                   e# [[c-a d-b] [e-a f-b]]
       ~           e# Dump the array on the stack.
                   e# [c-a d-b] [e-a f-b]
        (+         e# Shift and append. Rotates the second array.
                   e# [c-a d-b] [f-b e-a]
          .*       e# Vectorized product.
                   e# [(c-a)(f-b) (d-b)(e-a)]
            :-     e# Reduce by subtraction.
                   e# (c-a)(f-b) - (d-b)(e-a)
              z    e# Apply absolute value.
                   e# |(c-a)(f-b) - (d-b)(e-a)|
               .5* e# Multiply by 0.5.
                   e# |(c-a)(f-b) - (d-b)(e-a)| / 2

10

Mathematica, 27 байт

Area@Polygon@Partition[t,2]

17
Мне нравится, как это использует встроенный и все еще дольше, чем ответ CJAM.
карцигеникат

2
@Carcigenicate реальная проблема - это то Partition[t,2], что соответствует 2/CJam. ;)
Мартин Эндер

10

JavaScript (ES6) 42,44 .

Изменить Формат ввода изменен, я могу сохранить 2 байта

Анонимная функция, которая принимает массив в качестве параметра и возвращает вычисленное значение.

(a,b,c,d,e,f)=>(a*(d-f)+c*(f-b)+e*(b-d))/2

Попробуйте запустить приведенный ниже фрагмент в браузере, совместимом с EcmaScript 6.

f=(a,b,c,d,e,f)=>(a*(d-f)+c*(f-b)+e*(b-d))/2

function test()
{
  var v=I.value.match(/\d+/g)
  I.value = v
  R.innerHTML=f(...v)
}
<input id=I onchange="test()"><button onclick="test()">-></button><span id=R></span>


1
Не могли бы вы просто взять значения в качестве стандартных параметров и сэкономить 2 символа при создании массива?
Mwr247

@ Mwr247 говорит вызовThe input will be a vector with six base 10 positive integers.
edc65

Ага. Сначала я интерпретировал это как значение каждой пары, составляющей вектор координат (такой как пример Вольфрама), в отличие от того, что сам вход ограничен массивом, и поэтому может использовать другие форматы. Теперь имеет больше смысла.
Mwr247

@ Mwr247 теперь ты прав
edc65

8

Юлия, 32 байта

abs(det(t[1:2].-t[[3 5;4 6]]))/2

Создает матрицу соответствующих терминов кросс-произведения, использует detдля получения результирующего значения, принимает абсолютное значение для работы с негативами, а затем делит на 2, потому что это треугольник, а не параллелограмм.


7

Matlab / Octave, 26 байтов

Я не знал об этом встроенном до сих пор =)

polyarea(t(1:2:5),t(2:2:6))

6

Ява, 79 88 байт

float f(int[]a){return Math.abs(a[0]*(a[3]-a[5])+a[2]*(a[5]-a[1])+a[4]*(a[1]-a[3]))/2f;}

Просто использует основную формулу, ничего особенного.

Изменить: Забыл принять абсолютное значение :(


Вам не нужно сделать его работоспособным?
downrep_nation

3
В примере просто показан вызов функции, и это здесь относительно нормальное значение по умолчанию.
Geobits

2
На вопрос: • Вы можете предположить, что входные данные уже сохранены в переменной, такой как 't'. Итак, return(t[0]*(t[3]...должно хватить, нет?
AdmBorkBork

@TimmyD ощущению теневым это делать, но это было бы довести его до 62 байт. Хммм .... Я собираюсь оставить все как есть, пока, по крайней мере.
Geobits

5

Минколанг 0,8 , 34 байта

ndndn0g-n1g-n0g-n0g-1R*1R*-$~2$:N.

Кто-нибудь хочет яйца n0g?

объяснение

Очень просто. Использует формулу |(x2-x1)(y3-y1) - (x3-x1)(y2-y1)|/2.

nd      x1, x1
nd      x1, x1, y1, y1
n0g-    x1, y1, y1, x2-x1
n1g-    x1, y1, x2-x1, y2-y1
n0g-    y1, x2-x1, y2-y1, x3-x1
n0g-    x2-x1, y2-y1, x3-x1, y3-y1
1R*     y3-y1, x2-x1, (y2-y1)(x3-x1)
1R*     (y2-y1)(x3-x1), (y3-y1)(x2-x1)
-       (y2-y1)(x3-x1) - (y3-y1)(x2-x1)
$~      |(y2-y1)(x3-x1) - (y3-y1)(x2-x1)|
2$:     |(y2-y1)(x3-x1) - (y3-y1)(x2-x1)|/2 (float division)
N.      Output as integer and quit.

3

JayScript , 58 байт

Объявляет анонимную функцию:

function(a,b,c,d,e,f){return (a*(d-f)+c*(f-b)+e*(b-d))/2};

Пример:

var nFunct = function(a,b,c,d,e,f){return (a*(d-f)+c*(f-b)+e*(b-d))/2};
print(nFunct(1,2,4,2,3,7));

что г делает?
Уровень Река St

@ steveverrill Ничего, я просто идиот. Исправление ...
mınxomaτ


3

PHP - 68 88 89 байт

Спасибо Марджину за отличные советы!

<?=.5*abs(($t[1]-$t[5])*($t[4]-$t[2])-($t[1]-$t[3])*($t[6]-$t[2]))?>

Чтобы его использовать, создайте файл area.phpс этим содержимым, дополнительная строка соответствует предположению, что данные сохранены в переменнойt части спецификации, а символ ␍ в конце добавляет возврат каретки, так что вывод будет красивым и разделенным:

<?php $t = $argv; ?>
<?=.5*abs(($t[1]-$t[5])*($t[4]-$t[2])-($t[1]-$t[3])*($t[6]-$t[2]))?>
␍

Затем укажите координаты в командной строке x₁ y₁ x₂ y₂ x₃ y₃, например,

$ php area.php 1 2 4 2 3 7
7.5

«Вы можете предположить, что входные данные уже сохранены в такой переменной, как t$a-> $t, удалить $a=$argv;сохранение 9 байтов
Мартейн

After that, you can replace <?php echo with <?=, saving another 7 bytes
Martijn

You can say that this is PHP4.1, with register_globals=On in your php.ini file (default). Read more at php.net/manual/en/security.globals.php
Ismael Miguel


2

R, 37 bytes

cat(abs(det(rbind(matrix(t,2),1))/2))

Converts the vector of coordinates into a matrix and tacks on a row of 1's.
Calculates the determinant and divides by 2.
Returns the absolute result. If the order was always clockwise the abs would not be required.

> t = c(1,2,4,2,3,7)
> cat(det(rbind(matrix(t,2),1))/2)
7.5

2

Python 2, 48 47 50 bytes

Very simple; follows the standard equation:

lambda a,b,c,d,e,f:abs(a*(d-f)+c*(f-b)+e*(b-d))/2.

The other, similarly simple approaches are longer:

def a(a,b,c,d,e,f):return abs(a*(d-f)+c*(f-b)+e*(b-d))/2. # 57
lambda t:abs(t[0]*(t[3]-t[5])+t[2]*(t[5]-t[1])+t[4]*(t[1]-t[3]))/2. # 67
def a(t):return abs(t[0]*(t[3]-t[5])+t[2]*(t[5]-t[1])+t[4]*(t[1]-t[3]))/2. # 74

Python's access to a determinate function is through numpy.

Thanks to muddyfish for 1 byte and xnor for catching an error.


you can remove the 0 from 2.0 to leave 2.
Blue

Quite true, @muddyfish, thanks!
Celeo

Is this Python 2 or 3? Division works differently depending on the version...
mbomb007

Clarified, @mbomb007.
Celeo

1
You need an abs to make the answer positive.
xnor

2

PHP, 77

Based on @Yimin Rong's answer, I felt I could improve upon it by a few bytes by using list() rather than straight $argv to abbreviate some variables. Also echo doesn't need a space if there is delimiter between echo and the thing being echoed.

echo$variable;, echo(4+2);, and echo'some string'; are equally valid whereas echofunction($variable) confuses PHP.

On the other hand, I also added abs() to be mathematically accurate, since some combinations of vertices yielded "negative area"

list($t,$a,$b,$c,$d,$e,$f)=$argv;echo.5*abs(($a-$e)*($d-$b)-($a-$c)*($f-$b));

You can run it via CLI

php -r "list($t,$a,$b,$c,$d,$e,$f)=$argv;echo.5*abs(($a-$e)*($d-$b)-($a-$c)*($f-$b));" 1 2 4 2 3 7
7.5

2

AWK – 51 42 bytes

AWK has no built-in abs so using sqrt(x^2) to substitute.

{print sqrt((($1-$5)*($4-$2)-($1-$3)*($6-$2))^2)/2}

Save as area.awk and use as echo x₁ y₁ x₂ y₂ x₃ y₃ | awk -f area.awk, e.g.

$ echo 1 2 4 2 3 7 | awk -f area.awk
7.5

1

PowerShell, 70 Bytes

[math]::Abs(($t[0]-$t[4])*($t[3]-$t[1])-($t[0]-$t[2])*($t[5]-$t[1]))/2

Uses the same standard formula as other solutions. Per the question, assumes the array is pre-populated, e.g. $t=(1,2,4,2,3,7). But ooof, does the $ and [] syntax kill this one...


Your comment about the penalty from using $ and [] inspired me to try an AWK solution which, by length, is not uncompetitive!

1

dc, 52 bytes

Assumes the input is in register t as: x1 y1 x2 y2 x3 y3 with x1 at the top of t's stack.

1kLtLtsaLtsbLtdscLtltrlalclbltla-*sd-*se-*leld++2/p

1 2 4 2 3 7stStStStStSt #puts coordinates into register t (closest thing dc has to variables) 1kLtLtsaLtsbLtdscLtltrlalclbltla-*sd-*se-*leld++2/p 7.5

This uses the following formula for area:

(x1(y2-y3) + x2(y3-y1) + x3(y1 - y2))/2

And for a quick breakdown of the process:

  • 1k Lt Lt sa Lt sb Lt d sc Lt lt r: set decimal precision to 1 place, move parts of the stack in t to the main stack and move various parts of the main stack to other registers for storage (d duplicates the top of main stack, r reverses the top two elements of main stack, L/l move/copy from the given register to main, s moves top of main stack to the given register)

    Main: y3 x3 y2 x1

    a: y1, b: x2, c: y2, t: y3

  • la lc lb lt la: copy the top of the stacks in registers a, c, b, t, and a to the main stack in that order

    Main: y1 y3 x2 y2 y1 y3 x3 y2 x1

    a: y1, b: x2, c: y2, t: y3

  • - * sd: calculate ((y3-y1)*x2) and put result in d (registers a, b, c, and t are no longer used so I'll drop them from the list of stacks now)

    Main: y2 y1 y3 x3 y2 x1

    d:((y3-y1)*x2)

  • - * se - *: compute ((y1-y2)*y3) and ((y2-x3)*x1); store the former in e and leave the latter on the main stack

    Main: ((y2-x3)*x1)

    d:((y3-y1)*x2), e:((y1-y2)*y3)

  • le ld + +: copy top of register e and d to the main stack, calculate sum of top 2 stack values (pushing result back to main stack) twice

    Main: (((y3-y1)*x2)+((y1-y2)*y3)+((y2-x3)*x1))

    d:((y3-y1)*x2), e:((y1-y2)*y3)

  • 2 /: push 2 onto main stack, divide 2nd values on stack by the 1st (d and e are no longer used, dropping them from list of stacks)

    Main: (((y3-y1)*x2)+((y1-y2)*y3)+((y2-x3)*x1))/2

Rearranging the value on the stack we can see it's equivalent to the formula at the top of this explanation: (x1(y2-y3) + x2(y3-y1) + x3(y1 - y2))/2

  • p: Print top of main stack to output.
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.