首页 > 其他分享 >GCD辗转相除的经典套路

GCD辗转相除的经典套路

时间:2023-01-22 20:44:43浏览次数:63  
标签:GCD 套路 cout 相除 st int ans now gcd

 


  • 前两个移动很像辗转相除法(这个套路在 Codeforces 上已经出烂了)<br>
  • 后两个移动可以让 g 乘上 $2^k$

class Solution {
public:
    bool isReachable(int X, int Y) {
        while(!(X&1))X >>= 1;
        while(!(Y&1))Y >>= 1;
        return gcd(X,Y) == 1;
    }
};

 

 


  • https://codeforces.com/contest/1766/problem/D

题意:

找到连续的最长gcd(a+k,b+k) == 1(a < b,k = {0,1,2,...})


 

  • 思路:
  •  gcd(a+k,b+k) == gcd(a+k,b - a)
    1. a - b = 1时特判
    2. 可以推出$gcd(a+k,b+k) == gcd(a+k,b - a)$,具体证明见https://codeforces.com/blog/entry/110066
    3. 设两个的结果分别为g和h,显然a+k和b+k都是g的倍数,那么让其中一个大的倍数减掉一个小的倍数,显然还是g的倍数,充分性成立
    4. 反证法也同理,a+k和b-a肯定都是h的倍数,一个倍数加上另一个倍数,肯定还是h的倍数,必要性成立

  • 接下来只需要找到b-a的质因子判断a+k是否有重复因子
  • 一个一个加显然超时,我们可以先预处理出每个数字的最大质因子,然后利用传递性,得到关于b-a的所有质因子
  • 设当前质因子为temp,更新答案的方法就是$ans = min(ans,(temp - a%temp) % temp)$
#include<bits/stdc++.h>
#define debug1(a) cout<<#a<<'='<< a << endl;
#define debug2(a,b) cout<<#a<<" = "<<a<<"  "<<#b<<" = "<<b<<endl;
#define debug3(a,b,c) cout<<#a<<" = "<<a<<"  "<<#b<<" = "<<b<<"  "<<#c<<" = "<<c<<endl;
#define debug4(a,b,c,d) cout<<#a<<" = "<<a<<"  "<<#b<<" = "<<b<<"  "<<#c<<" = "<<c<<"  "<<#d<<" = "<<d<<endl;
#define debug5(a,b,c,d,e) cout<<#a<<" = "<<a<<"  "<<#b<<" = "<<b<<"  "<<#c<<" = "<<c<<"  "<<#d<<" = "<<d<<"  "<<#e<<" = "<<e<<endl;
#define debug0(x) cout << "debug0: " << x << endl
#define fr(t, i, n)for (long long i = t; i < n; i++)
#define YES cout<<"Yes"<<endl
#define nO cout<<"no"<<endl
#define fi first
#define se second
//#define int long long
using namespace std;

typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;
typedef pair<LL,LL> PLL;

//#pragma GCC optimize(3,"Ofast","inline")
//#pragma GCC optimize(2)

const int N = 2e5+10,M = 1e7+10;
int to[M];
void solve() 
{
    int x,y;scanf("%d%d",&x,&y);
    if(y == x+1)
    {
        printf("-1\n");
        return ;
    }
    
    int now = y - x;
    int ans = 1e9;
    for(;now != 1;now /= to[now])
        ans = min(ans, (now - x % to[now]) % to[now]);
    
    printf("%d\n",ans);
}

signed main()
{
    /*
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    */

    to[1] = 1;//每个数字到它最大的质因子
    for(int i = 2;i < M;i ++)if(!to[i])
    {
        for(int j = i;j < M;j += i)to[j] = i;
    }

    int T = 1;
    scanf("%d",&T);
    //cin >> T;
    while(T--){
        solve();
    }
    return 0;
}
  • https://codeforces.com/contest/1762/problem/D

题意:

给一个长度为n的permutation,每次一个询问,得到结果为gcd(i,j),请在2*n次之内找到那个是0(或者哪两个之中之一是0)

思路:

  • 三个指针i,j,k(i<j<k)

  • 令x=gcd(a[i],a[j]),y=gcd(a[i],a[k]);
    1. 如果x==y,显然a[i]>0     2. 如果x<y,可以证明a[j]>0     3. 如果x>y,可以证明a[k]>0
  • 这样就可以写出一个答案了
#include<bits/stdc++.h>
#define debug1(a) cout<<#a<<'='<< a << endl;
#define debug2(a,b) cout<<#a<<" = "<<a<<"  "<<#b<<" = "<<b<<endl;
#define debug3(a,b,c) cout<<#a<<" = "<<a<<"  "<<#b<<" = "<<b<<"  "<<#c<<" = "<<c<<endl;
#define debug4(a,b,c,d) cout<<#a<<" = "<<a<<"  "<<#b<<" = "<<b<<"  "<<#c<<" = "<<c<<"  "<<#d<<" = "<<d<<endl;
#define debug5(a,b,c,d,e) cout<<#a<<" = "<<a<<"  "<<#b<<" = "<<b<<"  "<<#c<<" = "<<c<<"  "<<#d<<" = "<<d<<"  "<<#e<<" = "<<e<<endl;
#define debug0(x) cout << "debug0: " << x << endl
#define fr(t, i, n)for (long long i = t; i < n; i++)
#define YES cout<<"Yes"<<endl
#define nO cout<<"no"<<endl
#define fi first
#define se second
// #define int long long
using namespace std;

