裴蜀定理简单应用「建议收藏」

裴蜀定理简单应用「建议收藏」裴蜀定理定理内容:设aaa,bbb是不全为000的整数,则存在整数xxx,yyy使得a⋅xa\cdotxa⋅x+++b⋅yb\cdotyb⋅y=gcd⁡(x,y)\gcd(x,y)gcd(x,y)。定理简单应用:例题:洛谷p4549https://www.luogu.com.cn/problem/P4549思路分析:给定一个序列,求一个SSS满足S=S=S=∑i=1n\sum\limits_{i=1}^ni=1∑n​Ai×XiA_i\timesX_iAi​×Xi​,而

大家好,又见面了,我是你们的朋友全栈君。如果您正在找激活码,请点击查看最新教程,关注关注公众号 “全栈程序员社区” 获取激活教程,可能之前旧版本教程已经失效.最新Idea2022.1教程亲测有效,一键激活。

Jetbrains全家桶1年46,售后保障稳定

裴蜀定理

定理内容:

  • a a a, b b b是不全为 0 0 0的整数,则存在整数 x x x y y y使得 a ⋅ x a\cdot x ax + + + b ⋅ y b\cdot y by = gcd ⁡ ( x , y ) \gcd(x,y) gcd(x,y)

定理简单应用:

例题:

洛谷p4549

https://www.luogu.com.cn/problem/P4549

思路分析:
  • 给定一个序列,求一个 S S S满足 S = S = S= ∑ i = 1 n \sum\limits_{i=1}^n i=1n A i × X i A_i\times X_i Ai×Xi,而且要求满足这个条件的S的最小值,这时我们想到对于任意不为 0 0 0的整数都有 a ⋅ x a\cdot x ax + + + b ⋅ y b\cdot y by = gcd ⁡ ( x , y ) \gcd(x,y) gcd(x,y)
  • 那么对这个定理的另一个解读就是总有 x x x y y y使得 a ⋅ x a\cdot x ax + + + b ⋅ y b\cdot y by = d d d,且 gcd ⁡ ( x , y ) \gcd(x,y) gcd(x,y) | d d d
  • 因此我们想到这个定理是否可以推广成 n n n个数呢?答案是肯定的。
  • 所以此题即要我们求这一个序列的 gcd ⁡ \gcd gcd即可。
代码如下:
#include <bits/stdc++.h>
using namespace std;
int a[100];
int gcd(int a, int b)
{ 
   
        if (b == 0)
                return a;
        else
                return gcd(b, a % b);
}
int main()
{ 
   
        int n;
        cin >> n;
        cin >> a[0];
        int ans = a[0];
        for (int i = 1; i < n; i++)
        { 
   
                cin >> a[i];
                ans = gcd(a[i], ans);
        }
        cout << abs(ans) << endl;
}

Jetbrains全家桶1年46,售后保障稳定

Codeforces Round #290 (Div. 2) D. Fox And Jumping

https://www.luogu.com.cn/problem/CF510D

思路分析:
  • 要到达每一格,那么我们就要使选上的数满足 gcd ⁡ \gcd gcd = 1。
  • 在这里我们用到了dp,也就是说要选上的数的价格(使 gcd ⁡ \gcd gcd = tmp)和之前的满足( gcd ⁡ \gcd gcd = tmp)的价格取最少即可。然后我们需要对每个数进行配对(不止是两两配对),所以我们想到了用一个map来储存下标为当前最大公因子数,value为所花价值来操作。
  • 具体细节见代码注释
