[OJ]信仰之跃、互质、求解同余方程、新世纪线段覆盖

最新推荐文章于 2023-02-16 10:55:34 发布
最新推荐文章于 2023-02-16 10:55:34 发布
阅读量1.9w 收藏 2
点赞数 1
分类专栏： C++编程 OJ 文章标签： C++ 经验分享 学习
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/mmmm0584/article/details/126183554
版权

1.信仰之跃

如图，一名刺客怀着信条，要进行他的信仰之跃。下面有个草垛，他必须跳到草垛里，不然会摔死。刺客可以看成一个质点，给出高度h和草垛距离前方刺客距离d的范围，请你计算他的初速度v（水平方向）的范围。（刺客在M星球，M星球的重力加速度g为18m/s^2）
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
图片来源：平抛运动8个公式是什么?
本题中只需要用到竖直位移和水平位移。对于本题来说，重要的是分数的表示。
分数的表示，参考：C++ 关于分数的处理与计算
其主要思想是用一个类表示分数，一个属性表示分子，另一个属性表示分母。再进行加减乘除。

/*
 * @Author: error: git config user.name && git config user.email & please set dead value or install git
 * @Date: 2022-08-05 15:15:22
 * @LastEditors: error: git config user.name && git config user.email & please set dead value or install git
 * @LastEditTime: 2022-08-05 18:30:35
 * @FilePath: \C++\testC++\JumpWall.cpp
 * @Description: 这是默认设置,请设置`customMade`, 打开koroFileHeader查看配置 进行设置: https://github.com/OBKoro1/koro1FileHeader/wiki/%E9%85%8D%E7%BD%AE
 */
#include<iostream>
#include<cmath>
#include<vector>
using namespace std;
const int  g=18;



class Assassin{
private:
    int x1;//水平方向区间的开始
    int x2;//区间结束
    int y;//垂直方向
public:
    Assassin(int x1,int x2,int y);
    int X1Get();
    void X1Set(int x1);
    int X2Get();
    void X2Set(int x2);
    int YGet();
    void YSet(int y);
    ~Assassin(){};
};

struct Fraction{
    int up;//分子
    int down;//分母
    Fraction(){};
    Fraction(int _up,int _down)
    {
        this->up = _up;
        this->down = _down;
    }
};

Fraction f1;
Fraction f2;
vector<Fraction> v1;//存放区间

//辗转相除法，递归，导致时间复杂度太大
int gcd(int a,int b)
{  
    if(b==0)
    {
        return a;
    }
    else
    {
        return gcd(b,a%b);
    }
    
}

//分数的化简
Fraction reduction(Fraction result)
{
    if(result.down < 0)
    {
        result.up *= (-1);
        result.down *= (-1);
    }
    if(result.up == 0)
        result.down = 1;
    int _gcd = gcd(abs(result.up),abs(result.down));
    result.up = result.up/_gcd;
    result.down = result.down/_gcd;
    return result;
}

Assassin::Assassin(int x1,int x2,int y)
{
    this->x1=x1;
    this->x2=x2;
    this->y=y;
}

int Assassin::X1Get()
{
    return x1;
}
    
void Assassin::X1Set(int x1)
{
    this->x1=x1;
}
int Assassin::X2Get()
{
    return x2;
}
    
void Assassin::X2Set(int x2)
{
    this->x2=x2;
}
int  Assassin::YGet()
{
    return y;
}

void Assassin::YSet(int y)
{
    this->y=y;
}

Fraction divide(Fraction f1 , Fraction f2)
{
    Fraction f3;
    if(f2.up == 0)
        return f2;//分母为0
    f3.up = f1.up*f2.down;
    f3.down = f1.down*f2.up;
    f3 = reduction(f3);
    return f3;
}

void print(Fraction f1,Fraction f2)
{
    if(f1.down==1 && f2.down==1)//如果分母是1,则当做整数输出
        {
            cout<<f1.up<<" "<<f2.up<<endl;
        }
        else if(f1.down==1 && f2.down>1)
        {
            cout<<f1.up<<" "<<f2.up<<"/"<<f2.down<<endl;
        }
        else if(f1.down>1 && f2.down==1)
        {
            cout<<f1.up<<"/"<<f1.down<<" "<<f2.up<<endl;
        }
        else
        {
            cout<<f1.up<<"/"<<f1.down<<" "<<f2.up<<"/"<<f2.down<<endl;
        }
}

