:: алгоритмы  и методы :: :: олимпиадные задачи :: :: связь :: :: о сайте ::
Путь: Математика » Комбинаторика и перебор » Минимальное количество монет
  У продавца и покупателя имеется неограниченное кол-во монет достоинством (1,2,5,10,20,50,100,200,500 к примеру). Покупатель купил товар на сумму n. Hужно найти минимальное кол-во монет, которые будут использованы при расплате. Деньги может давать как покупатель, так и продавец



Из переписки Ilia Kantor с Max Alekseyev

AD>>> У продавца и покупателя имеется неограниченное кол-во монет
AD>>> достоинством (1,2,5,10,20,50,100,200,500 к примеру). Покупатель
AD>>> купил товар на сумму n. Hужно найти минимальное кол-во монет,
AD>>> которые будут использованы при расплате. Деньги может давать как
AD>>> покупатель, так и продавец. Какие у олл соображения по поводу этой
AD>>> задачи?

MA>> Эта задача (относящаяся, кстати, к области целочисленного
MA>> программирования) красиво решается с помощью базисов Грёбнера. Я
MA>> приведу ниже свою реализацию этой идеи для случая, когда первая
MA>> монета имеет достоинство 1 (остальные монеты - любые). В общем
MA>> случае базисы Гребнера тоже прекрасно работают, но реализация будет
MA>> немного сложнее.

IK> Все очень интересно, но нельзя ли поподробнее ? Что есть базисы
IK> Гребнера

  Если вкратце, то теория такая:



Рассмотрим кольцо многочленов над некоторым полем. Пусть I некоторый идеал в этом кольце, порожденный многочленами g1, g2, ..., gn. При этом набор G={g1, g2, ..., gn} называется базисом (или системой образующих) идеала I (обозначение I=<G>).

Пусть на мономах задан линейный порядок (term order), удовлетворяющий условиям: * для любого монома a, отличного от 1, a>1; * если a<b, то для любого монома c имеем ac<bc. Примерами таких порядков являются: лексикографический (lex); общей степени, затем лексикографический (deglex); общей степени, затем обратный лексикографический (degrevlex).

У любого многочлена f однозначно определяется старший (в смысле заданного порядка) моном (leading term), который мы будем обозначать lt(f). Можно считать, что все рассматриваемые многочлены нормированы (коэффициент при старшем мономе равен 1). Если lt(f) делится на lt(g), то многочлен f можно редуцировать относительно многочлена g - результатом будет многочлен f-(lt(f)/lt(g))*g.

Многочлен f называется редуцированным относительно G, если lt(f) не делится на lt(g) ни для какого g из G. Любой многочлен f за конечное число редукций можно привести к редуцированному относительно G многочлену. Заметим, что результат такой редукции, вообще говоря, неоднозначен.

Определение Базис G идеала I называется _базисом Грёбнера_ (стандартным базисом) относительно порядка <, если в результате любой редукции элемента f идеала I к редуцированному относительно G многочлену всегда получается 0.

Каждый элемент g из G можно редуцировать относительно набора G-{g}. Результатом будет редуцированный базис, в котором старшие термы любых двух элементов не делятся один на другой.

Доказано, что
1) любой идеал обладает базисом Грёбнера относительно любого порядка;
2) редуцированный базис Грёбнера идеала I относительно порядка < определяется однозначно с точностью до порядка следования многочленов;
2') два идеала I и J равны между собой тогда и только тогда, когда их базисы Грёбнера совпадают.

Для построения базисов Грёбнера Бухбергер придумал алгоритм, который теперь носит его имя.

IK> и каким именно образом они используются ?

Пусть достоинства монет 1=c0, c1, ..., cm-1. Любую сумму вида (k0-km)*c0 + (k1-km+1)*c1 + ... + (km-1-k2m-1)*cm-1, где все ki, i=0..2m-1 неотрицательны, которые мы будем кодировать мономами x0^{k0}x1^{k1}...x2m-1^{k2m-1}.

Пусть I - идеал, порожденный мономами соответствующими нулевой сумме.

Утв 1. Базисом I является набор многочленов x0^{ci}-xi, x1^{ci}xm+i-1, i=0..m-1. Доказательство индукцией по m. [skipped]

Пусть G - базис Грёбнера идеала I относительно deglex. Пусть нам дана сумма s. Редуцируем многочлен x0^s относительно G.

Утв 2. Результатом этой редукции будет многочлен x0^{k0}x1^{k1}...x2m-1^{k2m-1}, причем в каждой паре (ki,km+i), i=0..m-1 не более одного числа отлично от нуля и s = (k0-km)*c0 + (k1-km+1)*c1 + ... + (km-1-k2m-1)*cm-1 и есть искомое минимальное по числу монет представление суммы s. Доказательство от противного. [skipped]

IK> Каков собственно алгоритм?

Алгоритм простой:
1) построить базис Грёбнера указанного идеала относительно deglex;
2) редуцировать многочлен x0^s;
3) интерпретировать результат.

