Создать Падовскую Спираль

34

Введение

Подобно последовательности Фибоначчи, последовательность Падована ( OEIS A000931 ) представляет собой последовательность чисел, которая получается путем добавления предыдущих терминов в последовательность. Начальные значения определяются как:

P(0) = P(1) = P(2) = 1

Все 0-й, 1-й и 2-й члены равны 1. Отношение повторения указано ниже:

P(n) = P(n - 2) + P(n - 3)

Таким образом, он дает следующую последовательность:

1, 1, 1, 2, 2, 3, 4, 5, 7, 9, 12, 16, 21, 28, 37, 49, 65, 86, 114, 151, 200, 265, 351, ...

Использование этих чисел в качестве длины сторон равносторонних треугольников дает хорошую спираль, когда вы размещаете их все вместе, во многом как спираль Фибоначчи:

введите описание изображения здесь

Изображение предоставлено Википедией


задача

Ваша задача - написать программу, воссоздающую эту спираль с помощью графического вывода, с вводом, соответствующим тому или иному термину.

правила

  • Ваша заявка должна быть в состоянии обработать по крайней мере до 10-го семестра (9)
  • Ваша заявка должна быть полной программой или функцией, которая принимает данные и отображает графический результат (выводит изображение или графики и т. Д.)
  • Вы должны показать подтверждение вашего графического вывода в вашем представлении
  • Вращения выходных данных допускаются в кратных 60 градусов с одинаковым представлением
  • Также возможно движение против часовой стрелки
  • Стандартные лазейки запрещены

Вы можете предположить, что ввод будет> 0, и будет указан правильный формат ввода.

счет

Это , поэтому выигрывает самый короткий код в байтах. Всех с новым годом!

Эндрю Ли
источник
Разрешается ли использовать пробел после строк?
Павел
@ Павел Да. Позвольте мне добавить это
Эндрю Ли
Должен ли результат быть идентичным примеру или допустимы отражения и повороты (кратные 60 градусам)?
Уровень River St
@LevelRiverSt Я бы позволил это. Позвольте мне уточнить это в посте.
Эндрю Ли
3
Не фанат того, чтобы использовать ASCII art и графический вывод в одной и той же задаче. Это очень разные задачи, и смешивание их вместе делает ответы на две разные возможности совершенно несопоставимыми. Было бы лучше иметь две отдельные задачи, одну для искусства ASCII и одну для графического вывода.
Мартин Эндер

Ответы:

12

Mathematica, 119 108 байт

Спасибо Мартину Эндеру за сохранение 11 байтов!

