Список, найти значение по формуле. Ошибка деления на нуль

{Программа решает следующую задачу:
Дан список из n целых неотрицательных чисел a1,a2,...,a n.
Найти: (a1+a2)div a n  +  (a2+a3)div a n-1  + ... +  (a n-1+a n)div a1}
 
program Teach;
 
  type dmas = ^rmas;
       rmas = record
         n    :integer;
         next :dmas;
         prev :dmas;
       end;
 
  var p1,pt,p2, pfirst, plast :dmas;
      e: integer;
      k: char;
      f: boolean;
 
{Создание массива}
{====================}
      procedure crmas;
          var i:integer;
        begin
          writeln;
          e:=random(20);
 
          for i:=1 to e do
            begin
              new(p2);
              p2^.n:=random(100)+1;
              p2^.next:=p1^.next;
              p2^.prev:=p1^.next^.prev;
              p1^.next^.prev:=p2;
              p1^.next:=p2;
            end;
        end;
{====================}
{Вывод массива}
{====================}
      procedure printlprya;      //прямой
        begin
          writeln;
          p1:=pfirst^.next;
          while p1<>pfirst do
             begin
              write(p1^.n:4);
              p1:=p1^.next;
             end;
          writeln;
        end;
{====================}
      procedure printlobra;      //обратный
        begin
          writeln;
          p2:=p1^.prev;
          while p2<>p1 do
            begin
              write(p2^.n:4);
              p2:=p2^.prev;
            end;
          writeln;
        end;
{====================}
{Условная сумма}
{====================}
       function usum: integer;      //обратный
       var ts,s:integer;
        begin
          ts:=0;
          p1:=pfirst^.next;
          p2:=p1^.prev;
          while (p1<>pfirst) do
             begin
              s:=p1^.n;
              p1:=p1^.next;
              s:=s+p1^.n;
              p1:=p1^.next;
 
              ts:=ts+(s div p2^.n);
              p2:=p2^.prev;
             end;
          usum:=ts;
        end;
{====================}
 
 
begin
       randomize;
       f:=true;
       new(p1);
       p1^.next:=p1;
       p1^.prev:=p1;
       pfirst:=p1;
 
       crmas;
       printlprya;
 
       writeln('Условная сумма:',usum:6);
       readln;
end.

       function usum: integer;      //обратный
       var ts,s:integer;
        begin
          ts:=0;
          p1:=pfirst^.next;
          p2:=p1^.prev^.prev; // <--- Раз
          while (p1<>pfirst) do
             begin
              s:=p1^.n;
              p1:=p1^.next;
              s:=s+p1^.n;
              // p1:=p1^.next; // <--- Два
               
              ts:=ts+(s div p2^.n);
              p2:=p2^.prev;
             end;
          usum:=ts;
        end;