SGU 106 The equation

原创 2012年03月23日 16:34:01

终于a了,这题还是比较麻烦的~wa的原因是没有考虑相除后是取它的上底还是下底。。。

#include<iostream>
#include<vector>
#include<algorithm>
#include<cstdio>
#include<queue>
#include<stack>
#include<string>
#include<map>
#include<set>
#include<cmath>
#include<cassert>
#include<cstring>
#include<iomanip>
using namespace std;

#ifdef _WIN32
#define i64 __int64
#define out64 "%I64d\n"
#define in64 "%I64d"
#else
#define i64 long long
#define out64 "%lld\n"
#define in64 "%lld"
#endif

#define FOR(i,a,b)      for( int i = (a) ; i <= (b) ; i ++)
#define FF(i,a)         for( int i = 0 ; i < (a) ; i ++)
#define FFD(i,a)        for( int i = (a)-1 ; i >= 0 ; i --)
#define S64(a)          scanf(in64,&a)
#define SS(a)           scanf("%d",&a)
#define LL(a)           ((a)<<1)
#define RR(a)           (((a)<<1)+1)
#define SZ(a)           ((int)a.size())
#define PP(n,m,a)       puts("---");FF(i,n){FF(j,m)cout << a[i][j] << ' ';puts("");}
#define pb              push_back
#define CL(Q)           while(!Q.empty())Q.pop()
#define MM(name,what)   memset(name,what,sizeof(name))
#define read            freopen("in.txt","r",stdin)
#define write           freopen("out.txt","w",stdout)

const int inf = 0x3f3f3f3f;
const i64 inf64 = 0x3f3f3f3f3f3f3f3fLL;
const double oo = 10e9;
const double eps = 10e-10;
const double pi = acos(-1.0);

i64 gcd(i64 _a, i64 _b)
{
    if (!_a || !_b)
    {
        return max(_a, _b);
    }
    i64 _t;
    while (_t = _a % _b)
    {
        _a = _b;
        _b = _t;
    }
    return _b;
};

i64 ext_gcd (i64 _a, i64 _b, i64 &_x, i64 &_y)
{
    if (!_b)
    {
        _x = 1;
        _y = 0;
        return _a;
    }
    i64 _d = ext_gcd (_b, _a % _b, _x, _y);
    i64 _t = _x;
    _x = _y;
    _y = _t - _a / _b * _y;
    return _d;
}

i64 exmod (i64 _a, i64 _p)
{
    i64 _ans, _y;
    ext_gcd (_a, _p, _ans, _y);
    _ans < 0 ? _ans += _p : 0;
    return _ans;
}

i64 a,b,c,x1,x2,y11,y22;
i64 gab;
i64 x,y;
i64 ans;


i64 findx(i64 q1,i64 q2,i64 w1,i64 w2)
{
    if(q1>q2)
    {
        swap(q1,q2);
    }
    if(w1>w2)
    {
        swap(w1,w2);
    }

    if(q1 > w2 || q2 < w1)
    {
        return 0;
    }
    else if(q1<=w1 && w2<=q2)
    {
        return w2-w1+1;
    }
    else if(q1>=w1 &&q1<=w2)
    {
        if(q2<=w2)
        {
            return q2-q1+1;
        }
        else
        {
            return w2 - q1+1;
        }
    }
    else if(q2>=w1 && q2 <= w2)
    {
        if(q1>=w1)
        {
            return q2 - q1 +1;
        }
        else
        {
            return q2 - w1 + 1;
        }
    }
}