void CalV(int &d1,int &d2,int &h)
{
    Assassin assassin(d1,d2,h);
    Fraction t;//表示t
    //根据y=1/2gt^2,求t,t可能为分数
    //int t=int(sqrt(2*assassin.YGet()/g));
    t.up=int(sqrt(assassin.YGet()));
    t.down=sqrt(g/2);
    t=reduction(t);//化简

    //根据x=vt,求v，v可能为分数
    
    //将f1(x1)看做一个分数
    f1.up=assassin.X1Get();
    f1.down=1;
    //v=x/t
    f1=divide(f1,t);

    f1=reduction(f1);//约分
    v1.push_back(f1);

    //f2(x2)
    f2.up=assassin.X2Get();
    f2.down=1;
    f2=divide(f2,t);
    f2=reduction(f2);//约分
    v1.push_back(f2);
}

int  main()
{
    int T;
    cin>>T;
    
    int h,d1,d2;
    for(int i=0;i<T;++i)
    {
        cin>>h>>d1>>d2;
        CalV(d1,d2,h);
    }
    //正向输出所有元素
    vector<Fraction>::iterator it;
    for(it=v1.begin();it!=v1.end();it+=2)
    {
        print(*it,*(it+1));
    }
    return 0;
}

在这里插入图片描述
在这里插入图片描述

2.互质

输入两个正整数n和k，求与n互质的第k个正整数。
在这里插入图片描述

/*
 * @Author: error: git config user.name && git config user.email & please set dead value or install git
 * @Date: 2022-08-04 16:47:37
 * @LastEditors: error: git config user.name && git config user.email & please set dead value or install git
 * @LastEditTime: 2022-08-04 23:28:24
 * @FilePath: \C++\testC++\PrimeNumber.cpp
 * @Description: 这是默认设置,请设置`customMade`, 打开koroFileHeader查看配置 进行设置: https://github.com/OBKoro1/koro1FileHeader/wiki/%E9%85%8D%E7%BD%AE
 */
#include<iostream>
using namespace std;
const int maxn=1000006;
int prime[maxn],num=0;

/*
采用欧拉函数进行求取
对正整数n，欧拉函数是小于等于n的正整数中与n互质的数的数目.
欧拉函数Euler(n)：求[2,n]中有多少个数与n互素

直接利用公式：φ(x)=x(1-1/p1)(1-1/p2)(1-1/p3)(1-1/p4)…..(1-1/pn)

其中：

pi为x的素因数
每个素因数只用一次
比如90=2 * 3^2 * 5
φ(90) = 90 * (1-1/2) * (1-1/3) * (1-1/5)
φ(1)=1（唯一和1互质的数(小于等于1)就是1本身）
如果不同的整数X和Y的最⼤公约数为1，则称X和Y为互质数。
Y的最⼤公约数为1，则称X和Y为互质数。显然，X和Y互质等价于Y和X互质。给出整数N，求1~N之间，有多少个数对
(X,Y)是互质数。
例如N=5时，(1,2) , (1,3), (1,4), (1,5), (2,3), (2,5), (3,4), (3,5), (4,5)共9对数是互质数。
*/

//辗转相除法，递归，导致时间复杂度太大
int gcd(int a,int b)
{  
    if(b==0)
    {
        return a;
    }
    else
    {
        return gcd(b,a%b);
    }
    
}

/*
由于当a与b互质时，a+b与b互质，可以求出a以内的互质的数，
由周期性推出剩下所有的与a互质的数
*/
int FindK(int n,int k)
{
    for(int i=1;i<n;i++)
    {
        if(gcd(n,i)==1)
        {
            prime[++num]=i;
        }
    }
   return k/num*n+prime[k%num];
}

int  main()
{
    int n,k;
    cin>>n>>k;
    int result=FindK(n,k);
    cout<<result;
    return 0;
}

在这里插入图片描述
在这里插入图片描述
参考：洛谷P1592 互质

3.求解同余方程

对于同余方程，求解它的最小整数解x0
在这里插入图片描述
方案：采用扩展欧几里得，即辗转相除法

