составить программу:

Составить программу нахождения среднего арифметического значения всех делителей заданного натурального числа N, кратных 3.
Вывести 0, если нет делителей.
(Само число и единица включаются в число делителей.)

Формат входных данных
Входной файл содержит число N (1≤N≤1000).

Формат выходных данных
Выведите искомую величину с точностью не менее 10^−9.

Примеры:
Вход
6
Выход
4.5000000000

1)
 program numOp;
  var a,b:longint;
 begin
  write('a='); readln(a);
  write('b='); readln(b); writeln;
  writeln('Среднее арифметическое: ',((a+b)/2):0:2);
  writeln('Сумма: ',a+b);
  writeln('Разность: ',a-b);
  writeln('Произведение: ',a*b); readln;
 end.
2)
 program triangleInf;
  var
   kat1,kat2:longint;
   gip,s:real;
  begin
   write('Катет1: '); readln(kat1);
   write('Катет2: '); readln(kat2);
   gip:=sqrt(sqr(kat1)+sqr(kat2)); s:=kat1*kat2/2;
   writeln('Гипотенуза: ',gip:0:2);
   writeln('Площадь: ',s:0:2); readln;
  end.

Var b,c,n,i,r:integer;
mass1:array[1..1000] of integer;
mass2:array[1..1000] of integer;
begin
writeln('Введите кол-во элементов массива');
read(n);
r:=1;
for i:=1 to n do
begin
mass1[i]:=random(19)-10;
if i mod 2=0 then begin
mass2[r]:=mass1[i];
r:=r+1;
end;
end;
for i:=1 to r do
if mass2[r]>0 then c:=c+1 else b:=b+1;
writeln('Исходный массив');
for i:=1 to n do
write(mass1[i]:4);
writeln;
for i:=1 to r do
write(mass2[i]:4);
writeln;
writeln('Кол-во положительных элементов массива равно ',c,' а отрицательных равно ',b);
end.