Diferencia del cuadrado de la suma


37

Encuentra la diferencia entre el cuadrado de las sumas y la suma de los cuadrados.

Esta es la representación matemática:

(n)2n2

Su programa / método debe tomar dos entradas, estos son sus límites inferior y superior del rango, e incluyen. Los límites serán enteros enteros por encima de 0.

Su programa / método debe devolver la respuesta.

Puede usar la base que desee, pero indique en su respuesta qué base ha utilizado.

Caso de prueba (Base 10)

5,9      970
91,123   12087152
1,10     2640

Este es el código de golf habitual, por lo que cuanto más corta sea la respuesta, mejor.


11
Me tomó un tiempo darme cuenta de que la entrada eran los puntos finales de un rango.
Brad Gilbert b2gills

@ BradGilbertb2gills editado para mayor claridad
george

Esto es más simple de lo que parece?
gato

@cat, ¿qué quieres decir con eso? Sí, las matemáticas son simples cosas de Alevel. Pero todo depende de cómo juegues golf
george

@george La pregunta y muchas de las respuestas que se vea como un montón de trabajo, pero es no
gato

Respuestas:


23

Python 2, 43 bytes

f=lambda a,b,s=0:b/a and 2*a*s+f(a+1,b,s+a)

Pruébalo en Ideone .

Cómo funciona

Llame a la función definida en la especificación g (a, b) . Tenemos eso

Defina la función f (x, y, s) de forma recursiva de la siguiente manera.

Al aplicar la relación de recurrencia de f (a, b, 0) un total de b - a veces, podemos mostrar eso.

Esta es la función f de la implementación. Mientrasb/a devuelve un número entero distinto de cero, andse ejecuta el siguiente código , implementando así la definición recursiva de f .

Una vez que b/allega a 0 , tenemos que b> a y la lambda devuelve False = 0 , implementando así el caso base de la definición de f .


Ah bien. ¿Podrías explicar tu método?
George

Lo haré, pero actualmente estoy tratando de jugar un poco más.
Dennis

Gracias por la fórmula. Supongo que nunca lo vi así porque no cubrimos sumas de series como esa en la escuela. ¡Aunque bastante interesante!
George

2
@george He terminado la explicación.
Dennis

¡Quiere contarnos un poco más sobre cómo en el mundo se le ocurrió la idea de definir f! ¡La motivacion! Estoy realmente interesado
Musa Al-hassy

15

MATL , 9 bytes

&:&*XRssE

Pruébalo en línea!

Explicación

&:   % Inclusive range between the two implicit inputs
&*   % Matrix of all pair-wise products
XR   % Upper triangular part of matrix, without the diagonal
ss   % Sum of all elements of the matrix
E    % Multiply by 2. Implicit display

Ejemplo

Estos son los resultados parciales de cada línea para entradas 5y 9:

  1. &:

    5 6 7 8 9
    
  2. &:&*

    25 30 35 40 45
    30 36 42 48 54
    35 42 49 56 63
    40 48 56 64 72
    45 54 63 72 81
    
  3. &:&*XR

    0 30 35 40 45
    0  0 42 48 54
    0  0  0 56 63
    0  0  0  0 72
    0  0  0  0  0
    
  4. &:&*XRss

    485
    
  5. &:&*XRssE

    970
    

77
Realmente me gusta ver los resultados parciales. Realmente ayudan a entender el programa. ¡Gracias por incluirlos!
DanTheMan

10

Jalea, 9 8 bytes

rµS²_²S$

Pruébalo en línea!

r         inclusive range from first input to second input
 µ        pass the range to a new monadic chain
  S       the sum
   ²      squared
    _     minus...
     ²S$  the squares summed

¡Gracias a FryAmTheEggman por un byte!


3
Por una vez, Jelly es realmente muy legible.
Adám

¿Puedo bifurcar esto a mi respuesta?
Leaky Nun

@LeakyNun, ¿qué significa eso?
Pomo de la puerta


66
Bonitos pendientes: S²_²S
Thomas Weller

10

Python 2, 45 bytes

lambda a,b:(a+~b)*(a-b)*(3*(a+b)**2+a-b-2)/12

Solución de formulario cerrado: no es la más corta, pero pensé que valdría la pena publicarla de todos modos.

Explicación

Vamos a p(n)ser el n º número piramidal cuadrado , y de t(n)ser el n º número triangular . Entonces, para n sobre el rango a , ..., b :

  • ∑n = t(b)-t(a-1), y
  • ∑n² = p(b) - p(a-1)
  • Entonces (∑n) ²-∑n² = (t(b)-t(a-1))² - (p(b) - p(a-1)).

