数学题代码模板

已于 2024-07-02 21:48:24 修改
阅读量613 收藏
点赞数
分类专栏： 个人学习 文章标签： 算法 数据结构
于 2021-01-19 09:48:56 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Cat_cat_cat__/article/details/112798964
版权

随我不断学习而持续更新

文章目录

1、质数的判定

试除法判定质数（只用循环sqrt（n）次）

bool is_prime(int n){
   if(n < 2)return false;
   for(int i = 2;i <= n/i; i++)
      if(n%i == 0)return false;
   return true;
}

2、质因数的分解

输出的是质因子及其个数

void divide(int n)
{
    for(int i = 2;i <= n/i;i ++ )
        if(n % i == 0)
        {
            int s=0;
            while(n%i==0)
            {
                n/=i;
                s++;
            }
            printf("%d %d",i,s);
        }
    if(n>1) printf("%d %d",n,1);

}

3、质数筛（埃式筛）

利用布尔值判定，再利用桶排序记录合数

int prime[];
bool st[];
void get_primes(int n)
{
    for(int i = 2;i <= n;i ++ )
    {
        if(!st[i])
        {
            primes[cnt ++] = i;//记录质数
            for(int j = i + i;j <= n;j += i) st[j] = true;//记录合数
        }
    }
}

4、质数（线性筛）//未学会

const int N=10000010;

int primes[N],cnt;
bool st[N];

void get_primes(int n)
{
    for(int i=2;i<=n;i++)
    {
        if(!st[i]) primes[cnt++]=i;
        for(int j=0;primes[j]<=n/i;j++)
        {
            st[primes[j]*i]=true;
            if(i%primes[j]==0) break;
        }
    }
}

5、求约数

试除法求约数

vector<int> get_divisors(int n)
{
    vector<int>res;
    for(int i=1;i<=n/i;i++)
        if(n%i==0)
        {
            res.push_back(i);
            if(i!=n/i) res.push_back(n/i);
        }
    sort(res.begin(),res.end());
    return res;
}

6、约数个数

#include <bits/stdc++.h>

using namespace std;

const int mod=1e9+5;
int main()
{
    int n;
    cin>>n;

    unordered_map<int,int>primes;//存底数和指数
    while(n--)
    {
        int x;
        cin>>x;
        for(int i=2;i<=x;i++)
            while(x%i==0)
            {
                x/=i;
                primes[i]++;
            }
        if(x>1) primes[x]++;
    }
    long long res =1;
    for(auto prime : primes) res = res * (prime.second +1)%mod;//循环用公式
    //如果 N = p1^c1 * p2^c2 * ... *pk^ck
    //约数个数： (c1 + 1) * (c2 + 1) * ... * (ck + 1)
    cout<<res<<endl;
    return 0;
}

7、合并果子

洛谷P1090

#include <bits/stdc++.h>

using namespace std;

int main()
{
    int n,c,ans,x,y;
    priority_queue<int,vector<int>,greater<int> >vis;
    cin>>n;
    ans=0;
    for(int i=1;i<=n;i++)
    {
        cin>>c;
        vis.push(c);
    }
        while(vis.size()>1)
    {
        x=vis.top(); vis.pop();
        y=vis.top(); vis.pop();
        ans+=x+y;
        vis.push(x+y);//注意不能直接将ans入堆，ans是最后的总和，而入堆的只是堆前两个数之和
    }
    cout<<ans<<endl;
    vis.pop();
    return 0;
}

8、约数的和

#include <bits/stdc++.h>

using namespace std;

const int mod=1e9+5;
int main()
{
    int n;
    cin>>n;

    unordered_map<int,int>primes;
    while(n--)
    {
        int x;
        cin>>x;
        for(int i=2;i<=x;i++)
            while(x%i==0)
            {
                x/=i;
                primes[i]++;
            }
        if(x>1) primes[x]++;
    }//求解约数结束
    long long res =1;
    for(auto prime : primes)//遍历求和开始
    {
        int p=prime.first,a=prime.second;
        long long t=1;
        while(a--)
            t=(t*p+1)%mod;//题干要求，答案取模
        res =res * t % mod;
    }
    cout<<res<<endl;

    return 0;
}

