Trova i massimi e i minimi locali


14

Definizione

I massimi e i minimi di una determinata funzione sono i valori più grandi e più piccoli della funzione all'interno di un determinato intervallo o comunque nell'intero dominio della funzione.

Sfida

La sfida è trovare i massimi e i minimi locali di una determinata funzione polinomiale usando qualsiasi metodo che ti piaccia . Non preoccuparti, farò del mio meglio per spiegare la sfida e renderla il più semplice possibile.

L'input conterrà tutti i coefficienti del singolo polinomio variabile in ordine decrescente o crescente di potenza (fino a te). Per esempio,

  • [3,-7,1] rappresenterà 3x2 - 7x + 1 = 0
  • [4,0,0,-3] rappresenterà 4x3-3=0.

Come risolvere (utilizzando i derivati)?

Ora, supponiamo che il nostro input sia [1,-12,45,8], che non è altro che la funzione .x3 - 12x2 + 45x + 8

  1. Il primo compito è trovare la derivata di quella funzione. Dal momento che è una funzione polinomiale, quindi è davvero un compito semplice da fare.

    Il derivato di è . Tutti i termini costanti presenti sono semplicemente moltiplicati. Inoltre, se ci sono termini aggiunti / sottratti, anche i loro derivati ​​vengono aggiunti o sottratti rispettivamente. Ricorda, la derivata di qualsiasi valore numerico costante è zero. Ecco alcuni esempi:xnn*xn-1xn

    • x3 -> 3x2
    • 9x4 -> 9*4*x3 = 36x3
    • -5x2 -> -5*2*x = - 10x
    • 2x3 - 3x2 + 7x -> 6x2 - 6x + 7
    • 4x2 - 3 -> 8x - 0 = 8x
  2. Ora risolvi l'equazione equiparando il nuovo polinomio a zero e ottieni solo i valori integrali di x.

  3. Inserire quei valori di x nella funzione originale e restituire i risultati. Questo dovrebbe essere l'output .

Esempio

Facciamo l'esempio che abbiamo menzionato prima, cioè [1,-12,45,8].

  • Ingresso: [1,-12,45,8]
  • Funzione: x3 - 12x2 + 45x + 8
  • Derivata -> 3x2 - 24x + 45 + 0 -> [3,-24,45]
  • Risolvendo l'equazione , otteniamo o .3x2 - 24x + 45 = 0x = 3x = 5
  • Ora inserendo x = 3e x = 5nella funzione, otteniamo i valori (62,58).
  • Uscita -> [62,58]

ipotesi

  1. Supponiamo che tutti i coefficienti di input siano numeri interi . Possono essere in ordine crescente o decrescente di potere.

  2. Supponiamo che l' input sia almeno un polinomio di 2 gradi . Se il polinomio non ha soluzioni intere, è possibile restituire qualsiasi cosa.

  3. Supponiamo che il risultato finale sarà solo numeri interi.

  4. È possibile stampare i risultati in qualsiasi ordine. Il grado del polinomio di input non sarebbe superiore a 5, in modo che il codice possa gestirlo.

  5. L'input sarà valido in modo che le soluzioni di x non siano punti di sella.

Inoltre, non sei obbligato a farlo con il metodo derivato. Puoi usare qualsiasi metodo tu voglia.

Esempio di input e output

[2,-8,0] -> (-8)
[2,3,-36,10] -> (91,-34)
[1,-8,22,-24,8] -> (-1,0,-1) 
[1,0,0] -> (0)

punteggio

Questo è quindi vince il codice più corto.


1
Se ho capito bene: nell'esempio il passaggio " Risolvere l'equazione " sarebbe in parte questa tua precedente sfida ? Inoltre, il passaggio " Ora inserendo x = 3 e x = 5 nella funzione " indica la funzione originale in " Funzione " e non la funzione in " Derivata ", giusto?
Kevin Cruijssen,

1
Per il campione I / O 3, sto ottenendo (-1, 0, 1), che credo sia la risposta corretta effettiva ... non sono sicuro però. Se non sei d'accordo con me, chiamami in chat.
HyperNeutrino,

