10

Votre tâche consiste à prendre un tableau de nombres et un nombre réel et à renvoyer la valeur à ce point du tableau. Les tableaux commencent à et sont comptés dans des intervalles . Le fait est que nous allons en fait interpoler entre les éléments étant donné l’indice. Par exemple: $\pi$ $\pi$

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

Parce que c'est , nous devons faire la trigonométrie obligatoire, nous allons donc utiliser l'interpolation en cosinus en utilisant la formule suivante: $\pi$

${\cos(i \mod \pi) + 1 \over 2} * (\alpha - \beta) + \beta$

où:

$i$ est l'entrée "index"
$\alpha$ est la valeur de l'élément immédiatement avant l '"index"
$\beta$ est la valeur de l'élément immédiatement après l '"index"
$\cos$ prend son angle en radians

Exemple

Étant donné [1.3, 3.7, 6.9], 5.3:

L'index 5.3 est compris entre $1\pi$ et $2\pi$ , donc 1.3 sera utilisé pour beforeet 3.7 sera utilisé pour after. En le mettant dans la formule, nous obtenons:

${\cos(5.3 \mod \pi) + 1 \over 2} * (1.3 - 3.7) + 3.7$

Ce qui sort à 3.165

Remarques

L'entrée et la sortie peuvent être dans n'importe quel format pratique
Vous pouvez supposer que le nombre d'entrée est supérieur à $\pi$ et inférieur à array length* $\pi$
Vous pouvez supposer que le tableau d'entrée comportera au moins 2 éléments.
Votre résultat doit avoir au moins deux décimales de précision, être précis à 0,05 près et prendre en charge jusqu'à 100 pour cette précision / exactitude. (les flotteurs à simple précision sont plus que suffisants pour répondre à cette exigence)

Bon golf!

code-golf array

— Beefster
source

8

Pour les golfeurs de FYI, il pourrait être plus court d'écrire réécrire en utilisant la formule demi-angle pour .

(\cos (x) + 1) / 2

$(\cos(x)+1)/2$

\cos (x / 2)^{2}

$\cos(x/2)^2$

\cos

$\cos$

— xnor

Puis-je prendre un dictionnaire avec des clés comme clés? Les doubles seront des nombres entiers, bien sûr.

— Incarnation de l'ignorance

@EmbodimentofIgnorance, bien sûr. Je doute que cela vous aide, mais c'est une représentation parfaitement raisonnable des tableaux puisque c'est ainsi que Lua le fait.

— Beefster

@KevinCruijssen Je ne vois pas pourquoi cela importerait. 3.7 se situe entre pi et 2pi.

— Beefster

5

R , 59 53 octets

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

Essayez-le en ligne!

Rien de trop intelligent ici - juste une version R de la formule de la question. Merci à @MickyT d'avoir sauvegardé un octet, et à @Giueseppe et indirectement @xnor pour deux autres, et merci à @RobinRyder d'avoir enregistré 3 autres.

— Nick Kennedy
source

Je pense que vous pouvez laisser tomber un octet avec...*(cos(i%%pi)+1)/2

— MickyT

@MickyT merci, j'avais à l'origine mis le +1 entre parenthèses, mais j'avais ajouté une paire de parenthèses redondante pour finir avec 60 octets

— Nick Kennedy

56 octets suite au commentaire de xnor sur la formule demi-angle.

— Giuseppe

1

53 octets

— Robin Ryder

4

Python 3.8 (pré-version) , 85 74 octets

-8 octets grâce à @xnor
-2 octets grâce à @Quintec

Cela profite du nouvel :=opérateur d'affectation de la pré-version de Python 3.8 . En dehors de cela, ce n'est vraiment que l'équation écrite en Python.

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

Usage:

>>> 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

Essayez-le en ligne!

— senox13
source

1

Vous pouvez simplement attribuer jle premier endroit mentionné - une partie du pouvoir des expressions d'affectation est qu'elles évaluent la valeur ainsi que l'assigner.

— xnor

1

Une autre sauvegarde d'octets: utilisez les identités trigonométriques pour convertir (cos(i%pi)+1)/2 encos(i%pi/2)**2

— xnor

@xnor Bon point. Je savais que j'utilisais mal

— senox13

1

Vous pouvez laisser tomber p=puisque les fonctions anonymes sont correctes

— Quintec

1

Vous avez oublié de mettre à jour bytecount :)

— Quintec

3

Gelée , 17 octets

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

Un programme complet acceptant et le tableau qui imprime la valeur interpolée. $i$

Essayez-le en ligne!

Comment?

Interpole entre tous les voisins en utilisant puis sélectionne la valeur appropriée. $\frac{\cos(i \mod \pi)+1}2$

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)

— Jonathan Allan
source

2

C # (Visual C # Interactive Compiler) , 69 octets

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

~~J'ai battu Python!~~ ~~Bon sang, Python m'a battu.~~ J'ai encore battu Python!

Essayez-le en ligne!

— Incarnation de l'ignorance
source

1

Regarder le bon côté des choses. Aucun de nous

— n'avait

@ senox13 Sauf peut-être Stax

— Embodiment of Ignorance

1

Röda , 51 octets

f a,i{j=i//PI;[(cos(i%PI)+1)/2*(a[j-1]-a[j])+a[j]]}

Essayez-le en ligne!

— fergusq
source

1

Stax , 17 octets

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

Exécuter et déboguer

Déballé, non golfé et commenté, cela ressemble à ceci.

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

Exécutez celui-ci

— récursif
source

1

Japt , 47 46 38 octets

u0;J=V/MP[A=UgJc) nUgJf))½McVuMP)Ä]×+A

À suivre ... (golf)

Essayez-le en ligne!

— Quintec
source

1

APL + WIN, 39 37 octets

2 octets économisés grâce à Adám

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

Essayez-le en ligne! Dyalog Classic

Explication:

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

— Graham
source

×.5×→ ÷2÷et m[2]→2⊃m

— Adám

1

Haskell , 65 octets

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

Essayez-le en ligne!

Remarque: le tableau est représenté sous forme de liste.

Merci à @xnor pour la pointe demi-angle.

— dfeuer
source

0

Gelée , 23 20 18 octets

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

Essayez-le en ligne!

÷Ø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
source

0

Attaché , 54 octets

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

Essayez-le en ligne!

Explication

${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

— Conor O'Brien
source

0

C (GCC) 99 79 octets

-20 octets plafondcat

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

Essayez-le en ligne!

Indicatif d'appel

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

notez qu'il avait besoin de l'indicateur du compilateur -lmpour se lier aux bibliothèques mathématiques, donc +3 octets si vous comptez cela.

— rtpax
source

0

05AB1E , 22 21 20 19 octets

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

Essayez-le en ligne ou vérifiez d'autres cas de test .

Explication:

ž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)

— Kevin Cruijssen
source

0

Rubis , 67 octets

->a,i{z=Math::PI;Math.cos(i%z/2)**2*(a[-1+j=(i/z).to_i]-a[j])+a[j]}

Essayez-le en ligne!

— Encre de valeur
source