Fraction Construction Problem

思路

$\frac{c}{d}?\frac{e}{f}=\frac{a}{b}$

$a<b\&f<b$

$1\le 1,e\le 4\times 10^{12}$

• 當(dāng)b = 1時(shí)，一定無解。

• gcd(a, b) != 1，能夠較為簡單地構(gòu)造出來。

• $b=p^k$，我們設(shè)$d=p^{k?x},f=p^x$，左邊方程通分得到$\frac{c{p}^x?e{p}^{k?x}}{p^k}$，則方程$c{p}^x?e{p}^{k?x}=a$有解，那么一定有$gcd(p^x,p^{k?x})\mid a$，

  只能$x=0$,所以$d=p^k=b$不符合要求。

• 所以這時(shí)$b$一定有至少兩個(gè)質(zhì)因，我們不妨設(shè)$d\times f=b$，$d,f$互質(zhì)，這樣我們只要通過拓展歐幾里德去求解即可。

代碼

/*Author : lifehappy */ #pragma GCC optimize(2) #pragma GCC optimize(3) #include <bits/stdc++.h>using namespace std;typedef long long ll;const int inf = 0x3f3f3f3f; const double eps = 1e-7;const int N = 1e4 + 10;int prime[N], cnt;ll c, d, e, f;bool st[N];void init() {for(int i = 2; i < N; i++) {if(!st[i]) {prime[cnt++] = i;}for(int j = 0; j < cnt && i * prime[j] < N; j++) {st[i * prime[j]] = 1;if(i % prime[j] == 0) break;}} }bool get_fac(ll n) {for(int i = 0; prime[i] * prime[i] <= n; i++) {if(n % prime[i] == 0) {d = 1;while(n % prime[i] == 0) {n /= prime[i];d *= prime[i];}if(n == 1) return false;return true;}}return false; }ll exgcd(ll a, ll b, ll & x, ll & y) {if(!b) {x = 1, y = 0;return a;}ll d = exgcd(b, a % b, x, y);ll temp = x;x = y;y = temp - a / b * y;return d; }int main() {// freopen("in.txt", "r", stdin);// freopen("out.txt", "w", stdout);// ios::sync_with_stdio(false), cin.tie(0), cout.tie(0);int T;init();scanf("%d", &T);while(T--) {ll a, b;scanf("%lld %lld", &a, &b);int gcd = __gcd(a, b);if(gcd != 1) {printf("%lld %lld %lld %lld\n", a / gcd + 1, b / gcd, 1, b / gcd);continue;}if(!get_fac(b)) {puts("-1 -1 -1 -1");continue;}f = b / d;if(f < d) swap(f, d);ll x, y;ll g = exgcd(f, d, x, y);while(y > 0) x += d, y -= f;printf("%lld %lld %lld %lld\n", x * a, d, -y * a, f);}return 0; }

總結(jié)

以上是生活随笔為你收集整理的Fraction Construction Problem（拓展欧几里德）的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。