Массивы начинаются с

10

Ваша задача - взять массив чисел и вещественное число и вернуть значение в этой точке массива. Массивы начинаются с и считаются через интервалы . Дело в том, что мы на самом деле собираемся интерполировать между элементами с учетом «индекса». Например:ππ

Index:    1π   2π   3π   4π   5π   6π
Array: [ 1.1, 1.3, 6.9, 4.2, 1.3, 3.7 ]

Поскольку это , мы должны сделать обязательную тригонометрию, поэтому мы будем использовать косинусную интерполяцию по следующей формуле:π

cos(imodπ)+12(αβ)+β

где:

  • i вход "индекс"
  • α - значение элемента непосредственно перед индексом
  • β - значение элемента сразу после «индекса»
  • cos принимает угол в радианах

пример

Учитывая [1.3, 3.7, 6.9], 5.3:

Индекс 5.3 находится между 1π и 2π , поэтому для 1.3 будет использоваться, а для before3.7 - after. Подставив его в формулу, получим:

cos(5.3modπ)+12(1.33.7)+3.7

Который выходит на 3.165

Ноты

  • Ввод и вывод могут быть в любом удобном формате
  • Вы можете предположить, что введенное число больше π и меньше array length* π
  • Вы можете предположить, что входной массив будет иметь длину не менее 2 элементов.
  • Ваш результат должен иметь как минимум две десятичные точки с точностью до 0,05 и поддерживать числа до 100 для этой точности / точности. (поплавки одинарной точности более чем достаточны для удовлетворения этого требования)

Счастливого гольфа!

Beefster
источник
8
Для игроков в гольф, было бы короче записать переписать как используя формулу угла для . (cos(x)+1)/2cos(x/2)2cos
xnor
Могу ли я взять в словаре с двойными в качестве ключей? Двойники, конечно, будут целыми числами.
Воплощение невежества
@EmbodimentofIgnorance, конечно. Я сомневаюсь, что это поможет вам, но это вполне разумное представление массивов, так как это делает Lua.
Beefster
@KevinCruijssen Я не понимаю, почему это важно. 3.7 находится между пи и 2pi.
Beefster

Ответы:

5

R , 59 53 байта

function(x,i)x[0:1+i%/%pi]%*%c(a<-cos(i%%pi/2)^2,1-a)

Попробуйте онлайн!

Ничего слишком умного здесь - просто R-версия формулы в вопросе. Спасибо @MickyT за сохранение байта, @Giueseppe и косвенно @xnor за еще два, и спасибо @RobinRyder за сохранение еще 3.

Ник Кеннеди
источник
Я думаю, что вы можете бросить байт с...*(cos(i%%pi)+1)/2
MickyT
@MickyT спасибо, я изначально поставил +1 в скобках, но добавил лишнюю пару скобок, заканчивая тем самым 60 байтами
Ник Кеннеди
56 байт после комментария xnor о формуле полуугольника.
Джузеппе
1
53 байта
Робин Райдер
4

Python 3.8 (предварительная версия) , 85 74 байта

-8 байт благодаря @xnor
-2 байт благодаря @Quintec

Это использует преимущества нового :=оператора присваивания предварительной версии Python 3.8 . Кроме этого, это действительно просто уравнение, написанное на Python.

import math
lambda l,i:cos(i%math.pi/2)**2*(l[(j:=int(i/pi))-1]-l[j])+l[j]

Применение:

>>> p=lambda l,i:cos(i%math.pi/2)**2*(l[(j:=int(i/pi))-1]-l[j])+l[j]
>>> print(p([1.3, 3.7, 6.9],5.3))
3.165249203414993

Попробуйте онлайн!

senox13
источник
1
Вы можете просто назначить jпервое упомянутое место - часть силы выражений присваивания заключается в том, что они оценивают значение так же, как и присваивают его.
xnor
1
Сохранение еще (cos(i%pi)+1)/2cos(i%pi/2)**2
xnor
@xnor Хорошая мысль. Я знал, что использовал это неправильно
senox13
1
Вы можете отказаться от p=анонимных функций в порядке
Quintec
1
Забыл обновить bytecount :)
Quintec
3

Желе , 17 байт

d©ØPṪÆẠ‘H×I_@Ḋ}®ị

Полная программа, принимающая и массив, который печатает интерполированное значение.i

Попробуйте онлайн!

Как?

Интерполирует между всеми соседями, используя затем выбирает соответствующее значение.cos(imodπ)+12

d©ØPṪÆẠ‘H×I_@Ḋ}®ị - Link: number, i; list of numbers, A
  ØP              - pi (ish) = 3.141592653589793
