poj 2115 C Looooops

扩展的欧几里得算法……

链接http://poj.org/problem?id=2115

#include<iostream>
#include<stdio.h>
#include<algorithm>
#include<vector>
#include<cmath>
#include<iomanip>
#include<string>
using namespace std;
__int64 extend_gcd(__int64 a,__int64 b,__int64 &x,__int64 &y)
{
    __int64 d;
    if(b==0)
    {
        x=1;
        y=0;
        return a;
    }
    else
    {
        d=extend_gcd(b,a%b,x,y);
        __int64 temp=x;
        x=y;
        y=temp-a/b*y;
    }
    return d;
}
int main()
{
    __int64 x,y,ans,m,n,l,i,a,b,k,c;
    scanf("%I64d%I64d%I64d%I64d%I64d",&x,&y,&m,&n,&l);
        k=n-m;c=x-y;
        if(k<0)
        {
            k=-k;
            c=-c;
        }
        i=extend_gcd(k,l,a,b);
        if((x-y)%i)
            cout<<"Impossible
"<<endl;
        else 
        {
            ans=(c)/i*a%l+l;
            ans=(ans)%(l/i);
            printf("%I64d
",ans);
        }
    return 0;
}