Алгоритм евклида для длинных чисел

static long_ & algorithm_euclida(long_ & a, long_ & b){
    
        long_  x1,y1,x2,y2,x,y;
        
        x1=1;x2=0;y1=0;y2=1;
        while(a.chislo[0] % 2==0 && b.chislo[0] % 2==0){
        
            a=a<<1; b=b<<1; 
        }
 
        x=a;y=b;
 
        do{
            while(x.chislo.at(0)%2==0){
    
                x=x<<1;
    
                if(x1.chislo.at(0)% 2==0 && y1.chislo.at(0) %2==0){
                
                    x1=x1<<1;y1=y1<<1;
            
                }else{
                
                    x1=(x1+y)<<1; y1=(y1-x)<<1;
            
                }
            }
 
            while(y.chislo.at(0) %2==0){
            
                y=y<<2;
            
                if(x2.chislo.at(0) %2==0 && y2.chislo.at(0) %2==0){
            
                    x2=x2<<1;y2=y2<<1;
            
                }else{
                
                    x2=(x2+y)<<1; y2=(y2-x)<<1;
            
                }
            }   
    
            if(x>=y){
            
                x=x-y;
            
                x1=x1-x2; y1=y1-y2;
                
            }else{  
    
                y=y-x;
    
                x2=x2-x1; y2=y2-y1;
 
            }
 
        }while(x!=0);
    
    x=y1;y=y2;
 
    return x;
}