Esta expresión se reduce a la del código.


Hola, ¿podrías explicar tu ecuación si es posible? Mi versión pitón es de 16 bytes más largo y no puedo averiguar cómo se derivó la ecuación
George

1
@george Sea p(n)el nnúmero piramidal cuadrado th y t(n)el nnúmero triangular th . Entonces esta es una versión simplificada de (t(b)-t(a-1))^2 - (p(b) - p(a-1)).
Martin Ender

@MartinEnder Así que esa es la fórmula exacta que he usado, pero Sp3000 la ha simplificado de una manera que no puedo entender. Mi script de Python es: (b * - ~ ba * ~ -a) ** 2 / 4- (b * - ~ b * (2 * b + 1) -a * ~ -a * (2 * a-1) ) / 6 si eso es de alguna utilidad. He jugado al golf tanto como puedo las dos fórmulas
George

@george A veces, con problemas como estos, la forma más fácil es hacer que Wolfram | Alpha haga la parte tediosa, luego verifique dos veces para asegurarse de que sea correcto. Para ser honesto, no creo que pudiera haber eliminado el (a-b-1)factor por (b*(b+1)*(2b+1)-a*(a-1)*(2a-1))/6mi cuenta.
Sp3000

@Sp3000 that's a great way to do it. I'll try that in future
george

6

05AB1E, 8 bytes

ŸDOnsnO-

Explained

ŸD       # range from a to b, duplicate
  On     # sum and square first range
    s    # swap top 2 elements
     nO  # square and sum 2nd range
       - # take difference

Try it online


Is 05AB1E a ROT13 version of Jelly maybe? Substitute r by Ÿ, µ by D, S by O, ² by n, _ by s and $ by -.
Thomas Weller

4
@ThomasWeller: Son bastante diferentes en realidad. Un desplazamiento común entre algunas "funciones" es muy probable que sea una coincidencia. Jelly es un lenguaje tácito sobre funciones de encadenamiento (afaik), mientras que 05AB1E es un lenguaje basado en pila.
Emigna

6

Mathematica, 21 bytes