1
The input will be valid so that the solutions of x are not saddle points, il caso [1,0,0,3]sembra dare un punto a sella.
JungHwan Min

1
@JungHwanMin ah quell'esempio è stato aggiunto prima che fosse fatta la regola. Rimosso ora.
Manish Kundu,

1
x^3 - 12x^2 + 45x + 8 = 0 , anche se personalmente preferisco scriverlo come f(x)=x^3-12x^2+45x+8senza il =0perché =0non ha senso dal momento che abbiamo a che fare con una funzione, non risolvere un'equazione.
Weijun Zhou

Risposte:


4

Gelatina , 20 byte

ASŒRḅ@Ðḟ
J’U×µṖÇḅ@€³

Provalo online!

Spiegazione

ASŒRḅ@Ðḟ     Helper Function; find all integer solutions to a polynomial
             All integer roots are within the symmetric range of the sum of the absolute values of the coefficients
A            Absolute Value (Of Each)
 S           Sum
  ŒR         Symmetric Range; `n -> [-n, n]`
      Ðḟ     Filter; keep elements where the result is falsy for:
    ḅ@       Base conversion, which acts like the application of the polynomial
J’U×µṖÇḅ@€³  Main Link
J                             Range of length
 ’                    Lowered
  U          Reversed
   ×         Multiplied with the original list (last value is 0)
    µ        Begin new monadic chain
     Ṗ       Pop; all but the last element
      Ç      Apply last link (get all integer solutions of the derivative)
       ḅ@€³  Base conversion of the polynomial into each of the solutions; apply polynomial to each solution of the derivative.

La funzione di aiuto in questo programma è stata presa dalla risposta del signor Xcoder qui che si basava sulla risposta di Luis qui


@JungHwanMin Lo segnalerò a OP. Questa è una violazione diretta dell'affermazione che non ci saranno punti di sella perché la derivata del polinomio 3è 0. modifica oh hai già fatto nvm appena votato il commento allora
HyperNeutrino

3

JavaScript (ES7), 129 120 byte

Prende i coefficienti in ordine crescente di potenza.

a=>(g=x=>x+k?(A=g(x-1),h=e=>a.reduce((s,n,i)=>s+n*(e||i&&i--)*x**i,0))()?A:[h(1),...A]:[])(k=Math.max(...a.map(n=>n*n)))

Casi test

Commentate

a => (                        // given the input array a[]
  g = x =>                    // g = recursive function checking whether x is a solution
    x + k ? (                 //   if x != -k:
      A = g(x - 1),           //     A[] = result of a recursive call with x - 1
      h = e =>                //     h = function evaluating the polynomial:
        a.reduce((s, n, i) => //       for each coefficient n at position i:
          s +                 //         add to s
          n                   //         the coefficient multiplied by
          * (e || i && i--)   //         either 1 (if e = 1) or i (if e is undefined)
          * x**i,             //         multiplied by x**i or x**(i-1)
          0                   //         initial value of s
        )                     //       end of reduce()
      )() ?                   //     if h() is non-zero:
        A                     //       just return A[]
      :                       //     else:
        [h(1), ...A]          //       prepend h(1) to A[]
    :                         //   else:
      []                      //     stop recursion
  )(k = Math.max(             // initial call to g() with x = k = maximum of
    ...a.map(n => n * n)      // the squared coefficients of the polynomial
  ))                          // (Math.abs would be more efficient, but longer)

1
non riesce per 0,0,1(x ^ 2 = 0)
betseg il

@betseg Grazie per aver segnalato questo. Fisso.
Arnauld,

3

Julia 0.6 (con Polynomialspacchetto), 57 byte

using Polynomials
x->map(Poly(x),roots(polyder(Poly(x))))

Provalo online!

Accetta i coefficienti in ordine crescente, ovvero il primo input è il termine costante.

Esempio di esecuzione:

julia> using Polynomials

