Singular

LIB "poly.lib";
proc DivAlgo(list L) {
   poly p = L[1];
   list a = list(0:size(L));
   L = delete(L, 1);
   while(p != 0) {
      poly LTp = lead(p);
      int flag = 0;   
      for(int i = 0; i < size(L); i++) {
         if(gcd(LTp, L[i]) != 1) {
            p = p - (LTp/lead(L[i])) * L[i];
            a[i] = a[i] + LTp/lead(L[i]);
            flag = 1;
            break;
         };
      };
      if(flag == 0) {
         p = p - LTp;
         a[size(a)] = a[size(a)] + LTp;
      };
   };
   return(a);
};