Tr[x=Range@##]^2-x.x&

Una función sin nombre que toma dos argumentos y devuelve la diferencia. Uso:

Tr[x=Range@##]^2-x.x&[91, 123]
(* 12087152 *)

Aquí hay tres trucos de golf pequeños (y bastante estándar):

  • ##representa ambos argumentos a la vez, por lo que podemos usar la notación de prefijo para Range. Range@##es una abreviatura para la Range[##]que se expande Range[a, b]y nos brinda un rango inclusivo según sea necesario.
  • Tres para rastrear, pero usarlo en un vector simplemente suma ese vector, ahorrando tres bytes Total.
  • x.xes un producto de punto que ahorra cuatro bytes Tr[x^2].

Ayudaría Variance?
Leaky Nun

@LeakyNun I don't see how, because one of the two terms in Variance is divided by n and the other by n^2 and I don't see an easy way to undo those separately.
Martin Ender

1
Tr@#^2-#.#&@*Range is only 18 bytes.
Misha Lavrov

@MishaLavrov neat! Feel free to make it a separate answer. :)
Martin Ender

5

Labyrinth, 28 24 bytes

?:?:}+=-:(:(#{:**+**#2/!

Try it online!

Explanation

Since loops tend to be expensive in Labyrinth, I figured the explicit formula should be shortest, as it can be expressed as linear code.

Cmd Explanation                 Stacks [ Main | Aux ]
?   Read M.                     [ M | ]
:   Duplicate.                  [ M M | ]
?   Read N.                     [ M M N | ]
:   Duplicate.                  [ M M N N | ]
}   Move copy to aux.           [ M M N | N ]
+   Add.                        [ M (M+N) | N ]
=   Swap tops of stacks.        [ M N | (M+N) ]
-   Subtract.                   [ (M-N) | (M+N) ]
:   Duplicate.                  [ (M-N) (M-N) | (M+N) ]
(   Decrement.                  [ (M-N) (M-N-1) | (M+N) ]
:   Duplicate.                  [ (M-N) (M-N-1) (M-N-1) | (M+N) ]
(   Decrement.                  [ (M-N) (M-N-1) (M-N-2) | (M+N) ]
#   Push stack depth.           [ (M-N) (M-N-1) (M-N-2) 3 | (M+N) ]
{   Pull (M+N) over from aux.   [ (M-N) (M-N-1) (M-N-2) 3 (M+N) | ]
:   Duplicate.                  [ (M-N) (M-N-1) (M-N-2) 3 (M+N) (M+N) | ]
*   Multiply.                   [ (M-N) (M-N-1) (M-N-2) 3 ((M+N)^2) | ]
*   Multiply.                   [ (M-N) (M-N-1) (M-N-2) (3*(M+N)^2) | ]
+   Add.                        [ (M-N) (M-N-1) (3*(M+N)^2 + M - N - 2) | ]
*   Multiply.                   [ (M-N) ((M-N-1)*(3*(M+N)^2 + M - N - 2)) | ]
*   Multiply.                   [ ((M-N)*(M-N-1)*(3*(M+N)^2 + M - N - 2)) | ]
#   Push stack depth.           [ ((M-N)*(M-N-1)*(3*(M+N)^2 + M - N - 2)) 1 | ]
2   Multiply by 10, add 2.      [ ((M-N)*(M-N-1)*(3*(M+N)^2 + M - N - 2)) 12 | ]
/   Divide.                     [ ((M-N)*(M-N-1)*(3*(M+N)^2 + M - N - 2)/12) | ]
!   Print.                      [ | ]

El puntero de instrucciones llega a un callejón sin salida y tiene que darse la vuelta. Cuando ahora se encuentra /, intenta una división por cero (ya que la parte inferior de la pila está implícitamente llena de ceros), lo que termina el programa.


4

Haskell, 34 bytes

a#b=sum[a..b]^2-sum(map(^2)[a..b])

Ejemplo de uso: 91 # 123->12087152.

Nada que explicar


3

Matlab, 30 29 28 bytes

Usar la idea de Suever normnos da 2 bytes menos

@(x,y)sum(x:y)^2-norm(x:y)^2

Versión antigua (simple):

@(x,y)sum(x:y)^2-sum((x:y).^2)

3

Octava, 27 23 bytes

@(x,y)sum(z=x:y)^2-z*z'

Crea una función anónima llamada ansque acepta dos entradas:ans(lower, upper)

Demo en línea

Explicación

Crea un vector de fila de xa y(inclusive) y lo almacena en z. Luego sumamos todos los elementos usando sumy lo elevamos al cuadrado ( ^2). Para calcular la suma de los cuadrados, realizamos una multiplicación matricial entre el vector fila y su transposición. Esto cuadrará efectivamente cada elemento y resumirá el resultado. Luego restamos los dos.


3

Java, 84 77 caracteres, 84 77 bytes

7 bytes más pequeños debido a Martin Ender y FryAmTheEggMan, gracias.

public int a(int b,int c){int e=0,f=0;for(;b<=c;e+=b,f+=b*b++);return e*e-f;}

Usando los tres casos de prueba en la publicación original: http://ideone.com/q9MZSZ

Sin golf:

public int g(int b, int c) {
    int e = 0, f = 0;
    for (; b <= c; e += b, f += b * b++);
    return e*e-f;
}

El proceso se explica por sí mismo. Declaré dos variables para representar el cuadrado de las sumas y la suma de los cuadrados y las incrementé repetidamente de manera apropiada. Finalmente, devuelvo la diferencia calculada.


Bienvenido a PPCG! Es probable que pueda guardar un byte que al poner ++en f+=b*b++(lo que puede dejar la tercera ranura del forvacío) y también no necesita de la plaza eantes de devolverlo (es decir, sólo lo hacen return e*e-f).
Martin Ender

Actually instead of leaving the third slot of the for empty, move the f+=b*b++ in there, so you can save on both a semicolon and the braces.
Martin Ender

Great catch @MartinEnder, thank you :)
Mario Ishac

También basado en lo que Martin tenía en mente, esto parece ser un poco más corto.
FryAmTheEggman

1
Aparentemente, mi último comentario fue incorrecto. En realidad, es una parte especial de la gramática de Java: la declaración final de un for es en realidad un tipo especial de declaración, que se denomina lista de expresiones de declaración. Esta declaración especial puede tener más de una declaración unida por una coma. Consulte 14.14.1 (tendrá que navegar allí usted mismo, no pude encontrar una manera de hacer un enlace más preciso) de la especificación del idioma.
FryAmTheEggman


3

JavaScript (ES6), 50 37 bytes

f=(n,m,s=0)=>n>m?0:2*n*s+f(n+1,m,n+s)

Ahora es un puerto de la solución Python de @ Dennis ♦.


Intente usarn=>m=>eval(`for(s=t=0;n<=m;t+=n++)s+=n*n;t*t-s`)
Mama Fun Roll

@MamaFunRoll Por otro lado, podría intentar portar la solución Python de Dennis ♦ ...
Neil

3

Factor, 48 bytes

[ [a,b] [ [ sq ] map sum ] [ sum sq ] bi - abs ]

Una función anónima.

[ 
  [a,b] ! a range from a to b 
  [ 
    [ sq ] map sum ! anonymous function: map sq over the range and sum the result 
  ] 
  [ sum sq ] ! the same thing, in reverse order
  bi - abs   ! apply both anon funcs to the range, subtract them and abs the result
]

3

Haskell, 36 bytes

m#n=sum[2*i*j|i<-[m..n],j<-[i+1..n]]

λ> m # n = sum [ 2*i*j | i <- [m..n], j <- [i+1..n] ]
λ> 5 # 9
970
λ> 91 # 123
12087152
λ> 1 # 10
2640

Tenga en cuenta que

(k=mnk)2k=mnk2==k1=mnk2=mk2k1nk1k2=k1=mnk2=k1+1n2k1k2

1
You don't need the parens around i+1.
Wheat Wizard

2
Also if you want to talk Haskell and Haskell golfing you can join us in the chat room.
Wheat Wizard

3

Perl 6,  36 32  31 bytes

{([+] $_=@_[0]..@_[1])²-[+] $_»²}
{([+] $_=$^a..$^b)²-[+] $_»²}
{[+]($_=$^a..$^b)²-[+] $_»²}

Test it

Explanation:

{ # bare block with placeholder parameters $a and $b

  [+](# reduce with &infix:<+>
      # create a range, and store it in $_
      $_ = $^a .. $^b
  
  -
  [+] # reduce with &infix:<+>
    # square each element of $_ ( possibly in parallel )
    $_»²
}

Test:

#! /usr/bin/env perl6
use v6.c;
use Test;

my @tests = (
  (5,9) => 970,
  (91,123) => 12087152,
  (1,10) => 2640,
);

plan +@tests;

my &diff-sq-of-sum = {[+]($_=$^a..$^b)²-[+] $_»²}

for @tests -> $_ ( :key(@input), :value($expected) ) {
  is diff-sq-of-sum(|@input), $expected, .gist
}
1..3
ok 1 - (5 9) => 970
ok 2 - (91 123) => 12087152
ok 3 - (1 10) => 2640

1
Save a byte moving the assignment and evading parens: {$_=$^a..$^b;.sum²-[+] $_»²}
Phil H

1
25 bytes: {.sum²-[+] $_»²}o&[..]
nwellnhof

2

Brachylog, 24 bytes

:efL:{:2^.}a+S,L+:2^:S-.

Expects the 2 numbers in Input as a list, e.g. [91:123].

Explanation

:efL                     Find the list L of all integers in the range given in Input
    :{:2^.}a             Apply squaring to each element of that list
            +S,          Unify S with the sum of the elements of that list
               L+:2^     Sum the elements of L, then square the result
                    :S-. Unify the Output with that number minus S

2

APL, 23 20 bytes

-/+/¨2*⍨{(+/⍵)⍵}⎕..⎕

Funciona en NARS2000.


2

MATL, 11 bytes

&:ts2^w2^s-

Pruébalo en línea!

Explicación:

&:           #Create a range from the input
  t          #Duplicate it
   s2^       #Sum it and square it
      w      #swap the two ranges
       2^s   #Square it and sum it
          -  #Take the difference

2

Pyth, 11 bytes

s*M-F#^}FQ2

Pruébalo en línea!

s*M-F#^}FQ2
       }FQ    Compute the range
      ^   2   Generate all pairs
   -F#        Remove those pairs who have identical elements
 *M           Product of all pairs
s             Sum.

Buen uso del filtro. Aunque ya hay una función s*M.P}FQ2
incorporada


1

CJam, 17 bytes

q~),>_:+2#\2f#:+-

Pruébalo aquí.

Explicación

q~       e# Read and evaluate input, dumping M and N on the stack.
),       e# Increment, create range [0 1 ... N].
>        e# Discard first M elements, yielding [M M+1 ... N].
_        e# Duplicate.
:+2#     e# Sum and square.
\2f#:+   e# Swap with other copy. Square and sum.
-        e# Subtract.

