zoukankan      html  css  js  c++  java
  • POJ 1061 青蛙的约会(拓展欧几里得算法求解模线性方程组详解)

    题目链接:

    BZOJ:

    https://www.lydsy.com/JudgeOnline/problem.php?id=1477

    POJ:

    https://cn.vjudge.net/problem/POJ-1061

    题目描述:

    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

    关于如何使用拓展欧几里得算法求解模线性方程组的证明请参考另一篇博客:

     1 /*
     2 题意描述
     3 计算满足(x+km)(mod l)=(y+kn)(mod l) 的k的最小正整数
     4 
     5 解题思路
     6 使用拓展欧几里得算法求解模线性方程组
     7 由题知,(x+km)≡(y+kn)(mod l) ,由它的充要条件可得
     8          (x+km)-(y+kn)=tl(其中t属于整数),整理可得
     9          k(m-n)-tl=y-x
    10          另a=m-n,b=-l,c=y-x
    11          可得ak+bl=c,即二元一次方程,利用拓展欧几里得算法求得
    12          a和b的最大公约数d以及满足方程的一组解(k0,t0)
    13          进而判断是否有整数解,有整数解后求出最小正整数解 
    14 */
    15 #include<cstdio>
    16 typedef long long LL;
    17 
    18 void extgcd(LL a, LL b, LL &d, LL &k0, LL &t0);
    19 
    20 int main()
    21 {
    22     LL x,y,m,n,l,a,b,c,d,k0,t0,k,t;
    23     while(scanf("%lld%lld%lld%lld%lld",&x,&y,&m,&n,&l) != EOF){
    24         a=m-n;
    25         b=l;
    26         c=y-x;
    27         if(a < 0){
    28             a = -a;
    29             c = -c;
    30         }
    31         extgcd(a,b,d,k0,t0);
    32         if(c%d != 0)
    33             printf("Impossible
    ");
    34         else
    35         {
    36             k=k0*c/d;
    37             l=l/d;
    38             if(k >= 0)
    39                 k %= l;
    40             else
    41                 k = k%l + l;
    42             printf("%lld
    ",k);
    43         }
    44     }    
    45     return 0;
    46 }
    47 
    48 void extgcd(LL a, LL b, LL &d, LL &k0, LL &t0)
    49 {
    50     if(b == 0){
    51         d=a;k0=1;t0=0;
    52     }
    53     else{
    54         extgcd(b,a%b,d,t0,k0);
    55         t0 -= k0*(a/b);    
    56     }
    57 }

     

  • 相关阅读:
    关于虚拟机链接本地磁盘文件的问题
    javaScript学习笔记
    html学习笔记
    eclipse svn插件安装
    python学习笔记一

    hive数据处理
    WordCount实验
    暑假第六周总结
    暑假第五周总结
  • 原文地址:https://www.cnblogs.com/wenzhixin/p/9343420.html
Copyright © 2011-2022 走看看