9、求两个数的最大公约数

#include <bits/stdc++.h>

using namespace std;

int gcd(int a,int b)
{
    return b ? gcd(b,a%b) : a;
}
int main()
{
    int n;
    scanf("%d",&n);
    while(n--)
    {
        int a,b;
        scanf("%d%d",&a,&b);
        printf("%d\n",gcd(a,b));
    }
    return 0;
}

10、快排模板

#include <bits/stdc++.h>

using namespace std;

const int N = 1e6 + 10;

int n;
int q[N];

void quick_sort(int q[],int l,int r)
{
    if(l >= r) return;//递归中止条件
    int x=q[(l+r)/2],i=l-1,j=r+1;//最先取中间值为断点，i为左指针，j为右指针，因为在循环里先移动指针，所以i，j要多初始化一个单位
    while(i<j)
    {
        do i++;while(q[i]<x);//当找到断点左边不小于x的值时循环结束
        do j--;while(q[j]>x);//当找到断点右边不大于x的值时循环结束
        if(i<j) swap(q[i],q[j]);//交换两个值
    }
    quick_sort(q,l,j);//递归继续排序
    quick_sort(q,j+1,r);
}
int main()
{
    scanf("%d",&n);
    for(int i = 0;i < n; i ++ ) scanf("%d",&q[i]);

    quick_sort(q,0,n-1);

    for(int i=0;i<n;i++) printf("%d ",q[i]);

    return 0;
}

11、归并排序 //未学会

#include <bits/stdc++.h>

using namespace std;

const int N = 1e6 + 10;

int n;
int q[N],tmp[N];

void merge_sort(int q[],int l,int r)
{
    if(l >= r) return;

    int mid = (l + r)>> 1;

    merge_sort(q,l,mid)，merge_sort(q,mid + 1,r);//谜一样的递归排序
    
    int k = 0,i = l,j = mid + 1;
    while(i <= mid && j <= r)
        if(q[i] <= q[j]) tmp[k ++ ]=q[i ++ ];
        else tmp[k ++ ]=q[j ++ ];
    while(i <= mid) tmp[k ++ ]=q[i ++];
    while(j <= r) tmp[k ++ ]=q[j ++];

    for(i =l,j=0;i<=r;i++,j++) q[i] = tmp[j];
}

int main()
{
    scanf("%d",&n);
    for(int i = 0;i < n; i ++ ) scanf("%d",&q[i]);

    merge_sort(q,0,n-1);

    for(int i=0;i<n;i++) printf("%d ",q[i]);

    return 0;
}

12、整数二分

来源于Acwing算法基础课-基础算法（一）

#include <bits/stdc++.h>

using namespace std;

const int N = 100010;

int n,m;
int q[N];

int main()
{
    scanf("%d%d",&n,&m);
    for(int i = 0;i < n; i ++ ) scanf("%d",&q[i]);

    while(m--)
    {
        int x;
        scanf("%d",&x);

        int l = 0,r = n - 1;
        while(l<r)
        {
            int mid = l+r >>1;
            if(q[mid]>=x) r=mid;//向左逼近答案
            else l=mid+1;
        }

        if(q[l]!=x) cout<<"-1 -1"<<endl;
        else
        {
            cout<<l<<' ';

            int l=0,r=n-1;
            while(l<r)
            {
                int mid=l+r+1 >>1;//需要加一，防止边界问题死循环
                if(q[mid]<=x) l=mid;//向右逼近答案
                else r=mid -1;
            }
            cout<<l<<endl;
        }
    }

    return 0;
}

13、浮点数二分

求一个数的平方根

#include <bits/stdc++.h>

using namespace std;

int main()
{
    double x;
    cin>>x;

    double l = 0,r=x;
    while(r-l>1e-6)
    {
        double mid = (l+r)/2;
        if(mid*mid>=x) r=mid;
        else l=mid;
    }
    printf("%lf\n",l);
    return 0;
}

14、高精度模板

