【DP】筷子

筷子

源程序名 　　　CHOP.??? (PAS,C,CPP)

可执行文件名   CHOP.EXE

输入文件名　　 CHOP.IN

输出文件名     CHOP.OUT

 

A先生有很多双筷子。确切的说应该是很多根，因为筷子的长度不一，很难判断出哪两根是一双的。这天，A先生家里来了K个客人，A先生留下他们吃晚饭。加上A先生，A夫人和他们的孩子小A，共K+3个人。每人需要用一双筷子。A先生只好清理了一下筷子，共N根，长度为T1,T2,T3,……,TN.现在他想用这些筷子组合成K+3双，使每双的筷子长度差的平方和最小。(怎么不是和最小？？这要去问A先生了，呵呵)

 

输入

输入文件共有两行，第一行为两个用空格隔开的整数，表示N,K(1≤N≤100, 0<K<50)，第二行共有N个用空格隔开的整数，为Ti.每个整数为1～50之间的数。

 

输出

输出文件仅一行。如果凑不齐K+3双，输出-1，否则输出长度差平方和的最小值。

 

样例

CHOP.IN

10 1

1 1 2 3 3 3 4 6 10 20

 

CHOP.OUT

5

 

说明

第一双 1 1

第二双 2 3

第三双 3 3

第四双 4 6

(1-1)^2+(2-3)^2+(3-3)^2+(4-6)^2=5

============================

=======================================

var
  n,k:longint;
  a:array[1..100]of longint;
  f:array[0..100,0..100]of longint;
procedure init;
begin
  assign(input,'chop.in');
  assign(output,'chop.out');
  reset(input); rewrite(output);
end;

procedure terminate;
begin
  close(input); close(output);
  halt;
end;

procedure qsort(s,t:longint);
var
  i,j:longint;
  x:longint;
  tem:longint;
begin
  i:=s; j:=t;
  x:=a[(s+t)shr 1];
  repeat
    while x<a[j] do dec(j);
    while a[i]<x do inc(i);
    if i<=j then
      begin
        tem:=a[i];
        a[i]:=a[j];
        a[j]:=tem;
        inc(i); dec(j);
      end;
  until i>j;
  if i<t then qsort(i,t);
  if s<j then qsort(s,j);
end;

function min(a,b:longint):longint;
begin
  if a>b then exit(b);
  exit(a);
end;

procedure main;
var
  i,j:longint;
begin
  readln(n,k);
  k:=k+3;
  //注意
  if n<2*k then
    begin
      writeln(-1);
      terminate;
    end;
  for i:=1 to n do read(a[i]);
  qsort(1,n);
  k:=k*2;
  fillchar(f,sizeof(f),$7);
  f[0,0]:=0;
  for i:=2 to n do
    for j:=2 to min(k,i) do
      if j and 1=0 then
      //如果是偶数则....
      begin
        f[i,j]:=f[i-2,j];
        //若不用第i-1个则直接传递...
        if f[i-1,j]<f[i,j] then f[i,j]:=f[i-1,j];
        //若不用第i个则直接传递.....
        f[i,j]:=min(f[i,j],f[i-2,j-2]+(a[i]-a[i-1])*(a[i]-a[i-1]));
      end;
  writeln(f[n,k]);
end;

begin
  init;
  main;
  terminate;
end.