d                 - divmod = [i//pi, i%pi]
 ©                - (copy to register for later)
    Ṫ             - tail (gets i%pi leaving register copy as [i//pi])  
     ÆẠ           - cosine = cos(i%pi)
       ‘          - increment
        H         - halve
         ×        - multiply by A (vectorises)
          I       - increments -- i.e. (cos(i%pi)+1)(r-l)/2 for neighbours [l,r]
             Ḋ}   - dequeue A
           _@     - swapped arg subtract (vectorises) -- i.e. r-(cos(i%pi)+1)(r-l)/2
                  -                                         = r+(cos(i%pi)+1)(l-r)/2
               ®  - recall value from the register
                ị - index into (vectorises) -- i.e. [β+(cos(i%pi)+1)(α-β)/2]
                  - implicit print of Jelly representation (only 1 entry so [] wont appear)
Джонатан Аллан
источник
2

C # (интерактивный компилятор Visual C #) , 69 байт

n=>m=>(Math.Cos(m%Math.PI)+1)/2*(n[m=(int)(m/Math.PI)-1]-n[++m])+n[m]

Я победил Питона! Черт возьми, Питон победил меня. Я снова победил Питона!

Попробуйте онлайн!

Воплощение невежества
источник
1
Смотреть на светлую сторону. Ни у
кого
@ senox13 За исключением, может быть, Stax
Воплощение Невежества
1

Stax , 17 байт

≈ëBü☺ÆssÅ¢â)KjjïΔ

Запустите и отладьте его

Распаковано, разглажено и прокомментировано это выглядит так.

VP|%    divmod with pi;  push div and mod results separately
|7^h    do (cos(modpart) + 1) / 2
sX      swap the original div result to top of stack, store it in the x register
v       decrement
;:-     pairwise differences of array
@       get element at index
N*      negate and multiply
;x@     get element from the original array at the x index, where x is the register
+       add

Запустите этот

рекурсивный
источник
1

APL + WIN, 39 37 байт

2 байта сохранены благодаря Адаму

2⊃m+(-/m←⎕[0 1+⌊n÷○1])÷2÷1+2○(○1)|n←⎕

Попробуйте онлайн! Dyalog Classic

Объяснение:

n←⎕ prompt for input of integer

2÷1+2○(○1)|n evaluate first term of formula

[0 1+⌊n÷○1] identify indices of alpha and beta

m←⎕[...] prompt for input of vector and select alpha and beta

-/m alpha-beta

2⊃m+ take result of adding beta to complete the equation 
Грэхем
источник
1

Haskell , 65 байт

v!i|(c,r)<-properFraction$i/pi=cos(r*pi/2)^2*(v!!(c-1)-v!!c)+v!!c

Попробуйте онлайн!

Примечание: массив представлен в виде списка.

Спасибо @xnor за половинный угол.

dfeuer
источник
0

Желе , 23 20 18 байт

³%ØPÆẠ×_++H
÷ØPịÇ/

Попробуйте онлайн!

÷ØPịṁؽµ³%ØPÆẠ×I_@SH    Dyadic link, arguments x (index) and Z (array):
֯P                     x/pi
   ị                    Index (into Z).
                        When x/pi is an integer, returns that elt of Z.
                        Otherwise returns 2 elements at floor and ceiling.
     ؽ                   [1,2] (generic 2 element array)
    ṁؽ                 Mold; shape like [1,2] to ensure we have 2 elements.
       µ                Start a new, monadic chain with the result [a,b]
        ³%ØPÆẠ×I_@SH    Monadic chain
        ³               x
         %ØP            x mod pi
            ÆẠ          Unarccosine; cos(x mod pi).
               I          Increment; b-a.
              ×I        (b-a) cos(x mod pi)
                  S       a+b
                _@S     a + b - (b-a) cos(x mod pi)
                   H    Halve; this is equivalent to our desired result.
lirtosiast
источник
0

Атташе , 54 байта

${Cos[y%PI/2]^2*&`-@(j:=x[1'-1*Floor[y'-y/PI]-1])+j@1}

Попробуйте онлайн!

объяснение

${Cos[y%PI/2]^2*&`-@(j:=x[1'-1*Floor[y'-y/PI]-1])+j@1}
${                                                   }  parameters: x, y
  Cos[y%PI/2]^2                                         the scaling function factor
               *                                        times
                     j:=                                set j to
                        x[                     ]        the element in x at
                          1'-1*Floor[y'-y/PI]-1         the closest indices scaled by PI
                &`-@(                           )       spread subtraction over bounds
                                                 +j@1   add the upper bound
Конор О'Брайен
источник
0

C (GCC) 99 79 байт

Потолок -20 байт

float P=3.141593;b;
#define f(i,a)(cos(fmod(i,P))+1)/2*(a[b=i/P-1]-a[++b])+a[b]

Попробуйте онлайн!

Телефонный код

int main() {
  float a[3] = {1.3,3.7,6.9};
  printf("%f\n", f(5.3,a));
}

обратите внимание, что -lmдля связи с математическими библиотеками необходим флаг компилятора , так что +3 байта, если считать это.

rtpax
источник
0

05AB1E , 22 21 20 19 байтов

žq‰`ž>;UÝèÐÁ-θX*-θ

Попробуйте онлайн или проверьте еще несколько тестов .

Объяснение:

žq        # Take the divmod PI of the (implicit) input-decimal
           # (part = input integer-divided by PI, remainder = input modulo-PI)
           #  i.e. 5.3 → [1, 2.158...]
   `       # Push both values separately to the stack
    ž     # Take the cosine of the remainder
           #  i.e. 2.158... → -0.554...
      >    # Increase it by 1
           #  i.e. -0.554... → 0.554...
       ;   # Halve it
           #  i.e. 0.554... → 0.222...
        U  # Pop and store it in variable `X`
    Ý      # Pop the part, and push a list in the range [0, part]
           #  i.e. 1 → [0, 1]
     è     # (0-based) index all of them into the (implicit) input-list
           #   i.e. [1.3, 3.7, 6.9] and [0, 1] → [1.3, 3.7]
Ð          # Triplicate this list
 Á         # Rotate the last copy once towards the right
           #  i.e. [1.3, 3.7] → [3.7, 1.3]
  -        # Subtract the values in the top two lists from one another
           #  i.e. [1.3, 3.7] and [3.7, 1.3] → [-2.4, 2.4]
   θ       # Pop and only leave the last value of this list
           #  i.e. [-2.4, 2.4] → 2.4
    X*     # Multiply it by `X`
           #  i.e. 2.4 * `X`=0.222... → 0.534...
     -     # Subtract it from each of the values in the list we triplicated
           #  i.e. [1.3, 3.7] - 0.534... → [0.765..., 3.165...]
      θ    # And only leave the last value of this list
           #  i.e. [0.765..., 3.165...] → 3.165...
           # (which is output implicitly as result)
Кевин Круйссен
источник