2020牛客NOIP赛前集训营-提高组(第一场)牛牛的方程式 (求解不定方程)

传送门
题意:
ax+by+cz=d中有没有至少存在一组{{x,y,z}}{x,y,z}的解,且x,y,z都为整数,使得方程式成立。如果成立yes 否者no,给你a,b,c,d,的值。

题解:看这个方程首先应该很容易想到扩展欧几里得,就是解不定方程的,那么首先我们可以写出ax+by=m; 这个方程要想有解必须ax+by=gcd(a,b)*k=m, 那么是不是现在我们就可以换元,
ax+by=gcd(a,b)*k, gcd(a,b)是已知,那么未知的k,我们就又可以把他看成是一个 gcd(a,b)*k+cz=d; 是不是就成一个标准的扩展欧几里得方程,要想此方程有解 那么我们就应该 找到 gcd(a,b,c),然后判断d%gcd(a,b,c)==0,gcd整除d,那么就有解,注意a,b,c,d种 0的时候;

篇外话:如果是找到x,y,z的值呢,那么就先算gcd(a,b)k+cz=d,我们可以找到K,和z的值,然后反代回去 计算ax+by=gcd(a,b)k,之后求出答案。

#pragma GCC optimize(2)
#pragma GCC optimize(3)
#include<bits/stdc++.h>
#define ll long long
#define ull unsigned long long
#define fio ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
#define mse(a,b) memset(a,b,sizeof a)
#define pb push_back
using namespace std;
const int mod=1e9+7;
const int maxx=1e6+10;
using namespace std;
const long double PI = 3.14159265358979323846;
#define fio ios_base::sync_with_stdio(0); cin.tie(0); cout.tie(0);
//inline int read(){ int x=0,f=1; char ch=getchar(); while (!isdigit(ch)) { if (ch=='-') f=-1; ch=getchar(); } while (isdigit(ch)) { x=x*10+ch-48;  ch=getchar();  } return x*f;}
//     ll cc = ((1ll << 62) - 1 + (1ll << 62));
int gcd(ll a,ll b)
{
    return b==0? a:gcd(b,a%b);
}
signed main()
{
    int n;
    cin>>n;
    while(n--){
    ll a,b,c,d;
     cin>>a>>b>>c>>d;
     ll dd=-0x7fffffffffffffff;
     if(a!=0){
     ll aa=gcd(gcd(a,b),c);
     dd=max(dd,aa);
     }
     if(a!=0){
     ll bb=gcd(gcd(a,c),b);
     dd=max(dd,bb);
     }
     if(b!=0){
     ll cc=gcd(gcd(b,c),a);
     dd=max(dd,cc);
     }
     //dd=max(aa,max(bb,cc));
     ///cout<<dd<<endl;
    if(d%dd==0) cout<<"YES"<<endl;
    else if(a==0&&b==0&&c==0&&d==0) cout<<"YES"<<endl;
    else cout<<"NO"<<endl;
    }
    return 0;
}


评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

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

抵扣说明:

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

余额充值