POJ1061 青蛙的约会 数论

                                            青蛙的约会
                            Time Limit: 1000MS      Memory Limit: 10000K
                            Total Submissions: 99911        Accepted: 19128

Description

两只青蛙在网上相识了，它们聊得很开心，于是觉得很有必要见一面。它们很高兴地发现它们住在同一条纬度线上，于是它们约定各自朝西跳，直到碰面为止。可是它们出发之前忘记了一件很重要的事情，既没有问清楚对方的特征，也没有约定见面的具体位置。不过青蛙们都是很乐观的，它们觉得只要一直朝着某个方向跳下去，总能碰到对方的。但是除非这两只青蛙在同一时间跳到同一点上，不然是永远都不可能碰面的。为了帮助这两只乐观的青蛙，你被要求写一个程序来判断这两只青蛙是否能够碰面，会在什么时候碰面。
我们把这两只青蛙分别叫做青蛙A和青蛙B，并且规定纬度线上东经0度处为原点，由东往西为正方向，单位长度1米，这样我们就得到了一条首尾相接的数轴。设青蛙A的出发点坐标是x，青蛙B的出发点坐标是y。青蛙A一次能跳m米，青蛙B一次能跳n米，两只青蛙跳一次所花费的时间相同。纬度线总长L米。现在要你求出它们跳了几次以后才会碰面。

Input

输入只包括一行5个整数x，y，m，n，L，其中x≠y < 2000000000，0 < m、n < 2000000000，0 < L < 2100000000。

Output

输出碰面所需要的跳跃次数，如果永远不可能碰面则输出一行”Impossible”

Sample Input

1 2 3 4 5

Sample Output

4

Source

浙江

题解:

自古浙江出神题……
这题是一道很显然的数论，读完题就应该能建立出如下公式：

sx + am = sy + an ( mod L) （a是要求的答案）
转化一下：a(n-m) = x-y (mod L)

接下来我们就进入漫长的求解过程：
通过阅读题目我们发现这题要求的是上述同余方程的最小非负整数解，那么无解情况就是方程无解的情况，下面我们来了解几个定理：

定理1： 在 gcd(a,b) == 1的情况下 ax = b(mod n) 在[0,n-1]上有唯一的解

证明：
设存在两个解x1,x2，则有：ax1 = b (mod n),ax2 = b (mod n)；
两式想减，得到a(x1-x2) = 0 (mod n)；
因为gcd(a,n) = 1，所以只有可能是(x1-x2)%n = 0；
因为x1,x2都属于[0,n-1]，所以他们的差值（取大于零）也属于[0,n-1]；
如果在[0,n-1]范围内一个数mod n为0，那么这个数只可能是0。（0%n = 0）；
则 x1 - x2 == 0 —–> x1 == x2；
这样就证明了解得唯一性~

定理2：在gcd(a,b) == d的情况下 ax=b(mod n) 在[0,n/d-1]上有唯一解

证明：
因为gcd(a,b) == d；
所以gcd(a/d,b/d) == 1；
将ax = b (mod n)展成ax + by = n;
则根据等式的基本性质 ： a/d*x + b/d*y = n/d;
即 a/d * x = b/d (mod n/d);
这两个式子是等价的。
得到下面的式子后我们发现gcd(a/d,n/d) == 1，满足定理1，即该式在，[0,b/d-1]上有唯一解。

呼~总算证明完毕了QAQ
那样我们就可以用扩展欧几里得定理求出一组解，然后通过上述定理的变换求出最小解咯~
我直接套的我以前整理的模板~

ll exgcd(ll a,ll b,ll &x,ll &y)
{
    if(a == 0 && b == 0)    return -1;
    if(b == 0){x = 1;y = 0;return a;}
    ll d = exgcd(b,a%b,x,y),t = x;
    x = y,y = t-a/b*y;  return d;
}

ll get(ll a,ll b,ll n)  // ax = b (mod n)
{
    ll x,y,d = exgcd(a,n,x,y);
    if(b % d)   return -1;  //b不能整除a和n的最大公约数 
    x = (x % n + n) % n;
    return (x*(b/d))%(n/d);//x*b/d仍然是一组解，上面证明时说过，然后因为唯一解再取模一下，让答案保证在[0,b/d-1]的范围内~
}

当然还得注意一下细节：负数全部直接转化成取模意义下的正数就可以
全部代码~

#include<iostream>
#include<cstring>
#include<cstdio>
#include<algorithm>
#include<queue>
#include<stack>
using namespace std;
typedef long long ll;

ll exgcd(ll a,ll b,ll &x,ll &y)
{
    if(a == 0 && b == 0)    return -1;
    if(b == 0){x = 1;y = 0;return a;}
    ll d = exgcd(b,a%b,x,y),t = x;
    x = y,y = t-a/b*y;  return d;
}

ll get(ll a,ll b,ll n)  // ax = b (mod n)
{
    ll x,y,d = exgcd(a,n,x,y);
    if(b % d)   return -1;  //b不能整除a和n的最大公约数 
    //上面那东西写成d > 1也行，d > 1时无逆元
    x = (x % n + n) % n;
    return (x*(b/d))%(n/d);
}

ll sx,sy,m,n,l;

int main()
{
    scanf("%lld%lld%lld%lld%lld",&sx,&sy,&m,&n,&l);
    ll qaq = sx - sy;
    if(qaq < 0 )    qaq = (qaq + l) % l;
    ll qwq = n - m;
    if(qwq < 0)     qwq = (qwq + l) % l;
    ll ans = get(qwq,qaq,l);
    if(ans == -1)   puts("Impossible");
    else    printf("%lld",ans);
 //   system("pause");
    return 0;
}
// 5 3

pause叫的时候一定要去掉=-= NOIP要是交上pause就是1=和3=的区别了，一定要注意！