poj 1061 + poj 2115 + poj 2142 （exgcd）

poj 1061：

题意：两只青蛙，在地球的同一纬度上，各自相对起点的初始坐标为x，y，现在他们从起点开始向同一个方向开始跳，第一只每次可以跳m，第二只可以跳n。

         给这个圆的长度l，问最小第几次时，两只青蛙能相遇。

        

解析：1. 设当t次跳跃时两娃相遇，则第一只青蛙的坐标为：x + m * t，第二只为：y + n * t 。

              所以当他们相遇时： (x + m * t) - (y + n * t ) = k * l  （0 <= k && k为整数）（因为整个路径是圆，等价与坐标相同，即距离为0）。

              整理得：( n - m ) * t  + k * l = x - y.

              即：          a * x   + b * y = c.    其中 a -> (n - m)   ,   x -> t,     b -> k,       y -> l,          c -> (x - y) .

          2. 用拓展gcd求出当 a * x + b * y = gcd(a, b) 时的一个解与d = gcd（a，b）

          3. 当c是d的倍数时，原方程的一组解为(x*c/d, y*c/d). 当c不是d的倍数时无解

          4. 方程ax≡b (mod n)的最小整数解为：( x % r + r ) % r.   r = b / d.

代码：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <stack>
#include <vector>
#include <queue>
#include <map>
#include <climits>
#include <cassert>
#define LL long long

using namespace std;
const int maxn = 1e6;
const int inf = 0x3f3f3f3f;
const double eps = 1e-8;
const double pi = 4 * atan(1.0);
const double ee = exp(1.0);

void exgcd(LL a, LL b, LL& d, LL& x, LL& y)//ax + by = d, d = gcd(a, b)
{
    if (b == 0)
    {
        d = a;
        x = 1;
        y = 0;
    }
    else
    {
        exgcd(b, a % b, d, y, x);
        y -= x * (a / b);
    }
}

int main()
{
    #ifdef LOCAL
    freopen("in.txt", "r", stdin);
    #endif // LOCAL
    LL x, y, m, n, l;
    while (scanf("%lld%lld%lld%lld%lld", &x, &y, &m, &n, &l) !=  EOF)
    {
        LL a = n - m;
        LL b = l;
        LL c = x - y;
        LL X, Y, D;
        exgcd(a, b, D, X, Y);
        if (c % D != 0)
        {
            printf("Impossible\n");
        }
        else
        {
            LL r = b / D;
            X = c / D * X;
            X = (X % r + r) % r;
            printf("%lld\n", X);
        }
    }
    return 0;
}

poj 2115 解法同上，不同的就是几个参数。

poj 2142：

题意：

给a， b， c， 求：

ax + by = c

的解x，y，要求x + y最小以及|ax| + |by|最小。

解析：

因为题目保证有解，所以先将整个方程两边同时除以（a，b），这样可以保证求解出的个数x和y是最小的。

在保证秤砣的个数时最小之后，只要保证和最小就行了，和最小就在两种方案间产生，一种是最小的a×x，一种时最小的b×y。

见代码。

代码：

#include <iostream>
#include <cstdio>
#include <cstdlib>
#include <algorithm>
#include <cstring>
#include <cmath>
#include <stack>
#include <vector>
#include <queue>
#include <map>
#include <climits>
#include <cassert>
#define LL long long

using namespace std;
const int inf = 0x3f3f3f3f;
const double eps = 1e-8;
const double pi = 4 * atan(1.0);
const double ee = exp(1.0);

const int maxn = 30 + 10;

LL gcd(LL a, LL b)
{
    return b ? gcd(b, a % b) : a;
}

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

int main()
{
#ifdef LOCAL
    freopen("in.txt", "r", stdin);
#endif // LOCAL
    LL a, b, c;
    while (scanf("%lld%lld%lld", &a, &b, &c) && a && b && c)
    {
        ///保证个数最小
        LL d = gcd(a, b);
        a /= d;
        b /= d;
        c /= d;

        LL x, y;
        exgcd(a, b, d, x, y);

        LL x1, y1, x2, y2;
        ///计算两种情况
        ///x 最小时
        x1 = x * c;///   x' * c / d
        x1 = (x1 % b + b) % b;
        y1 = (c - a * x1) / b;
        ///y 最小时
        y2 = y * c;
        y2 = (y2 % a + a) % a;
        x2 = (c - b * y2) / a;

        if (abs(x1 + 0.0) + abs(y1 + 0.0) < abs(x2 + 0.0) + abs(y2 + 0.0))
            x = abs(x1 + 0.0), y = abs(y1 + 0.0);
        else
            x = abs(x2 + 0.0), y = abs(y2 + 0.0);
        printf("%lld %lld\n", x, y);
    }
    return 0;
}