[7.22NOIP模拟测试7]方程的解 题解(扩展欧几里得)

Orz

送分比较慷慨的一道题，疯狂特判能拿不少分。

对于$a>0,b>0$的情况：

用exgcd求出方程通解，然后通过操作得到最小正整数解和最大正整数解

他们以及他们之间的解满足等差数列性质，小学数奥求项数即可

（其实就是(末项-首项)/公差+1）

其他情况特判掉或者转化为可处理情况即可（比如全负），不多说，代码里写的还是比较清晰的

//#define XR
#include<cstdio>
#include<iostream>
#include<cstring>
using namespace std;
int T,a,b,c,x,y;
//void exgcd(int a,int b,int &x)

int exgcd(int a,int b,int &x,int &y,int c)
{
    if(!b)
    {
        x=c/a;
        y=;
        return a;
    }
    int g=exgcd(b,a%b,y,x,c);
    y-=a/b*x;
    return g;
}
void work()
{
    scanf("%d%d%d",&a,&b,&c);x=y=;
    if(a==&&b==)
    {
        if(c==)
        {
            puts("ZenMeZheMeDuo");
            return ;
        }
        else
        {
            puts("");
            return ;
        }
    }
    if(a==||b==)
    {
        long long now=a+b;
        if(c==||(c%now==&&(long long)now*c>))
        {
            puts("ZenMeZheMeDuo");
            return ;
        }
        else
        {
            puts("");
            return ;
        }
    }
    if((a>&&b>&&c<=)||(a<&&b<&&c>=)||(a>&&b>&&a+b>c)||(a<&&b<&&a+b<c))
    {
        puts("");
        return ;
    }
    if(a==b&&a==)
    {
        if(c>)puts("ZenMeZheMeDuo");
        else if(c<=)puts("");
        else cout<<c-<<endl;
        return ;
    }
    if(a+b==c)
    {
        puts("");
        return ;
    }
    if(a<&&b<)a=-a,b=-b,c=-c;
    int GCD=exgcd(a,b,x,y,c);//cout<<GCD<<endl;
    if(c%GCD!=)
    {
        puts("");
        return ;
    }

    if((long long)a*b<)
    {

        puts("ZenMeZheMeDuo");
        return ;
    }
    a/=GCD;b/=GCD;c/=GCD;x%=b;
    while(x<=)x+=b;
    y=(c-a*x)/b;
    int ym=y%a;
    while(ym<=)ym+=a;int ans;
    if(ym>y)ans=;
    else ans=(y-ym)/a+;
    if(ans>)puts("ZenMeZheMeDuo");
    else cout<<ans<<endl;

}
void test()
{
    scanf("%d%d",&a,&b);
    exgcd(a,b,x,y,c);
    cout<<x<<' '<<y<<endl;
}
int main()
{
//    cout<<(18%(-5))<<endl;
    //while(1)test();
#ifdef XR
    freopen("data.in","r",stdin);
    freopen("data.out","w",stdout);
#endif
    scanf("%d",&T);
    while(T--)work();
    return ;
}