/*
 * @Author: error: git config user.name && git config user.email & please set dead value or install git
 * @Date: 2022-08-05 09:42:32
 * @LastEditors: error: git config user.name && git config user.email & please set dead value or install git
 * @LastEditTime: 2022-08-05 10:31:03
 * @FilePath: \C++\testC++\congruentEquation.cpp
 * @Description: 这是默认设置,请设置`customMade`, 打开koroFileHeader查看配置 进行设置: https://github.com/OBKoro1/koro1FileHeader/wiki/%E9%85%8D%E7%BD%AE
 */
#include<iostream>
using namespace std;

/*
 扩展欧几里得是用来求这样的一组解的：ax + by = gcd(a, b)，求x和y，
 
    （求出x后自然知道了y，所以算是求一个x）
    逆元呢是求这样的一个解：ax ≡ 1 (mod b)，（为了方便理解，改了一下变量名），
    求x，貌似并没有一点点相似处，那么我们变一下，
    ax+by = gcd(a,b)，变成 ax+by = c；
    ax ≡ 1 (mod b)，变成 ax-by = 1；如果将y看成负的，ax+by = 1；
    完全一样,所以直接套用扩展欧几里得求
*/

int exgcd(int a,int b,int &x,int &y)
{
    //gcd(a,b)表示a和b的最大公约数,0与一个数的最大公约数是这个数本身
    //此时，毫无疑问gcd(a,b) == a,此时存在解x = 1，y = 0;
    if(b==0)
    {
        y=0;
        x=1;
        return a;
    }
    int r=exgcd(b,a%b,x,y);
    /*
       '''
    b*x1 + (a%b)*y1 = gcd
    由 a%b = a - (a/b)*b（这里的 “/” 指的是整除，例如 5/2=2 , 1/3=0）可得
    gcd = b*x1 + (a-(a/b)*b)*y1
        = a*y1 + b*(x1 – a/b*y1)
    所以
    x = y1
    y = x1 – a/b*y1
    '''
    */
    int temp=x;
    x=y;
    y=temp-a/b*y;
   
    return r;

}

int  main()
{
    int a,b,x,y;
    cin>>a>>b;
    int num=exgcd(a,b,x,y);
    cout<<(x%b+b)%b<<endl;
    return 0;
}

在这里插入图片描述
在这里插入图片描述

4.新世纪线段覆盖

现给定一根坐标轴，坐标轴上有n条线段，每条线段从1开始编号，若线段编号为k，则该线段覆盖的区间是：[（k-1）×m，（k-1）×m+l]。
有一个质点从原点开始向x轴正半轴移动，每次移动d距离。问第一个不被线段覆盖的点的位置在哪里？
在这里插入图片描述

