GCD Path CH Round #53 -【Nescafé 32】杯NOIP模拟赛
描述
给定一张N个点的有向图,点i到点j有一条长度为 i/(gcd(i,j))的边。有Q个询问,每个询问包含两个数x和y,求x到y的最短距离。
输入格式
第一行包含两个用空格隔开的整数,N和Q。
接下来Q行,每行两个数x和y。
输出格式
输出Q行整数,表示从x到y的最短距离。
样例输入
6 2 4 6 2 5
样例输出
2 2
数据范围与约定
- 对于30%的数据,1<=N<=100。
- 对于70%的数据,1<=N<=10^5。
- 对于100%的数据,1<=N<=10^7,1<=x,y<=N,Q<=10^5。
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++
这道题刚开始看用Floyd写了一遍,可是只过30……然后看别人的题解,发现原来用数学的知识…
假设求x->y的最短路,则直接走 这条路 长度为x/gcd(x,y)设这个长度的质因数分解为a1*a2*a3*a4……(两项可以相等)
然后要用到一个结论:
若a>=2,b>=2,则a+b<=a*b
移项就是 (1-a)*(1-b)>=1 这是显然的。
所以我们不妨把这个长度分开来走,
因为每拆一项都会使答案减小或不边,那我们不妨直接将该数全部分解为质数,一个一个质数来走。
举个例子
100-1,则100/gcd(100,1)的质因数分解为2*2*5*5
我们不妨使每次走的长度为2 2 5 5,而2+2+5+5=14<100 这样使长度之和达到最小。
所以我们可以这样走 100->50->25->5->1->1
所以就有了下面的代码……可是发现会有错误答案……
经过仔细的研究我发现有一种情况:自己到自己……按题目的要求应该是1,可是很明显答案是0……
然后我又加了一行if,最后过了……
哎……题目有歧义啊!!
+++++++++++++++++++++++++++++++++++++++++++++++++++++++++
#include <iostream>
#include <cstdio>
#include <cmath>
#include <cstring>
using namespace std;
int n,q,x,y;
int is_prime(int x)
{
for(int i=2;i<=(int)sqrt(x);i++)
if (x%i==0)return 0;
return 1;
}
int pr(int x)
{
int sum=0;
if (is_prime(x))return x;
int a=2;
while(a*a<=x)
{
while(x%a==0)
{
sum+=a;
x=x/a;
}
a++;
}
if(x>1)sum+=x;
return sum;
}
int i_gcd(int x,int y)
{
int ii=x;
if (!x||!y)return x>y?x:y;
for(int i;i=x%y;x=y,y=i);
return pr(ii/y);
}
int main()
{
scanf("%d%d",&n,&q);
for (int i=1;i<=q;i++)
{
scanf("%d%d",&x,&y);
if(x==y)printf("0\n");
else printf("%d\n",i_gcd(x,y));
}
return 0;
}