模拟计算的过程
将数据存到数组里
加法：

#include <bits/stdc++.h>

using namespace std;

const int N = 1e6 + 10;

vector<int> add(vector<int> &A,vector<int> &B)
{
    vector<int>C;

    int t=0;
    for(int i=0;i<A.size()||i<B.size();i++)
    {
        if(i<A.size()) t+=A[i];
        if(i<B.size()) t+=B[i];
        C.push_back(t%10);
        t/=10;
    }

    if(t) C.push_back(1);
    return C;
}

int main()
{
    string a,b;
    vector<int>A,B;

    cin>>a>>b;
    for(int i=a.size()-1;i>=0;i--) A.push_back(a[i]-'0');
    for(int i=b.size()-1;i>=0;i--) B.push_back(b[i]-'0');

    auto C = add(A,B);

    for(int i = C.size()- 1;i>=0;i--) printf("%d",C[i]);
    return 0;
}

减法：

#include <bits/stdc++.h>

using namespace std;

const int N = 1e6 + 10;

bool cmp(vector<int> &A,vector<int> &B)
{
    if(A.size()!=B.size()) return A.size() > B.size();
    for(int i=A.size()-1;i>=0;i--)
        if(A[i] != B[i])
            return A[i] > B[i];
    return true;
}
vector<int> sub(vector<int>&A,vector<int>&B)
{
    vector<int>C;
    for(int i = 0,t=0;i<A.size();i++)
    {
        t=A[i]-t;
        if(i<B.size()) t-=B[i];
        C.push_back((t+10)%10);
        if(t<0) t=1;
        else t=0;
    }
    while(C.size() > 1&& C.back() == 0) C.pop_back();
    return C;
}

int main()
{
    string a,b;
    vector<int>A,B;
    cin>>a>>b;
    for(int i=a.size()-1;i>=0;i--) A.push_back(a[i]-'0');
    for(int i=b.size()-1;i>=0;i--) B.push_back(b[i]-'0');

    if(cmp(A,B))
    {
        auto C = sub(A,B);
        for(int i = C.size()- 1;i>=0;i--) printf("%d",C[i]);
    }
    else
    {
        auto C = sub(A,B);
        printf("-");
        for(int i = C.size()- 1;i>=0;i--) printf("%d",C[i]);
    }
    return 0;
}

乘法：

#include <bits/stdc++.h>

using namespace std;

vector<int> mul(vector<int> &A,int b)
{
    vector<int>C;
    int t=0;
    for(int i = 0;i<A.size()||t;i++)
    {
        if(i<A.size()) t += A[i]*b;
        C.push_back(t%10);
        t/=10;
    }
    return C;
}

int main()
{
    string a;
    vector<int>A;
    int b;
    cin>>a>>b;
    for(int i=a.size()-1;i>=0;i--) A.push_back(a[i]-'0');

    auto C = mul(A,b);

    for(int i=C.size() -1;i>=0;i--) printf("%d",C[i]);

    return 0;
}

除法：

#include <bits/stdc++.h>

using namespace std;

vector<int> div(vector<int> &A,int b,int &r)
{
    vector<int>C;
    r=0;
    for(int i= A.size()-1;i>=0;i--)
    {
        r=r*10+A[i];
        C.push_back(r/b);
        r%=b;
    }
    reverse(C.begin(),C.end());
    while(C.size()>1&&C.back()==0) C.pop_back();
    return C;
}
int main()
{
    string a;
    vector<int>A;
    int b;
    cin>>a>>b;
    for(int i=a.size()-1;i>=0;i--) A.push_back(a[i]-'0');
    int r;
    auto C = div(A,b,r);
    for(int i=C.size() -1;i>=0;i--) printf("%d",C[i]);
    cout<<endl<<r<<endl;
    return 0;
}

15、前缀和

#include <bits/stdc++.h>

using namespace std;

const int N = 100010;

int n,m;
int a[N],s[N];

int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);

    for(int i=1;i<=n;i++) s[i] = s[i-1] + a[i];

    while(m--)
    {
        int l,r;
        scanf("%d%d",&l,&r);
        printf("%d\n",s[r] - s[l-1]);
    }
    return 0;
}

