È una sequenza aritmetico-geometrica?


11

Una sequenza aritmetico-geometrica è il prodotto elementally di una sequenza aritmetica e una sequenza geometrica. Ad esempio, 1 -4 12 -32è il prodotto della sequenza aritmetica 1 2 3 4e della sequenza geometrica 1 -2 4 -8. L'ennesimo termine di una sequenza aritmetico-geometrica intera può essere espresso come

an=rn(a0+nd)

per qualche numero reale , diverso da zero reale e intero . Nota che e non sono necessariamente numeri interi.dra0rd

Ad esempio, la sequenza 2 11 36 100 256 624 1472 3392ha , e .a0=2r=2d=3.5

Ingresso

Un elenco ordinato di numeri interi come input in qualsiasi formato ragionevole. Poiché alcune definizioni di sequenza geometrica consentono e definiscono , se un input è una sequenza aritmetico-geometrica non dipenderà dal fatto che sia autorizzato a 0. Ad esempio, non si verificherà come input.n2r=000=1r123 0 0 0 0

Produzione

Che si tratti di una sequenza aritmetico-geometrica. Emette un valore di verità / falsa o due diversi valori coerenti.

Casi test

Vero:

1 -4 12 -32
0 0 0
-192 0 432 -1296 2916 -5832 10935 -19683
2 11 36 100 256 624 1472 3392
-4374 729 972 567 270 117 48 19
24601 1337 42
0 -2718
-1 -1 0 4 16
2 4 8 16 32 64
2 3 4 5 6 7
0 2 8 24

falso:

4 8 15 16 23 42
3 1 4 1
24601 42 1337
0 0 0 1
0 0 1 0 0
1 -1 0 4 16

1
Cordiali saluti, è possibile utilizzare la modalità matematica in linea con \$per scrivere cose come . a0
FryAmTheEggman,

Sono davvero possibili input a due termini? Non ce ne sono nei casi di test.
xnor

@xnor Trivialmente puoi impostare o quindi le sequenze non sono uniche in quel caso, ma l'output dovrebbe essere sempre veritierod = 0r=1d=0
Giuseppe

1
Suggerisci una prova 0 2 8 24, 0 0 1, 0 0 0 1
tsh

1
1 -1 0 4 16sarebbe un falso caso utile, poiché condivide quattro elementi consecutivi con ciascuno dei casi True 1 -1 0 4 -16e -1 -1 0 4 16.
Anders Kaseorg,

Risposte:


2

Perl 6 , 184 128 135 byte

{3>$_||->\x,\y,\z{?grep ->\r{min (x,{r&&r*$_+(y/r -x)*($×=r)}...*)Z==$_},x??map (y+*×sqrt(y²-x*z).narrow)/x,1,-1!!y&&z/y/2}(|.[^3])}

Provalo online!

Calcola e dai primi tre elementi e verifica se la sequenza risultante corrisponde all'input. Sfortunatamente, Rakudo genera un'eccezione quando si divide per zero, anche quando si usano numeri in virgola mobile, con un costo di ~ 9 byte.rd

Enumera la sequenza usando .an=ran1+rnd

Alcuni miglioramenti sono ispirati alla risposta JavaScript di Arnauld.

Spiegazione

3>$_||  # Return true if there are less than three elements

->\x,\y,\z{ ... }(|.[^3])}  # Bind x,y,z to first three elements

# Candidates for r
x  # If x != 0
??map (y+*×sqrt(y²-x*z).narrow)/x,1,-1  # then solutions of quadratic equation
!!y&&z/y/2  # else solution of linear equation or 0 if y==0

?grep ->\r{ ... },  # Is there an r for which the following is true?

    ( ,                         ...*)  # Create infinite sequence
     x  # Start with x
       {                       }  # Compute next term
        r&&  # 0 if r==0
                (y/r -x)  # d
           r*$_  # r*a(n-1)
                          ($×=r)  # r^n
                +        *  # r*a(n-1)+d*r^n
                                     Z==$_  # Compare with each element of input
min  # All elements are equal?

2

JavaScript (ES7), 135 127 byte

a=>!([x,y,z]=a,1/z)|!a.some(n=>n)|[y/x+(d=(y*y-x*z)**.5/x),y/x-d,z/y/2].some(r=>a.every((v,n)=>(v-(x+n*y/r-n*x)*r**n)**2<1e-9))

Provalo online!

Come?

rd<109

Caso speciale n. 1: meno di 3 termini

Se ci sono meno di 3 termini, è sempre possibile trovare una sequenza corrispondente. Quindi forziamo un valore veritiero.

Caso speciale n. 2: solo zeri

0a0=0d=0r0

a0=0

a0=0

an=rn×n×d

Che dà:

a1=r×da2=2r2×d

d0a10

r=a22a1

a00

an+1an

an+1=r.an+rn+1d

an+2

an+2=r.an+1+rn+2d=r(r.an+rn+1d)+rn+2d=r2an+2r.rn+1d=r2an+2r(an+1r.an)=r2an+2r.an+1

Abbiamo in particolare:

a2=r2a0+2r.a1

Portando alla seguente quadratica:

r2a02r.a1+a2=0

Le cui radici sono:

r0=a1+a12a0a2a0r1=a1a12a0a2a0


Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.