julia> g = x -> map(Poly(x), roots(polyder(Poly(x))))
(::#1) (generic function with 1 method)

julia> g([8,45,-12,1])
2-element Array{Float64,1}:
 58.0
 62.0

3

Java 8, 364 239 227 226 218 byte

a->{int l=a.length,A[]=a.clone(),s=0,i,r,f=l,p;for(;f>0;A[--f]*=f);for(int n:A)s+=n<0?-n:n;for(r=~s;r++<s;){for(p=0,i=f=1;i<l;f*=r)p+=A[i++]*f;if(p==0){for(f=i=0;i<l;f+=a[i++]*Math.pow(r,p++));System.out.println(f);}}}

Utilizza la stessa funzionalità da questa mia risposta.

-8 byte grazie a @ OlivierGrégoire portando l'array in ordine inverso.

Spiegazione:

Provalo online.

a->{                  // Method with integer-varargs parameter and integer return-type
  int l=a.length,     //  The length of the input-array
      A[]=a.clone(),  //  Copy of the input-array
      s=0,            //  Sum-integer, starting at 0
      i,              //  Index-integer
      r,              //  Range-integer
      f=l,            //  Factor-integer, starting at `l`
      p;              //  Polynomial-integer
  for(;f>0;           //  Loop over the copy-array
    A[--f]*=f);       //   And multiply each value with it's index
                      //   (i.e. [8,45,-12,1] becomes [0,45,-24,3])
  for(int n:A)        //  Loop over this copy-array:
    s+=n<0?-n:n;      //   And sum their absolute values
  for(r=~s;r++<s;){   //  Loop `r` from `-s` up to `s` (inclusive) (where `s` is the sum)
    for(p=0,          //   Reset `p` to 0
        i=f=1;        //   and `f` to 1
                      //   (`i` is 1 to skip the first item in the copy-array)
        i<l;          //   Inner loop over the input again, this time with index (`i`)
        f*=r)         //     After every iteration: multiply `f` with the current `r`
      p+=             //    Sum the Polynomial-integer `p` with:
         A[i++]       //     The value of the input at index `i`,
               *f;}   //     multiplied with the current factor `f`
    if(p==0){         //   If `p` is now 0:
      for(f=i=0;      //    Use `f` as sum, and reset it to 0
          i<l;        //    Loop over the input-array
        f+=a[i++]*Math.pow(r,p++));
                      //     Fill in `r` in the parts of the input-function
      System.out.println(f);}}}
                      //    And print the sum

2
non riesce per 1,0,0(x ^ 2 = 0)
betseg il

@betseg Grazie! Fisso e golf.
Kevin Cruijssen,

1
Si dovrebbe accettare l'ingresso in ordine inverso (è esplicitamente consentito) per ridurre il valore in questo modo: int... ,i, ...; for(;f>0;)A[--f]*=f;. A meno che non mi sbagli, questo dovrebbe farti risparmiare almeno 4 byte. In tal caso, assicurarsi di invertire tutti gli accessi all'input.
Olivier Grégoire,

@ OlivierGrégoire Grazie, 8 byte salvati!
Kevin Cruijssen,


2

Haskell , 89 byte

-3 byte grazie a Laikoni.

r#l=foldr1((.(r*)).(+))l
r l|_:d<-zipWith(*)[0..]l,t<-sum$abs<$>d=[i#l|i<-[-t..t],i#d==0]

Provalo online!

Accetta i coefficienti invertiti.


2
d<-tail$può essere abbreviato in _:d<-.
Laikoni,

1

Python 3 , 156 byte

def f(p,E=enumerate):a=lambda p,v:sum(e*v**i for i,e in E(p));d=[e*i for i,e in E(p)][1:];r=sum(map(abs,d));return[a(p,x)for x in range(-r,r+1)if a(d,x)==0]

Provalo online!

-2 byte grazie a Mr. Xcoder
-22 byte grazie a ovs


1

Python + numpy, 91

  • 1 byte salvato grazie a @KevinCruijssen
from numpy import*
p=poly1d(input())
print map(lambda x:int(polyval(p,x)),roots(p.deriv()))

Provalo online .


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.