Если нужно решить несколько задач для одного и того же набора номиналов, но разных сумм, то базис Грёбнера нужно построить только один раз (это самая трудоемкая операция), и затем очень быстро проводить редукцию для каждой из сумм.

{$B-}
const m=9;
      c:array[1..m] of integer=(1,2,5,10,20,50,100,200,500);

      m2=2*m;

type  TTerm=array[1..m2] of integer;
      PTermList=^TTermList;
      TTermList=record
        u,v:TTerm;
        next:PTermList;
      end;

var   Basis:TTermList;
      p,q,r,t:PTermList;
      i,j,k:integer;

function cmp(var x,y:TTerm):integer;
var t,i,dx,dy:integer;
begin
  dx:=0; dy:=0; t:=0;
  for i:=1 to m2 do
  begin
    inc(dx,x[i]); inc(dy,y[i]);
    if t=0 then
    begin
      if x[i]>y[i] then t:=1;
      if y[i]>x[i] then t:=-1;
    end;
  end;
  if dx=dy then cmp:=t else if dx>dy then cmp:=1 else cmp:=-1;
end;

function S(var a,b,c:PTermList):boolean;
var x,y:TTerm;
    i,t:integer;
begin
  S:=False;
  for i:=1 to m2 do
  begin
    if a^.u[i]>b^.u[i] then t:=a^.u[i] else t:=b^.u[i];
    x[i]:=t-a^.u[i]+a^.v[i]; y[i]:=t-b^.u[i]+b^.v[i];
    if x[i]>y[i] then t:=y[i] else t:=x[i];
    dec(x[i],t); dec(y[i],t);
  end;
  case cmp(x,y) of
   1: begin c^.u:=x; c^.v:=y; end;
  -1: begin c^.u:=y; c^.v:=x; end;
  else S:=True; 
  end;
end;

function Reduce(var a:PTermList):boolean;
var p:PTermList;
    x:TTerm;
    i:integer;
begin
  p:=Basis.next;
  repeat
  if a<>p then
  begin
    i:=1; while (i<=m2) and (a^.u[i]>=p^.u[i]) do inc(i);
    if i>m2 then
    begin
      if S(a,p,a) then begin Reduce:=true; exit; end;
      p:=@Basis;
    end;
  end;
  p:=p^.next;
  until (p=nil);
  Reduce:=false;
end;

procedure ReduceBasis;
var p,q,t:PTermList;
begin
  p:=@Basis;
  repeat
    q:=p; p:=p^.next;
    while (p<>nil) and Reduce(p) do
    begin
      q^.next:=p^.next; dispose(p); p:=q^.next;
    end;
    if p=nil then break;
    t:=p;
    while (t^.next<>nil) and (cmp(t^.next^.u,p^.u)=1) do t:=t^.next;
    if t<>p then
    begin
      q^.next:=p^.next; p^.next:=t^.next; t^.next:=p;
    end;
    p:=q^.next;
  until p=nil;
end;

begin
  Basis.next:=nil;
  for i:=1 to m do
  begin
    if i>1 then
    begin
      new(p); with p^ do
      begin
        FillChar(u,SizeOf(u),0); u[1]:=c[i];
        FillChar(v,SizeOf(v),0); v[i]:=1;
        next:=Basis.next;
      end;
      Basis.next:=p;
    end;
    new(p); with p^ do
    begin
      FillChar(u,SizeOf(u),0); u[1]:=c[i]; u[m+i]:=1;
      FillChar(v,SizeOf(v),0);
      next:=Basis.next;
    end;
    Basis.next:=p;
  end;

  write('Construct Groebner basis');
  p:=@Basis; new(r);
  repeat
    q:=p^.next;
    while (q<>nil) do
    begin
      if (not S(p,q,r)) and (not Reduce(r)) then
      begin
        t:=@Basis;
        while (t^.next<>nil) and (cmp(t^.next^.u,r^.u)=1) do t:=t^.next;
        r^.next:=t^.next; t^.next:=r; new(r);
        write('.');
        ReduceBasis;
        p:=@Basis; break;
      end;
      q:=q^.next;
    end;
    p:=p^.next;
  until p=nil;
  writeln(' Done');

  with r^ do
  repeat
    FillChar(u,SizeOf(u),0); FillChar(v,SizeOf(v),0);
    write('Amount of money (0 - exit): '); readln(u[1]);
    if u[1]=0 then break;    
    Reduce(r);
    write('Coins: ');
    j:=0;
    for i:=1 to m2 do
    begin
      for k:=1 to u[i] do
        if i<=m then write('+',c[i]) else write('-',c[i-m]);
      inc(j,u[i]);
    end;
    writeln; writeln('Total number: ',j);
  until false; 
  
  dispose(r);
  p:=Basis.next;
  while p<>nil do
  begin
    q:=p; p:=p^.next; dispose(q);  end;end.