经典算法与数据结构实现-CSDN博客

本文涵盖多种核心算法实现，包括ST表、最大公约数、最小公倍数、逆元等数学算法，以及高精度计算、扩展欧几里得算法、SPFA最短路径算法、LCA算法、状压DP、树状数组、背包问题等多种数据结构和算法的应用。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

ST表

for(int j=1;1<<j<=n;j++)
	    for(int i=1;i+(1<<j)-1<=n;i++){
	        f1[i][j]=max(f1[i][j-1],f1[i+(1<<j-1)][j-1]);                                            
	        f2[i][j]=min(f2[i][j-1],f2[i+(1<<j-1)][j-1]);
	    }

最大公约数

int gcd(int n,int m){
    if(n%m==0) return m;
    else return gcd(m,n%m);
}

最小公倍数=a*b/gcd(a,b)

逆元

long long exgcd(long long a,long long b,long long &x,long long &y){
    if(b==0){
        x=1;
        y=0;
        return a;
    }
    LL gcd=exgcd(b,a%b,y,x);
    y-=x*(a/b);
    return gcd;
}

二元一次不定方程求解，乘法

高精度

读入：

void read(int x[]){
	string s;
	cin>>s;
	x[0]=s.size();
	for(int i=0;i<x[0];i++)
		x[x[0]-i]=s[i]-'0';
}

加：

int add(int z[],int x[],int y[]){
	int g=0;
	int len=max(x[0],y[0]);
	for(int i=1;i<=len;i++){
		z[i]=x[i]+y[i]+g;
		g=z[i]/10;
		z[i]%=10;
	}
	if(g>0)
		z[++len]=g;
	z[0]=len;
}

减：

int cut(int z[],int x[],int y[]){
	int len=x[0];
	for(int i=1;i<=len;i++){
		if(x[i]<y[i]){
			x[i]+=10;
			x[i+1]--;
		}
		z[i]=x[i]-y[i];
	}
	while(len>1&&z[len]==0)
		len--;
	z[0]=len;
}

乘：

void mulbbign(int z[],int x[],int y[]){
	int lenz=x[0]+y[0]-1;
	for(int i=1;i<=x[0];i++){
		for(int j=1;j<=y[0];j++){
			z[i+j-1]+=x[i]*y[j];
			z[i+j]+=z[i+j-1]/10;
			z[i+j-1]%=10;
		}
	}
	if(z[lenz+1]!=0)
		lenz++;
	while(lenz>1&&z[lenz]==0)
		lenz--;
	z[0]=lenz;
}

扩展欧几里得

void exgcd(int a,int b){
    if(b==0){
        x=1;
        y=0;
        return;
    }
    exgcd(b,a%b);
    k=x;
    x=y;
    y=k-a/b*y;
    return;
}

SPFA

#include<bits/stdc++.h>
const int INF = 0x3f3f3f3f;
const int N = 10000+10;

using namespace std;

struct Node{
	int v,w;
};
vector<Node> E[N]; 
int n,m,s,t;
int dis[N];
bool vis[N];

void SPFA(int s){
	for(int i = 1;i <= n; ++i) 
		vis[i] = false,dis[i] = INF;
	queue<int> que;
	que.push(s);
	dis[s] = 0,vis[s] = true;
	while(!que.empty()){
		int t = que.front();
		que.pop();
		vis[t] = false;
		for(int i = 0,l = E[t].size();i < l; ++i) {
			int j = E[t][i].v;
			int k = E[t][i].w;
			if(dis[j] > dis[t] + k){
				dis[j] = dis[t] + k;
				if(!vis[j]){
					vis[j] = true;
					que.push(j);
				}
			}
		}
	}
}
int main()
{
	cin>>n>>m>>s>>t;
	int u,v,w;
	for(int i = 1;i <= m; ++i) {
		cin>>u>>v>>w;
		E[u].push_back({v,w});
		E[v].push_back({u,w});
	}
	
	SPFA(s);
	if(dis[t] >= INF / 2) cout<<"-1"<<endl;
	else cout<<dis[t]<<endl;
	return 0;
}

LCA

状压DP

树状数组

