luogu P1516 青蛙的约会

最近想刷一刷数论的题，就先找了道简单的入手。

 

这一看就是exgcd求不定方程。

首先可以列出来：x + m * k ≡ y + n * k (mod L)

　　　　　化简：x + m * k - y - n * k = p * L

　　　　　　　　(m - n) * k - p * L = y - x

令a = m - n, b = L, c = y - x,于是

　　　　　　　　a * k - b * p = c

首先要保证a,b是正的。

代入exgcd求得ak + bp = 1的解k', p'，在gcd(a, b) | c的情况下，得出解k = k' * c / gcd(a, b)。

最小整数解就是(k + m) % m。

 1 #include<cstdio>
 2 #include<iostream>
 3 #include<cmath>
 4 #include<algorithm>
 5 #include<cstring>
 6 #include<cstdlib>
 7 #include<cctype>
 8 #include<vector>
 9 #include<stack>
10 #include<queue>
11 using namespace std;
12 #define enter puts("") 
13 #define space putchar(' ')
14 #define Mem(a, x) memset(a, x, sizeof(a))
15 #define rg register
16 typedef long long ll;
17 typedef double db;
18 const int INF = 0x3f3f3f3f;
19 const db eps = 1e-8;
20 const int maxn = 105;
21 inline ll read()
22 {
23     ll ans = 0;
24     char ch = getchar(), last = ' ';
25     while(!isdigit(ch)) {last = ch; ch = getchar();}
26     while(isdigit(ch)) {ans = ans * 10 + ch - '0'; ch = getchar();}
27     if(last == '-') ans = -ans;
28     return ans;
29 }
30 inline void write(ll x)
31 {
32     if(x < 0) x = -x, putchar('-');
33     if(x >= 10) write(x / 10);
34     putchar(x % 10 + '0');
35 }
36 
37 ll n, m, x, y, L;
38 
39 ll gcd(ll a, ll b) 
40 {
41     return b ? gcd(b, a % b) : a;
42 }
43 void exgcd(ll a, ll b, ll &x, ll &y)
44 {
45     if(!b) x = 1, y = 0;
46     else exgcd(b, a % b, y, x), y -= a / b * x;
47 }
48 
49 int main()
50 {
51     x = read(); y = read(); m = read(); n = read(); L = read();
52     if(m - n < 0) swap(n, m), swap(x, y);
53     ll a = m - n, b = L, c = y - x, k, p;
54     exgcd(a, b, k, p);
55     ll g = gcd(a, b), m = b / g;
56     if(c % g) printf("Impossible\n");
57     else write((k * (c / g) % m + m) % m);
58     return 0;
59 }

View Code

 

转载于:https://www.cnblogs.com/mrclr/p/9751728.html