算法模板（未完）

衣谷金锋

已于 2022-02-04 23:10:42 修改

阅读量470

点赞数

分类专栏：算法日记文章标签：算法

于 2022-01-09 20:00:50 首次发布

本文链接：https://blog.csdn.net/qq_52760858/article/details/122398093

版权

算法日记专栏收录该内容

1 篇文章 0 订阅

订阅专栏

文章目录

比赛前复制粘贴
最大公约数函数
质数筛法
分解质因数
约数个数以及约数之和
快速幂
快速乘
大数相乘取余运算（结合快速幂和快速和）
求逆元（快速幂）
堆（优先队列）
二分
BFS
并查集
前缀和
差分
归并排序（应用：逆序对的数量）
逆序对的数量
有用的函数

比赛前复制粘贴

#include <bits/stdc++.h>
using namespace std;
#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0)
typedef long long ll;
typedef pair<int,int> PII;
int dx[] = {-1,0,1,0},dy[] = {0,1,0,-1};

最大公约数函数

int gcd(int a, int b)
{
    return b?gcd(b, a%b):a;
}

质数筛法

求解1~n中有多少个素数
//埃氏筛法 （略）
//线性筛法 （因为比较叼就决定用他了）
const int N = 1e6+10;
int cnt,prime[N];
bool st[N];

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

//时间复杂度O(n), cnt为1~n中的质数个数，大约为n/ln n。
//prime数组存放的是素数

分解质因数

输入：
6
输出：
2 1
3 1

void get_divisors(int x)
{
    for(int i = 2; i <= x/i; i ++)
    {
        if(x % i == 0)
        {
            int s = 0;
            while(x % i == 0)
            {
                x /= i;
                s ++;
            }
            cout << i << " " << s << endl;
        }
    }
    if(x > 1) cout << x << " " << 1 << endl;     //输出大于根号x的约数
    cout << endl;
}

约数个数以及约数之和

//应用了上述的分解质因数 此处的s更新为哈希表存下了每一个i所对应的s = mp[i]


unordered_map<int,int> mp;
void get_divisors(int x)
{
    for(int i = 2; i <= x/i; i ++)
    {
        if(x % i == 0)
        {
            while(x % i == 0)
            {
                x /= i;
                mp[i] ++;
            }
        }
    }
    if(x > 1) mp[x] ++;
}

//约数个数 将每一个质因数的个数+1 相乘累积在一起
ll ans = 1;
for(auto i = mp.begin(); i != mp.end(); i ++)
{
    ans = ans*(i->second + 1)%mod;
}
// ans为该数的约数个数之和


//约数之和 6 = 2^1 * 3^1 约数之和1+2+3+6=12 算法： (2^0+2^1)*(3^0+3^1) = 12
ll ans = 1;
for(auto i = mp.begin(); i != mp.end(); i ++)
{
    ll a = i->first, b = i->second, t = 1;
    while(b --) t = (t*a+1)%mod;
    ans = ans*t%mod;
}
// ans为该数的约数之和

快速幂

ll quickpow(ll a,ll b) //a^b % p 利用a^b = (a^2)^b/2减少幂的次数
{
   ll ans = 1, base = a;
    while(b)
    {
        if(b & 1) ans = ans * base % mod;
        base = base * base % mod;
        b >>= 1;
    }
    return ans;
}

快速乘

ll quickmul(ll a,ll b) //(a+b) % p 用a*b = 2a*b/2减少加的次数
{
   ll ans = 0, base = a;
    while(b)
    {
        if(b & 1) ans = (ans + base) % mod;
        base = (base + base) % mod;
        b >>= 1;
    }
    return ans;
}

大数相乘取余运算（结合快速幂和快速和）

计算72356825451438的53727354327038次方 % 27335634901889

#include <iostream>
using namespace std;
typedef long long ll;

const ll mod = 27335634901889;

ll quickadd(ll a,ll b,ll mod)
{
    ll ans = 0, base = a;
    while(b)
    {
        if(b & 1)   ans = (ans+base)%mod;
        base = (base+base)%mod;
        b >>= 1;
    }
    return ans;
}

ll quickpow(ll a,ll b,ll mod)
{
    ll ans = 1, base = a;
    while(b)
    {
        if(b & 1) ans = quickadd(ans,base,mod); //此处 base*bae会溢出 使用快速和不会溢出
        base = quickadd(base,base,mod);
        b >>= 1;

    }
    return ans;
}
int main()
{
    cout << quickpow(72356825451438,53727354327038,mod) << endl;
    return 0;
}

求逆元（快速幂）

//定义：满足a/b ≡ a*x(mod p)式子时，x = b^-1; 则称x为b的模p逆元
/*
b存在乘法逆元的充要条件是 b 与模数 p 互质。否则不存在乘法逆元。
且当模数 m 为质数时，b^m−2 即为 b 的乘法逆元。
*/
求 b模p的乘法逆元就是 (b^(p-2)) % p
  使用快速幂即可

堆（优先队列）

#include <bits/stdc++.h>
using namespace std;

priority_queue <int, vector<int>, greater<int>> q; //优先队列，从小到大排序
//priority_queue<int> 默认为priority_queue<int, vector<int>, less<int>> 从大到小排序

