UVA 11768 Lattice Point or Not(扩展欧几里得)

题意:给出两个点,求两个点之间的线段之间的整数点的个数。

思路:这题跟SGU的一个题差不多,就是利用扩展欧几里得来寻找整数点,只不过这个点是小数,刚开始想的麻烦了,搞的都快傻了、、、、细节的地方感觉不是很好写,uva太坑了,Judging Error 了一下午。

代码:

#include<iostream>
#include<cstdio>
#include<cstring>
#include<string>
#include<algorithm>
#include<map>
#include<queue>
#include<stack>
#include<set>
#include<cmath>
#include<vector>
#define inf 0x3f3f3f3f
#define Inf 0x3FFFFFFFFFFFFFFFLL
#define eps 1e-9
#define pi acos(-1.0)
using namespace std;
typedef long long ll;
ll ext_gcd(ll a,ll b,ll &x,ll &y)
{
    if(b==0)
    {
        x=1;y=0;
        return a;
    }
    ll d=ext_gcd(b,a%b,y,x);
    y-=a/b*x;
    return d;
}
int dcmp(double x)
{
    if(fabs(x)<eps) return 0;
    return x<0?-1:1;
}
ll solve(double xx1,double yy1,double xx2,double yy2)
{
    ll x1,y1,x2,y2;
    x1=(xx1)*10;
    y1=(yy1)*10;
    x2=(xx2)*10;
    y2=(yy2)*10;
    if(x1==x2)
    {
        if(x1%10) return 0;
        if(yy2<yy1) swap(yy1,yy2);
        y1=ceil(yy1);y2=floor(yy2);
        return max(y2-y1+1,0LL);
    }
    if(y1==y2)
    {
        if(y1%10) return 0;
        if(xx2<xx1) swap(xx1,xx2);
        x1=ceil(xx1);x2=floor(xx2);
        return max(x2-x1+1,0LL);
    }
    ll a=(y2-y1)*10,b=(x1-x2)*10,c=x1*y2-x2*y1;
    ll x,y,d,k1=0;
    if(yy2<yy1) swap(yy1,yy2);
        y1=ceil(yy1);y2=floor(yy2);
    if(xx2<xx1) swap(xx1,xx2);
        x1=ceil(xx1);x2=floor(xx2);
    if(y1>y2||x1>x2) return 0;
    d=ext_gcd(a,b,x,y);
    if(c%d!=0) return 0;
    a/=d;b/=d;
    x*=(c/d);y*=(c/d);
    if(b<0) b=-b;
    x=x-(x-x1)/b*b;
    x-=b;
    while(x<x1) x+=b;
    k1=(x2-x)/b;
    while(x+k1*b<=x2) k1++;
    return k1;
}
int main()
{
    //freopen("in.txt","r",stdin);
    //freopen("out.txt","w",stdout);
    double x1,y1,x2,y2;
    int t;
    scanf("%d",&t);
    while(t--)
    {
        scanf("%lf%lf%lf%lf",&x1,&y1,&x2,&y2);
        ll ans=solve(x1,y1,x2,y2);
        printf("%lld\n",ans);
    }
    return 0;
}


  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值