Alternativamente, uno puede sumar los productos de todos los pares distintos (básicamente multiplicando el cuadrado de la suma y eliminando los cuadrados), pero eso es un byte más largo:

q~),>2m*{)-},::*:+

1

PowerShell v2 +, 47 bytes

Dos variaciones

param($n,$m)$n..$m|%{$o+=$_;$p+=$_*$_};$o*$o-$p

$args-join'..'|iex|%{$o+=$_;$p+=$_*$_};$o*$o-$p

En ambos casos, estamos generando un rango con el ..operador, canalizándolo a un bucle |%{...}. Cada iteración, estamos acumulando $oy $pcomo la suma o la suma de los cuadrados. Luego calculamos el cuadrado de las sumas con $o*$oy restamos $p. La salida se deja en la tubería y la impresión es implícita.


1

JavaScript (ES6), 67 bytes

a=>b=>([s=q=0,...Array(b-a)].map((_,i)=>q+=(s+=(n=i+a),n*n)),s*s-q)

Banco de pruebas

f=a=>b=>([s=q=0,...Array(b-a)].map((_,i)=>q+=(s+=(n=i+a),n*n)),s*s-q)
e=s=>`${s} => ${eval(s[0])}` // template tag format for tests
console.log(e`f(5)(9)`)
console.log(e`f(91)(123)`)
console.log(e`f(1)(10)`)


