先贴SPFA算法代码:(和我原来写的dijsktra算法非常相似)
bool queue<int> q,int n) //配套spfa函数用以检查queue中是否存在元素n
{
while(!q.empty())
{
int a = q.front();
q.pop();
if(a==n) return 1;
}
return 0;
}
vector<int> spfa(vector<vector<int>> G,int source)
{
int n = G.size();
vector<int> road(n,inf);
road[source] = 0;
queue<int> q;
q.push(source);
while(!q.empty())
{
int front = q.front();
q.pop();
for(int i=0;i<n;i++)
if(road[front]+G[front][i]<road[i] && front!=i)
{
road[i] = road[front]+G[front][i];
if(!spfa_check(q,i)) q.push(i);
}
}
return road;
}
再贴一个欧几里得算法(coding学习自b站大牛“大雪菜”)
int gcd(int a,int b) //欧几里得算法,辗转相除求最大公约数
{
return b?gcd(b,a%b):a;
}
最后贴一个完整解决该问题main函数代码:
int main()
{
vector<vector<int>> G(2021);
for(int i=1;i<=2021;i++)
{
for(int j=1;j<=2021;j++)
{
if(fabs(j-i)<=21){
int m = gcd(j,i);
int mn = i*j/m;
G[i-1].push_back(mn);
}
else
{
G[i-1].push_back(inf);
}
}
}
vector<int> a = spfa(G, 0);
cout<<a[2020]<<endl;
return 0;
}
最后算出答案为:10266837
dijsktra算法模版见:https://blog.csdn.net/qq_39743377/article/details/123992565?spm=1001.2014.3001.5501
同样适用于该题,但是dijkstra只是用没有负权值的最短路
该题不做讲解, 没啥好讲的, 纯纯考最短路写代码
讲一讲SPFA算法:
SPFA有“辗转”的味道, 先说基本思路
从原点开始让一个点先进入队列 (queue队列是先进先出), 如果这个点连接的任何一个点能让现存的路径值更小, 然后让所有满足上述要求的点进入这个队列**(如果队列里不存在这个点的话)**, 并且更新权值. 周而复始, 直到队列里没有任何一个元素为止,说明此时没有权值可以更新 (一定收敛)
SPFA算法适用于单源负权值最短路搜索(没有负环)
更详细的SPFA算法来日再更新,可以先看大佬的文章:
http://www.360doc.com/content/13/1208/22/14357424_335569176.shtml