±n_:=If[n<4,1,±(n-2)+±(n-3)];Graphics@Line@ReIm@Accumulate@Flatten@{0,z=I^(2/3),±# z^(#+{2,4,1})&~Array~#}&@

Безымянная функция принимает положительный целочисленный аргумент (1-indexed) и возвращает графический вывод. Пример вывода для ввода 16:

введите описание изображения здесь

Разработан одновременно с ответом Flawr на Matlab, но со многими сходствами в дизайне - даже включая определение I^(2/3)шестого корня единства! Более легкая для чтения версия:

1  (±n_:=If[n<4,1,±(n-2)+±(n-3)];
2   Graphics@Line@ReIm@
3   Accumulate@Flatten@
4   {0,z=I^(2/3),±# z^(#+{2,4,1})&~Array~#}
5  ])&

Строка 1 определяет последовательность Падована ±n = P(n). Строка 4 создает вложенный массив комплексных чисел, определяющих zпо пути; последняя часть ±# z^(#+{2,4,1})&~Array~#генерирует множество троек, каждый из которых соответствует векторам, которые мы должны нарисовать, чтобы завершить соответствующий треугольник ( ±#контролирует длину, а z^(#+{2,4,1})контролирует направления). Строка 3 избавляет от вложенности списка, а затем вычисляет промежуточные суммы комплексных чисел для преобразования из векторов в чистые координаты; Затем строка 2 преобразует комплексные числа в упорядоченные пары действительных чисел и выводит соответствующую многоугольную линию.

Грег Мартин
источник
1
не обращайте внимания, что эта часть была только мной, будучи глупым.
Мартин Эндер
9

Matlab, 202 190 байт

N=input('');e=i^(2/3);f=1/e;s=[0,e,1,f,-e,e-2];l=[1,1,1,2];M=N+9;T=[l,2:M-3;2:M+1;3:M+2];for k=5:N;l(k)=l(k-2)+l(k-3);s(k+2)=s(k+1)+e*l(k);e=e*f;end;T=[T;T(1,:)];plot(s(T(:,1:N)));axis equal

Выход для N=19(индексация на основе 1):

введите описание изображения здесь

объяснение

Грубая идея в основном работает с комплексными числами. Тогда края треугольников всегда указывают в направлении шестого корня из единства.

N=input('');                         % Fetch input
e=i^(2/3);                           % 6th root of unity
f=1/e;                               %  "
s=[0,e,1,f,-e,e-2];                  % "s" is a list of vertices in the order as the spiral is defined
l=[1,1,1,2];                         % "l" is a list of edge-lengths of the triangles
for k=5:N;                           % if we need more values in "l"/"s" we calculate those
    l(k)=l(k-2)+l(k-3);
    s(k+2)=s(k+1)+e*l(k);
    e=e*f;
end;
M=N+9;
T=[[1,1,1,2,2:M-3];2:M+1;3:M+2]';    % this matrix describes which vertices from s are needed for each triangle (the cannonical way how meshes of triangles are stored)
trimesh(T(1:N,:),real(s),imag(s));   % plotting the mesh, according to "T"
axis equal
flawr
источник
Хорошая работа! Есть ли возможность объяснения?
Андрей Ли
объяснение добавлено!
flawr
очень нравится использование комплексных чисел здесь.
Дон
7

PHP + SVG, 738 байт

<?php
$a=[1,1,1];
for($i=0;$i<99;)$a[]=$a[$i]+$a[++$i];
$d=$e=$f=$g=$x=$y=0;
$c=[333,999];
$z="";
foreach($a as$k=>$v){
if($k==$_GET['n'])break;
$h=$v/2*sqrt(3);
if($k%6<1){$r=$x+$v/2;$s=$y+$h;$t=$r-$v;$u=$s;}
if($k%6==1){$r=$x-$v/2;$s=$y+$h;$t=$x-$v;$u=$y;}
if($k%6==2){$r=$x-$v;$s=$y;$t=$r+$v/2;$u=$y-$h;}
if($k%6==3){$r=$x-$v/2;$s=$y-$h;$t=$r+$v;$u=$s;}
if($k%6==4){$r=$x+$v/2;$s=$y-$h;$t=$r+$v/2;$u=$y;}
if($k%6>4){$r=$x+$v;$s=$y;$t=$r-$v/2;$u=$y+$h;}
$d=min([$d,$r,$t]);
$e=max([$e,$r,$t]);
$f=min([$f,$s,$u]);
$g=max([$g,$s,$u]); 
$p="M$x,{$y}L$r,{$s}L$t,{$u}Z";
$z.="<path d=$p fill=#{$c[$k%2]} />";
$x=$r;
$y=$s;
}
?>
<svg viewBox=<?="$d,$f,".($e-$d).",".($g-$f)?> width=100% height=100%>
<?=$z?>
</svg>

Выход за 16

<svg viewBox=-53,-12.124355652982,75.5,42.435244785437 width=100% height=100%>
<path d=M0,0L0.5,0.86602540378444L-0.5,0.86602540378444Z fill=#333 /><path d=M0.5,0.86602540378444L0,1.7320508075689L-0.5,0.86602540378444Z fill=#999 /><path d=M0,1.7320508075689L-1,1.7320508075689L-0.5,0.86602540378444Z fill=#333 /><path d=M-1,1.7320508075689L-2,0L0,0Z fill=#999 /><path d=M-2,0L-1,-1.7320508075689L0,0Z fill=#333 /><path d=M-1,-1.7320508075689L2,-1.7320508075689L0.5,0.86602540378444Z fill=#999 /><path d=M2,-1.7320508075689L4,1.7320508075689L0,1.7320508075689Z fill=#333 /><path d=M4,1.7320508075689L1.5,6.0621778264911L-1,1.7320508075689Z fill=#999 /><path d=M1.5,6.0621778264911L-5.5,6.0621778264911L-2,-8.8817841970013E-16Z fill=#333 /><path d=M-5.5,6.0621778264911L-10,-1.7320508075689L-1,-1.7320508075689Z fill=#999 /><path d=M-10,-1.7320508075689L-4,-12.124355652982L2,-1.7320508075689Z fill=#333 /><path d=M-4,-12.124355652982L12,-12.124355652982L4,1.7320508075689Z fill=#999 /><path d=M12,-12.124355652982L22.5,6.0621778264911L1.5,6.0621778264911Z fill=#333 /><path d=M22.5,6.0621778264911L8.5,30.310889132455L-5.5,6.0621778264911Z fill=#999 /><path d=M8.5,30.310889132455L-28.5,30.310889132455L-10,-1.7320508075689Z fill=#333 /><path d=M-28.5,30.310889132455L-53,-12.124355652982L-4,-12.124355652982Z fill=#999 /></svg>

Йорг Хюльсерманн
источник
1
Две маленькие вещи в гольф: $k%6==0может быть $k%6<1и $k%6==5может быть $k%6>4.
Кевин Круйссен,
4

Python 3, 280 , 262 байта

18 байтов сохранено благодаря ovs

Golfed:

import turtle
P=lambda n:n<4or P(n-3)+P(n-2)
N=int(input())
M=9
t=turtle.Turtle()
Q=range
R=t.right
L=t.left
F=t.forward
S=[P(x)*M for x in Q(N,0,-1)]
A=S[0]
F(A)
R(120)
F(A)
R(120)
F(A)
L(120)
i=1
while i<N:
 A=S[i]
 for j in Q(3):F(A);L(120)
 F(A)
 L(60)
 i+=1

То же самое с некоторыми комментариями:

import turtle

# P(n) returns nth term in the sequence
P=lambda n:n<4or P(n-3)+P(n-2)

# M: scales the triangle side-length
M=9
# N: show triangles from 1 to (and including) N from sequence
N=int(input())
t=turtle.Turtle()
Q=range
R=t.right # R(a) -> turn right "a" degrees
L=t.left  # L(a) -> turn left "a" degrees
F=t.forward # F(l) -> move forward "l" units

# S: M*P(N),M*P(N-1), ... M*P(1)
S=[P(x)*M for x in Q(N,0,-1)]

# draw the largest triangle
A=S[0]
F(A)
R(120)
F(A)
R(120)
F(A)
L(120)
i=1

# draw the next N-1 smaller triangles
while i<N:
 A=S[i]
 for j in Q(3):F(A);L(120)
 F(A)
 L(60)
 i+=1

Скриншот для N=9:

N = 9

Bobas_Pett
источник
2

Двиттер 151

s=(n)=>{P=(N)=>N<3||P(N-3)+P(N-2)
for(a=i=0,X=Y=500,x.moveTo(X,Y);i<n*4;i++)k=9*P(i/4),x.lineTo(X+=C(a)
*k,Y+=S(a)*k),x.stroke(),a+=i%4>2?1.047:2.094}

можно проверить на http://dwitter.net (используйте полноэкранный режим)

введите описание изображения здесь

Основная идея - логотип черепаха, игра в гольф. украл функцию P () сверху!

Я полагаю, что больше может быть в гольф рекурсией, но это не плохо.

не яркий
источник
1

LOGO, 119 байт

to s:n
make"x 10
make"y:x
make"z:y
bk:z
repeat:n[lt 60
fw:z
rt 120
fw:z
bk:z
make"w:y+:z
make"z:y
make"y:x
make"x:w]end

Чтобы использовать, сделайте что-то вроде этого :

reset
lt 150
s 12

Пример вывода (не может быть встроен, потому что это не HTTPS и не удалось загрузить в imgur)

Нил
источник