int main()
{
    i64 temp;
    while(cin>>a>>b>>c)
    {
        cin>>x1>>x2;
        cin>>y11>>y22;
        if(a==0 && b==0)
        {
            if(c==0)
            {
                temp = (x2-x1+1)*(y22-y11+1);
                cout<<temp<<endl;
                continue;
            }
            else
            {
                cout<<"0"<<endl;
                continue;
            }
        }
        else if(a==0)
        {
            if(c%b==0)
            {
                temp = -c/b;
                if(temp<=y22 && temp >= y11)
                {
                    temp = x2-x1+1;
                    cout<<temp<<endl;
                    continue;
                }
            }
            else
            {
                cout<<"0"<<endl;
                continue;
            }
        }
        else if(b==0)
        {
            if(c%a==0)
            {
                temp = -c/a;
                if(temp <= x2 && temp >= x1)
                {
                    temp = y22 - y11 +1;
                    cout<<temp<<endl;
                    continue;
                }
            }
            else
            {
                cout<<"0"<<endl;
                continue;
            }
        }
        else
        {
            i64 t1,t2;
            if(c>0)
            {
                c =- c;
                a = -a;
                b = -b;
            }
            c = -c;
            if(a<0)
            {
                a = -a;
                x1=-x1;
                x2=-x2;
                if(x1>x2)
                swap(x1,x2);
            }
            if(b<0)
            {
                b=-b;
                y11=-y11;
                y22=-y22;
                if(y11>y22)
                swap(y11,y22);
            }
            gab = gcd(a,b);
            a /= gab;
            b /= gab;
            if(c%gab!=0)
            {
                cout<<"0"<<endl;
                continue;
            }

            x = exmod(a,b);
            x *= c;
            x %= b;
            y = (c-a*x)/b;
            i64 lix,rix;
            i64 liy,riy;
            lix = (x1-x)/b;
            if(x1-x > 0 && (x1-x)%b!=0 )
            {
                lix ++;
            }
            rix = (x2-x)/b;
            if(x2-x<0 && (x2-x)%b!=0)
            {
                rix--;
            }
            liy = (y-y22)/a;
            if( y - y22 >0 && (y-y22)%a!=0)
            {
                liy++;
            }
            riy = (y-y11)/a;
            if(y-y11<0 && (y-y11)%a != 0 )
            {
                riy--;
            }
            cout<<findx(lix,rix,liy,riy)<<endl;
        }
    }
    return 0;
}





SGU106 The equation[扩展欧几里德算法]

A - The equation Time Limit:250MS     Memory Limit:4096KB     64bit IO Format:%I64d & %I64u Submit ...
  • ControlBear
  • ControlBear
  • 2016年07月24日 18:55
  • 314

SGU 106 The equation【扩展欧几里得】

先放一张搞笑图。。 不知这是哪位大神,我一直wa2,他一直wa9。。。这样搞笑的局面持续了一个晚上。。。最后各wa了10发才A。。。题目链接:http://acm.hust.edu.cn/vju...
  • Yukizzz
  • Yukizzz
  • 2016年04月13日 13:55
  • 224

SGU 106 The equation(拓展欧几里得)

通解形式 然后利用x,y去求出范围,就能得到解的个数 注意特判a和b都为0的情况 代码: #include #include #include #include using namespa...
  • u011217342
  • u011217342
  • 2015年04月21日 21:06
  • 320

SGU 106 The equation(扩展欧几里德)

题目链接: SGU 106 The equation 题意: 给出a,b,c,x1,x2,y1,y2求满足ax + by + c = 0且x1 ...
  • Ramay7
  • Ramay7
  • 2016年05月27日 02:10
  • 196

SGU 106 The equation ----扩展欧几里得

http://hi.baidu.com/lycanlancelot/item/8ad98b4aefeef8a4de2a9f53 方程ax+by+c=0;    问有多少对x,y. 满足 x1...
  • u010489389
  • u010489389
  • 2014年02月23日 16:54
  • 428

SGU 106 The equation

H - The equation Time Limit:250MS     Memory Limit:4096KB     64bit IO Format:%I64d & %I64u Submit...
  • qq_33328072
  • qq_33328072
  • 2016年05月24日 10:12
  • 103

SGU 106 The equation(扩展欧几里德)

没想起来怎么做 参考:https://www.cnblogs.com/zjbztianya/archive/2013/03/12/2956835.html思路还是很简单的,ax+by=c, gcd(...
  • gyhguoge01234
  • gyhguoge01234
  • 2017年11月14日 19:19
  • 106

H - The equation——(SGU 106)

传送门 Password:nefu 题目大意: 给出a,b,c,x1,x2,y1,y2,求满足ax+by+c=0,且x∈[x1,x2],y∈[y1,y2]的整数解个数。解题思路: 首先想到这...
  • qingshui23
  • qingshui23
  • 2016年03月19日 19:20
  • 525

SGU 106 The equation(扩展欧几里得)

Description 求方程ax+by+c=0的满足x1
  • V5ZSQ
  • V5ZSQ
  • 2016年04月14日 09:41
  • 243

SGU - 106 - The equation (扩展欧几里得)

题目传送:106. The equation 106. The equation time limit per test: 0.25 sec.  memory l...
  • u014355480
  • u014355480
  • 2015年04月26日 13:22
  • 691
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:SGU 106 The equation
举报原因:
原因补充:

(最多只允许输入30个字)