代码如下:
#include <bits/stdc++.h>
using namespace std;
map<int, int> dp;
int l[301];
int c[301];
int gcd(int a, int b)
{ 
   
        if (b == 0)
                return a;
        else
                return gcd(b, a % b);
}
//求最大公因子
int main()
{ 
   
        int n;
        scanf("%d", &n);
        dp.clear();
        //容器清空
        for (int i = 1; i <= n; i++)
        { 
   
                scanf("%d", &l[i]);
        }
        //读入数
        for (int i = 1; i <= n; i++)
        { 
   
                scanf("%d", &c[i]);
        }
        //读入选择该数的费用
        dp[0] = 0;
        //初始化
        for (int i = 1; i <= n; i++)
        { 
   

                map<int, int>::iterator it = dp.begin();
                //用迭代器
                for (; it != dp.end(); it++)
                { 
   
                        int tmp = gcd(it->first, l[i]);
                        //tmp即为选上的这个数和map里已经选好的数进行gcd运算
                        if (dp.count(tmp))
                        { 
   
                                dp[tmp] = min(dp[tmp], it->second + c[i]);
                                //就是说要取到gcd为tmp时的最小花费
                        }
                        else
                                dp[tmp] = it->second + c[i];
                                //之前数配对时没有出现过的gcd,因此直接储存
                }
        }
        if (dp.count(1))
                cout << dp[1] << endl;
                //最后由裴蜀定理可知我们要的是gcd = 1 的最小花费
        else
                cout << -1 << endl;
                //没有就代表无论怎么选gcd都不为1,那么总有格子跳不到
        return 0;
}
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 举报,一经查实,本站将立刻删除。

发布者:全栈程序员-用户IM,转载请注明出处:https://javaforall.cn/230606.html原文链接:https://javaforall.cn

【正版授权,激活自己账号】: Jetbrains全家桶Ide使用,1年售后保障,每天仅需1毛

【官方授权 正版激活】: 官方授权 正版激活 支持Jetbrains家族下所有IDE 使用个人JB账号...

(0)


相关推荐

  • 为什么0xffffffff是-1?(计算机对整型的存储)[通俗易懂]

    为什么0xffffffff是-1?(计算机对整型的存储)[通俗易懂]一个数字在计算机中都是以二进制补码的形式存储的。先了解这句核心。。。我们认为中的int整型数值顺序java中int类型是4个字节,也就是32位,其中第一位是符号位,int数值的存储结构我们利用System.out.println(Integer.toBinaryString(Integer.MAX_VALUE));拿到int的最大值,是1111111111111111111111111111111,31个1,首位是0(代表正数,省略了)那我们给int的最大值+1,会发生什么呢?Sys

  • ubuntu20.04清华源_ubuntu20.04更换国内源

    ubuntu20.04清华源_ubuntu20.04更换国内源Ubuntu22.04的稳定版计划于2022年4月21日发布。开发工作已经在紧锣密鼓地进行,它将遵循如下发布时间表:2022年2月24日:功能冻结2022年3月17日:用户界面冻结2022年3月31日:测试版发布2022年4月14日:候选版本2022年4月21日:最终稳定版本Ubuntu22.04仍在积极开发中。您不应该在生产机器或主系统上使用它。如果你想在备用机器或虚拟机上测试它,你可以从Ubuntu的网站下载每日

  • opkg list 报错

    opkg list 报错opkglistCollectederrors:*opkg_conf_load:Couldnotlock/var/lock/opkg.lock:Resourcetemporarilyunavail echo”nameserver114.114.114.114″>/tmp/resolv.conf rm-f/var/lock/opkg.lock opkgupdate

  • cavas的使用汇总[通俗易懂]

    cavas的使用汇总[通俗易懂]canvas图像处理汇总一、canvas的情况canvas自从出来了之后,在前端的图像处理上面提供了各种各样的遍历,虽然很多的操作其实都是要应用到算法的,但是这个也给前端提供了很多的可能性,其中最

  • pycharm激活码2021年3月最新在线激活[通俗易懂]

    pycharm激活码2021年3月最新在线激活,https://javaforall.cn/100143.html。详细ieda激活码不妨到全栈程序员必看教程网一起来了解一下吧!

  • 语义分割的定义_语义分割模型

    语义分割的定义_语义分割模型点击上方“AI算法修炼营”,选择加星标或“置顶”标题以下,全是干货引言语义分割结合了图像分类、目标检测和图像分割,通过一定的方法将图像分割成具有一定语义含义的区域块,并识别出每个区域块的…

发表回复

您的电子邮箱地址不会被公开。

关注全栈程序员社区公众号