#include <stdio.h>
#include <stdlib.h>

int GCD(int a,int b)
{
ldiv_t lx;
int c;

if(a==0 || b==0) return 0;
if(a==b) return a;
if(a<0) a= -1*a;
if(b<0) b= -1*b;

  if(a<b)
    {
    c=b;
    b=a;
    a=c;
    }
  else
  ;

for(;;)
  {
  lx =ldiv(a,b);

  if(lx.rem== 1)  /* relative  prim numbers */
     return 0;
  else
    if(lx.rem==0)
      return(b);
    else
      {
       a=b;
       b=lx.rem;
      }
  }



}