二维前缀和:

#include <bits/stdc++.h>

using namespace std;

const int N = 1010;

int n,m,q;
int a[N][N],s[N][N];

int main()
{
    scanf("%d%d%d",&n,&m,&q);
    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            scanf("%d",&a[i][j]);

    for(int i=1;i<=n;i++)
        for(int j=1;j<=m;j++)
            s[i][j] = s[i-1][j] + s[i][j -1]-s[i-1][j-1] +a[i][j];
    while(q--)
    {
        int x1,y1,x2,y2;
        scanf("%d%d%d%d",&x1,&y1,&x2,&y2);
        printf("%d\n",s[x2][y2] - s[x1-1][y2] - s[x2][y1- 1]+s[x1-1][y1-1]);
    }
    return 0;
}

16、差分//未学会

#include <bits/stdc++.h>

using namespace std;

const int N = 100010;

int n,m;
int a[N],b[N];

void insert(int l,int r,int c)
{
    b[l]+=c;
    b[r+1]-=c;
}

int main()
{
    scanf("%d%d",&n,&m);
    for(int i=1;i<=n;i++) scanf("%d",&a[i]);

    for(int i=1;i<=n;i++) insert(i,i,a[i]);

    while(m--)
    {
        int l,r,c;
        scanf("%d%d%d",&l,&r,&c);
        insert(l,r,c);
    }
    for(int i=1;i<=n;i++) b[i] += b[i-1];

    for(int i=1;i<=n;i++) printf("%d",b[i]);

    return 0;
}

17、双指针算法

for(i=0,j=0;i<n;i++)
{
   while(j<i&&check(i,j))j++;
   //每道题目的具体逻辑
}

18、离散化//不太会

Acwing区间和

#include <bits/stdc++.h>

using namespace std;

typedef pair<int,int>PII;

const int N = 300010;

int n,m;
int a[N],s[N];

vector<int> alls;
vector<PII> add,query;

int find(int x)//一个二分查找的过程
{
    int l=0,r=alls.size()-1;
    while(l<r)
    {
        int mid = l+r >> 1;
        if(alls[mid]>=x) r=mid;
        else l=mid + 1;
    }
    return r+1;
}
int main()
{
    cin>>n>>m;
    for(int i=0;i<n;i++)//记录插入数据
    {
        int x,c;
        cin>>x>>c;
        add.push_back({x,c});
        alls.push_back(x);
    }
    for(int i=0;i<m;i++)//记录询问数据
    {
        int l,r;
        cin>>l>>r;
        query.push_back({l,r});
        
        alls.push_back(l);
        alls.push_back(r);
    }
    //去重
    sort(alls.begin(),alls.end());
    alls.erase(unique(alls.begin(),alls.end()),alls.end());
    //处理插入
    for(auto item : add)
    {
        int x=find(item.first);
        a[x] += item.second;
    }
    //处理前缀和
    for(int i=1;i<=alls.size();i++) s[i] = s[i - 1]+a[i];
    //处理询问
    for(auto item : query)
    {
        int l = find(item.first),r = find(item.second);
        cout<<s[r] - s[l - 1]<<endl;
    }
    return 0;
}

19、单调栈

AcWing-单调栈

#include <iostream>

using namespace std;

const int N = 100010;

int n;
int stk[N],tt;

