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

inline LL exgcd(LL a,LL b,LL &x,LL &y) {
        LL r;
        if (b == 0) {
                x = 1;
                y = 0;
                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 0;
    
}
复制代码

 

posted @   evenbao  阅读(144)  评论(0编辑  收藏  举报
编辑推荐:
· 如何编写易于单元测试的代码
· 10年+ .NET Coder 心语,封装的思维:从隐藏、稳定开始理解其本质意义
· .NET Core 中如何实现缓存的预热?
· 从 HTTP 原因短语缺失研究 HTTP/2 和 HTTP/3 的设计差异
· AI与.NET技术实操系列:向量存储与相似性搜索在 .NET 中的实现
阅读排行:
· 10年+ .NET Coder 心语 ── 封装的思维:从隐藏、稳定开始理解其本质意义
· 地球OL攻略 —— 某应届生求职总结
· 提示词工程——AI应用必不可少的技术
· Open-Sora 2.0 重磅开源!
· 周边上新:园子的第一款马克杯温暖上架
点击右上角即可分享
微信分享提示