int main()
{
    int n; cin >> n;
    for(int i = 1; i <= n; i ++)
    {
        int x; cin >> x;
        q.push(x);
    }
    int ans = 0;
    while(q.size()!=1)
    {
        int a = q.top(); q.pop();
        int b = q.top(); q.pop();
        ans += (a+b);
        q.push(a+b);
    }
    cout << ans << endl;
    return 0;
}

二分

二分后 满足左边 去掉右边区域
while(l < r)
{
    int mid = l + r >> 1;
    if(check()) r = mid;
    else l = mid + 1;
}

二分后 满足右边 去掉左边区域
while(l < r)
{
    int mid = l + r + 1 >> 1;
    if(a[mid] <= k) l = mid;
    else r = mid-1;
}

BFS

采用队列来做
1.起点和终点 不一定是按题目来
2.如果起点有多个x,y，先把起点全部先push进去，bfs里面再直接while（q.size()){...}

#include <iostream>
#include <cstring>
#include <queue>
using namespace std;

typedef pair<int,int> PII;
const int N = 110;
int a[N][N],d[N][N];
int n,m;

int bfs()
{
    queue<PII>q;
    q.push({0,0});
    memset(d,-1,sizeof(d));
    d[0][0] = 0;
    int dx[]={-1,0,1,0},dy[]={0,1,0,-1};
    while(q.size())
    {
        PII t = q.front();
        q.pop();
        for(int i = 0; i < 4; i ++)
        {
            int x = t.first+dx[i],y = t.second+dy[i];
            if(x>=0 && x<=n-1 && y>=0 && y<=m-1 && d[x][y]==-1 && a[x][y]==0)
            {
                d[x][y] = d[t.first][t.second]+1;
                q.push({x,y});
            }
        }
    }
    return d[n-1][m-1];
}
int main()
{
    cin >> n >> m;
    for(int i = 0; i < n; i ++)
        for(int j = 0; j < m; j ++)
            cin >> a[i][j];
    cout << bfs() << endl;
    return 0;
}

并查集

int find(int x)
{
    if(p[x]!=x) p[x]=find(p[x]);
    return p[x];
}

前缀和

1.一维数组a[n]的[l,r]前缀和
//增加sum[n]数组

for(int i = 1; i < n; i ++) sum[i] = sum[i-1] + a[i];
ans = sum[r] - sum[l-1];

2.二维数组a[n][n]的[(x1,y1),(x2,y2)]前缀和
//增加sum[n][n]数组
for(int i = 1; i < n; i ++)
    for(int j = 1; j < n; j ++)
        sum[i][j] = sum[i-1][j] + sum[i][j-1] - sum[i-1][j-1] + a[i][j];
ans = sum[x1][y1] - sum[x1-1][y2] - sum[x2][y1-1] + sum[x1-1][y1-1];

差分

1.一维数组a[n]在[l,r]上每个数加上常数c
//使用一个新的数组b[n],使得b[n]的前缀和为a[n],这样将b[l]+=c,b[r+1]-=c达到目的

void insert(int l, int r, int c)
{
    b[l] += c;
    b[r+1] -= c;
}
for(int i = 1; i <= n; i ++) insert(i,i,a[i]); //这样b[n]构造完成 b[n]的前缀和就是a[n]
//如果在[l,r]插入常数c,  insert(l,r,c);

2.差分矩阵
//未看

归并排序（应用：逆序对的数量）

const int N = 1e5+10;
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 = 1, 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(int i = l,j = 1; i <= r; i ++,j ++)
        q[i] = tmp[j];
}

merge_sort(q,1,n);
//即进行归并排序 此处是从小到大排序，将以下代码修改 则是从大到小
while(i <= mid && j <= r)
    if(q[i] >= q[j]) tmp[k ++] = q[i ++];
     else tmp[k ++] = q[j ++];

逆序对的数量

const int N = 1e5+10;
int q[N],tmp[N];
ll ans;    //记录逆序对数量，采取边归并边记录数量的方式

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 = 1, i = l,j = mid+1;
    while(i <= mid && j <= r)
        if(q[i] <= q[j]) tmp[k ++] = q[i ++];
        else{
            ans += mid-i+1;
            tmp[k ++] = q[j ++];
        }
    while(i <= mid) tmp[k ++] = q[i ++];
    while(j <= r) tmp[k ++] = q[j ++];

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

有用的函数

ceil(x)     --向上取整
floor(x)     --向下取整
acos(x)     --可用于求pi = acos（-1）
lowbit(s)    --求一个数二进制最低的1代表的数，6对应的110 对应的就是10 也就是2。lowbit（6）= 2
/*
#define lowbit(x) (x & -x) //函数实现

cin >> x;
while(x)
{
    if(x == lowbit(x)) high = lowbit(x); //为最高位1所代表的数
    x-= lowbit(x);
    ans ++;
}
ans的值即为x二进制数1的个数

二进制最后一位取反 等价于 ^ 1 与1异或
判断奇偶 等价于& 1  与1
*/

衣谷金锋

关注

0
点赞
踩
0

收藏

觉得还不错? 一键收藏
打赏
0
评论
算法模板（未完）

文件前复制粘贴#include <bits/stdc++.h>using namespace std;#define rep(i,a,n) for(int i = a; i <= n; i ++)#define IOS ios::sync_with_stdio(false);cin.tie(0);cout.tie(0)typedef long long ll;typedef pair<int,int> PII;文章目录文件前复制粘贴最大公约数函数质数筛法快速幂快
复制链接

扫一扫