【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=; char c = getchar(); x=;
        for (; !isdigit(c); c = getchar()) { if (c=='-') f=-; }
        for (; isdigit(c); c = getchar()) x=x*+c-'';
        x*=f;
}

inline LL exgcd(LL a,LL b,LL &x,LL &y) {
        LL r;
        if (b == ) {
                x = ;
                y = ;
                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 ;

}