#include<bits/stdc++.h>
#define lowbit(x) (x&(-x))
typedef long long ll; 
using namespace std;
int c[2000010];
int n,m;
ll ans;
int add_dandian(int x,int k)
{
	for(int i=x;i<=n;i+=lowbit(i))
	c[i]+=k;
}
int search(int begin,int end)
{
	for(int i=end;i;i-=lowbit(i))
	ans+=c[i];
	for(int i=begin-1;i;i-=lowbit(i))
	ans-=c[i];
	return 0;
}
int main()
{
	scanf("%d %d",&n,&m);
	for(int i=1;i<=n;i++)
	{
		int number;
		scanf("%d",&number);
		add_dandian(i,number);
	}
	for(int i=1;i<=m;i++)
	{
		int choice,x,y;
		scanf("%d %d %d",&choice,&x,&y);
		if(choice==1) add_dandian(x,y);
		else
		{
			ans=0;
			search(x,y);
			printf("%lld\n",ans);
		}
	}
	return 0;
}

#include<bits/stdc++.h>
#define lowbit(x) (x&(-x)) 
using namespace std;
typedef long long ll;
const int Maxn=1e6+10;
int a[500010];
int d[Maxn]={0};
ll c[Maxn]; 
int n,m;
int update(int pos,int k)
{
	for(int i=pos;i<=n;i+=lowbit(i))
	c[i]+=k;
	return 0;
}
ll ask_qujian(int pos)
{
	ll ans=0;
	for(int i=pos;i;i-=lowbit(i)) ans+=c[i];
	return ans;
} 
int main()
{
	memset(c,0,sizeof(c));
	scanf("%d %d",&n,&m);
	a[0]=0;
	for(int i=1;i<=n;i++)
	{
		scanf("%d",&a[i]);
		d[i]=a[i]-a[i-1];
		update(i,d[i]);
	}
	for(int i=1;i<=m;i++)
	{
		int choice,x,y,k;
		scanf("%d",&choice);
		if(choice==1)
		{
			scanf("%d %d %d",&x,&y,&k);
			update(x,k);
			update(y+1,-k);
		}
		else
		{
			scanf("%d",&x);
			printf("%lld\n",ask_qujian(x));
		}
	}
	return 0;
}

最短路

dijkstra

void Dijkstra(int x){
	memset(dist,INF,sizeof(dist));
	dist[x]=0;
	memset(visit,0,sizeof(visit));
	for(int i=1;i<=n;i++){
		int j=0;
		for(int k=1;k<=n;k++){
			if(!visit[k]&&dist[k]<=dist[j]){
				j=k;
			}
		}
		visit[j]=1;
		for(int k=1;k<=n;k++){
			dist[k]=min(dist[k],dist[j]+map[j][k]);
		}
	}
}

floyd

#include<bits/stdc++.h>
using namespace std;
const int inf = 0x7fffff - 1;
int e[1000][1000];
int n, m;
int main() {
	cin >> n >> m;
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			if (i == j)e[i][j] = 0;
			else e[i][j] = inf;
		}
	}
	int src, dst, val;
	for (int i = 0; i < m; i++) {
		cin >> src >> dst >> val;
		e[src][dst] = val;
	}
	for (int k = 0; k < n; k++) {
		for (int i = 0; i < n; i++) {
			for (int j = 0; j < n; j++) {
				if (e[i][k] + e[k][j] < e[i][j]) {
					e[i][j] = e[i][k] + e[k][j];
				}
			}
		}
	}
	for (int i = 0; i < n; i++) {
		for (int j = 0; j < n; j++) {
			printf("%5d", e[i][j]);
		}
		cout << endl;
	}
	return 0;
}

背包

无限

#include <bits/stdc++.h>
using namespace std;
long long n,m,a[100010],b[100010],f[10000010];
int main(){
	cin>>m>>n;
	for(int i=1;i<=n;i++)
		cin>>a[i]>>b[i];
	for(int i=1;i<=n;i++)
		for(int j=a[i];j<=m;j++)
			f[j]=max(f[j],f[j-a[i]]+b[i]);
	cout<<f[m];
	return 0;
}

#include <bits/stdc++.h>
using namespace std;
long long n,m,a[100010],b[100010],f[10000010];
int main(){
	cin>>m>>n;
	for(int i=1;i<=n;i++)
		cin>>a[i]>>b[i];
	for(int i=1;i<=n;i++)
		for(int j=m;j>=a[i];j--)
			f[j]=max(f[j],f[j-a[i]]+b[i]);
	cout<<f[m];
	return 0;
}