1

J, 29 bytes

Port of Doorknob's Jelly respuesta .

[:(+/@(^&2)-~2^~+/)[}.[:i.1+]

Uso

>> f = [:(+/@(^&2)-~2^~+/)[}.[:i.1+]
>> 91 f 123x
<< 12087152

Donde >>está STDIN, <<es STDOUT y xes para una precisión extendida.


1

Pyke, 11 bytes

h1:Ds]MXXs-

Pruébalo aquí!

h1:         - inclusive_range(input)
   Ds]      -     [^, sum(^)]
      MX    -    deep_map(^, <--**2)
         s  -   ^[1] = sum(^[1])
          - -  ^[0]-^[1]

1

Julia, 25 bytes

f(a,b,x=a:b)=sum(x)^2-x'x

This is a function that accepts two integers and returns a 1x1 integer array.

The approach is simple: Construct a UnitRange from the endpoints a and b and call it x, then sum x, square it, and subtract its norm, which is computed as transpose(x) * x.

Pruébalo en línea! (incluye todos los casos de prueba)


1
a\b=-(x=a:b)'x+sum(x)^2Guarda algunos bytes.
Dennis

1

TI-BASIC, 19 bytes

Prompt N,M
randIntNoRep(N,M
sum(Ans)2-sum(Ans2

randIntNoRepobtiene el rango (barajado). El resto es bastante autoexplicativo.


1

Fith , 52 bytes

{ 1 + range dup sum 2 pow swap { 2 pow } map sum - }

Esta es una función anónima que toma los dos números en la pila y deja un solo número.

Explicación:

{
    1 + range dup      2 ranges from a to b inclusive
    sum 2 pow          Sum one and square it
    swap               Bring a fresh range to the top
    { 2 pow } map sum  Square every element and sum the list
    -                  Subtract
}

1
If you like postfix, point-free and stack-based functional prorgamming you might like Factor :D
cat

1

GeoGebra, 91 bytes

a(x)=(x²+x)/2
b(x)=x³/3+x²/2+x/6
c(x,y)=(a(y)-a(x))²
d(x,y)=b(y)-b(x)
c(x-1,y)-d(x-1,y)

Define una función (probablemente e(x,y)) que calcula la diferencia deseada.
a(x)calcula la suma de números naturales entre 0y x.
b(x)calcula la suma de los cuadrados de los números naturales entre 0y x.
c(x,y)primero calcula la suma de los números naturales entre xy y, luego cuadra esa suma.
d(x,y)calcula la suma de cuadrados entre b(x)y b(y).
La última línea define una función multivariable que finaliza el cálculo. A la función se le asigna automáticamente un nombre, guardando unos pocos bytes.


Hi, how do I call the function that this defines? I was able to figure out the input at geogebra.org/classic#cas , but couldn't figure out how to find or call the final function.
sundar - Reinstate Monica

@sundar: La última línea es una expresión en x e y. Podríamos e(x,y)=pretender darle un nombre, pero para guardar bytes, no lo hacemos aquí. GeoGebra asigna automáticamente un nombre a la expresión (probablemente e, ya que esa es la siguiente letra disponible). No tengo el entorno disponible en este momento, pero no usaría el panel CAS. El panel de álgebra y la barra de entrada deberían hacer el trabajo correctamente. (Ha pasado un tiempo desde que usé GGb en línea; mi imagen mental puede estar desactualizada.)
Joe
Al usar nuestro sitio, usted reconoce que ha leído y comprende nuestra Política de Cookies y Política de Privacidad.
Licensed under cc by-sa 3.0 with attribution required.