int main()
{
    scanf("%d",&n);

    for(int i=0;i<n;i++)
    {
        int x;
        scanf("%d",&x);
        while(tt && stk[tt]>=x) tt --;
        if(tt)printf("%d ",stk[tt]);
        else printf("-1");

        stk[++tt]=x;
    }
    return 0;
}
nefu_red_stone
关注
专栏目录
AcWing算法题常用代码模板
Lyf_Ah的博客
11-18 6805
一、基础算法 1、常用代码模板 快速排序算法模板 —— 模板题 AcWing 785. 快速排序 算法步骤： ​ ①确定分界点 ​ ②调整区间 ​ ③递归处理左右两段 Tips： ​ ①如果输入数据量比较大的话，用scanf输入更快，不要使用cin。 ​ ②只要背诵一种万能模板即可避开所有边界问题。 ​ ③考试时一般不会用到快排，面试时喜欢考察。 ​ ④ #include<bits/stdc++.h>包含了目前c++所包含的所有头文件 。 ​ ⑤分治一般用不到，只要掌握快速排序和归并排序就行。
Python10个与数学有关的简单实例代码
weixin_59071011的博客
07-12 1557
注意：我用的python2.7，大家如果用Python3.0以上的版本，请记得在print()函数哦！如果因为版本问题评论的，不做回复哦！！！ 1.题目：有1、2、3、4个数字，能组成多少个互不相同且无重复数字的三位数？都是多少？ 程序分析：可填在百位、十位、个位的数字都是1、2、3、4。组成所有的排列后再去 掉不满足条件的排列。 程序源代码： # -*- coding: UTF-8 -*- for i in range(1,5): for j in range(1,5): for k in.
微信小程序数学题库代码
04-27
https://blog.csdn.net/saywhat_sayhello/article/details/80110735 配合文章食用
用代码解决小学数学题
m0_45443233的博客
02-20 1078
笔主分享有趣的小学题代码 ##求1-100能被7而不能被5整除所有整数 for i in range(1,101): if i%7==0 and i%5!=0: print(i) #求200以内的能被17整除的最大整数 a=[] for i in range(1,201): if i%7==0: a.append(i) print(max(a)) ...
算法题代码模板2
crazyliu的博客
05-13 179
// 静态单链表 int head, e[N], ne[N], idx; void init() { head = -1; idx = 0; } // 头插 void insert(int x) { e[idx] = x, ne[idx] = head, head = idx++; } void remove() { head = ne[head]; } // 栈 int stk[N], top = -1; stk[++top] = x; // 入栈 --top;.
一个简单的数学题
一个和代码打交道の小孩
08-02 511
题目描述 zyc最近迷上了数学，一天，dj想出了一道数学题来难住他。算出1/n,但zyc一时答不上来希望大家能编程帮助他。 输入 第一行整数T，表示测试组数。后面T行，每行一个整数 n (1&lt;=|n|&lt;=10^5). 输出 输出1/n. (是循环小数的,只输出第一个循环节). 样例输入 4 2 3 7 168 样例输出 0.5 0.3 0.142857 0....
数据标注平台后端代码
A_Coding_Man
12-31 494
数据标注平台后端代码1. 文件目录设置csv文件的内容aa_lable_data_server.pyab_db_info.pyac_create_tables.pyad_config.py 1. 文件目录设置 csv文件的内容 aa_lable_data_server.py from ad_config import get_param,param_error, app, db_execute, db_execute_error, \ resultmsg,
基于JavaSwing离散数学题库管理系统（源代码+SQL脚本+论文+外文翻译）
06-19
适合人群：工作中项目需要、模板框架实战学习、大学毕业设计参考、找工作项目实战积累等等 题目名称 离散数学题库管理系统 毕设地点 xx学院 为帮助教师解决出题难,特别是出一套质量高同时又适合不同层次学生的...
2024年五月最新五一数学建模常用代码模板以及论文
05-02
2024五一数学建模，最新五一数学建模常用代码模板以及论文 为了进一步提升学生实践创新能力、促进数学建模教育事业的发展、为经济社会发展中的一些实际问题提供解决方案，组委会研究决定举办2024年第二十一届五一...
JAVA+SQL离散数学题库管理系统(源代码+论文模板+外文翻译).zip
06-20
标题中有的内容，压缩包中就有，适用于各位学子在完成大作业、课题、毕设时进行参考。 这些项目来自国内外大神提供，可以作为毕业设计的项目参考，也可以作为大作业项目参考，如果需要对项目进行修改，需要具备一定...
JAVA+SQL离散数学题库管理系统(源代码+毕设文档+外文翻译).zip
04-19
JAVA+SQL离散数学题库管理系统(源代码+毕设文档+外文翻译).zip JAVA+SQL离散数学题库管理系统(源代码+毕设文档+外文翻译).zip JAVA+SQL离散数学题库管理系统(源代码+毕设文档+外文翻译).zip JAVA+SQL离散数学题库...
简单数学编程题
辉小歌的博客
03-19 2367
目录1019 数字黑洞问题 A: 守形数问题 B: 反序数问题 C: 百鸡问题问题 D: abc问题 E: 众数 1019 数字黑洞 https://pintia.cn/problem-sets/994805260223102976/problems/994805302786899968 思路: 用数组保存每位的数字，并用bool数组标记一下。 根据bool数组 判断各个数字是不是完全的相同， 相同根据题目输出 不相同用sort排序后 计算最大值 和最小值。 用差接着重复执行上面的步骤直到差为6174。
4.如何用代码写计算数学问题
m0_68521667的博客
07-14 1657
在上一节课中，我给大家讲述了如何用代码计算任意两个数相乘，这节课我将继续深入教学。 数学题一直都是一个难搞的问题，像从一加到100，利息啥的问题，我相信大家列出公式应该是轻而易举的，如果不是当我没说。但是这些公式计算起来极为复杂，这时候我们就可以利用电脑的算力将问题的答案解答出来。那么这些数学题的代码应该如何实现呢？听我细细道来。 我们就从1加到100的数学问题来解答，代码如下 我们从上面的代码可以看到，大都是我们前面学的内容，只有一个for语句不是。老规矩一定要写框架，框架是 #include .
【算法编程】小学数学题难倒博士
青峰碧陋室
05-27 6439
昨天在科学网上得知这样一个新闻《越南小学数学题难倒博士》，据悉题目来自越南保禄小学三年班，不过报道称该题难倒了上至博士下至家长，未免也太言过其实了。 题目描述 学生需要在下图表格中按由上至下、从左到右的顺序，填入1~9的数字，可重复填写，并按先乘除后加减（图中冒号代表除法）的运算法则，完成整条算式。 解题方法 显然，这题对于我们这种程序员来说完全不是问题，只要在大一上过C语言的
数学类编程题
qq_43003252的博客
07-23 1302
1. 求数值中的整数次方 给定一个double类型的浮点数base和int类型的整数exponent。求base的exponent次方。保证base和exponent不同时为0 思路：方法1，直接求，注意任意数的0次方等于1，负数的次方等于它的倒数（时间复杂度O(n)）。方法2，利用平方的性质，an=a(n/2) * a(n/2)或an=a^(n-1/2) * a^(n-1/2)*a（分别为偶数、奇数的情况）（时间复杂度O(logn)） 代码一: class Solution { public: d
标注 画框 转写,亲亲，你做错了没？
weixin_42456278的博客
03-27 3599
基本要求和说明：1.1 印刷体、表格和图形的不框；1.2 画框层级最多只有3层；第一级：1 (必框)第二级：1-1 、1-2、1-3等 (注意：以题号区分，和文字块的阅读顺序 )第三级：1-1-1、1-1-2等 (注意：第二级下面的每一行，一个框)1.3 框属性的继承标注方法：画好第1个框--->先标好属性(手写文本)--->开启界面右下角的继承属性--->选中框...
编程题_数学问题
RuanMeng的博客
03-14 1640
1.百万富翁问题 2.球半径和体积问题 3.计算糖果 4.跟奥巴马一起编程 5.另类加法
dp算法练习【7】
MADAO的博客
09-08 1372
给你一个按排序的整数数组nums，返回组成的新数组，要求也按排序。平方后，数组变为 [16,1,0,9,100]排序后，数组变为 [0,1,9,16,100]非常简单，思路平方后排序。优化：插入法。
数学建模c题matlab代码
最新发布
09-09
由于问题的多样性，没有统一的模板可以适用于所有C题的MATLAB代码。但是，我可以提供一个大致的框架，以及如何用MATLAB解决数学建模问题的一般步骤： 1. 问题理解：首先仔细阅读题目，理解要解决的问题以及所给定的...
评论 1
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值