// Magma program for generating a b-file for OEIS A367566
// Jon E. Schoenfield, Nov 24 2023

for n in [1..10000] do
   if IsEven(n) then
      prod:=1;
   else
      prod:=2;
   end if;
   if (n gt 1) and (n mod 3 ne 0) then
      prod*:=3;
   end if;
   P:=[];
   for j in [-1,1] do
      t:=n+j;
      if t gt 0 then
         F:=Factorization(t);
         for f in F do
            p:=f[1];
            if p gt 3 then
               B:=p*(p-1) div 2;
               if n mod B in [p-1,B-(p-1)] then
                  prod*:=p;
               end if;
            end if;
         end for ;
      end if;
   end for;
   n, prod;
end for;