博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
UVA 11768 - Lattice Point or Not(数论)
阅读量:4629 次
发布时间:2019-06-09

本文共 2237 字,大约阅读时间需要 7 分钟。

UVA 11768 - Lattice Point or Not

题意:给定两个点,构成一条线段。这些点都是十分位形式的,求落在这个直线上的正数点。

思路:先把直线表达成a x + b y = c的形式,a,b, c都化为整数表示。然后利用扩展gcd求出x和y的通解,然后已知min(x1, x2) <= x <= max(x1, x2), min(y1, y2) <= y <= max(y1, y2)。这样一来就能够求出通解中t的范围,t能取的整数就是整数解。就得到答案。

值得注意的是。直线为平行坐标系的情况。要特殊推断一下

代码:

#include 
#include
#include
#include
using namespace std;const long long INF = 0x3f3f3f3f3f3f3f;int t;long long xx1, yy1, xx2, yy2;long long a, b, c;long long read(){ double t; scanf("%lf", &t); return (long long)(10 * (t + 0.05));}long long gcd(long long a, long long b) { if (!b) return a; return gcd(b, a % b);}long long exgcd(long long a, long long b, long long &x, long long &y) { if (!b) {x = 1; y = 0; return a;} long long d = exgcd(b, a % b, y, x); y -= a / b * x; return d;}void build() { a = (yy2 - yy1) * 10; b = (xx1 - xx2) * 10; c = (yy2 - yy1) * xx1 + (xx1 - xx2) * yy1; long long t = gcd(gcd(a, b), c); a /= t; b /= t; c /= t;}long long solve() { long long ans = 0; long long x, y; long long d = exgcd(a, b, x, y); long long up = INF, down = -INF; if (xx1 > xx2) swap(xx1, xx2); if (yy1 > yy2) swap(yy1, yy2); if (c % d) return ans; if (b / d > 0) { down = max(down, (long long)ceil((xx1 * d * 1.0 / 10 - x * c * 1.0) / b)); up = min(up, (long long)floor((xx2 * d * 1.0 / 10 - x * c * 1.0) / b)); } else if (b / d < 0) { up = min(up, (long long)floor((xx1 * d * 1.0 / 10 - x * c * 1.0) / b)); down = max(down, (long long)ceil((xx2 * d * 1.0 / 10 - x * c * 1.0) / b)); } else if (xx1 % 10) return ans; if (a / d > 0) { down = max(down, (long long)ceil((y * c * 1.0 - d * yy2 * 1.0 / 10) / a)); up = min(up, (long long)floor((y * c * 1.0 - d * yy1 * 1.0 / 10) / a)); } else if (a / d < 0) { up = min(up, (long long)floor((y * c * 1.0 - d * yy2 * 1.0 / 10) / a)); down = max(down, (long long)ceil((y * c * 1.0 - d * yy1 * 1.0 / 10) / a)); } else if (yy1 % 10) return ans; if (down <= up) ans += up - down + 1; return ans;}int main() { scanf("%d", &t); while (t--) { xx1 = read(); yy1 = read(); xx2 = read(); yy2 = read(); build(); printf("%lld\n", solve()); } return 0;}

转载于:https://www.cnblogs.com/mengfanrong/p/5164221.html

你可能感兴趣的文章
Android通过ksoap2调用.net(c#)的webservice
查看>>
你真的以为了解java.io吗 呕心沥血 绝对干货 别把我移出首页了
查看>>
【Henu ACM Round#15 A】 A and B and Chess
查看>>
【BZOJ 3879】SvT
查看>>
C#制作安装包
查看>>
NodeJS 模块
查看>>
数据图表之圆柱图
查看>>
SQL 学习笔记
查看>>
在CentOS上把MySQL从5.5升级到5.6
查看>>
SurfaceView 间取得焦点
查看>>
软件工程网络15结对编程作业
查看>>
nginx配置ssl证书
查看>>
Java学习笔记07--日期操作类
查看>>
3.commonjs模块
查看>>
case when then else end
查看>>
利用反射对应数据库字段
查看>>
Linux常用命令
查看>>
结对项目开发电梯调度 - 整体设计
查看>>
CSS3边框背景-边框背景(-border-image)
查看>>
CSS3关于过渡效果的问题
查看>>