牛客练习赛111-D

我们可以用扩欧求出来一组特解

也就是说如果有解必然是

a*x+b*y=c;

那么我们求出一组特解之后呢？

也就是说我们如何让a的系数变到[l,r]的范围呢

我们要想如果要变化a的系数同时b也可以满足系数变化使得答案有解那么所移动的步数必然是

h, h=`lcm(a,b)/a 不然不会满足b的系数

这里再稍作解释

由于我们已经找到了一组解

如果说 要对a的系数变化  那么他变化的值一定来自于b的系数

也就是要维持有解的话

x*a变化的值也一定是 b的倍数

那么最小就是lcm(a,b);

对应到a的系数变化上就是lcm(a,b)/a
 

要让x到  [l,r] 第一个大于l的解

那么就是这个步骤  最后判断跳了之后是否小于r是的话就是有解 不然没解

if(x<l)
    {
        
      if((l-x)%tt==0) //l到x的距离是系数变化的倍数,即刚好可以到l
      {
          cout<<"YES"<<endl;
          return ;
      }
    
        LL ll=(l-x)/tt*tt+tt+x;// 
        if(ll<=r)
        {
            cout<<"YES"<<endl;
            return ;
        }
    }

这是ac代码：

// 数学公式要变形
// 莫急莫急先读题
#include <bits/stdc++.h>
using namespace std;
#define lowbit(x) (x&(-x))
#define endl "\n"
#define ios ios::sync_with_stdio(0); cin.tie(0),cout.tie(0);
#define LF(x)   fixed<<setprecision(x)// c++ 保留小数
typedef long long LL;
typedef pair<int,int> PII;
typedef unsigned long long ULL;
const int N=1000010,M=1010,F=2*N,INF=0x3f3f3f3f,pp=133331;
const double pai=acos(-1.0);// pai
map<int,int> q;
LL t,n,m;
LL gcd(LL a,LL b)
{
    return b?gcd(b,a%b):a;
}
 
LL exgcd(LL a,LL b,LL &x,LL &y)
{
    if(!b)
    {
        x=1,y=0;
        return a;
    }
    LL d=exgcd(b,a%b,y,x);
    y-=a/b*x;
    return d;
}
LL lcm(LL a,LL b)
{
    return a/gcd(a,b)*b;
}
LL a,b,l,r;
void solve()
{
    cin>>a>>b>>n>>l>>r;
    LL ans=gcd(a,b),res=lcm(a,b);
    if(n%ans){
        cout<<"NO"<<endl;
        return ;
    }
    LL x,y;
    LL tt=res/a;
    exgcd(a,b,x,y);
    x*=n/ans,y*=n/ans;
    
    if(x>=l&&x<=r){
        cout<<"YES"<<endl;
        return ;
    }
    if(x<l)
    {
        
      if((l-x)%tt==0)
      {
          cout<<"YES"<<endl;
          return ;
      }
    
        LL ll=(l-x)/tt*tt+tt+x;
        if(ll<=r)
        {
            cout<<"YES"<<endl;
            return ;
        }
    }
    if(x>r)
    {
        
       if((x-r)%tt==0)
      {
          cout<<"YES"<<endl;
          return ;
      }
      LL rr=x-(x-r)/tt*tt-tt;
        if(l<=rr)
        {
            cout<<"YES"<<endl;
            return ;
        }
    }
    cout<<"NO"<<endl;
        return ;
}
 
int main ()
{
    ios// 不能有printf puts scanf
    cin>>t;
    while(t--)
    solve();
    return 0;
}