POJ

来源:互联网 发布:淘宝韩妆正品店铺推荐 编辑:程序博客网 时间:2024/06/05 19:54

题目链接:http://poj.org/problem?id=1061

中文题0.0

思路:扩展欧几里得算法的模板题,通过把等式变形得到a=m-n,b=l,c=y-x,然后如果c%gcd不等于0的话,那显然是不行的,如果等于0,那么a,b,c都同除以gcd,然后x,y还要同除以c/gcd,这时才满足ax`+by`=1的模式,此时x`=x/(c/gcd),那么x=x`*(c/gcd),再保证都是正数就行了。


代码:

#include <cstdio>
#include <cmath>
#include <iostream>
#include <cstring>
#include <algorithm>
#include <queue>
#include <stack>
#include <vector>
#include <map>
#include <numeric>
#include <set>
#include <string>
#include <cctype>
#include <sstream>
#define INF 0x3f3f3f3f
using namespace std;
typedef long long LL;
typedef pair<int, int> P;
const int maxn = 5e5 + 5;

LL extgcd(LL a,LL b,LL &t,LL &k){
    LL d=a;
    if (b!=0){
        d=extgcd(b,a%b,k,t);
        k-=(a/b)*t;
    }
    else {
        t=1;k=0;
    }
    return d;
}
int main () {
    //freopen ("in.txt", "r", stdin);
    LL x,y,m,n,l,t,k;
    cin>>x>>y>>m>>n>>l;
    LL a=m-n,b=l,c=y-x;
    LL g=extgcd(a,b,t,k);
    if (c%g) cout<<"Impossible"<<endl;
    else {
        b/=g;
        c/=g;
        t=t*c;
        if (b<0) b=-b;
        t%=b;
        if (t<=0) t+=b;
        cout<<t<<endl;
    }
    return 0;
}