typedef long long LL;
typedef unsigned long long ULL;
typedef pair<int,int> PII;
typedef pair<LL,LL> PLL;

//#pragma GCC optimize(3,"Ofast","inline")
//#pragma GCC optimize(2)

const int N = 2e5+10,mod = 998244353;
bool st[N];

int ask(int a, int b) {
    cout << "? " << a << " " << b << endl;
    int ans = 0;
    cin >> ans;
    return ans;
}

void solve() 
{
    memset(st,0,sizeof st);
    int n;cin >> n;

    int a[3] = {1,2,3};
    for(;;)
    {
        sort(a,a+3);
        if(a[2] > n)break;
        int x = ask(a[0],a[1]),y = ask(a[0],a[2]);
        if(x == y)
        {
            st[a[0]] = 1;
            a[0] = a[2] + 1;
        }else if(x < y)
        {
            st[a[1]] = 1;
            a[1] = a[2] + 1;
        }else{
            st[a[2]] = 1;
            a[2] = a[2] + 1;
        }
    }

    cout << "! " << a[0] << " " << a[1] << endl;
    int t;cin >> t;
}

signed main()
{
    /*
    ios::sync_with_stdio(false);
    cin.tie(0);
    cout.tie(0);
    */
    int T = 1;cin >> T;
    while(T--){
        solve();
    }
    return 0;
}

 

标签:GCD,套路,cout,相除,st,int,ans,now,gcd
From: https://www.cnblogs.com/cfddfc/p/17064646.html

相关文章

  • webFlux routers 嵌套路由解析规则
    有下列路由publicRouterFunction<ServerResponse>doctorRoutes(DoctorHandlerhandler){returnRouterFunctions.route().path("/doctors"......
  • 「学习笔记」莫比乌斯反演套路总结
    基本套路Crash的数字表格/JZPTAB求:\[\sum_{i=1}^n\sum_{j=1}^m\mathrm{lcm}(i,j)\]\[\begin{aligned}&\sum_{i=1}^n\sum_{j=1}^m\mathrm{lcm}(i,j)\\=&\sum_{i=......
  • exgcd & 线性同余方程
    前置芝士裴蜀定理exgcdexgcd即扩展欧几里得定理,常用来求解\(ax+by=gcd(a,b)\)的可行解问题推导过程:考虑我们有:​ \(ax+by=gcd(a,b)\)——裴蜀定理​ \(a_......
  • leetcode 每日一题 gcd+枚举
    1819.序列中不同最大公约数的数目给你一个由正整数组成的数组nums。数字序列的最大公约数定义为序列中所有整数的共有约数中的最大整数。  例如,序列[4,6,16]......
  • abc254 F - Rectangle GCD
    题意:给定等长的正整数组\(a[],b[]\),它们确定了一个矩阵\(A_{i,j}=a_i+b_j\)。\(q\)次询问,回答矩阵中一个矩形区域内所有数的\(\gcd\)\(n,q\le2e5\)思路:差分,绝对......
  • Apple开发_使用 GCDAsyncUdpSocket 发送组播消息报错"No route to host"的解决
    1、问题描述苹果手机升级到ios14.5系统后,使用GCDAsyncUdpSocke发送组播消息的时候,发现报错了,ErrorDomain=NSPOSIXErrorDomainCode=65"Noroutetohost"UserInfo=......
  • 扩展欧几里得(exgcd(a,b))
    回顾       由贝祖定理可知:ax+by=gcd(a,b)       (a,y)为一组解      (x+kb,y-ka)也为解       所以有无数组解一、扩展欧几里得算法......
  • 学习开源项目的几个实用套路
    记得我的leader之前说过,很多人工作之后就丧失了钻研技术的热情,这个确实,我发现自己多少也有这个问题。转眼已经毕业一年多了,回想这一年,有些惭愧,感觉不仅技术能力上并没有......
  • gcd(a, b, c) = gcd(gcd(a, b), c)
    某一天,我正苦逼的刷题看题解,看到下面的代码inttmp=0; for(inti=1;i<=n;++i){ scanf("%d",&a[i]); tmp=gcd(tmp,a[i]); }​ 我心中一惊:wc,这就能求gcd(a1,a2......
  • P2398 GCD SUM——欧拉函数
    此题可以拓展为\(\sum\limits^n_{i=1}\sum\limits^m_{j=1}\gcd(i,j)\)结论是\(\sum\limits^{\min(n,m)}_{d=1}\varphi(d)\lfloor\dfrac{n}{d}\rfloor\lfloor\dfrac{m}{......