【POJ 1061】 青蛙的约会

【题目链接】

           点击打开链接

【算法】

         列出同余方程，然后用exgcd求解

【代码】

          

#include <algorithm>
#include <bitset>
#include <cctype>
#include <cerrno>
#include <clocale>
#include <cmath>
#include <complex>
#include <cstdio>
#include <cstdlib>
#include <cstring>
#include <ctime>
#include <deque>
#include <exception>
#include <fstream>
#include <functional>
#include <limits>
#include <list>
#include <map>
#include <iomanip>
#include <ios>
#include <iosfwd>
#include <iostream>
#include <istream>
#include <ostream>
#include <queue>
#include <set>
#include <sstream>
#include <stdexcept>
#include <streambuf>
#include <string>
#include <utility>
#include <vector>
#include <cwchar>
#include <cwctype>
#include <stack>
#include <limits.h>
using namespace std;

typedef long long LL;

LL p,q,x,y,n,m,L,r;

template <typename T> void read(T &x) {
        LL f=1; char c = getchar(); x=0;
        for (; !isdigit(c); c = getchar()) { if (c=='-') f=-1; }
        for (; isdigit(c); c = getchar()) x=x*10+c-'0';
        x*=f;
}

inline LL exgcd(LL a,LL b,LL &x,LL &y) {
        LL r;
        if (b == 0) {
                x = 1;
                y = 0;
                return a;
        }else {
                r = exgcd(b,a%b,y,x);
                y -= a / b * x;
                return r;
        }
}

int main() {
        

        read(x); read(y); read(m); read(n); read(L);
        
        r = exgcd(L,n-m,p,q);
        
        if ((x - y) % r) puts("Impossible");    
        else cout<< (((q * (x - y) / r) % L) + L) % L << endl;
        
        return 0;
    
}