第十二届蓝桥杯省赛 E:路径

分析：这道题目直接按照题目所给条件进行建图，然后从1号点跑一个最短路即可。可能有小伙伴不知道怎么求解两个数a和b的最小公倍数，下面我着重讲解一下最小公倍数和最大公约数之间的关系：

对a和b分别进行质因子分解有：

a = ( p1^x1 ) * ( p2^x2 ) * ( p3^x3 ) * …… * ( pn^xn )

b = ( p1^y1 ) * ( p2^y2 ) * ( p3^y3 ) * …… * ( pn^yn )

则 lcm( a , b ) = ( p1^max (x1 , y1) ) * ( p2^max (x2 , y2) ) * ( p3^ max (x3 , y3) * …… * pn^max (xn , yn) )

而gcd( a, b ) = ( p1^min (x1 , y1) ) * ( p2^min (x2 , y2) ) * ( p3^ min (x3 , y3) * …… * pn^min (xn , yn) )

那么我们可以发现对于任意一个i均有max (xi , yi)+min (xi , yi)=xi+yi，那么我们可以得到lcm( a , b )*gcd( a, b )=p1^(x1+y1) * p2^(x2+y2) *p3^(x3+y3)* …… *pn^(xn+yn)我们可以发现这就是a*b，那么我们求解lcm(a,b)时可以用a*b/gcd(a,b)进行求解。

至于gcd(a,b)如何求解我们是利用gcd(a,b)=gcd(b%a,a)来递归求解的，为什么这个等式会成立呢？下面来给出证明：

首先证明gcd(a,b)是gcd(b%a,a)的因子，假设r=gcd(a,b)

那么r肯定是a的因子，所以我们只需要证明r是b%a的因子即可，因为b%a=b-b/a*a（这里的/是下取整的），那么我们可以发现b是r的倍数，而且a也是r的倍数，那么b和a的线性组合也是r的倍数，所以就得到了gcd(a,b)是gcd(b%a,a)的因子

下面我们再来证明gcd(b%a,a)是gcd(a,b)的因子

假设r=gcd(b%a,a)，那么r肯定是a的因子，接下来只需要证明r是b的因子即可

由于b%a=b-b/a*a,由于r是a的因子，所以r肯定也是b/a*a的因子，又因为r是b%a的因子，所以r肯定是b%a+b/a*a的因子，也就是b的因子

两个数互为因子，那么这两个数肯定相同，证毕！

细节见代码：

#include<cstdio>
#include<iostream>
#include<algorithm>
#include<cstring>
#include<map>
#include<queue>
#include<vector>
#include<cmath>
using namespace std;
const int N=1e7+10;
bool vis[N];
int e[N],ne[N],h[N],w[N],idx;
int d[N];
typedef pair<int,int> PII;
void add(int x,int y,int z)
{
	e[idx]=y;
	w[idx]=z;
	ne[idx]=h[x];
	h[x]=idx++;
}
int gcd(int x,int y)
{
	if(x==0) return y;
	return gcd(y%x,x);
}
int lcm(int x,int y)
{
	return x*y/gcd(x,y);
}
void dijkstra(int x)
{
	priority_queue<PII,vector<PII>,greater<PII> >q;
	memset(d,0x3f,sizeof d);
	d[x]=0;
	q.push({d[x],x});
	while(!q.empty())
	{
		int t=q.top().second;
		q.pop();
		if(vis[t]) continue;
		vis[t]=true;
		for(int i=h[t];i!=-1;i=ne[i])
		{
			int j=e[i];
			if(d[j]>d[t]+w[i])
			{
				d[j]=d[t]+w[i];
				q.push({d[j],j});
			}
		}
	}
}
int main()
{
	for(int i=1;i<=2021;i++)
		h[i]=-1;
	for(int i=1;i<=2021;i++)
	for(int j=i+1;j<=min(i+21,2021);j++)
	{
		int t=lcm(i,j);
		add(i,j,t);
		add(j,i,t);
	}
	dijkstra(1);
	printf("%d",d[2021]);
	return 0;
}