[BZOJ2712][[Violet 2]棒球][类欧几里得算法] 类似于下面这道题吧,只要把小数转换成分数就好了。
http://blog.csdn.net/g1n0st/article/details/62044709
代码:
#include
using namespace std;
typedef long long ll;
typedef pair abcd;
inline ll Min(const ll &a, const ll &b) {
return a < b ? a : b;
}
inline ll gcd(const ll &a, const ll &b) {
if (!b) return a;
return gcd(b, a % b);
}
inline abcd solve(ll p1, ll q1, ll p2, ll q2) {
ll l = p1 / q1 + 1;
abcd ret(0,0);
if (l * q2 < p2) return abcd(l, 1);
if (p1 == 0) return abcd(1, q2 / p2 + 1);
if (p1 <= q1 && p2 <= q2) {
ret = solve(q2, p2, q1, p1);
return abcd(ret.second, ret.first);
}
ll t = p1 / q1;
ret = solve(p1 - q1 * t, q1, p2 - q2 * t, q2);
ret.first += ret.second * t;
return ret;
}
int main(void) {
//freopen("in.txt", "r", stdin);
ll p1, q1, p2, q2, r;
int n;
while (~scanf("%d 0.%lld", &n, &r)) {
if (r == 0) {
printf("1\n");
continue;
}
ll x = 10;
while (n--) x *= 10;
p1 = r * 10 - 5, q1 = x, p2 = r * 10 + 5, q2 = x;
ll g = gcd(p1, q1);
p1 /= g, q1 /= g;
g = gcd(p2, q2);
p2 /= g, q2 /= g;
abcd ans = solve(p1, q1, p2, q2);
printf("%lld\n", Min(ans.second, q1));
}
return 0;
}
完。
【类欧几里得算法|[BZOJ2712][[Violet 2]棒球][类欧几里得算法]】By g1n0st
推荐阅读
- 2017/7/16现代诗
- 2017/11/04|2017/11/04 秀峰 日志
- HDU 5528【2015长春现场赛 B】 Count a * b
- 类欧几里得算法|[类欧几里得算法 数论] BZOJ 2987 Earthquake
- [数论] Codeforces 819D R #421 D.Mister B and Astronomers & 516E R #292 E. Drazil and His Happy Friends
- 模板 poj2947 Widget Factory 高斯消元
- 【扩展欧几里得】练习题
- 扩展欧几里得【数论
- 数论|hdu 5322 Hope(分治+NTT)
- bzoj|Bzoj3817:Sum