对于该题来说，暴力解法的时间复杂度相当高，因此采用区间求解。
以样例1为例：此时n=5,d=4,m=11,l=8;
对于1到n的区间内，从1开始遍历到n，
当k=1时，[(k-1)*m,(k-1)*m+l]为[0,8];
当k=2时，[(k-1)*m,(k-1)*m+l]为[11,19];
当k=3时，[(k-1)*m,(k-1)*m+l]为[22,30];
…
因此，只需要判断(8,11）（19,22）之间的数是否是d的倍数（注意，这里是开区间），仅且判断一次，如果是，这直接跳出循环，返回数字。
当超出n的时候，例如n=2,d=4,m=11,l=8;
则用(n-1)*m+l来对d取余，设j=（(n-1)*m+l）%d,那么还差d-j个数便是下一个d的倍数。
代码如下：

/*
 * @Author: error: git config user.name && git config user.email & please set dead value or install git
 * @Date: 2022-08-05 19:04:25
 * @LastEditors: error: git config user.name && git config user.email & please set dead value or install git
 * @LastEditTime: 2022-08-05 20:05:51
 * @FilePath: \C++\testC++\LineCoverage.cpp
 * @Description: 这是默认设置,请设置`customMade`, 打开koroFileHeader查看配置 进行设置: https://github.com/OBKoro1/koro1FileHeader/wiki/%E9%85%8D%E7%BD%AE
 */
#include<iostream>
using namespace std;

/*
每次移动d,记录x当前的位置
如果采用暴力法，时间复杂度很高
*/
int LineCoverage(int n,int d,int m,int l)
{
    if(n==1)
    {
         //如果超出n条线段的最大区间
        int num=((n-1)*m+l)%d;
        return ((n-1)*m+l)+(d-num);
    }
    else if(n>1)
    {
        //当n>1时
        for (int i=1;i<n;i++)
        {  
            //在区间((i-1)*m+l)--((i+1-1)*m)中的数，是否是d的倍数
            //如果是，则寻找超出((i-1)*m+l)的第一个数
            for(int k=((i-1)*m+l)+1;k<((i+1-1)*m);k++)
            {
                if(k%d==0)
                {
                    return k;
                }
            }     
            
        }
        //如果超出n条线段的最大区间
        int num=((n-1)*m+l)%d;
        return ((n-1)*m+l)+(d-num);
    }
    return 0;
}


int  main()
{
    int n,d,m,l;
    cin>>n>>d>>m>>l;
    int result=LineCoverage(n,d,m,l);
    cout<<result<<endl;
    return 0;
}

在这里插入图片描述
在提交时，有一处错误，应该是有种情况没有考虑到
在这里插入图片描述
经过反复检查发现当n,d,m,l取1e6时，超出int范围，导致错误发生，因此修改变量类型，改为long long ，代码如下：

/*
 * @Author: error: git config user.name && git config user.email & please set dead value or install git
 * @Date: 2022-08-05 19:04:25
 * @LastEditors: error: git config user.name && git config user.email & please set dead value or install git
 * @LastEditTime: 2022-08-05 22:13:14
 * @FilePath: \C++\testC++\LineCoverage.cpp
 * @Description: 这是默认设置,请设置`customMade`, 打开koroFileHeader查看配置 进行设置: https://github.com/OBKoro1/koro1FileHeader/wiki/%E9%85%8D%E7%BD%AE
 */
#include<iostream>
using namespace std;

/*
每次移动d,记录x当前的位置
如果采用暴力法，时间复杂度很高
*/
long long LineCoverage(long long n,long long d,long long m,long long l)
{
    if(n==1)
    {
         //如果超出n条线段的最大区间
        int num=((n-1)*m+l)%d;
        return ((n-1)*m+l)+(d-num);
    }
    else if(n>1)
    {
        //当n>1时
        for (int i=1;i<n;i++)
        {  
            //在区间((i-1)*m+l)--((i+1-1)*m)中的数，是否是d的倍数
            //如果是，则寻找超出((i-1)*m+l)的第一个数
            for(long long k=((i-1)*m+l)+1;k<((i+1-1)*m);k++)
            {
                if(k%d==0)
                {
                    return k;
                }
            }     
            
        }
        //如果超出n条线段的最大区间
        int num=((n-1)*m+l)%d;
        return ((n-1)*m+l)+(d-num);
    }
    return 0;
}


int  main()
{
    long long n,d,m,l;
    cin>>n>>d>>m>>l;
    long long result=LineCoverage(n,d,m,l);
    cout<<result<<endl;
    return 0;
}

在这里插入图片描述

臧初之
关注
专栏目录
东北大学OJ-1539: 实验3-10：解方程
亓官劼的博客
06-04 1386
东北大学OJ-1539: 实验3-10：解方程   大家好，我叫亓官劼（qí guān jié ），在CSDN中记录学习的点滴历程，时光荏苒，未来可期，加油~博客地址为：亓官劼的博客，B站昵称为：亓官劼，地址为亓官劼的B站 本文原创为亓官劼，请大家支持原创，部分平台一直在盗取博主的文章！！！ 本专栏为东北大学C语言课程题库的题解，内含盖课程要求的50题的题解，请各位支持原创，每个都是一字一字亲手玛出来的代码，目前仅在CSDN发布。 本栏目的文章导航为：东北大学C语言课程题库题解专栏目录 如果订阅
oj做题c语言
最新发布
Hello大家好，我是Dream，如果帮得到你，那我深感荣幸！交流学习、商务合作：https://bbs.csdn.net/topics/614347534
08-18 4万+
输入两个整数A,B，用空格隔开输出一个整数，表示这两个数的和数据范围0≤A,B≤108。
今日学习在线编程题：新世纪线段覆盖
matiji2008的博客
04-09 423
今日学习在线编程题：新世纪线段覆盖
数论-中国剩余定理（模数互质的同余方程组）
qq_53269459的博客
08-10 1693
求同余方程组（中国剩余定理） 在《孙子兵法》中有一这样问题：“今有物不知其数，三三数之剩二（除以3余2），五五剩值三（除以5余3），七七数之剩二（除以7余2），问物几何？”这个问题称为“孙子问题”，该问题解法称为中国剩余定理。具体分三步： 1.找出三个数：3和5公倍数中找出被7除余1的最小数为15，从3和7的公倍数中找出被5除余1的最小数为21，最后从5和7的公倍数中找出除3余1的最小数70。 2.用15乘2（得到30的是除7余2的最终结果），用21乘3（得到的63是除5余2的最终结果），用70乘2（140
杭电oj1014 互质就可以了
西柚茶茶的博客
07-28 369
#include #include using namespace std; int a(int step, int mod) { if (mod == 0) return step; else a(mod, step%mod); } int main() { int step, mod; while (cin >> step >> mod) { i
《C语言学习笔记10》指针入门
qq_49594899的博客
07-22 863
常见的指针问题及疑惑的解决博文
vue2版本《后台管理模式》（中）
老余的博客
02-16 816
登录页需要先设token值 看看token 是否存在 如果存在则跳转页面 不存在留着继续输入。还需要输入框的正则方式 格式不对无法提交账号信息 接下来就是需要登录接口了。看着挺多的但是你仔细捋一下思路 都是有规律的 好比上学 需要语文课本 就需要领语文课本，而语文课本是有编写者写的。你要是有足够时间你也可以自己配置一本语文书 然后想学就拿！！！回到代码上 登录页需要什么 你配置什么最后引入就行。
洛谷P1803 凌乱的yyy / 线段覆盖 第6个用例 输入
10-11
洛谷P1803 凌乱的yyy / 线段覆盖 第6个用例 输入 比较长所以这里就不直接输入了，大家可以自行下载当作自己的测试用例。 这个用例应该可以算是比较具有普遍性的用例了。 输入格式： 第一行是一个整数 n，接下来 n 行...
【如皋OJ】C++判断互质
格其非心的博客 C++学习
12-24 3275
1075: ex423：判断互质 时间限制: 1 Sec  内存限制: 128 MB 提交: 162  解决: 78 [提交] [状态] [讨论版] [命题人:外部导入] 题目描述 输入两个正整数m和n，判断m和n是否互质（即最大公约数为1）。   输入 只有1行2个正整数m和n（1&lt;=n,m&lt;231）。   输出 如果2个整数互质，就输出Yes；否则输出No，并在第...
【求解三对角矩阵线性方程组】第七次OJ的总结
SOAR
12-05 1637
问题描述 三对角矩阵是指除主对角线和其相邻的上下两条对角线之外，其他所有元素都为0的矩阵。 本次实验给定一个三对角矩阵A和右端矩阵B，求解矩阵方程AX=B。 可以认为B=[b1,b2,...,bn]为一些列向量的组合，从而解出的矩阵X=[x1,x2,...,xn]为一些解向量的组合。 本次实验给出的矩阵A均为三对角非奇异方阵，解的准确程度由残差值 ||Axi - bi||2给出，这里i从1...
【C语言——递归】迷宫问题
weixin_44348253的博客
10-03 614
【代码】【C语言——递归】迷宫问题。
【828】
weixin_62369651的博客
08-28 392
自己刷题
Vue框架的学习(Vue的基础指令操作二 V-For)第二课
qq_56248592的博客
11-11 786
V-For
1.2继承性
qq_62865634的博客
09-11 226
本文仅仅简单介绍了继承（子承父业）、继承行高的关系。
Vue框架的学习(Vue的基础指令操作一)第二课
qq_56248592的博客
11-11 985
Vue基础指令一
go 验证字符串中是否包含中文或英文
是鲤鱼啊
09-16 1608
【代码】go 验证字符串中是否包含中文或英文。
git 提交代码报错，提示邮箱格式不正确
qq_41572507的博客
03-19 4161
首先查看email和name是否是公司要求正确格式 查看name： git config user.name 查看email：git config user.email 如果email不正确进行修改 修改email：git config user.email "Author Email" 修改name：git config user.name "Author Name" 修改之后进行再次进行...
git用户名配置不显示user.name及user.email报错及解决
jinyeran的博客
12-28 1万+
解决方案：仔细核对输入格式 注意邮箱没有“”，name和email后有空格
oj求一元二次方程的解
10-28
oj可以通过以下代码求解一元二次方程的解： ``` #include #include using namespace std; int main() { double a, b, c; cin >> a >> b >> c; double delta = b * b - 4 * a * c; if (delta >= 0) { double...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值