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;
}





[SGU]106. The Equation

Analysis     使用扩展欧几里德算法……其实是第一次听说这算法,百度了一下才有了思路……程序里有几个是否非零的判断。 Accepted Code var     x1,x2,x,...

sgu106 The Equation 扩展欧几里德

题意:求ax+by+c=0

X - The equation SGU oj106

Description There is an equation ax + by + c = 0. Given a,b,c,x1,x2,y1,y2 you must determ...

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

A - The equation Time Limit:250MS     Memory Limit:4096KB     64bit IO Format:%I64d & %I64u Submit ...

SGU 106. The equation 解题报告(模线性方程)

106. The equation time limit per test: 0.25 sec.  memory limit per test: 4096 KB There is an ...
  • kbdwo
  • kbdwo
  • 2014年04月24日 16:16
  • 487

SGU 106 The Equation 扩展欧几里得应用

Sol:线性不定方程+不等式求解 证明的去搜下别人的证明就好了。。。数学题。 #include #include #include using namespace std; long l...
  • imutzcy
  • imutzcy
  • 2013年12月15日 23:59
  • 630

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
  • 177

SGU_106 The equation 同余方程

http://acm.sgu.ru/problem.php?contest=0&problem=106 题意: 给你一个 ax+by+c = 0的方程, 然后分别给你x、y的区间[x1 , x2]...

【扩展欧几里得+解不等式】sgu106The equation

题目链接 题目大意:对于一个不定方程ax+by+c=0ax+by+c=0,其中x∈[x1,x2]x∈[x1,x2]和y∈[y1,y2]y∈[y1,y2]。求有多少组解。 数据范围:每一个数的绝对...
  • cqbztsy
  • cqbztsy
  • 2015年09月10日 16:50
  • 408

SGU106 The equation(数论)

给定一个二元一次方程ax+by+c=0,求该方程在x1 扩展欧几里得的应用。 先移项,c=-c,得到ax+by=c 特判一下a,b是否为0(貌似这道题的数据不用特判也没有关系) 令d=...
  • cqbzwja
  • cqbzwja
  • 2015年09月10日 20:32
  • 242
内容举报
返回顶部
收藏助手
不良信息举报
您举报文章:SGU 106 The equation
举报原因:
原因补充:

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